CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Main CFD Forum

Linear Algebra Package for general sparse matrices

Register Blogs Members List Search Today's Posts Mark Forums Read

Reply
 
LinkBack Thread Tools Display Modes
Old   December 18, 2012, 01:42
Question Linear Algebra Package for general sparse matrices
  #1
Senior Member
 
Prapanch Nair
Join Date: Mar 2009
Location: Bangalore, India
Posts: 105
Rep Power: 8
prapanj is on a distinguished road
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.
prapanj is offline   Reply With Quote

Old   December 18, 2012, 03:45
Default
  #2
Super Moderator
 
praveen's Avatar
 
Praveen. C
Join Date: Mar 2009
Location: Bangalore
Posts: 251
Blog Entries: 6
Rep Power: 9
praveen is on a distinguished road
PETSC and TRILINOS are the most comprehensive linear algebra packages available today.
praveen is offline   Reply With Quote

Old   December 18, 2012, 08:25
Default
  #3
Senior Member
 
Chris DeGroot
Join Date: Nov 2011
Location: Canada
Posts: 387
Rep Power: 6
cdegroot is on a distinguished road
I use PETSc in my code and I have been happy with it. It has many solver options and is pretty straightforward to use.
cdegroot is offline   Reply With Quote

Old   January 2, 2013, 07:58
Default
  #4
New Member
 
Join Date: Jan 2013
Posts: 3
Rep Power: 4
ignatius is on a distinguished road
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.
ignatius is offline   Reply With Quote

Reply

Tags
bicgstab, lapack, linear algebra, poisson equation, solver

Thread Tools
Display Modes

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are On
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
A turbulent test case for rhoCentralFoam immortality OpenFOAM Running, Solving & CFD 13 April 20, 2014 06: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 01:23
Errors running allwmake in OpenFOAM141dev with WM_COMPILE_OPTION%3ddebug unoder OpenFOAM Installation 11 January 30, 2008 21:30


All times are GMT -4. The time now is 21:06.