you can try openNL , there is also a commercial cula package that has some basic materials for free. there is no magic created for sparse matrix solvers, it depends a lot on the sparseness matrix template and how you want to use it, do you have a solution? matrix dose remains unchanged? In my case, we ended up writing our own solver for the reasons mentioned above.
source
share