Linear Algebra Package for general sparse matrices
 

I am programming a particle based method for flow problems. It uses linked list data structure. It solves a poisson equation and since the particles can have arbitrary positions, the matrix need not be banded. So far I have been using BiCGSTAB which I programmed myself. Since the data structure is linked list, I find my code inefficient. I am considering temporarily writing to arrays to solve the poisson equations and also use a Linear Algebra package. I would like to know what Linear Algebra package may be apt. I searched through LAPACK and similar websites and they mention that they exclude general sparse matrices in their features. Considering this, please suggest me some package that would be easy to use and doesn't have restrictions. I program in C and would like to know about packages that are in C.

PETSC and TRILINOS are the most comprehensive linear algebra packages available today.

I use PETSc in my code and I have been happy with it. It has many solver options and is pretty straightforward to use.

I have used both, Petsc and Trilinos, petsc is easier to install. On the other hand I know trilinos has some features to handle dense matrices.

With respect to petsc, here there is a thread about it, in some other website
http://scicomp.stackexchange.com/que...dense-matrices
Anyway, how did you get a dense matrix if you are solving a poisson equation?

I recommend you to first try with petsc and if it is not good enough, make an extra effort and install trilinos.

regards,
I.