Hi,
I am running a parameter estimation on PEtab and I get the following error when I use InterpolatingAdjoint(autojacvec = ReverseDiffVJP(true)) or
InterpolatingAdjoint(autojacvec = EnzymeVJP()))
type FakeIntegrator has no field sol
Stacktrace:
[1] getproperty(fi::SciMLSensitivity.FakeIntegrator{Vector{…}, Vector{…}, ReverseDiff.TrackedReal{…}, Float64}, s::Symbol)
@ SciMLSensitivity C:\Users\memmo\.julia\packages\SciMLSensitivity\HRhwU\src\callback_tracking.jl:174
[2] macro expansion
@ .\none:5 [inlined]
[3] macro expansion
@ C:\Users\memmo\.julia\packages\RuntimeGeneratedFunctions\M9ZX8\src\RuntimeGeneratedFunctions.jl:163 [inlined]
[4] macro expansion
@ .\none:0 [inlined]
[5] generated_callfunc
@ .\none:0 [inlined]
[6] (::RuntimeGeneratedFunctions.RuntimeGeneratedFunction{…})(args::SciMLSensitivity.FakeIntegrator{…})
@ RuntimeGeneratedFunctions C:\Users\memmo\.julia\packages\RuntimeGeneratedFunctions\M9ZX8\src\RuntimeGeneratedFunctions.jl:150
[7] (::SciMLSensitivity.var"#276#278"{…})(du::Vector{…}, u::ReverseDiff.TrackedArray{…}, p::ReverseDiff.TrackedArray{…}, t::ReverseDiff.TrackedReal{…})
@ SciMLSensitivity C:\Users\memmo\.julia\packages\SciMLSensitivity\HRhwU\src\callback_tracking.jl:389
[8] (::SciMLSensitivity.var"#126#130"{…})(u::ReverseDiff.TrackedArray{…}, p::ReverseDiff.TrackedArray{…}, t::ReverseDiff.TrackedArray{…})
@ SciMLSensitivity C:\Users\memmo\.julia\packages\SciMLSensitivity\HRhwU\src\adjoint_common.jl:385
[9] ReverseDiff.GradientTape(f::Function, input::Tuple{…}, cfg::ReverseDiff.GradientConfig{…})
@ ReverseDiff C:\Users\memmo\.julia\packages\ReverseDiff\p1MzG\src\api\tape.jl:207
[10] ReverseDiff.GradientTape(f::Function, input::Tuple{Vector{Float64}, Vector{Float64}, Vector{Float64}})
@ ReverseDiff C:\Users\memmo\.julia\packages\ReverseDiff\p1MzG\src\api\tape.jl:204
[11] get_paramjac_config(autojacvec::ReverseDiffVJP{…}, p::Vector{…}, f::SciMLSensitivity.var"#276#278"{…}, y::Vector{…}, _p::Vector{…}, _t::Float64; numindvar::Int64, alg::Nothing, isinplace::Bool, isRODE::Bool, _W::Nothing)
@ SciMLSensitivity C:\Users\memmo\.julia\packages\SciMLSensitivity\HRhwU\src\adjoint_common.jl:381
[12] get_cb_diffcaches(cb::DiscreteCallback{…}, autojacvec::ReverseDiffVJP{…})
@ SciMLSensitivity C:\Users\memmo\.julia\packages\SciMLSensitivity\HRhwU\src\callback_tracking.jl:468
[13] _setup_reverse_callbacks(cb::DiscreteCallback{…}, affect::SciMLSensitivity.TrackedAffect{…}, sensealg::InterpolatingAdjoint{…}, dgdu::Function, dgdp::Nothing, loss_ref::Base.RefValue{…}, terminated::Bool)
@ SciMLSensitivity C:\Users\memmo\.julia\packages\SciMLSensitivity\HRhwU\src\callback_tracking.jl:255
[14] _setup_reverse_callbacks
@ C:\Users\memmo\.julia\packages\SciMLSensitivity\HRhwU\src\callback_tracking.jl:219 [inlined]
[15] _broadcast_getindex_evalf
@ .\broadcast.jl:709 [inlined]
[16] _broadcast_getindex
@ .\broadcast.jl:682 [inlined]
[17] (::Base.Broadcast.var"#31#32"{Base.Broadcast.Broadcasted{…}})(k::Int64)
@ Base.Broadcast .\broadcast.jl:1118
[18] ntuple
@ .\ntuple.jl:48 [inlined]
[19] copy
@ .\broadcast.jl:1118 [inlined]
[20] materialize
@ .\broadcast.jl:903 [inlined]
[21] setup_reverse_callbacks(cb::CallbackSet{…}, sensealg::InterpolatingAdjoint{…}, dgdu::Function, dgdp::Nothing, cur_time::Base.RefValue{…}, terminated::Bool)
@ SciMLSensitivity C:\Users\memmo\.julia\packages\SciMLSensitivity\HRhwU\src\callback_tracking.jl:200
[22] generate_callbacks(sensefun::SciMLSensitivity.ODEInterpolatingAdjointSensitivityFunction{…}, dgdu::Function, dgdp::Nothing, λ::Vector{…}, t::Vector{…}, t0::Float64, callback::CallbackSet{…}, init_cb::Bool, terminated::Bool)
@ SciMLSensitivity C:\Users\memmo\.julia\packages\SciMLSensitivity\HRhwU\src\adjoint_common.jl:617
[23] ODEAdjointProblem(sol::ODESolution{…}, sensealg::InterpolatingAdjoint{…}, alg::TRBDF2{…}, t::Vector{…}, dgdu_discrete::PEtab.var"#135#143"{…}, dgdp_discrete::Nothing, dgdu_continuous::Nothing, dgdp_continuous::Nothing, g::Nothing, ::Val{…}; checkpoints::Vector{…}, callback::CallbackSet{…}, reltol::Float64, abstol::Float64, kwargs::@Kwargs{})
@ SciMLSensitivity C:\Users\memmo\.julia\packages\SciMLSensitivity\HRhwU\src\interpolating_adjoint.jl:357
[24] __adjoint_sensitivities!(_du::Vector{…}, _dp::Vector{…}, sol::ODESolution{…}, sensealg::InterpolatingAdjoint{…}, t::Vector{…}, solver::TRBDF2{…}, abstol::Float64, reltol::Float64, callback::CallbackSet{…}, compute_∂G∂u::PEtab.var"#135#143"{…}; kwargs::@Kwargs{…})
@ PEtabSciMLSensitivityExtension C:\Users\memmo\.julia\packages\PEtab\MalGN\ext\PEtabSciMLSensitivityExtension\adjoint.jl:219
[25] _grad_adjoint_cond!(grad::Vector{…}, xdynamic::Vector{…}, xnoise::Vector{…}, xobservable::Vector{…}, xnondynamic::Vector{…}, icid::Int64, probinfo::PEtab.PEtabODEProblemInfo{…}, model_info::PEtab.ModelInfo, vjp_ss_cid::typeof(identity))
@ PEtabSciMLSensitivityExtension C:\Users\memmo\.julia\packages\PEtab\MalGN\ext\PEtabSciMLSensitivityExtension\adjoint.jl:174
[26] _grad_adjoint_xdynamic!(grad::Vector{…}, probinfo::PEtab.PEtabODEProblemInfo{…}, model_info::PEtab.ModelInfo; cids::Vector{…})
@ PEtabSciMLSensitivityExtension C:\Users\memmo\.julia\packages\PEtab\MalGN\ext\PEtabSciMLSensitivityExtension\adjoint.jl:63
[27] _grad_adjoint_xdynamic!
@ C:\Users\memmo\.julia\packages\PEtab\MalGN\ext\PEtabSciMLSensitivityExtension\adjoint.jl:26 [inlined]
[28] grad_adjoint!(grad::Vector{…}, x::Vector{…}, _nllh_not_solveode!::Function, probinfo::PEtab.PEtabODEProblemInfo{…}, model_info::PEtab.ModelInfo; cids::Vector{…})
@ PEtabSciMLSensitivityExtension C:\Users\memmo\.julia\packages\PEtab\MalGN\ext\PEtabSciMLSensitivityExtension\adjoint.jl:10
[29] grad_adjoint!
@ C:\Users\memmo\.julia\packages\PEtab\MalGN\ext\PEtabSciMLSensitivityExtension\adjoint.jl:2 [inlined]
[30] #12
@ C:\Users\memmo\.julia\packages\PEtab\MalGN\ext\PEtabSciMLSensitivityExtension\create_adjoint.jl:13 [inlined]
[31] (::PEtabSciMLSensitivityExtension.var"#13#18"{…})(g::ComponentArrays.ComponentVector{…}, x::ComponentArrays.ComponentVector{…}; prior::Bool)
@ PEtabSciMLSensitivityExtension C:\Users\memmo\.julia\packages\PEtab\MalGN\ext\PEtabSciMLSensitivityExtension\create_adjoint.jl:20
[32] (::PEtabSciMLSensitivityExtension.var"#13#18"{…})(g::ComponentArrays.ComponentVector{…}, x::ComponentArrays.ComponentVector{…})
@ PEtabSciMLSensitivityExtension C:\Users\memmo\.julia\packages\PEtab\MalGN\ext\PEtabSciMLSensitivityExtension\create_adjoint.jl:17
[33] top-level scope
I was not getting the following error in the previous PEtab version