Hello and welcome to the OpenModelica forum,

What kind of linear solver are you planning to integrate?

We use linear solvers for several purposes, therefore we have different modules at different places.

On one hand there is the implementation of SUNDIALS/IDA, which has it's own interface to integrate
different linear solvers. Here the user can choose between internal linear solver, LAPACK, KLU,
SUPERLU and some more iterative solvers. Here is OpenModelica the user.
Currently we interfacing from IDA: KLU and LAPACK .

On the other Hand we need to solve, so called linear "algebraic loops" to obtain the RHS for the integrator.
Here we are also interfacing different external solver e.g. KLU, UMFPACK, LAPACK, but we have also
some own implementations like the totalpivot solver.

As far as I understand you want to integrate your linear solver into the IDA integrator, then it's enough
to add it under 3rdParty/sundials/idas and interfacing it in SimulationRuntime/c/simulation/ida_solver.c,
that you do not need to use "linearSystem" module.

Do you have your code on github? Then you could probably start a Pull-Request and we can further see
what you have done how it would fit into the project.

so long.
WIlli

Hello

thanks for the quick answer.

I am currently working on implementing one (or more) CPU-parallelized solver into IDA to see if I can gain any speedups in circuit simulation. However, I am still working on getting test models large enough to receive speedups over KLU or LAPACK. I would first like to test the security and robustness of my implementation before I upload anything to GitHub.

Best Regards,
Lennart