Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Error testing GXBeam #68

Closed
Kevin-Mattheus-Moerman opened this issue Feb 4, 2022 · 3 comments
Closed

Error testing GXBeam #68

Kevin-Mattheus-Moerman opened this issue Feb 4, 2022 · 3 comments

Comments

@Kevin-Mattheus-Moerman
Copy link

I ran ] test GXBeam and got the following error message:

┌ Warning: dt <= dtmin. Aborting. There is either an error in your model specification or the true solution is unstable.
└ @ SciMLBase ~/.julia/packages/SciMLBase/L7Nun/src/integrator_interface.jl:345
DifferentialEquations: Test Failed at /home/kevin/.julia/packages/GXBeam/S8Yc6/test/runtests.jl:1200
  Expression: sol.t[end] == 2.0
   Evaluated: 0.4756178300699102 == 2.0
Stacktrace:
 [1] macro expansion
   @ ~/.julia/packages/GXBeam/S8Yc6/test/runtests.jl:1200 [inlined]
 [2] macro expansion
   @ /buildworker/worker/package_linux64/build/usr/share/julia/stdlib/v1.6/Test/src/Test.jl:1151 [inlined]
 [3] top-level scope
   @ ~/.julia/packages/GXBeam/S8Yc6/test/runtests.jl:1141
