I am implementing a code for solving a DAE system using sundials and IDA (through the DifferentialEquation.jl interface). However, as my simulation requires the solution of several equations > 500. It is taking a long time for ending each simulation. So, I am wondering, if it is possible to use the native sundials capabilities to achieve a parallel solution. It is possible from the Julia sundials or DifferentialEquation.jl wrappers?
Hi! For now I am getting some parallel capabilities after using linear_solver=:LapackDense. However, I am wondering if there is a way to increase the number of threads. Because, I have the impression that the solver didn’t use all the computer threads. So there is some kind of configuration key that I can use through the DifferentialEquation.jl wrapper?