coherentstructures / coherentstructures.jl Goto Github PK
View Code? Open in Web Editor NEWTools for computing (Lagrangian) coherent structures objectively in Julia
License: Other
coherentstructures.jl's People
Contributors
![github-actions[bot] avatar](https://avatars.githubusercontent.com/in/15368?v=4 "github-actions[bot]")
Stargazers
Watchers
Forkers
hijiker edmondztt m20190649 mahermann chanjeunlam moredread gaioguy adrilow blaszm jamesthesnakecoherentstructures.jl's Issues
Cannot index SymmetricTensor
julia> VERSION
v"1.7.1"
And CoherentStructures is up to date.
In this example https://coherentstructures.github.io/CoherentStructures.jl/stable/generated/bickley/#Bickley-jet
a tensor D is defined for a non-identity diffusion tensor, which causes this error
C̅ = pmap(mCG_tensor, P; batch_size=ceil(Int, length(P)/nprocs()^2))
ERROR: On worker 2:
MethodError: objects of type SymmetricTensor{2, 2, Float64, 3} are not callable
Use square brackets [] for indexing an Array.
If I do not pass D by doing
mCG_tensor = u -> av_weighted_CG_tensor(bickley, u, tspan, δ;
tolerance=1e-6, solver=Tsit5())
instead of
mCG_tensor = u -> av_weighted_CG_tensor(bickley, u, tspan, δ;
D=D, tolerance=1e-6, solver=Tsit5())
the problem disappears.
OECS results in intersecting eddies
I did the example for OECS shown here but with the following changes:
xmin, xmax, ymin, ymax = minimum(xs), maximum(xs), minimum(ys), maximum(ys)
and (note I am using the fifth time instead of the first as in the example):
const V = scale(interpolate(SVector{2}.(us[:,:,5], vs[:,:,5]), BSpline(Quadratic(Free(OnGrid())))), xs, ys)
I get two eddies that intersect each other between 25 and 30 South and around 6 West (figure below), making it impossible for these "eddy" boundaries to be transport barriers. What could be the problem?
Thanks.
Replace Arpack.jl by a pure-Julia package
Options:
- KrylovKit.jl
- ArnoldiMethod.jl
- IterativeSolvers.jl
Support for Hexahedral elements
As discussed via mail, extending pointlocator and gridcontext objects to support 3d Hexahedron elements would be very helpful.
nested task error: BoundsError: attempt to access Tuple{Float64, Float64} at index [3]
@time res = CoherentStructures.plot_ftle( # uv_tricubic,p, my_interp, p2, [times[1], times[end]], LL, UR, 100, 100; tolerance=1e-6, aspect_ratio=1, color=:rainbow, pass_on_errors=true );
TaskFailedException
nested task error: BoundsError: attempt to access Tuple{Float64, Float64} at index [3]
Stacktrace:
[1] getindex
@ .\tuple.jl:29 [inlined]
[2] getindex
@ C:\Users\51342\.julia\packages\StaticArrays\LJQEe\src\SVector.jl:39 [inlined]
[3] my_interp(u::SVector{2, Float64}, p::Tuple{Interpolations.Extrapolation{Float64, 3, ScaledInterpolation{Float64, 3, Interpolations.BSplineInterpolation{Float64, 3, OffsetArrays.OffsetArray{Float64, 3, Array{Float64, 3}}, BSpline{Cubic{Free{OnGrid}}}, Tuple{Base.OneTo{Int64}, Base.OneTo{Int64}, Base.OneTo{Int64}}}, BSpline{Cubic{Free{OnGrid}}}, Tuple{StepRangeLen{Float64, Base.TwicePrecision{Float64}, Base.TwicePrecision{Float64}}, StepRangeLen{Float64, Base.TwicePrecision{Float64}, Base.TwicePrecision{Float64}}, StepRangeLen{Float64, Base.TwicePrecision{Float64}, Base.TwicePrecision{Float64}}}}, BSpline{Cubic{Free{OnGrid}}}, Tuple{Periodic{Nothing}, Periodic{Nothing}, Flat{Nothing}}}, Interpolations.Extrapolation{Float64, 3, ScaledInterpolation{Float64, 3, Interpolations.BSplineInterpolation{Float64, 3, OffsetArrays.OffsetArray{Float64, 3, Array{Float64, 3}}, BSpline{Cubic{Free{OnGrid}}}, Tuple{Base.OneTo{Int64}, Base.OneTo{Int64}, Base.OneTo{Int64}}}, BSpline{Cubic{Free{OnGrid}}}, Tuple{StepRangeLen{Float64, Base.TwicePrecision{Float64}, Base.TwicePrecision{Float64}}, StepRangeLen{Float64, Base.TwicePrecision{Float64}, Base.TwicePrecision{Float64}}, StepRangeLen{Float64, Base.TwicePrecision{Float64}, Base.TwicePrecision{Float64}}}}, BSpline{Cubic{Free{OnGrid}}}, Tuple{Periodic{Nothing}, Periodic{Nothing}, Flat{Nothing}}}}, t::Float64)
@ Main .\In[1]:42
[4] (::SciMLBase.ODEFunction{false, typeof(my_interp), LinearAlgebra.UniformScaling{Bool}, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, typeof(SciMLBase.DEFAULT_OBSERVED), Nothing})(::SVector{2, Float64}, ::Vararg{Any, N} where N)
@ SciMLBase C:\Users\51342\.julia\packages\SciMLBase\h4Gxc\src\scimlfunctions.jl:334
[5] arraymap2(u::SVector{10, Float64}, p::Tuple{Interpolations.Extrapolation{Float64, 3, ScaledInterpolation{Float64, 3, Interpolations.BSplineInterpolation{Float64, 3, OffsetArrays.OffsetArray{Float64, 3, Array{Float64, 3}}, BSpline{Cubic{Free{OnGrid}}}, Tuple{Base.OneTo{Int64}, Base.OneTo{Int64}, Base.OneTo{Int64}}}, BSpline{Cubic{Free{OnGrid}}}, Tuple{StepRangeLen{Float64, Base.TwicePrecision{Float64}, Base.TwicePrecision{Float64}}, StepRangeLen{Float64, Base.TwicePrecision{Float64}, Base.TwicePrecision{Float64}}, StepRangeLen{Float64, Base.TwicePrecision{Float64}, Base.TwicePrecision{Float64}}}}, BSpline{Cubic{Free{OnGrid}}}, Tuple{Periodic{Nothing}, Periodic{Nothing}, Flat{Nothing}}}, Interpolations.Extrapolation{Float64, 3, ScaledInterpolation{Float64, 3, Interpolations.BSplineInterpolation{Float64, 3, OffsetArrays.OffsetArray{Float64, 3, Array{Float64, 3}}, BSpline{Cubic{Free{OnGrid}}}, Tuple{Base.OneTo{Int64}, Base.OneTo{Int64}, Base.OneTo{Int64}}}, BSpline{Cubic{Free{OnGrid}}}, Tuple{StepRangeLen{Float64, Base.TwicePrecision{Float64}, Base.TwicePrecision{Float64}}, StepRangeLen{Float64, Base.TwicePrecision{Float64}, Base.TwicePrecision{Float64}}, StepRangeLen{Float64, Base.TwicePrecision{Float64}, Base.TwicePrecision{Float64}}}}, BSpline{Cubic{Free{OnGrid}}}, Tuple{Periodic{Nothing}, Periodic{Nothing}, Flat{Nothing}}}}, t::Float64, odefun::SciMLBase.ODEFunction{false, typeof(my_interp), LinearAlgebra.UniformScaling{Bool}, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, typeof(SciMLBase.DEFAULT_OBSERVED), Nothing})
@ CoherentStructures C:\Users\51342\.julia\packages\CoherentStructures\hfDt1\src\util.jl:84
[6] (::CoherentStructures.var"#56#66"{SciMLBase.ODEFunction{false, typeof(my_interp), LinearAlgebra.UniformScaling{Bool}, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, typeof(SciMLBase.DEFAULT_OBSERVED), Nothing}})(u::SVector{10, Float64}, p::Tuple{Interpolations.Extrapolation{Float64, 3, ScaledInterpolation{Float64, 3, Interpolations.BSplineInterpolation{Float64, 3, OffsetArrays.OffsetArray{Float64, 3, Array{Float64, 3}}, BSpline{Cubic{Free{OnGrid}}}, Tuple{Base.OneTo{Int64}, Base.OneTo{Int64}, Base.OneTo{Int64}}}, BSpline{Cubic{Free{OnGrid}}}, Tuple{StepRangeLen{Float64, Base.TwicePrecision{Float64}, Base.TwicePrecision{Float64}}, StepRangeLen{Float64, Base.TwicePrecision{Float64}, Base.TwicePrecision{Float64}}, StepRangeLen{Float64, Base.TwicePrecision{Float64}, Base.TwicePrecision{Float64}}}}, BSpline{Cubic{Free{OnGrid}}}, Tuple{Periodic{Nothing}, Periodic{Nothing}, Flat{Nothing}}}, Interpolations.Extrapolation{Float64, 3, ScaledInterpolation{Float64, 3, Interpolations.BSplineInterpolation{Float64, 3, OffsetArrays.OffsetArray{Float64, 3, Array{Float64, 3}}, BSpline{Cubic{Free{OnGrid}}}, Tuple{Base.OneTo{Int64}, Base.OneTo{Int64}, Base.OneTo{Int64}}}, BSpline{Cubic{Free{OnGrid}}}, Tuple{StepRangeLen{Float64, Base.TwicePrecision{Float64}, Base.TwicePrecision{Float64}}, StepRangeLen{Float64, Base.TwicePrecision{Float64}, Base.TwicePrecision{Float64}}, StepRangeLen{Float64, Base.TwicePrecision{Float64}, Base.TwicePrecision{Float64}}}}, BSpline{Cubic{Free{OnGrid}}}, Tuple{Periodic{Nothing}, Periodic{Nothing}, Flat{Nothing}}}}, t::Float64)
@ CoherentStructures C:\Users\51342\.julia\packages\CoherentStructures\hfDt1\src\pullbacktensors.jl:145
[7] (::SciMLBase.ODEFunction{false, CoherentStructures.var"#56#66"{SciMLBase.ODEFunction{false, typeof(my_interp), LinearAlgebra.UniformScaling{Bool}, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, typeof(SciMLBase.DEFAULT_OBSERVED), Nothing}}, LinearAlgebra.UniformScaling{Bool}, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, typeof(SciMLBase.DEFAULT_OBSERVED), Nothing})(::SVector{10, Float64}, ::Vararg{Any, N} where N)
@ SciMLBase C:\Users\51342\.julia\packages\SciMLBase\h4Gxc\src\scimlfunctions.jl:334
[8] initialize!(integrator::OrdinaryDiffEq.ODEIntegrator{OrdinaryDiffEq.BS5, false, SVector{10, Float64}, Nothing, Float64, Tuple{Interpolations.Extrapolation{Float64, 3, ScaledInterpolation{Float64, 3, Interpolations.BSplineInterpolation{Float64, 3, OffsetArrays.OffsetArray{Float64, 3, Array{Float64, 3}}, BSpline{Cubic{Free{OnGrid}}}, Tuple{Base.OneTo{Int64}, Base.OneTo{Int64}, Base.OneTo{Int64}}}, BSpline{Cubic{Free{OnGrid}}}, Tuple{StepRangeLen{Float64, Base.TwicePrecision{Float64}, Base.TwicePrecision{Float64}}, StepRangeLen{Float64, Base.TwicePrecision{Float64}, Base.TwicePrecision{Float64}}, StepRangeLen{Float64, Base.TwicePrecision{Float64}, Base.TwicePrecision{Float64}}}}, BSpline{Cubic{Free{OnGrid}}}, Tuple{Periodic{Nothing}, Periodic{Nothing}, Flat{Nothing}}}, Interpolations.Extrapolation{Float64, 3, ScaledInterpolation{Float64, 3, Interpolations.BSplineInterpolation{Float64, 3, OffsetArrays.OffsetArray{Float64, 3, Array{Float64, 3}}, BSpline{Cubic{Free{OnGrid}}}, Tuple{Base.OneTo{Int64}, Base.OneTo{Int64}, Base.OneTo{Int64}}}, BSpline{Cubic{Free{OnGrid}}}, Tuple{StepRangeLen{Float64, Base.TwicePrecision{Float64}, Base.TwicePrecision{Float64}}, StepRangeLen{Float64, Base.TwicePrecision{Float64}, Base.TwicePrecision{Float64}}, StepRangeLen{Float64, Base.TwicePrecision{Float64}, Base.TwicePrecision{Float64}}}}, BSpline{Cubic{Free{OnGrid}}}, Tuple{Periodic{Nothing}, Periodic{Nothing}, Flat{Nothing}}}}, Float64, Float64, Float64, Float64, Vector{SVector{10, Float64}}, SciMLBase.ODESolution{Float64, 2, Vector{SVector{10, Float64}}, Nothing, Nothing, Vector{Float64}, Vector{Vector{SVector{10, Float64}}}, SciMLBase.ODEProblem{SVector{10, Float64}, Tuple{Float64, Float64}, false, Tuple{Interpolations.Extrapolation{Float64, 3, ScaledInterpolation{Float64, 3, Interpolations.BSplineInterpolation{Float64, 3, OffsetArrays.OffsetArray{Float64, 3, Array{Float64, 3}}, BSpline{Cubic{Free{OnGrid}}}, Tuple{Base.OneTo{Int64}, Base.OneTo{Int64}, Base.OneTo{Int64}}}, BSpline{Cubic{Free{OnGrid}}}, Tuple{StepRangeLen{Float64, Base.TwicePrecision{Float64}, Base.TwicePrecision{Float64}}, StepRangeLen{Float64, Base.TwicePrecision{Float64}, Base.TwicePrecision{Float64}}, StepRangeLen{Float64, Base.TwicePrecision{Float64}, Base.TwicePrecision{Float64}}}}, BSpline{Cubic{Free{OnGrid}}}, Tuple{Periodic{Nothing}, Periodic{Nothing}, Flat{Nothing}}}, Interpolations.Extrapolation{Float64, 3, ScaledInterpolation{Float64, 3, Interpolations.BSplineInterpolation{Float64, 3, OffsetArrays.OffsetArray{Float64, 3, Array{Float64, 3}}, BSpline{Cubic{Free{OnGrid}}}, Tuple{Base.OneTo{Int64}, Base.OneTo{Int64}, Base.OneTo{Int64}}}, BSpline{Cubic{Free{OnGrid}}}, Tuple{StepRangeLen{Float64, Base.TwicePrecision{Float64}, Base.TwicePrecision{Float64}}, StepRangeLen{Float64, Base.TwicePrecision{Float64}, Base.TwicePrecision{Float64}}, StepRangeLen{Float64, Base.TwicePrecision{Float64}, Base.TwicePrecision{Float64}}}}, BSpline{Cubic{Free{OnGrid}}}, Tuple{Periodic{Nothing}, Periodic{Nothing}, Flat{Nothing}}}}, SciMLBase.ODEFunction{false, CoherentStructures.var"#56#66"{SciMLBase.ODEFunction{false, typeof(my_interp), LinearAlgebra.UniformScaling{Bool}, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, typeof(SciMLBase.DEFAULT_OBSERVED), Nothing}}, LinearAlgebra.UniformScaling{Bool}, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, typeof(SciMLBase.DEFAULT_OBSERVED), Nothing}, Base.Iterators.Pairs{Union{}, Union{}, Tuple{}, NamedTuple{(), Tuple{}}}, SciMLBase.StandardODEProblem}, OrdinaryDiffEq.BS5, OrdinaryDiffEq.InterpolationData{SciMLBase.ODEFunction{false, CoherentStructures.var"#56#66"{SciMLBase.ODEFunction{false, typeof(my_interp), LinearAlgebra.UniformScaling{Bool}, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, typeof(SciMLBase.DEFAULT_OBSERVED), Nothing}}, LinearAlgebra.UniformScaling{Bool}, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, typeof(SciMLBase.DEFAULT_OBSERVED), Nothing}, Vector{SVector{10, Float64}}, Vector{Float64}, Vector{Vector{SVector{10, Float64}}}, OrdinaryDiffEq.BS5ConstantCache{Float64, Float64}}, DiffEqBase.DEStats}, SciMLBase.ODEFunction{false, CoherentStructures.var"#56#66"{SciMLBase.ODEFunction{false, typeof(my_interp), LinearAlgebra.UniformScaling{Bool}, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, typeof(SciMLBase.DEFAULT_OBSERVED), Nothing}}, LinearAlgebra.UniformScaling{Bool}, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, typeof(SciMLBase.DEFAULT_OBSERVED), Nothing}, OrdinaryDiffEq.BS5ConstantCache{Float64, Float64}, OrdinaryDiffEq.DEOptions{Float64, Float64, Float64, Float64, OrdinaryDiffEq.PIController{Rational{Int64}}, typeof(DiffEqBase.ODE_DEFAULT_NORM), typeof(LinearAlgebra.opnorm), Nothing, DiffEqBase.CallbackSet{Tuple{}, Tuple{}}, typeof(DiffEqBase.ODE_DEFAULT_ISOUTOFDOMAIN), typeof(DiffEqBase.ODE_DEFAULT_PROG_MESSAGE), typeof(DiffEqBase.ODE_DEFAULT_UNSTABLE_CHECK), DataStructures.BinaryHeap{Float64, DataStructures.FasterForward}, DataStructures.BinaryHeap{Float64, DataStructures.FasterForward}, Nothing, Nothing, Int64, Tuple{}, Vector{Float64}, Tuple{}}, SVector{10, Float64}, Float64, Nothing, OrdinaryDiffEq.DefaultInit}, cache::OrdinaryDiffEq.BS5ConstantCache{Float64, Float64})
@ OrdinaryDiffEq C:\Users\51342\.julia\packages\OrdinaryDiffEq\PIjOZ\src\perform_step\low_order_rk_perform_step.jl:370
[9] __init(prob::SciMLBase.ODEProblem{SVector{10, Float64}, Tuple{Float64, Float64}, false, Tuple{Interpolations.Extrapolation{Float64, 3, ScaledInterpolation{Float64, 3, Interpolations.BSplineInterpolation{Float64, 3, OffsetArrays.OffsetArray{Float64, 3, Array{Float64, 3}}, BSpline{Cubic{Free{OnGrid}}}, Tuple{Base.OneTo{Int64}, Base.OneTo{Int64}, Base.OneTo{Int64}}}, BSpline{Cubic{Free{OnGrid}}}, Tuple{StepRangeLen{Float64, Base.TwicePrecision{Float64}, Base.TwicePrecision{Float64}}, StepRangeLen{Float64, Base.TwicePrecision{Float64}, Base.TwicePrecision{Float64}}, StepRangeLen{Float64, Base.TwicePrecision{Float64}, Base.TwicePrecision{Float64}}}}, BSpline{Cubic{Free{OnGrid}}}, Tuple{Periodic{Nothing}, Periodic{Nothing}, Flat{Nothing}}}, Interpolations.Extrapolation{Float64, 3, ScaledInterpolation{Float64, 3, Interpolations.BSplineInterpolation{Float64, 3, OffsetArrays.OffsetArray{Float64, 3, Array{Float64, 3}}, BSpline{Cubic{Free{OnGrid}}}, Tuple{Base.OneTo{Int64}, Base.OneTo{Int64}, Base.OneTo{Int64}}}, BSpline{Cubic{Free{OnGrid}}}, Tuple{StepRangeLen{Float64, Base.TwicePrecision{Float64}, Base.TwicePrecision{Float64}}, StepRangeLen{Float64, Base.TwicePrecision{Float64}, Base.TwicePrecision{Float64}}, StepRangeLen{Float64, Base.TwicePrecision{Float64}, Base.TwicePrecision{Float64}}}}, BSpline{Cubic{Free{OnGrid}}}, Tuple{Periodic{Nothing}, Periodic{Nothing}, Flat{Nothing}}}}, SciMLBase.ODEFunction{false, CoherentStructures.var"#56#66"{SciMLBase.ODEFunction{false, typeof(my_interp), LinearAlgebra.UniformScaling{Bool}, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, typeof(SciMLBase.DEFAULT_OBSERVED), Nothing}}, LinearAlgebra.UniformScaling{Bool}, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, typeof(SciMLBase.DEFAULT_OBSERVED), Nothing}, Base.Iterators.Pairs{Union{}, Union{}, Tuple{}, NamedTuple{(), Tuple{}}}, SciMLBase.StandardODEProblem}, alg::OrdinaryDiffEq.BS5, timeseries_init::Tuple{}, ts_init::Tuple{}, ks_init::Tuple{}, recompile::Type{Val{true}}; saveat::Vector{Float64}, tstops::Tuple{}, d_discontinuities::Tuple{}, save_idxs::Nothing, save_everystep::Bool, save_on::Bool, save_start::Bool, save_end::Nothing, callback::Nothing, dense::Bool, calck::Bool, dt::Float64, dtmin::Nothing, dtmax::Float64, force_dtmin::Bool, adaptive::Bool, gamma::Rational{Int64}, abstol::Float64, reltol::Float64, qmin::Rational{Int64}, qmax::Int64, qsteady_min::Int64, qsteady_max::Int64, beta1::Nothing, beta2::Nothing, qoldinit::Rational{Int64}, controller::Nothing, fullnormalize::Bool, failfactor::Int64, maxiters::Int64, internalnorm::typeof(DiffEqBase.ODE_DEFAULT_NORM), internalopnorm::typeof(LinearAlgebra.opnorm), isoutofdomain::typeof(DiffEqBase.ODE_DEFAULT_ISOUTOFDOMAIN), unstable_check::typeof(DiffEqBase.ODE_DEFAULT_UNSTABLE_CHECK), verbose::Bool, timeseries_errors::Bool, dense_errors::Bool, advance_to_tstop::Bool, stop_at_next_tstop::Bool, initialize_save::Bool, progress::Bool, progress_steps::Int64, progress_name::String, progress_message::typeof(DiffEqBase.ODE_DEFAULT_PROG_MESSAGE), userdata::Nothing, allow_extrapolation::Bool, initialize_integrator::Bool, alias_u0::Bool, alias_du0::Bool, initializealg::OrdinaryDiffEq.DefaultInit, kwargs::Base.Iterators.Pairs{Union{}, Union{}, Tuple{}, NamedTuple{(), Tuple{}}})
@ OrdinaryDiffEq C:\Users\51342\.julia\packages\OrdinaryDiffEq\PIjOZ\src\solve.jl:456
[10] #__solve#465
@ C:\Users\51342\.julia\packages\OrdinaryDiffEq\PIjOZ\src\solve.jl:4 [inlined]
[11] #solve_call#56
@ C:\Users\51342\.julia\packages\DiffEqBase\ge0vq\src\solve.jl:61 [inlined]
[12] solve_up(prob::SciMLBase.ODEProblem{SVector{10, Float64}, Tuple{Float64, Float64}, false, Tuple{Interpolations.Extrapolation{Float64, 3, ScaledInterpolation{Float64, 3, Interpolations.BSplineInterpolation{Float64, 3, OffsetArrays.OffsetArray{Float64, 3, Array{Float64, 3}}, BSpline{Cubic{Free{OnGrid}}}, Tuple{Base.OneTo{Int64}, Base.OneTo{Int64}, Base.OneTo{Int64}}}, BSpline{Cubic{Free{OnGrid}}}, Tuple{StepRangeLen{Float64, Base.TwicePrecision{Float64}, Base.TwicePrecision{Float64}}, StepRangeLen{Float64, Base.TwicePrecision{Float64}, Base.TwicePrecision{Float64}}, StepRangeLen{Float64, Base.TwicePrecision{Float64}, Base.TwicePrecision{Float64}}}}, BSpline{Cubic{Free{OnGrid}}}, Tuple{Periodic{Nothing}, Periodic{Nothing}, Flat{Nothing}}}, Interpolations.Extrapolation{Float64, 3, ScaledInterpolation{Float64, 3, Interpolations.BSplineInterpolation{Float64, 3, OffsetArrays.OffsetArray{Float64, 3, Array{Float64, 3}}, BSpline{Cubic{Free{OnGrid}}}, Tuple{Base.OneTo{Int64}, Base.OneTo{Int64}, Base.OneTo{Int64}}}, BSpline{Cubic{Free{OnGrid}}}, Tuple{StepRangeLen{Float64, Base.TwicePrecision{Float64}, Base.TwicePrecision{Float64}}, StepRangeLen{Float64, Base.TwicePrecision{Float64}, Base.TwicePrecision{Float64}}, StepRangeLen{Float64, Base.TwicePrecision{Float64}, Base.TwicePrecision{Float64}}}}, BSpline{Cubic{Free{OnGrid}}}, Tuple{Periodic{Nothing}, Periodic{Nothing}, Flat{Nothing}}}}, SciMLBase.ODEFunction{false, CoherentStructures.var"#56#66"{SciMLBase.ODEFunction{false, typeof(my_interp), LinearAlgebra.UniformScaling{Bool}, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, typeof(SciMLBase.DEFAULT_OBSERVED), Nothing}}, LinearAlgebra.UniformScaling{Bool}, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, typeof(SciMLBase.DEFAULT_OBSERVED), Nothing}, Base.Iterators.Pairs{Union{}, Union{}, Tuple{}, NamedTuple{(), Tuple{}}}, SciMLBase.StandardODEProblem}, sensealg::Nothing, u0::SVector{10, Float64}, p::Tuple{Interpolations.Extrapolation{Float64, 3, ScaledInterpolation{Float64, 3, Interpolations.BSplineInterpolation{Float64, 3, OffsetArrays.OffsetArray{Float64, 3, Array{Float64, 3}}, BSpline{Cubic{Free{OnGrid}}}, Tuple{Base.OneTo{Int64}, Base.OneTo{Int64}, Base.OneTo{Int64}}}, BSpline{Cubic{Free{OnGrid}}}, Tuple{StepRangeLen{Float64, Base.TwicePrecision{Float64}, Base.TwicePrecision{Float64}}, StepRangeLen{Float64, Base.TwicePrecision{Float64}, Base.TwicePrecision{Float64}}, StepRangeLen{Float64, Base.TwicePrecision{Float64}, Base.TwicePrecision{Float64}}}}, BSpline{Cubic{Free{OnGrid}}}, Tuple{Periodic{Nothing}, Periodic{Nothing}, Flat{Nothing}}}, Interpolations.Extrapolation{Float64, 3, ScaledInterpolation{Float64, 3, Interpolations.BSplineInterpolation{Float64, 3, OffsetArrays.OffsetArray{Float64, 3, Array{Float64, 3}}, BSpline{Cubic{Free{OnGrid}}}, Tuple{Base.OneTo{Int64}, Base.OneTo{Int64}, Base.OneTo{Int64}}}, BSpline{Cubic{Free{OnGrid}}}, Tuple{StepRangeLen{Float64, Base.TwicePrecision{Float64}, Base.TwicePrecision{Float64}}, StepRangeLen{Float64, Base.TwicePrecision{Float64}, Base.TwicePrecision{Float64}}, StepRangeLen{Float64, Base.TwicePrecision{Float64}, Base.TwicePrecision{Float64}}}}, BSpline{Cubic{Free{OnGrid}}}, Tuple{Periodic{Nothing}, Periodic{Nothing}, Flat{Nothing}}}}, args::OrdinaryDiffEq.BS5; kwargs::Base.Iterators.Pairs{Symbol, Any, NTuple{6, Symbol}, NamedTuple{(:saveat, :save_everystep, :dense, :reltol, :abstol, :force_dtmin), Tuple{Vector{Float64}, Bool, Bool, Float64, Float64, Bool}}})
@ DiffEqBase C:\Users\51342\.julia\packages\DiffEqBase\ge0vq\src\solve.jl:85
[13] #solve#57
@ C:\Users\51342\.julia\packages\DiffEqBase\ge0vq\src\solve.jl:73 [inlined]
[14] #_flow#48
@ C:\Users\51342\.julia\packages\CoherentStructures\hfDt1\src\pullbacktensors.jl:69 [inlined]
[15] linearized_flow(odefun::SciMLBase.ODEFunction{false, typeof(my_interp), LinearAlgebra.UniformScaling{Bool}, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, typeof(SciMLBase.DEFAULT_OBSERVED), Nothing}, x::SVector{2, Float64}, tspan::Vector{Float64}, δ::Float64; tolerance::Float64, solver::OrdinaryDiffEq.BS5, p::Tuple{Interpolations.Extrapolation{Float64, 3, ScaledInterpolation{Float64, 3, Interpolations.BSplineInterpolation{Float64, 3, OffsetArrays.OffsetArray{Float64, 3, Array{Float64, 3}}, BSpline{Cubic{Free{OnGrid}}}, Tuple{Base.OneTo{Int64}, Base.OneTo{Int64}, Base.OneTo{Int64}}}, BSpline{Cubic{Free{OnGrid}}}, Tuple{StepRangeLen{Float64, Base.TwicePrecision{Float64}, Base.TwicePrecision{Float64}}, StepRangeLen{Float64, Base.TwicePrecision{Float64}, Base.TwicePrecision{Float64}}, StepRangeLen{Float64, Base.TwicePrecision{Float64}, Base.TwicePrecision{Float64}}}}, BSpline{Cubic{Free{OnGrid}}}, Tuple{Periodic{Nothing}, Periodic{Nothing}, Flat{Nothing}}}, Interpolations.Extrapolation{Float64, 3, ScaledInterpolation{Float64, 3, Interpolations.BSplineInterpolation{Float64, 3, OffsetArrays.OffsetArray{Float64, 3, Array{Float64, 3}}, BSpline{Cubic{Free{OnGrid}}}, Tuple{Base.OneTo{Int64}, Base.OneTo{Int64}, Base.OneTo{Int64}}}, BSpline{Cubic{Free{OnGrid}}}, Tuple{StepRangeLen{Float64, Base.TwicePrecision{Float64}, Base.TwicePrecision{Float64}}, StepRangeLen{Float64, Base.TwicePrecision{Float64}, Base.TwicePrecision{Float64}}, StepRangeLen{Float64, Base.TwicePrecision{Float64}, Base.TwicePrecision{Float64}}}}, BSpline{Cubic{Free{OnGrid}}}, Tuple{Periodic{Nothing}, Periodic{Nothing}, Flat{Nothing}}}})
@ CoherentStructures C:\Users\51342\.julia\packages\CoherentStructures\hfDt1\src\pullbacktensors.jl:146
[16] linearized_flow(odefun::Function, x::Vector{Float64}, tspan::Vector{Float64}, δ::Float64; kwargs::Base.Iterators.Pairs{Symbol, Any, Tuple{Symbol, Symbol, Symbol}, NamedTuple{(:tolerance, :p, :solver), Tuple{Float64, Tuple{Interpolations.Extrapolation{Float64, 3, ScaledInterpolation{Float64, 3, Interpolations.BSplineInterpolation{Float64, 3, OffsetArrays.OffsetArray{Float64, 3, Array{Float64, 3}}, BSpline{Cubic{Free{OnGrid}}}, Tuple{Base.OneTo{Int64}, Base.OneTo{Int64}, Base.OneTo{Int64}}}, BSpline{Cubic{Free{OnGrid}}}, Tuple{StepRangeLen{Float64, Base.TwicePrecision{Float64}, Base.TwicePrecision{Float64}}, StepRangeLen{Float64, Base.TwicePrecision{Float64}, Base.TwicePrecision{Float64}}, StepRangeLen{Float64, Base.TwicePrecision{Float64}, Base.TwicePrecision{Float64}}}}, BSpline{Cubic{Free{OnGrid}}}, Tuple{Periodic{Nothing}, Periodic{Nothing}, Flat{Nothing}}}, Interpolations.Extrapolation{Float64, 3, ScaledInterpolation{Float64, 3, Interpolations.BSplineInterpolation{Float64, 3, OffsetArrays.OffsetArray{Float64, 3, Array{Float64, 3}}, BSpline{Cubic{Free{OnGrid}}}, Tuple{Base.OneTo{Int64}, Base.OneTo{Int64}, Base.OneTo{Int64}}}, BSpline{Cubic{Free{OnGrid}}}, Tuple{StepRangeLen{Float64, Base.TwicePrecision{Float64}, Base.TwicePrecision{Float64}}, StepRangeLen{Float64, Base.TwicePrecision{Float64}, Base.TwicePrecision{Float64}}, StepRangeLen{Float64, Base.TwicePrecision{Float64}, Base.TwicePrecision{Float64}}}}, BSpline{Cubic{Free{OnGrid}}}, Tuple{Periodic{Nothing}, Periodic{Nothing}, Flat{Nothing}}}}, OrdinaryDiffEq.BS5}}})
@ CoherentStructures C:\Users\51342\.julia\packages\CoherentStructures\hfDt1\src\pullbacktensors.jl:93
[17] CG_tensor(odefun::Function, u::Vector{Float64}, tspan::Vector{Float64}, δ::Float64; kwargs::Base.Iterators.Pairs{Symbol, Any, Tuple{Symbol, Symbol, Symbol}, NamedTuple{(:tolerance, :p, :solver), Tuple{Float64, Tuple{Interpolations.Extrapolation{Float64, 3, ScaledInterpolation{Float64, 3, Interpolations.BSplineInterpolation{Float64, 3, OffsetArrays.OffsetArray{Float64, 3, Array{Float64, 3}}, BSpline{Cubic{Free{OnGrid}}}, Tuple{Base.OneTo{Int64}, Base.OneTo{Int64}, Base.OneTo{Int64}}}, BSpline{Cubic{Free{OnGrid}}}, Tuple{StepRangeLen{Float64, Base.TwicePrecision{Float64}, Base.TwicePrecision{Float64}}, StepRangeLen{Float64, Base.TwicePrecision{Float64}, Base.TwicePrecision{Float64}}, StepRangeLen{Float64, Base.TwicePrecision{Float64}, Base.TwicePrecision{Float64}}}}, BSpline{Cubic{Free{OnGrid}}}, Tuple{Periodic{Nothing}, Periodic{Nothing}, Flat{Nothing}}}, Interpolations.Extrapolation{Float64, 3, ScaledInterpolation{Float64, 3, Interpolations.BSplineInterpolation{Float64, 3, OffsetArrays.OffsetArray{Float64, 3, Array{Float64, 3}}, BSpline{Cubic{Free{OnGrid}}}, Tuple{Base.OneTo{Int64}, Base.OneTo{Int64}, Base.OneTo{Int64}}}, BSpline{Cubic{Free{OnGrid}}}, Tuple{StepRangeLen{Float64, Base.TwicePrecision{Float64}, Base.TwicePrecision{Float64}}, StepRangeLen{Float64, Base.TwicePrecision{Float64}, Base.TwicePrecision{Float64}}, StepRangeLen{Float64, Base.TwicePrecision{Float64}, Base.TwicePrecision{Float64}}}}, BSpline{Cubic{Free{OnGrid}}}, Tuple{Periodic{Nothing}, Periodic{Nothing}, Flat{Nothing}}}}, OrdinaryDiffEq.BS5}}})
@ CoherentStructures C:\Users\51342\.julia\packages\CoherentStructures\hfDt1\src\pullbacktensors.jl:314
[18] macro expansion
@ C:\Users\51342\.julia\packages\CoherentStructures\hfDt1\src\plotting.jl:528 [inlined]
[19] (::CoherentStructures.var"#559#561"{SharedArrays.SharedMatrix{Float64}, Vector{Float64}, Float64, Int64, Vector{Float64}, Tuple{Interpolations.Extrapolation{Float64, 3, ScaledInterpolation{Float64, 3, Interpolations.BSplineInterpolation{Float64, 3, OffsetArrays.OffsetArray{Float64, 3, Array{Float64, 3}}, BSpline{Cubic{Free{OnGrid}}}, Tuple{Base.OneTo{Int64}, Base.OneTo{Int64}, Base.OneTo{Int64}}}, BSpline{Cubic{Free{OnGrid}}}, Tuple{StepRangeLen{Float64, Base.TwicePrecision{Float64}, Base.TwicePrecision{Float64}}, StepRangeLen{Float64, Base.TwicePrecision{Float64}, Base.TwicePrecision{Float64}}, StepRangeLen{Float64, Base.TwicePrecision{Float64}, Base.TwicePrecision{Float64}}}}, BSpline{Cubic{Free{OnGrid}}}, Tuple{Periodic{Nothing}, Periodic{Nothing}, Flat{Nothing}}}, Interpolations.Extrapolation{Float64, 3, ScaledInterpolation{Float64, 3, Interpolations.BSplineInterpolation{Float64, 3, OffsetArrays.OffsetArray{Float64, 3, Array{Float64, 3}}, BSpline{Cubic{Free{OnGrid}}}, Tuple{Base.OneTo{Int64}, Base.OneTo{Int64}, Base.OneTo{Int64}}}, BSpline{Cubic{Free{OnGrid}}}, Tuple{StepRangeLen{Float64, Base.TwicePrecision{Float64}, Base.TwicePrecision{Float64}}, StepRangeLen{Float64, Base.TwicePrecision{Float64}, Base.TwicePrecision{Float64}}, StepRangeLen{Float64, Base.TwicePrecision{Float64}, Base.TwicePrecision{Float64}}}}, BSpline{Cubic{Free{OnGrid}}}, Tuple{Periodic{Nothing}, Periodic{Nothing}, Flat{Nothing}}}}, typeof(my_interp), Bool, typeof(CoherentStructures.always_true), OrdinaryDiffEq.BS5, Float64})(reducer::CoherentStructures.var"#558#560", R::UnitRange{Int64}, lo::Int64, hi::Int64)
@ CoherentStructures C:\buildbot\worker\package_win64\build\usr\share\julia\stdlib\v1.6\Distributed\src\macros.jl:291
[20] #137
@ C:\buildbot\worker\package_win64\build\usr\share\julia\stdlib\v1.6\Distributed\src\remotecall.jl:354 [inlined]
[21] run_work_thunk(thunk::Distributed.var"#137#138"{CoherentStructures.var"#559#561"{SharedArrays.SharedMatrix{Float64}, Vector{Float64}, Float64, Int64, Vector{Float64}, Tuple{Interpolations.Extrapolation{Float64, 3, ScaledInterpolation{Float64, 3, Interpolations.BSplineInterpolation{Float64, 3, OffsetArrays.OffsetArray{Float64, 3, Array{Float64, 3}}, BSpline{Cubic{Free{OnGrid}}}, Tuple{Base.OneTo{Int64}, Base.OneTo{Int64}, Base.OneTo{Int64}}}, BSpline{Cubic{Free{OnGrid}}}, Tuple{StepRangeLen{Float64, Base.TwicePrecision{Float64}, Base.TwicePrecision{Float64}}, StepRangeLen{Float64, Base.TwicePrecision{Float64}, Base.TwicePrecision{Float64}}, StepRangeLen{Float64, Base.TwicePrecision{Float64}, Base.TwicePrecision{Float64}}}}, BSpline{Cubic{Free{OnGrid}}}, Tuple{Periodic{Nothing}, Periodic{Nothing}, Flat{Nothing}}}, Interpolations.Extrapolation{Float64, 3, ScaledInterpolation{Float64, 3, Interpolations.BSplineInterpolation{Float64, 3, OffsetArrays.OffsetArray{Float64, 3, Array{Float64, 3}}, BSpline{Cubic{Free{OnGrid}}}, Tuple{Base.OneTo{Int64}, Base.OneTo{Int64}, Base.OneTo{Int64}}}, BSpline{Cubic{Free{OnGrid}}}, Tuple{StepRangeLen{Float64, Base.TwicePrecision{Float64}, Base.TwicePrecision{Float64}}, StepRangeLen{Float64, Base.TwicePrecision{Float64}, Base.TwicePrecision{Float64}}, StepRangeLen{Float64, Base.TwicePrecision{Float64}, Base.TwicePrecision{Float64}}}}, BSpline{Cubic{Free{OnGrid}}}, Tuple{Periodic{Nothing}, Periodic{Nothing}, Flat{Nothing}}}}, typeof(my_interp), Bool, typeof(CoherentStructures.always_true), OrdinaryDiffEq.BS5, Float64}, Tuple{CoherentStructures.var"#558#560", UnitRange{Int64}, Int64, Int64}, Base.Iterators.Pairs{Union{}, Union{}, Tuple{}, NamedTuple{(), Tuple{}}}}, print_error::Bool)
@ Distributed C:\buildbot\worker\package_win64\build\usr\share\julia\stdlib\v1.6\Distributed\src\process_messages.jl:63
[22] remotecall_fetch(::Function, ::Distributed.LocalProcess, ::Function, ::Vararg{Any, N} where N; kwargs::Base.Iterators.Pairs{Union{}, Union{}, Tuple{}, NamedTuple{(), Tuple{}}})
@ Distributed C:\buildbot\worker\package_win64\build\usr\share\julia\stdlib\v1.6\Distributed\src\remotecall.jl:379
[23] remotecall_fetch(::Function, ::Distributed.LocalProcess, ::Function, ::Vararg{Any, N} where N)
@ Distributed C:\buildbot\worker\package_win64\build\usr\share\julia\stdlib\v1.6\Distributed\src\remotecall.jl:379
[24] remotecall_fetch(::Function, ::Int64, ::Function, ::Vararg{Any, N} where N; kwargs::Base.Iterators.Pairs{Union{}, Union{}, Tuple{}, NamedTuple{(), Tuple{}}})
@ Distributed C:\buildbot\worker\package_win64\build\usr\share\julia\stdlib\v1.6\Distributed\src\remotecall.jl:421
[25] remotecall_fetch
@ C:\buildbot\worker\package_win64\build\usr\share\julia\stdlib\v1.6\Distributed\src\remotecall.jl:421 [inlined]
[26] (::Distributed.var"#157#158"{CoherentStructures.var"#558#560", CoherentStructures.var"#559#561"{SharedArrays.SharedMatrix{Float64}, Vector{Float64}, Float64, Int64, Vector{Float64}, Tuple{Interpolations.Extrapolation{Float64, 3, ScaledInterpolation{Float64, 3, Interpolations.BSplineInterpolation{Float64, 3, OffsetArrays.OffsetArray{Float64, 3, Array{Float64, 3}}, BSpline{Cubic{Free{OnGrid}}}, Tuple{Base.OneTo{Int64}, Base.OneTo{Int64}, Base.OneTo{Int64}}}, BSpline{Cubic{Free{OnGrid}}}, Tuple{StepRangeLen{Float64, Base.TwicePrecision{Float64}, Base.TwicePrecision{Float64}}, StepRangeLen{Float64, Base.TwicePrecision{Float64}, Base.TwicePrecision{Float64}}, StepRangeLen{Float64, Base.TwicePrecision{Float64}, Base.TwicePrecision{Float64}}}}, BSpline{Cubic{Free{OnGrid}}}, Tuple{Periodic{Nothing}, Periodic{Nothing}, Flat{Nothing}}}, Interpolations.Extrapolation{Float64, 3, ScaledInterpolation{Float64, 3, Interpolations.BSplineInterpolation{Float64, 3, OffsetArrays.OffsetArray{Float64, 3, Array{Float64, 3}}, BSpline{Cubic{Free{OnGrid}}}, Tuple{Base.OneTo{Int64}, Base.OneTo{Int64}, Base.OneTo{Int64}}}, BSpline{Cubic{Free{OnGrid}}}, Tuple{StepRangeLen{Float64, Base.TwicePrecision{Float64}, Base.TwicePrecision{Float64}}, StepRangeLen{Float64, Base.TwicePrecision{Float64}, Base.TwicePrecision{Float64}}, StepRangeLen{Float64, Base.TwicePrecision{Float64}, Base.TwicePrecision{Float64}}}}, BSpline{Cubic{Free{OnGrid}}}, Tuple{Periodic{Nothing}, Periodic{Nothing}, Flat{Nothing}}}}, typeof(my_interp), Bool, typeof(CoherentStructures.always_true), OrdinaryDiffEq.BS5, Float64}, UnitRange{Int64}, Vector{UnitRange{Int64}}, Int64, Int64})()
@ Distributed C:\buildbot\worker\package_win64\build\usr\share\julia\stdlib\v1.6\Distributed\src\macros.jl:274
Stacktrace:
[1] remotecall_fetch(::Function, ::Distributed.LocalProcess, ::Function, ::Vararg{Any, N} where N; kwargs::Base.Iterators.Pairs{Union{}, Union{}, Tuple{}, NamedTuple{(), Tuple{}}})
@ Distributed C:\buildbot\worker\package_win64\build\usr\share\julia\stdlib\v1.6\Distributed\src\remotecall.jl:380
[2] remotecall_fetch(::Function, ::Distributed.LocalProcess, ::Function, ::Vararg{Any, N} where N)
@ Distributed C:\buildbot\worker\package_win64\build\usr\share\julia\stdlib\v1.6\Distributed\src\remotecall.jl:379
[3] remotecall_fetch(::Function, ::Int64, ::Function, ::Vararg{Any, N} where N; kwargs::Base.Iterators.Pairs{Union{}, Union{}, Tuple{}, NamedTuple{(), Tuple{}}})
@ Distributed C:\buildbot\worker\package_win64\build\usr\share\julia\stdlib\v1.6\Distributed\src\remotecall.jl:421
[4] remotecall_fetch
@ C:\buildbot\worker\package_win64\build\usr\share\julia\stdlib\v1.6\Distributed\src\remotecall.jl:421 [inlined]
[5] (::Distributed.var"#157#158"{CoherentStructures.var"#558#560", CoherentStructures.var"#559#561"{SharedArrays.SharedMatrix{Float64}, Vector{Float64}, Float64, Int64, Vector{Float64}, Tuple{Interpolations.Extrapolation{Float64, 3, ScaledInterpolation{Float64, 3, Interpolations.BSplineInterpolation{Float64, 3, OffsetArrays.OffsetArray{Float64, 3, Array{Float64, 3}}, BSpline{Cubic{Free{OnGrid}}}, Tuple{Base.OneTo{Int64}, Base.OneTo{Int64}, Base.OneTo{Int64}}}, BSpline{Cubic{Free{OnGrid}}}, Tuple{StepRangeLen{Float64, Base.TwicePrecision{Float64}, Base.TwicePrecision{Float64}}, StepRangeLen{Float64, Base.TwicePrecision{Float64}, Base.TwicePrecision{Float64}}, StepRangeLen{Float64, Base.TwicePrecision{Float64}, Base.TwicePrecision{Float64}}}}, BSpline{Cubic{Free{OnGrid}}}, Tuple{Periodic{Nothing}, Periodic{Nothing}, Flat{Nothing}}}, Interpolations.Extrapolation{Float64, 3, ScaledInterpolation{Float64, 3, Interpolations.BSplineInterpolation{Float64, 3, OffsetArrays.OffsetArray{Float64, 3, Array{Float64, 3}}, BSpline{Cubic{Free{OnGrid}}}, Tuple{Base.OneTo{Int64}, Base.OneTo{Int64}, Base.OneTo{Int64}}}, BSpline{Cubic{Free{OnGrid}}}, Tuple{StepRangeLen{Float64, Base.TwicePrecision{Float64}, Base.TwicePrecision{Float64}}, StepRangeLen{Float64, Base.TwicePrecision{Float64}, Base.TwicePrecision{Float64}}, StepRangeLen{Float64, Base.TwicePrecision{Float64}, Base.TwicePrecision{Float64}}}}, BSpline{Cubic{Free{OnGrid}}}, Tuple{Periodic{Nothing}, Periodic{Nothing}, Flat{Nothing}}}}, typeof(my_interp), Bool, typeof(CoherentStructures.always_true), OrdinaryDiffEq.BS5, Float64}, UnitRange{Int64}, Vector{UnitRange{Int64}}, Int64, Int64})()
@ Distributed C:\buildbot\worker\package_win64\build\usr\share\julia\stdlib\v1.6\Distributed\src\macros.jl:274
Stacktrace:
[1] wait
@ .\task.jl:322 [inlined]
[2] fetch
@ .\task.jl:337 [inlined]
[3] preduce(reducer::Function, f::Function, R::UnitRange{Int64})
@ Distributed C:\buildbot\worker\package_win64\build\usr\share\julia\stdlib\v1.6\Distributed\src\macros.jl:278
[4] macro expansion
@ .\task.jl:387 [inlined]
[5] macro expansion
@ C:\Users\51342.julia\packages\CoherentStructures\hfDt1\src\plotting.jl:521 [inlined]
[6] apply_recipe(plotattributes::AbstractDict{Symbol, Any}, as::CoherentStructures.Plot_FTLE)
@ CoherentStructures C:\Users\51342.julia\packages\RecipesBase\3fzVq\src\RecipesBase.jl:283
[7] _process_userrecipes!(plt::Any, plotattributes::Any, args::Any)
@ RecipesPipeline C:\Users\51342.julia\packages\RecipesPipeline\Bxu2O\src\user_recipe.jl:36
[8] recipe_pipeline!(plt::Any, plotattributes::Any, args::Any)
@ RecipesPipeline C:\Users\51342.julia\packages\RecipesPipeline\Bxu2O\src\RecipesPipeline.jl:70
[9] _plot!(plt::Plots.Plot, plotattributes::Any, args::Any)
@ Plots C:\Users\51342.julia\packages\Plots\5kcBO\src\plot.jl:208
[10] #plot#157
@ C:\Users\51342.julia\packages\Plots\5kcBO\src\plot.jl:91 [inlined]
[11] plot_ftle(::Function, ::Vararg{Any, N} where N; kw::Base.Iterators.Pairs{Symbol, Any, NTuple{4, Symbol}, NamedTuple{(:tolerance, :aspect_ratio, :color, :pass_on_errors), Tuple{Float64, Int64, Symbol, Bool}}})
@ CoherentStructures C:\Users\51342.julia\packages\RecipesBase\3fzVq\src\RecipesBase.jl:358
[12] top-level scope
@ .\timing.jl:210 [inlined]
[13] top-level scope
@ .\In[2]:0
[14] eval
@ .\boot.jl:360 [inlined]
[15] include_string(mapexpr::typeof(REPL.softscope), mod::Module, code::String, filename::String)
@ Base .\loading.jl:1116
DISPLAY_PLOT not defined
When I test the examples "bickley.jl", an error occurs as follows: "ERROR: LoadError: UndefVarError: DISPLAY_PLOT not defined", Is other packages needed?
Cannot install due to package version incompatibility
JuAFEM might be necessary in [compat]?
(v1.3) pkg> add https://github.com/CoherentStructures/CoherentStructures.jl.git
Updating registry at `C:\Users\islent\.julia\registries\General`
Updating git-repo `https://github.com/JuliaRegistries/General.git`
Cloning git-repo `https://github.com/CoherentStructures/CoherentStructures.jl.git`
Updating git-repo `https://github.com/CoherentStructures/CoherentStructures.jl.git`
Updating git-repo `https://github.com/CoherentStructures/CoherentStructures.jl.git`
Resolving package versions...
ERROR: Unsatisfiable requirements detected for package JuAFEM [30d91d44]:
JuAFEM [30d91d44] log:
├─JuAFEM [30d91d44] has no known versions!
└─restricted to versions * by CoherentStructures [0c1513b4] — no versions left
└─CoherentStructures [0c1513b4] log:
├─possible versions are: 0.0.0 or uninstalled
└─CoherentStructures [0c1513b4] is fixed to version 0.0.0Improve documentation
We have some undocumented, but effective functionality, that we should expose to our growing user community. ;-) This issue is to track progress, so better tick off items rather than close this issue until we are completely happy with our documentation. Feel free to add items as they occur:
- advection-diffusion equation in Lagrangian coordinates/time-dependent diffusion equation solvers
- plotting functionality
Ocean flow example in documentation shows only 3 vortices
Towards v1.0
I think the package is in pretty good shape by now, so we could think about releasing a version 1.0. Before that, we should get the following things done:
- finish #12
- consciously handle the orientation issue in the closed orbit computation in
constrainedLCS(should be easy) - include linear ODE integrators and showcase solving the advection-diffusion equation in Lagrangian coordinates (that's been done by Philipp Rohrmüller in #16)
- maybe tests of the FEM part?
- would be good, perhaps, to also turn on coverage testing (well, that should be done after we have some FEM tests, otherwise coverage will just look disastrous
😄 ).
Anything else to include?
Installation error due to LinearImplicitEuler
I've tried to install the package via the package manager, but it gives me the following error:
1 dependency errored. To see a full report either run import Pkg; Pkg.precompile() or load the package
When I run Pkg.precompile(), it says there is invalid subtyping in definition of LinearImplicitEuler:
I attach the full error message, as I'm not sure what might be relevant. This was on a freshly installed Julia-1.6.2, on MacOS 11.4.
Error compiling in v1.7.0
Hello!
julia> VERSION v"1.7.0"
julia> Sys.iswindows() true
julia> import Pkg; Pkg.precompile()
Precompiling project...
✗ CoherentStructures
0 dependencies successfully precompiled in 16 seconds (249 already precompiled)
ERROR: The following 1 direct dependency failed to precompile:
CoherentStructures [0c1513b4-3a13-56f1-9cd2-8312eaec88c4]
Failed to precompile CoherentStructures [0c1513b4-3a13-56f1-9cd2-8312eaec88c4] to C:\Users\Rodrigo.julia\compiled\v1.7\CoherentStructures\jl_8DFF.tmp.
ERROR: LoadError: invalid subtyping in definition of LinearImplicitEuler
Stacktrace:
[1] top-level scope
@ C:\Users\Rodrigo.julia\packages\CoherentStructures\CxxKh\src\odesolvers.jl:7
[2] include(mod::Module, _path::String)
@ Base .\Base.jl:418
[3] include(x::String)
@ CoherentStructures C:\Users\Rodrigo.julia\packages\CoherentStructures\CxxKh\src\CoherentStructures.jl:1
[4] top-level scope
@ C:\Users\Rodrigo.julia\packages\CoherentStructures\CxxKh\src\CoherentStructures.jl:104
[5] include
@ .\Base.jl:418 [inlined]
[6] include_package_for_output(pkg::Base.PkgId, input::String, depot_path::Vector{String}, dl_load_path::Vector{String}, load_path::Vector{String}, concrete_deps::Vector{Pair{Base.PkgId, UInt64}}, source::Nothing)
@ Base .\loading.jl:1318
[7] top-level scope
@ none:1
[8] eval
@ .\boot.jl:373 [inlined]
[9] eval(x::Expr)
@ Base.MainInclude .\client.jl:453
[10] top-level scope
@ none:1
in expression starting at C:\Users\Rodrigo.julia\packages\CoherentStructures\CxxKh\src\odesolvers.jl:7
in expression starting at C:\Users\Rodrigo.julia\packages\CoherentStructures\CxxKh\src\CoherentStructures.jl:1
Stacktrace:
[1] pkgerror(msg::String)
@ Pkg.Types C:\Users\Rodrigo\AppData\Local\Programs\Julia-1.7.0\share\julia\stdlib\v1.7\Pkg\src\Types.jl:68
[2] precompile(ctx::Pkg.Types.Context; internal_call::Bool, strict::Bool, warn_loaded::Bool, kwargs::Base.Pairs{Symbol, Union{}, Tuple{}, NamedTuple{(), Tuple{}}})
@ Pkg.API C:\Users\Rodrigo\AppData\Local\Programs\Julia-1.7.0\share\julia\stdlib\v1.7\Pkg\src\API.jl:1362
[3] precompile
@ C:\Users\Rodrigo\AppData\Local\Programs\Julia-1.7.0\share\julia\stdlib\v1.7\Pkg\src\API.jl:1013 [inlined]
[4] #precompile#220
@ C:\Users\Rodrigo\AppData\Local\Programs\Julia-1.7.0\share\julia\stdlib\v1.7\Pkg\src\API.jl:1011 [inlined]
[5] precompile()
@ Pkg.API C:\Users\Rodrigo\AppData\Local\Programs\Julia-1.7.0\share\julia\stdlib\v1.7\Pkg\src\API.jl:1011
[6] top-level scope
No elliptic barrier found
Hi everyone,
I have used the CoherentStructure code in the past to find Lagrangian coherent contours from an altimetry-derived velocity field, and it worked fine. However, I updated all libraries last week; the code ran fine, but no elliptic barrier was identified for the same dataset.
I then tried to run the ocean_flow.jl (example), and the problem persists: no elliptic barrier is identified. Please, see below the results for both my analysis and the ocean_flow.jl example.
I am using Julia v1.8.5 and CoherentStructures v0.4.11.
Do you have insights on what could be happening?
Thank you!
Some trouble with DISPLAY_PLOT
Hi guys,
thanks for providing the toolbox! I have been trying to run the examples and I guess I have included all the necessary packages. Nevertheless I am always getting an error whenever the function DISPLAY_PLOT is called. In the case of the rot_double_gyre.jl example I have the following error:
ERROR: LoadError: UndefVarError: DISPLAY_PLOT not defined
Stacktrace:
[1] top-level scope at none:0
[2] include at ./boot.jl:317 [inlined]
[3] include_relative(::Module, ::String) at ./loading.jl:1041
[4] include(::Module, ::String) at ./sysimg.jl:29
[5] include(::String) at ./client.jl:388
[6] top-level scope at none:0
I already added PyPlot and Plots. Is there any Julia plot package which I am missing?
Compiling problem
Hello,
I just did ]update and I get a 1 dependency error, CoherentStructures won't precompile:
julia> using CoherentStructures
[ Info: Precompiling CoherentStructures [0c1513b4-3a13-56f1-9cd2-8312eaec88c4]
ERROR: LoadError: invalid subtyping in definition of LinearImplicitEuler
Stacktrace:
[1] top-level scope
@ C:\Users\Rodrigo.julia\packages\CoherentStructures\NPGM6\src\odesolvers.jl:7
[2] include(mod::Module, _path::String)
@ Base .\Base.jl:418
[3] include(x::String)
@ CoherentStructures C:\Users\Rodrigo.julia\packages\CoherentStructures\NPGM6\src\CoherentStructures.jl:1
[4] top-level scope
@ C:\Users\Rodrigo.julia\packages\CoherentStructures\NPGM6\src\CoherentStructures.jl:104
[5] include
@ .\Base.jl:418 [inlined]
[6] include_package_for_output(pkg::Base.PkgId, input::String, depot_path::Vector{String}, dl_load_path::Vector{String}, load_path::Vector{String}, concrete_deps::Vector{Pair{Base.PkgId, UInt64}}, source::Nothing)
@ Base .\loading.jl:1318
[7] top-level scope
@ none:1
[8] eval
@ .\boot.jl:373 [inlined]
[9] eval(x::Expr)
@ Base.MainInclude .\client.jl:453
[10] top-level scope
@ none:1
in expression starting at C:\Users\Rodrigo.julia\packages\CoherentStructures\NPGM6\src\odesolvers.jl:7
in expression starting at C:\Users\Rodrigo.julia\packages\CoherentStructures\NPGM6\src\CoherentStructures.jl:1
ERROR: Failed to precompile CoherentStructures [0c1513b4-3a13-56f1-9cd2-8312eaec88c4] to C:\Users\Rodrigo.julia\compiled\v1.7\CoherentStructures\jl_2785.tmp.
Stacktrace:
[1] error(s::String)
@ Base .\error.jl:33
[2] compilecache(pkg::Base.PkgId, path::String, internal_stderr::IO, internal_stdout::IO, ignore_loaded_modules::Bool)
@ Base .\loading.jl:1466
[3] compilecache(pkg::Base.PkgId, path::String)
@ Base .\loading.jl:1410
[4] _require(pkg::Base.PkgId)
@ Base .\loading.jl:1120
[5] require(uuidkey::Base.PkgId)
@ Base .\loading.jl:1013
[6] require(into::Module, mod::Symbol)
Recommend Projects
-
React
A declarative, efficient, and flexible JavaScript library for building user interfaces.
-
Vue.js
🖖 Vue.js is a progressive, incrementally-adoptable JavaScript framework for building UI on the web.
-
Typescript
TypeScript is a superset of JavaScript that compiles to clean JavaScript output.
-
TensorFlow
An Open Source Machine Learning Framework for Everyone
-
Django
The Web framework for perfectionists with deadlines.
-
Laravel
A PHP framework for web artisans
-
D3
Bring data to life with SVG, Canvas and HTML. 📊📈🎉
-
Recommend Topics
-
javascript
JavaScript (JS) is a lightweight interpreted programming language with first-class functions.
-
web
Some thing interesting about web. New door for the world.
-
server
A server is a program made to process requests and deliver data to clients.
-
Machine learning
Machine learning is a way of modeling and interpreting data that allows a piece of software to respond intelligently.
-
Visualization
Some thing interesting about visualization, use data art
-
Game
Some thing interesting about game, make everyone happy.
Recommend Org
-
Facebook
We are working to build community through open source technology. NB: members must have two-factor auth.
-
Microsoft
Open source projects and samples from Microsoft.
-
Google
Google ❤️ Open Source for everyone.
-
Alibaba
Alibaba Open Source for everyone
-
D3
Data-Driven Documents codes.
-
Tencent
China tencent open source team.