DifferentialEquations: Error During Test at /home/kevin/.julia/packages/GXBeam/S8Yc6/test/runtests.jl:1139
  Got exception outside of a @test
  type Nothing has no field α
  Stacktrace:
    [1] setproperty!(x::Nothing, f::Symbol, v::Float64)
      @ Base ./Base.jl:34
    [2] calc_W!
      @ ~/.julia/packages/OrdinaryDiffEq/Op0Oq/src/derivative_utils.jl:549 [inlined]
    [3] update_W!
      @ ~/.julia/packages/OrdinaryDiffEq/Op0Oq/src/derivative_utils.jl:653 [inlined]
    [4] nlsolve!(nlsolver::OrdinaryDiffEq.NLSolver{NLNewton{Rational{Int64}, Rational{Int64}, Rational{Int64}}, true, Vector{Float64}, Rational{Int64}, Nothing, Float64, OrdinaryDiffEq.NLNewtonCache{Vector{Float64}, Float64, Float64, Vector{Float64}, Matrix{Float64}, Matrix{Float64}, Nothing, Nothing, LinearSolve.LinearCache{Matrix{Float64}, Vector{Float64}, Vector{Float64}, SciMLBase.NullParameters, Nothing, LU{Float64, Matrix{Float64}}, LinearSolve.InvPreconditioner{Diagonal{Float64, Vector{Float64}}}, Diagonal{Float64, Vector{Float64}}, Float64}}}, integrator::OrdinaryDiffEq.ODEIntegrator{DABDF2{12, true, Nothing, NLNewton{Rational{Int64}, Rational{Int64}, Rational{Int64}}, typeof(OrdinaryDiffEq.DEFAULT_PRECS), Val{:forward}, true, nothing}, true, Vector{Float64}, Vector{Float64}, Float64, Tuple{var"#103#105"{Int64}, Dict{Int64, DistributedLoads{Float64}}, Dict{Int64, Vector{PointMass{Float64}}}, StaticArrays.SVector{3, Float64}, StaticArrays.SVector{3, Float64}, StaticArrays.SVector{3, Float64}, StaticArrays.SVector{3, Float64}, StaticArrays.SVector{3, Float64}, StaticArrays.SVector{3, Float64}}, Float64, Float64, Float64, Float64, Vector{Vector{Float64}}, DAESolution{Float64, 2, Vector{Vector{Float64}}, Nothing, Nothing, Nothing, Vector{Float64}, DAEProblem{Vector{Float64}, Vector{Float64}, Tuple{Float64, Float64}, true, Tuple{var"#103#105"{Int64}, Dict{Int64, DistributedLoads{Float64}}, Dict{Int64, Vector{PointMass{Float64}}}, StaticArrays.SVector{3, Float64}, StaticArrays.SVector{3, Float64}, StaticArrays.SVector{3, Float64}, StaticArrays.SVector{3, Float64}, StaticArrays.SVector{3, Float64}, StaticArrays.SVector{3, Float64}}, DAEFunction{true, GXBeam.var"#299#301"{Assembly{Float64, Vector{Vector{Float64}}, UnitRange{Int64}, Vector{GXBeam.Element{Float64}}}, Float64, Vector{Int64}, Vector{Int64}, Vector{Int64}, Vector{Int64}, Vector{Int64}, Vector{Int64}}, Nothing, Nothing, GXBeam.var"#300#302"{Assembly{Float64, Vector{Vector{Float64}}, UnitRange{Int64}, Vector{GXBeam.Element{Float64}}}, Float64, Vector{Int64}, Vector{Int64}, Vector{Int64}, Vector{Int64}, Vector{Int64}, Vector{Int64}}, Nothing, Nothing, Nothing, SparseArrays.SparseMatrixCSC{Bool, Int64}, Nothing, Nothing, Nothing, Nothing, Nothing}, Base.Iterators.Pairs{Union{}, Union{}, Tuple{}, NamedTuple{(), Tuple{}}}, Vector{Bool}}, DABDF2{12, true, Nothing, NLNewton{Rational{Int64}, Rational{Int64}, Rational{Int64}}, typeof(OrdinaryDiffEq.DEFAULT_PRECS), Val{:forward}, true, nothing}, OrdinaryDiffEq.InterpolationData{DAEFunction{true, GXBeam.var"#299#301"{Assembly{Float64, Vector{Vector{Float64}}, UnitRange{Int64}, Vector{GXBeam.Element{Float64}}}, Float64, Vector{Int64}, Vector{Int64}, Vector{Int64}, Vector{Int64}, Vector{Int64}, Vector{Int64}}, Nothing, Nothing, GXBeam.var"#300#302"{Assembly{Float64, Vector{Vector{Float64}}, UnitRange{Int64}, Vector{GXBeam.Element{Float64}}}, Float64, Vector{Int64}, Vector{Int64}, Vector{Int64}, Vector{Int64}, Vector{Int64}, Vector{Int64}}, Nothing, Nothing, Nothing, SparseArrays.SparseMatrixCSC{Bool, Int64}, Nothing, Nothing, Nothing, Nothing, Nothing}, Vector{Vector{Float64}}, Vector{Float64}, Vector{Vector{Vector{Float64}}}, OrdinaryDiffEq.DABDF2Cache{Vector{Float64}, Vector{Float64}, Vector{Float64}, OrdinaryDiffEq.NLSolver{NLNewton{Rational{Int64}, Rational{Int64}, Rational{Int64}}, true, Vector{Float64}, Rational{Int64}, Nothing, Float64, OrdinaryDiffEq.NLNewtonCache{Vector{Float64}, Float64, Float64, Vector{Float64}, Matrix{Float64}, Matrix{Float64}, Nothing, Nothing, LinearSolve.LinearCache{Matrix{Float64}, Vector{Float64}, Vector{Float64}, SciMLBase.NullParameters, Nothing, LU{Float64, Matrix{Float64}}, LinearSolve.InvPreconditioner{Diagonal{Float64, Vector{Float64}}}, Diagonal{Float64, Vector{Float64}}, Float64}}}, Float64}}, DiffEqBase.DEStats}, DAEFunction{true, GXBeam.var"#299#301"{Assembly{Float64, Vector{Vector{Float64}}, UnitRange{Int64}, Vector{GXBeam.Element{Float64}}}, Float64, Vector{Int64}, Vector{Int64}, Vector{Int64}, Vector{Int64}, Vector{Int64}, Vector{Int64}}, Nothing, Nothing, GXBeam.var"#300#302"{Assembly{Float64, Vector{Vector{Float64}}, UnitRange{Int64}, Vector{GXBeam.Element{Float64}}}, Float64, Vector{Int64}, Vector{Int64}, Vector{Int64}, Vector{Int64}, Vector{Int64}, Vector{Int64}}, Nothing, Nothing, Nothing, SparseArrays.SparseMatrixCSC{Bool, Int64}, Nothing, Nothing, Nothing, Nothing, Nothing}, OrdinaryDiffEq.DABDF2Cache{Vector{Float64}, Vector{Float64}, Vector{Float64}, OrdinaryDiffEq.NLSolver{NLNewton{Rational{Int64}, Rational{Int64}, Rational{Int64}}, true, Vector{Float64}, Rational{Int64}, Nothing, Float64, OrdinaryDiffEq.NLNewtonCache{Vector{Float64}, Float64, Float64, Vector{Float64}, Matrix{Float64}, Matrix{Float64}, Nothing, Nothing, LinearSolve.LinearCache{Matrix{Float64}, Vector{Float64}, Vector{Float64}, SciMLBase.NullParameters, Nothing, LU{Float64, Matrix{Float64}}, LinearSolve.InvPreconditioner{Diagonal{Float64, Vector{Float64}}}, Diagonal{Float64, Vector{Float64}}, Float64}}}, Float64}, OrdinaryDiffEq.DEOptions{Float64, Float64, Float64, Float64, PIController{Rational{Int64}}, typeof(DiffEqBase.ODE_DEFAULT_NORM), typeof(opnorm), Nothing, 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{}, Tuple{}, Tuple{}}, Vector{Float64}, Float64, Nothing, OrdinaryDiffEq.DefaultInit}, cache::OrdinaryDiffEq.DImplicitEulerCache{Vector{Float64}, Vector{Float64}, Vector{Float64}, OrdinaryDiffEq.NLSolver{NLNewton{Rational{Int64}, Rational{Int64}, Rational{Int64}}, true, Vector{Float64}, Rational{Int64}, Nothing, Float64, OrdinaryDiffEq.NLNewtonCache{Vector{Float64}, Float64, Float64, Vector{Float64}, Matrix{Float64}, Matrix{Float64}, Nothing, Nothing, LinearSolve.LinearCache{Matrix{Float64}, Vector{Float64}, Vector{Float64}, SciMLBase.NullParameters, Nothing, LU{Float64, Matrix{Float64}}, LinearSolve.InvPreconditioner{Diagonal{Float64, Vector{Float64}}}, Diagonal{Float64, Vector{Float64}}, Float64}}}}, repeat_step::Bool)
      @ OrdinaryDiffEq ~/.julia/packages/OrdinaryDiffEq/Op0Oq/src/nlsolve/nlsolve.jl:21
    [5] perform_step!(integrator::OrdinaryDiffEq.ODEIntegrator{DABDF2{12, true, Nothing, NLNewton{Rational{Int64}, Rational{Int64}, Rational{Int64}}, typeof(OrdinaryDiffEq.DEFAULT_PRECS), Val{:forward}, true, nothing}, true, Vector{Float64}, Vector{Float64}, Float64, Tuple{var"#103#105"{Int64}, Dict{Int64, DistributedLoads{Float64}}, Dict{Int64, Vector{PointMass{Float64}}}, StaticArrays.SVector{3, Float64}, StaticArrays.SVector{3, Float64}, StaticArrays.SVector{3, Float64}, StaticArrays.SVector{3, Float64}, StaticArrays.SVector{3, Float64}, StaticArrays.SVector{3, Float64}}, Float64, Float64, Float64, Float64, Vector{Vector{Float64}}, DAESolution{Float64, 2, Vector{Vector{Float64}}, Nothing, Nothing, Nothing, Vector{Float64}, DAEProblem{Vector{Float64}, Vector{Float64}, Tuple{Float64, Float64}, true, Tuple{var"#103#105"{Int64}, Dict{Int64, DistributedLoads{Float64}}, Dict{Int64, Vector{PointMass{Float64}}}, StaticArrays.SVector{3, Float64}, StaticArrays.SVector{3, Float64}, StaticArrays.SVector{3, Float64}, StaticArrays.SVector{3, Float64}, StaticArrays.SVector{3, Float64}, StaticArrays.SVector{3, Float64}}, DAEFunction{true, GXBeam.var"#299#301"{Assembly{Float64, Vector{Vector{Float64}}, UnitRange{Int64}, Vector{GXBeam.Element{Float64}}}, Float64, Vector{Int64}, Vector{Int64}, Vector{Int64}, Vector{Int64}, Vector{Int64}, Vector{Int64}}, Nothing, Nothing, GXBeam.var"#300#302"{Assembly{Float64, Vector{Vector{Float64}}, UnitRange{Int64}, Vector{GXBeam.Element{Float64}}}, Float64, Vector{Int64}, Vector{Int64}, Vector{Int64}, Vector{Int64}, Vector{Int64}, Vector{Int64}}, Nothing, Nothing, Nothing, SparseArrays.SparseMatrixCSC{Bool, Int64}, Nothing, Nothing, Nothing, Nothing, Nothing}, Base.Iterators.Pairs{Union{}, Union{}, Tuple{}, NamedTuple{(), Tuple{}}}, Vector{Bool}}, DABDF2{12, true, Nothing, NLNewton{Rational{Int64}, Rational{Int64}, Rational{Int64}}, typeof(OrdinaryDiffEq.DEFAULT_PRECS), Val{:forward}, true, nothing}, OrdinaryDiffEq.InterpolationData{DAEFunction{true, GXBeam.var"#299#301"{Assembly{Float64, Vector{Vector{Float64}}, UnitRange{Int64}, Vector{GXBeam.Element{Float64}}}, Float64, Vector{Int64}, Vector{Int64}, Vector{Int64}, Vector{Int64}, Vector{Int64}, Vector{Int64}}, Nothing, Nothing, GXBeam.var"#300#302"{Assembly{Float64, Vector{Vector{Float64}}, UnitRange{Int64}, Vector{GXBeam.Element{Float64}}}, Float64, Vector{Int64}, Vector{Int64}, Vector{Int64}, Vector{Int64}, Vector{Int64}, Vector{Int64}}, Nothing, Nothing, Nothing, SparseArrays.SparseMatrixCSC{Bool, Int64}, Nothing, Nothing, Nothing, Nothing, Nothing}, Vector{Vector{Float64}}, Vector{Float64}, Vector{Vector{Vector{Float64}}}, OrdinaryDiffEq.DABDF2Cache{Vector{Float64}, Vector{Float64}, Vector{Float64}, OrdinaryDiffEq.NLSolver{NLNewton{Rational{Int64}, Rational{Int64}, Rational{Int64}}, true, Vector{Float64}, Rational{Int64}, Nothing, Float64, OrdinaryDiffEq.NLNewtonCache{Vector{Float64}, Float64, Float64, Vector{Float64}, Matrix{Float64}, Matrix{Float64}, Nothing, Nothing, LinearSolve.LinearCache{Matrix{Float64}, Vector{Float64}, Vector{Float64}, SciMLBase.NullParameters, Nothing, LU{Float64, Matrix{Float64}}, LinearSolve.InvPreconditioner{Diagonal{Float64, Vector{Float64}}}, Diagonal{Float64, Vector{Float64}}, Float64}}}, Float64}}, DiffEqBase.DEStats}, DAEFunction{true, GXBeam.var"#299#301"{Assembly{Float64, Vector{Vector{Float64}}, UnitRange{Int64}, Vector{GXBeam.Element{Float64}}}, Float64, Vector{Int64}, Vector{Int64}, Vector{Int64}, Vector{Int64}, Vector{Int64}, Vector{Int64}}, Nothing, Nothing, GXBeam.var"#300#302"{Assembly{Float64, Vector{Vector{Float64}}, UnitRange{Int64}, Vector{GXBeam.Element{Float64}}}, Float64, Vector{Int64}, Vector{Int64}, Vector{Int64}, Vector{Int64}, Vector{Int64}, Vector{Int64}}, Nothing, Nothing, Nothing, SparseArrays.SparseMatrixCSC{Bool, Int64}, Nothing, Nothing, Nothing, Nothing, Nothing}, OrdinaryDiffEq.DABDF2Cache{Vector{Float64}, Vector{Float64}, Vector{Float64}, OrdinaryDiffEq.NLSolver{NLNewton{Rational{Int64}, Rational{Int64}, Rational{Int64}}, true, Vector{Float64}, Rational{Int64}, Nothing, Float64, OrdinaryDiffEq.NLNewtonCache{Vector{Float64}, Float64, Float64, Vector{Float64}, Matrix{Float64}, Matrix{Float64}, Nothing, Nothing, LinearSolve.LinearCache{Matrix{Float64}, Vector{Float64}, Vector{Float64}, SciMLBase.NullParameters, Nothing, LU{Float64, Matrix{Float64}}, LinearSolve.InvPreconditioner{Diagonal{Float64, Vector{Float64}}}, Diagonal{Float64, Vector{Float64}}, Float64}}}, Float64}, OrdinaryDiffEq.DEOptions{Float64, Float64, Float64, Float64, PIController{Rational{Int64}}, typeof(DiffEqBase.ODE_DEFAULT_NORM), typeof(opnorm), Nothing, 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{}, Tuple{}, Tuple{}}, Vector{Float64}, Float64, Nothing, OrdinaryDiffEq.DefaultInit}, cache::OrdinaryDiffEq.DImplicitEulerCache{Vector{Float64}, Vector{Float64}, Vector{Float64}, OrdinaryDiffEq.NLSolver{NLNewton{Rational{Int64}, Rational{Int64}, Rational{Int64}}, true, Vector{Float64}, Rational{Int64}, Nothing, Float64, OrdinaryDiffEq.NLNewtonCache{Vector{Float64}, Float64, Float64, Vector{Float64}, Matrix{Float64}, Matrix{Float64}, Nothing, Nothing, LinearSolve.LinearCache{Matrix{Float64}, Vector{Float64}, Vector{Float64}, SciMLBase.NullParameters, Nothing, LU{Float64, Matrix{Float64}}, LinearSolve.InvPreconditioner{Diagonal{Float64, Vector{Float64}}}, Diagonal{Float64, Vector{Float64}}, Float64}}}}, repeat_step::Bool)
      @ OrdinaryDiffEq ~/.julia/packages/OrdinaryDiffEq/Op0Oq/src/perform_step/dae_perform_step.jl:67
    [6] perform_step!(integrator::OrdinaryDiffEq.ODEIntegrator{DABDF2{12, true, Nothing, NLNewton{Rational{Int64}, Rational{Int64}, Rational{Int64}}, typeof(OrdinaryDiffEq.DEFAULT_PRECS), Val{:forward}, true, nothing}, true, Vector{Float64}, Vector{Float64}, Float64, Tuple{var"#103#105"{Int64}, Dict{Int64, DistributedLoads{Float64}}, Dict{Int64, Vector{PointMass{Float64}}}, StaticArrays.SVector{3, Float64}, StaticArrays.SVector{3, Float64}, StaticArrays.SVector{3, Float64}, StaticArrays.SVector{3, Float64}, StaticArrays.SVector{3, Float64}, StaticArrays.SVector{3, Float64}}, Float64, Float64, Float64, Float64, Vector{Vector{Float64}}, DAESolution{Float64, 2, Vector{Vector{Float64}}, Nothing, Nothing, Nothing, Vector{Float64}, DAEProblem{Vector{Float64}, Vector{Float64}, Tuple{Float64, Float64}, true, Tuple{var"#103#105"{Int64}, Dict{Int64, DistributedLoads{Float64}}, Dict{Int64, Vector{PointMass{Float64}}}, StaticArrays.SVector{3, Float64}, StaticArrays.SVector{3, Float64}, StaticArrays.SVector{3, Float64}, StaticArrays.SVector{3, Float64}, StaticArrays.SVector{3, Float64}, StaticArrays.SVector{3, Float64}}, DAEFunction{true, GXBeam.var"#299#301"{Assembly{Float64, Vector{Vector{Float64}}, UnitRange{Int64}, Vector{GXBeam.Element{Float64}}}, Float64, Vector{Int64}, Vector{Int64}, Vector{Int64}, Vector{Int64}, Vector{Int64}, Vector{Int64}}, Nothing, Nothing, GXBeam.var"#300#302"{Assembly{Float64, Vector{Vector{Float64}}, UnitRange{Int64}, Vector{GXBeam.Element{Float64}}}, Float64, Vector{Int64}, Vector{Int64}, Vector{Int64}, Vector{Int64}, Vector{Int64}, Vector{Int64}}, Nothing, Nothing, Nothing, SparseArrays.SparseMatrixCSC{Bool, Int64}, Nothing, Nothing, Nothing, Nothing, Nothing}, Base.Iterators.Pairs{Union{}, Union{}, Tuple{}, NamedTuple{(), Tuple{}}}, Vector{Bool}}, DABDF2{12, true, Nothing, NLNewton{Rational{Int64}, Rational{Int64}, Rational{Int64}}, typeof(OrdinaryDiffEq.DEFAULT_PRECS), Val{:forward}, true, nothing}, OrdinaryDiffEq.InterpolationData{DAEFunction{true, GXBeam.var"#299#301"{Assembly{Float64, Vector{Vector{Float64}}, UnitRange{Int64}, Vector{GXBeam.Element{Float64}}}, Float64, Vector{Int64}, Vector{Int64}, Vector{Int64}, Vector{Int64}, Vector{Int64}, Vector{Int64}}, Nothing, Nothing, GXBeam.var"#300#302"{Assembly{Float64, Vector{Vector{Float64}}, UnitRange{Int64}, Vector{GXBeam.Element{Float64}}}, Float64, Vector{Int64}, Vector{Int64}, Vector{Int64}, Vector{Int64}, Vector{Int64}, Vector{Int64}}, Nothing, Nothing, Nothing, SparseArrays.SparseMatrixCSC{Bool, Int64}, Nothing, Nothing, Nothing, Nothing, Nothing}, Vector{Vector{Float64}}, Vector{Float64}, Vector{Vector{Vector{Float64}}}, OrdinaryDiffEq.DABDF2Cache{Vector{Float64}, Vector{Float64}, Vector{Float64}, OrdinaryDiffEq.NLSolver{NLNewton{Rational{Int64}, Rational{Int64}, Rational{Int64}}, true, Vector{Float64}, Rational{Int64}, Nothing, Float64, OrdinaryDiffEq.NLNewtonCache{Vector{Float64}, Float64, Float64, Vector{Float64}, Matrix{Float64}, Matrix{Float64}, Nothing, Nothing, LinearSolve.LinearCache{Matrix{Float64}, Vector{Float64}, Vector{Float64}, SciMLBase.NullParameters, Nothing, LU{Float64, Matrix{Float64}}, LinearSolve.InvPreconditioner{Diagonal{Float64, Vector{Float64}}}, Diagonal{Float64, Vector{Float64}}, Float64}}}, Float64}}, DiffEqBase.DEStats}, DAEFunction{true, GXBeam.var"#299#301"{Assembly{Float64, Vector{Vector{Float64}}, UnitRange{Int64}, Vector{GXBeam.Element{Float64}}}, Float64, Vector{Int64}, Vector{Int64}, Vector{Int64}, Vector{Int64}, Vector{Int64}, Vector{Int64}}, Nothing, Nothing, GXBeam.var"#300#302"{Assembly{Float64, Vector{Vector{Float64}}, UnitRange{Int64}, Vector{GXBeam.Element{Float64}}}, Float64, Vector{Int64}, Vector{Int64}, Vector{Int64}, Vector{Int64}, Vector{Int64}, Vector{Int64}}, Nothing, Nothing, Nothing, SparseArrays.SparseMatrixCSC{Bool, Int64}, Nothing, Nothing, Nothing, Nothing, Nothing}, OrdinaryDiffEq.DABDF2Cache{Vector{Float64}, Vector{Float64}, Vector{Float64}, OrdinaryDiffEq.NLSolver{NLNewton{Rational{Int64}, Rational{Int64}, Rational{Int64}}, true, Vector{Float64}, Rational{Int64}, Nothing, Float64, OrdinaryDiffEq.NLNewtonCache{Vector{Float64}, Float64, Float64, Vector{Float64}, Matrix{Float64}, Matrix{Float64}, Nothing, Nothing, LinearSolve.LinearCache{Matrix{Float64}, Vector{Float64}, Vector{Float64}, SciMLBase.NullParameters, Nothing, LU{Float64, Matrix{Float64}}, LinearSolve.InvPreconditioner{Diagonal{Float64, Vector{Float64}}}, Diagonal{Float64, Vector{Float64}}, Float64}}}, Float64}, OrdinaryDiffEq.DEOptions{Float64, Float64, Float64, Float64, PIController{Rational{Int64}}, typeof(DiffEqBase.ODE_DEFAULT_NORM), typeof(opnorm), Nothing, 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{}, Tuple{}, Tuple{}}, Vector{Float64}, Float64, Nothing, OrdinaryDiffEq.DefaultInit}, cache::OrdinaryDiffEq.DABDF2Cache{Vector{Float64}, Vector{Float64}, Vector{Float64}, OrdinaryDiffEq.NLSolver{NLNewton{Rational{Int64}, Rational{Int64}, Rational{Int64}}, true, Vector{Float64}, Rational{Int64}, Nothing, Float64, OrdinaryDiffEq.NLNewtonCache{Vector{Float64}, Float64, Float64, Vector{Float64}, Matrix{Float64}, Matrix{Float64}, Nothing, Nothing, LinearSolve.LinearCache{Matrix{Float64}, Vector{Float64}, Vector{Float64}, SciMLBase.NullParameters, Nothing, LU{Float64, Matrix{Float64}}, LinearSolve.InvPreconditioner{Diagonal{Float64, Vector{Float64}}}, Diagonal{Float64, Vector{Float64}}, Float64}}}, Float64}, repeat_step::Bool)
      @ OrdinaryDiffEq ~/.julia/packages/OrdinaryDiffEq/Op0Oq/src/perform_step/dae_perform_step.jl:179
    [7] perform_step!
      @ ~/.julia/packages/OrdinaryDiffEq/Op0Oq/src/perform_step/dae_perform_step.jl:171 [inlined]
    [8] solve!(integrator::OrdinaryDiffEq.ODEIntegrator{DABDF2{12, true, Nothing, NLNewton{Rational{Int64}, Rational{Int64}, Rational{Int64}}, typeof(OrdinaryDiffEq.DEFAULT_PRECS), Val{:forward}, true, nothing}, true, Vector{Float64}, Vector{Float64}, Float64, Tuple{var"#103#105"{Int64}, Dict{Int64, DistributedLoads{Float64}}, Dict{Int64, Vector{PointMass{Float64}}}, StaticArrays.SVector{3, Float64}, StaticArrays.SVector{3, Float64}, StaticArrays.SVector{3, Float64}, StaticArrays.SVector{3, Float64}, StaticArrays.SVector{3, Float64}, StaticArrays.SVector{3, Float64}}, Float64, Float64, Float64, Float64, Vector{Vector{Float64}}, DAESolution{Float64, 2, Vector{Vector{Float64}}, Nothing, Nothing, Nothing, Vector{Float64}, DAEProblem{Vector{Float64}, Vector{Float64}, Tuple{Float64, Float64}, true, Tuple{var"#103#105"{Int64}, Dict{Int64, DistributedLoads{Float64}}, Dict{Int64, Vector{PointMass{Float64}}}, StaticArrays.SVector{3, Float64}, StaticArrays.SVector{3, Float64}, StaticArrays.SVector{3, Float64}, StaticArrays.SVector{3, Float64}, StaticArrays.SVector{3, Float64}, StaticArrays.SVector{3, Float64}}, DAEFunction{true, GXBeam.var"#299#301"{Assembly{Float64, Vector{Vector{Float64}}, UnitRange{Int64}, Vector{GXBeam.Element{Float64}}}, Float64, Vector{Int64}, Vector{Int64}, Vector{Int64}, Vector{Int64}, Vector{Int64}, Vector{Int64}}, Nothing, Nothing, GXBeam.var"#300#302"{Assembly{Float64, Vector{Vector{Float64}}, UnitRange{Int64}, Vector{GXBeam.Element{Float64}}}, Float64, Vector{Int64}, Vector{Int64}, Vector{Int64}, Vector{Int64}, Vector{Int64}, Vector{Int64}}, Nothing, Nothing, Nothing, SparseArrays.SparseMatrixCSC{Bool, Int64}, Nothing, Nothing, Nothing, Nothing, Nothing}, Base.Iterators.Pairs{Union{}, Union{}, Tuple{}, NamedTuple{(), Tuple{}}}, Vector{Bool}}, DABDF2{12, true, Nothing, NLNewton{Rational{Int64}, Rational{Int64}, Rational{Int64}}, typeof(OrdinaryDiffEq.DEFAULT_PRECS), Val{:forward}, true, nothing}, OrdinaryDiffEq.InterpolationData{DAEFunction{true, GXBeam.var"#299#301"{Assembly{Float64, Vector{Vector{Float64}}, UnitRange{Int64}, Vector{GXBeam.Element{Float64}}}, Float64, Vector{Int64}, Vector{Int64}, Vector{Int64}, Vector{Int64}, Vector{Int64}, Vector{Int64}}, Nothing, Nothing, GXBeam.var"#300#302"{Assembly{Float64, Vector{Vector{Float64}}, UnitRange{Int64}, Vector{GXBeam.Element{Float64}}}, Float64, Vector{Int64}, Vector{Int64}, Vector{Int64}, Vector{Int64}, Vector{Int64}, Vector{Int64}}, Nothing, Nothing, Nothing, SparseArrays.SparseMatrixCSC{Bool, Int64}, Nothing, Nothing, Nothing, Nothing, Nothing}, Vector{Vector{Float64}}, Vector{Float64}, Vector{Vector{Vector{Float64}}}, OrdinaryDiffEq.DABDF2Cache{Vector{Float64}, Vector{Float64}, Vector{Float64}, OrdinaryDiffEq.NLSolver{NLNewton{Rational{Int64}, Rational{Int64}, Rational{Int64}}, true, Vector{Float64}, Rational{Int64}, Nothing, Float64, OrdinaryDiffEq.NLNewtonCache{Vector{Float64}, Float64, Float64, Vector{Float64}, Matrix{Float64}, Matrix{Float64}, Nothing, Nothing, LinearSolve.LinearCache{Matrix{Float64}, Vector{Float64}, Vector{Float64}, SciMLBase.NullParameters, Nothing, LU{Float64, Matrix{Float64}}, LinearSolve.InvPreconditioner{Diagonal{Float64, Vector{Float64}}}, Diagonal{Float64, Vector{Float64}}, Float64}}}, Float64}}, DiffEqBase.DEStats}, DAEFunction{true, GXBeam.var"#299#301"{Assembly{Float64, Vector{Vector{Float64}}, UnitRange{Int64}, Vector{GXBeam.Element{Float64}}}, Float64, Vector{Int64}, Vector{Int64}, Vector{Int64}, Vector{Int64}, Vector{Int64}, Vector{Int64}}, Nothing, Nothing, GXBeam.var"#300#302"{Assembly{Float64, Vector{Vector{Float64}}, UnitRange{Int64}, Vector{GXBeam.Element{Float64}}}, Float64, Vector{Int64}, Vector{Int64}, Vector{Int64}, Vector{Int64}, Vector{Int64}, Vector{Int64}}, Nothing, Nothing, Nothing, SparseArrays.SparseMatrixCSC{Bool, Int64}, Nothing, Nothing, Nothing, Nothing, Nothing}, OrdinaryDiffEq.DABDF2Cache{Vector{Float64}, Vector{Float64}, Vector{Float64}, OrdinaryDiffEq.NLSolver{NLNewton{Rational{Int64}, Rational{Int64}, Rational{Int64}}, true, Vector{Float64}, Rational{Int64}, Nothing, Float64, OrdinaryDiffEq.NLNewtonCache{Vector{Float64}, Float64, Float64, Vector{Float64}, Matrix{Float64}, Matrix{Float64}, Nothing, Nothing, LinearSolve.LinearCache{Matrix{Float64}, Vector{Float64}, Vector{Float64}, SciMLBase.NullParameters, Nothing, LU{Float64, Matrix{Float64}}, LinearSolve.InvPreconditioner{Diagonal{Float64, Vector{Float64}}}, Diagonal{Float64, Vector{Float64}}, Float64}}}, Float64}, OrdinaryDiffEq.DEOptions{Float64, Float64, Float64, Float64, PIController{Rational{Int64}}, typeof(DiffEqBase.ODE_DEFAULT_NORM), typeof(opnorm), Nothing, 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{}, Tuple{}, Tuple{}}, Vector{Float64}, Float64, Nothing, OrdinaryDiffEq.DefaultInit})
      @ OrdinaryDiffEq ~/.julia/packages/OrdinaryDiffEq/Op0Oq/src/solve.jl:478
    [9] #__solve#495
      @ ~/.julia/packages/OrdinaryDiffEq/Op0Oq/src/solve.jl:5 [inlined]
   [10] __solve
      @ ~/.julia/packages/OrdinaryDiffEq/Op0Oq/src/solve.jl:4 [inlined]
   [11] #solve_call#37
      @ ~/.julia/packages/DiffEqBase/1V2xg/src/solve.jl:61 [inlined]
   [12] solve_call(_prob::DAEProblem{Vector{Float64}, Vector{Float64}, Tuple{Float64, Float64}, true, Tuple{var"#103#105"{Int64}, Dict{Int64, DistributedLoads{Float64}}, Dict{Int64, Vector{PointMass{Float64}}}, StaticArrays.SVector{3, Float64}, StaticArrays.SVector{3, Float64}, StaticArrays.SVector{3, Float64}, StaticArrays.SVector{3, Float64}, StaticArrays.SVector{3, Float64}, StaticArrays.SVector{3, Float64}}, DAEFunction{true, GXBeam.var"#299#301"{Assembly{Float64, Vector{Vector{Float64}}, UnitRange{Int64}, Vector{GXBeam.Element{Float64}}}, Float64, Vector{Int64}, Vector{Int64}, Vector{Int64}, Vector{Int64}, Vector{Int64}, Vector{Int64}}, Nothing, Nothing, GXBeam.var"#300#302"{Assembly{Float64, Vector{Vector{Float64}}, UnitRange{Int64}, Vector{GXBeam.Element{Float64}}}, Float64, Vector{Int64}, Vector{Int64}, Vector{Int64}, Vector{Int64}, Vector{Int64}, Vector{Int64}}, Nothing, Nothing, Nothing, SparseArrays.SparseMatrixCSC{Bool, Int64}, Nothing, Nothing, Nothing, Nothing, Nothing}, Base.Iterators.Pairs{Union{}, Union{}, Tuple{}, NamedTuple{(), Tuple{}}}, Vector{Bool}}, args::DABDF2{12, true, Nothing, NLNewton{Rational{Int64}, Rational{Int64}, Rational{Int64}}, typeof(OrdinaryDiffEq.DEFAULT_PRECS), Val{:forward}, true, nothing})
      @ DiffEqBase ~/.julia/packages/DiffEqBase/1V2xg/src/solve.jl:48
   [13] solve_up(prob::DAEProblem{Vector{Float64}, Vector{Float64}, Tuple{Float64, Float64}, true, Tuple{var"#103#105"{Int64}, Dict{Int64, DistributedLoads{Float64}}, Dict{Int64, Vector{PointMass{Float64}}}, StaticArrays.SVector{3, Float64}, StaticArrays.SVector{3, Float64}, StaticArrays.SVector{3, Float64}, StaticArrays.SVector{3, Float64}, StaticArrays.SVector{3, Float64}, StaticArrays.SVector{3, Float64}}, DAEFunction{true, GXBeam.var"#299#301"{Assembly{Float64, Vector{Vector{Float64}}, UnitRange{Int64}, Vector{GXBeam.Element{Float64}}}, Float64, Vector{Int64}, Vector{Int64}, Vector{Int64}, Vector{Int64}, Vector{Int64}, Vector{Int64}}, Nothing, Nothing, GXBeam.var"#300#302"{Assembly{Float64, Vector{Vector{Float64}}, UnitRange{Int64}, Vector{GXBeam.Element{Float64}}}, Float64, Vector{Int64}, Vector{Int64}, Vector{Int64}, Vector{Int64}, Vector{Int64}, Vector{Int64}}, Nothing, Nothing, Nothing, SparseArrays.SparseMatrixCSC{Bool, Int64}, Nothing, Nothing, Nothing, Nothing, Nothing}, Base.Iterators.Pairs{Union{}, Union{}, Tuple{}, NamedTuple{(), Tuple{}}}, Vector{Bool}}, sensealg::Nothing, u0::Vector{Float64}, p::Tuple{var"#103#105"{Int64}, Dict{Int64, DistributedLoads{Float64}}, Dict{Int64, Vector{PointMass{Float64}}}, StaticArrays.SVector{3, Float64}, StaticArrays.SVector{3, Float64}, StaticArrays.SVector{3, Float64}, StaticArrays.SVector{3, Float64}, StaticArrays.SVector{3, Float64}, StaticArrays.SVector{3, Float64}}, args::DABDF2{0, true, Nothing, NLNewton{Rational{Int64}, Rational{Int64}, Rational{Int64}}, typeof(OrdinaryDiffEq.DEFAULT_PRECS), Val{:forward}, true, nothing}; kwargs::Base.Iterators.Pairs{Union{}, Union{}, Tuple{}, NamedTuple{(), Tuple{}}})
      @ DiffEqBase ~/.julia/packages/DiffEqBase/1V2xg/src/solve.jl:87
   [14] solve_up
      @ ~/.julia/packages/DiffEqBase/1V2xg/src/solve.jl:78 [inlined]
   [15] #solve#38
      @ ~/.julia/packages/DiffEqBase/1V2xg/src/solve.jl:73 [inlined]
   [16] solve(prob::DAEProblem{Vector{Float64}, Vector{Float64}, Tuple{Float64, Float64}, true, Tuple{var"#103#105"{Int64}, Dict{Int64, DistributedLoads{Float64}}, Dict{Int64, Vector{PointMass{Float64}}}, StaticArrays.SVector{3, Float64}, StaticArrays.SVector{3, Float64}, StaticArrays.SVector{3, Float64}, StaticArrays.SVector{3, Float64}, StaticArrays.SVector{3, Float64}, StaticArrays.SVector{3, Float64}}, DAEFunction{true, GXBeam.var"#299#301"{Assembly{Float64, Vector{Vector{Float64}}, UnitRange{Int64}, Vector{GXBeam.Element{Float64}}}, Float64, Vector{Int64}, Vector{Int64}, Vector{Int64}, Vector{Int64}, Vector{Int64}, Vector{Int64}}, Nothing, Nothing, GXBeam.var"#300#302"{Assembly{Float64, Vector{Vector{Float64}}, UnitRange{Int64}, Vector{GXBeam.Element{Float64}}}, Float64, Vector{Int64}, Vector{Int64}, Vector{Int64}, Vector{Int64}, Vector{Int64}, Vector{Int64}}, Nothing, Nothing, Nothing, SparseArrays.SparseMatrixCSC{Bool, Int64}, Nothing, Nothing, Nothing, Nothing, Nothing}, Base.Iterators.Pairs{Union{}, Union{}, Tuple{}, NamedTuple{(), Tuple{}}}, Vector{Bool}}, args::DABDF2{0, true, Nothing, NLNewton{Rational{Int64}, Rational{Int64}, Rational{Int64}}, typeof(OrdinaryDiffEq.DEFAULT_PRECS), Val{:forward}, true, nothing})
      @ DiffEqBase ~/.julia/packages/DiffEqBase/1V2xg/src/solve.jl:68
   [17] macro expansion
      @ ~/.julia/packages/GXBeam/S8Yc6/test/runtests.jl:1206 [inlined]
   [18] macro expansion
      @ /buildworker/worker/package_linux64/build/usr/share/julia/stdlib/v1.6/Test/src/Test.jl:1151 [inlined]
   [19] top-level scope
      @ ~/.julia/packages/GXBeam/S8Yc6/test/runtests.jl:1141
   [20] include(fname::String)
      @ Base.MainInclude ./client.jl:444
   [21] top-level scope
      @ none:6
   [22] eval
      @ ./boot.jl:360 [inlined]
   [23] exec_options(opts::Base.JLOptions)
      @ Base ./client.jl:261
   [24] _start()
      @ Base ./client.jl:485
Test Summary:         | Fail  Error  Total
DifferentialEquations |    1      1      2
ERROR: LoadError: Some tests did not pass: 0 passed, 1 failed, 1 errored, 0 broken.
in expression starting at /home/kevin/.julia/packages/GXBeam/S8Yc6/test/runtests.jl:1139
ERROR: Package GXBeam errored during testing
@taylormcd
Copy link
Collaborator

taylormcd commented Feb 4, 2022

This is odd considering the CI tests pass, though I also recently observed this behavior in the development branch's documentation. It's probably just a poorly posed DAE problem, at least when used with the DifferentialEquations solvers. I'll look into this.

@taylormcd
Copy link
Collaborator

Looks like the issue happens when using the latest version of OrdinaryDiffEq. I'll probably just cap the version until I can get this resolved.

@taylormcd
Copy link
Collaborator

taylormcd commented Feb 7, 2022

I removed the DAE custom jacobians and will add them back in when SciML/DifferentialEquations.jl#846 is resolved. I also explicitly marked the ODE functions as experimental (but not the DAE functions) and removed them from the tests and documentation since non-constant mass matrices are not directly supported by DifferentialEquations.

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

No branches or pull requests

2 participants