rschiemann1

Hi, I appreciate the offered support in using MethodOfLines.jl for automatic discretization of PDEs. Infact, this is part of my use case and so I definitely see myself coming back...

Bump. If there is any interest from the main MTK developers to realize such a documentation, I am willing to help.

Thanks. Okay, so how should we get this started?

> @rschiemann1 I have submitted a model DAE/ODE Problem specified in MTK with domains and BCs #50 to MethodOfLines as a prototype for a tutorial I just checked your model...

Then how does everyone else specify their initial values for array variables? I don't think I am the first one to try this.

Thank you two for looking into this. Just for the record, this is (in my view) the same issue as in #1501.

Oh. Of course I have almost no insight into the inner workings here, but @ValentinKaisermayer reported a running code with `ODEProblem(sys, u0map, tspan, param_values)` and reported the same error as...

> ```julia > function build_problem(sys; tspan = (0.0, 1.5), param_values = [100e2, 200, 0.5]) > u0map = [ > (collect(sys.Tgas) .=> fill(100.0, N))... > (collect(sys.Tsol) .=> fill(20.0, N))... > (collect(sys.cR)...

Okay, thank you for the clarification. The workaround with collect errors on the solve on my side as well - both with and without jac and sparse in ODEProblem construction....

With yesterday's release of Symbolics (closing this [issue](https://github.com/JuliaSymbolics/Symbolics.jl/issues/567)), constructing the ODEProblem with jac=true, sparse = true works, which is one step closer. Unfortunately, solving the constructed ODEProblem fails with a...