Yingbo Ma
Yingbo Ma
Yeah, this solution sounds good.
In MTK, we had to add an explicit `dropzero!` because of this change https://github.com/SciML/ModelingToolkit.jl/blob/05525b540820569c01d2249096dbbc5880d943d0/src/structural_transformation/bareiss.jl#L146-L148 ```julia julia> using ModelingToolkit, BenchmarkTools, SparseArrays julia> A = BigInt.(sprand(Int, 100, 100, 0.01)) + I; julia>...
Fine Dulmage-Mendelsohn decompositions are just another name for BLT or SCC.
```julia @named fol_separate = ODESystem([ RHS ~ (1 - x)/τ, D(x) ~ RHS ]) prob = ODEProblem(structural_simplify(fol_separate), [x=>0.0, RHS=>0.0], (0.0,10.0), [τ => 3.0]) ```
Is there an MWE?
Yeah. I think it's fine. It won't be breaking either.
I am not sure about the cause, and I don't have a Windows machine. Does it crash Julia every time when you run it?
You can make a PR if you want this feature.
Should we fix this inaccuracy by inserting a sleep call in the dynamic broadcasting branch?
Yes, Shashi and I are looking into IRTools.jl and source-to-source transformation. But before that, we want to solve #33 first to make ForwardDiff2 functional.