Linear Algebra Package for general sparse matrices

prapanj · December 18, 2012, 01:42

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.

praveen · December 18, 2012, 03:45

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

cdegroot · December 18, 2012, 08:25

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

ignatius · January 2, 2013, 07:58

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.

December 18, 2012, 01:42	Linear Algebra Package for general sparse matrices	#1
prapanj Senior Member Prapanch Nair Join Date: Mar 2009 Location: Bangalore, India Posts: 105 Rep Power: 17	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.

Similar Threads
Thread	Thread Starter	Forum	Replies	Last Post
A turbulent test case for rhoCentralFoam	immortality	OpenFOAM Running, Solving & CFD	13	April 20, 2014 07:32
Extend Project 1.6-ext Binary Release for Mac OS X	hjasak	OpenFOAM Announcements from Other Sources	26	November 5, 2013 17:50
Problem? with Installation of OpenFOAM 2.1.0 on to Ubuntu 10.4.4 VirtualBox	Stubby	OpenFOAM Installation	8	March 5, 2012 16:12
solution diverges when linear upwind interpolation scheme is used	subash	OpenFOAM	0	May 29, 2010 02:23
Errors running allwmake in OpenFOAM141dev with WM_COMPILE_OPTION%3ddebug	unoder	OpenFOAM Installation	11	January 30, 2008 21:30

December 18, 2012, 03:45		#2
praveen Super Moderator Praveen. C Join Date: Mar 2009 Location: Bangalore Posts: 343 Blog Entries: 6 Rep Power: 18	PETSC and TRILINOS are the most comprehensive linear algebra packages available today.

December 18, 2012, 08:25		#3
cdegroot Senior Member Chris DeGroot Join Date: Nov 2011 Location: Canada Posts: 414 Rep Power: 18	I use PETSc in my code and I have been happy with it. It has many solver options and is pretty straightforward to use.

January 2, 2013, 07:58		#4
ignatius New Member Join Date: Jan 2013 Posts: 3 Rep Power: 13	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.