|
[Sponsors] |
May 17, 2000, 17:58 |
Good library to solve huge linear system
|
#1 |
Guest
Posts: n/a
|
Hello everybody,
I am developing some unsteady, viscous, unstructured triangular mesh codes to simulate the filling of a cavity with plastic material. I am in need of a good, reliable, free, well documented library of functions to solve linear systems. Iīve tried to work whit a C library called LASPack (not LAPACK), and itīs well documented and free, but has some difficulties to solve systems of more than 5000 variables. Could you advise me about a good library to solve this problem. I would prefer a C library. Thank you. |
|
May 18, 2000, 06:12 |
Re: Good library to solve huge linear system
|
#2 |
Guest
Posts: n/a
|
I use aztec from sandia labs http://www.cs.sandia.gov/CRF/aztec1.html but i have not yet made experiences with more than 5000 variables
cheers |
|
May 18, 2000, 10:08 |
Re: Good library to solve huge linear system
|
#3 |
Guest
Posts: n/a
|
Try SuperLU, http://www.nersc.gov/~xiaoye/SuperLU/
|
|
May 19, 2000, 07:24 |
Re: Good library to solve huge linear system
|
#4 |
Guest
Posts: n/a
|
The PETSc library from Argonne does a good job. It has C and fortran interfaces. You can use it in a uni-processor and parallel environment.
It is designed for sparse matrizes. 5000 unknowns is no problem. My problems have about 100000 unknowns, which is still a small problem. PETSc has got GMRES, BiCGSTAB ..various preconditioners (ILU, jacobi etc), and all kinds of stuff what one needs for solving large, sparse systems of linear equations. You can also define a matrix free GMRES context, include your own preconditioners etc... I don't have the address at hand right now, but if you look for "petsc", you should easily find it. The manual is a bit short sometimes (lots of functions), but the basic usage is not too difficult and there are many examples one should have a look at. The online manual is also extensive. I guess it's similar to the atztec solver from sandia. I have the feeling that more people use petsc, though. Hope this is useful. Frank |
|
Thread Tools | Search this Thread |
Display Modes | |
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
Submarine with SimpleFoam | alex_rubel | OpenFOAM | 20 | February 25, 2011 07:34 |
New linear system solvers | chegdan | OpenFOAM Running, Solving & CFD | 11 | April 30, 2010 11:22 |
Linearized NS euqations: how to solve them?(problem with Matrix operations..) | matteoL | OpenFOAM Running, Solving & CFD | 0 | November 18, 2009 07:58 |
OpenFOAM141dev linking error on IBM AIX 52 | matthias | OpenFOAM Installation | 24 | April 28, 2008 16:49 |
Vibrating system analysis - a good introduction? | Bren | Main CFD Forum | 0 | November 20, 2007 13:40 |