Isaac Wheeler
Isaac Wheeler
Using `JET.@report_opt`, the problem construction still doesn't have instabilities, but solving does: ```julia prob = ODEProblem(cbase) @report_opt solve(prob) ``` ``` ┌ Warning: Tuple{DiffEqBase.var"##solve#42", Any, Any, Any, Any, Base.Pairs{Symbol, V, NTuple{N,...
Exploring this MWE with Cthulhu, it looks like maybe this is related to #2613 ? Inside of `solve_up`, it fails to infer the full type of the `alg`:  which...
I checked this MWE again after #2617 merged, and this still occurs, so something else seems to be in play here.
Ruling out the instability from the ForwardDiff chunksize, I think this function might be at least one culprit: https://github.com/SciML/OrdinaryDiffEq.jl/blob/d227c902cff8f1da22458c06b7a4334541efa7a6/lib/OrdinaryDiffEqCore/src/misc_utils.jl#L108-L130 which gets used in `__init`: https://github.com/SciML/OrdinaryDiffEq.jl/blob/d227c902cff8f1da22458c06b7a4334541efa7a6/lib/OrdinaryDiffEqCore/src/solve.jl#L480-L502 Building on the above MWE:...
It seems like the behavior of `get_differential_vars` is informed partly by the preference that, in the event all the variables are differential, the result should be `nothing` rather than a...
> `|| all(!iszero, mm) ` I don't think that should go to the `nothing` branch. It should just be a `BitVector`. `mm isa UniformScaling` is a good special case but...
> [@Ickaser](https://github.com/Ickaser) is your test case FullSpecialize? Yep, which I didn't deliberately choose. When I construct the ODEFunction with mass matrix: ```julia dae_fc = ODEFunction{true}(dae!, mass_matrix=Diagonal([1.0, 0.0])) ``` the resulting...
On further inspection, the `ODEFunction` constructor defaults to `FullSpecialize` in general, whether or not the mass matrix is specified, and when it has `FullSpecialize` the chunksize is not inferred. In...
An MRE for the extension-loading issue: ``` julia> using Unitful julia> using Plots julia> Base.get_extension(PlotsBase, :UnitfulExt) julia> using Latexify julia> Base.get_extension(PlotsBase, :UnitfulExt) julia> using UnitfulLatexify julia> Base.get_extension(PlotsBase, :UnitfulExt) UnitfulExt ```
This isn't fully resolved yet