Yingbo Ma

I misunderstood the question as well. So you want to find conservation laws from ODEs and convert them into DAEs, so in this case, naturally, the algebraic equations take precedence....

I think by using the matching algorithm, one can find what differential equations to replace with added algebraic equations. I am not sure if it is fully correct. I have...

I cannot reproduce the stack overflow error.

If we are gonna make a breaking change, maybe we should change `max_iter` to `maxiters` to be consistant? > Maybe it doesn't have any impact on performance and e.g. the...

Hey, could you also provide the code that generates this model? Also, I think this should be fixed if you update your packages.

MTK tries to transform the system so that it is an almost identity matrix, i.e. the semi-explicit form. I tried it on the latest MTK, but initialization of the DAE...

The printing now looks like ```julia julia> constructTsitouras5() Tsitouras5{explicit, FSAL, stages=7, order=5, embedded order=4}: A: ⋅ ⋅ ⋅ ⋅ ⋅ ⋅ ⋅ 0.161 ⋅ ⋅ ⋅ ⋅ ⋅ ⋅ -0.008480655492356989...

With ```julia @generated function constructRadauIIA5(::Type{T} = Float64) where T sq6 = sqrt(convert(T, 6)) A = @SArray [11//45-7sq6/360 37//225-169sq6/1800 -2//225+sq6/75 37//225+169sq6/1800 11//45+7sq6/360 -2//225-sq6/75 4//9-sq6/36 4//9+sq6/36 1//9] c = @SArray [2//5-sq6/10;2/5+sq6/10;1] α...

Because it automatically detects the explicitness and FSAL on construction and puts the result in the type. We will need `@generated` either way when we use `BigFloat` inside.

The problem is caused by ```julia julia> dt = 4e-4; julia> sdt = 0.0; for i in 1:25 sdt += dt (0:dt:1)[i+1]