Open avik-pal opened 10 months ago
From the error message, it's probably trying to construct the u0
from the unsimplified system by mistake.
julia> s = symbolic_discretize(pdesys, discretization)[1];
JuliaSimCompiler: generate_system
julia> ss = structural_simplify(s);
julia> NonlinearProblem(structural_simplify(s), rand(length(states(ss))))
NonlinearProblem with uType Vector{Float64}. In-place: true
u0: 12-element Vector{Float64}:
0.5658693983716057
0.730452210930952
0.7999093647009754
0.34068068695778175
0.6066626175654133
0.49832174589693123
0.985735570069928
0.10313229637404808
0.5761788103850009
0.019531876564822248
0.4600618070405691
0.5048377501152175
julia> length(states(s))
121
Afraid that this isn't the case:
simpsys = structural_simplify(sys)
if tspan === nothing
add_metadata!(get_metadata(sys), sys)
return prob = NonlinearProblem(simpsys, ones(length(simpsys.states));
discretization.kwargs..., kwargs...)
Digging a bit further 8.71.2
for Modeling Toolkit works
I thought you said you were on latest?
I meant 8.71.2
works, but the latest 8.72
doesn't
https://docs.sciml.ai/MethodOfLines/dev/tutorials/heatss/ seems to be broken with the latest versions of the packages.
The last working versions we could figure out were:
cc @ChrisRackauckas