Re: LU factorization for sparse matrices

[Bill Allombert <Bill.Allombert@math.u-bordeaux1.fr>]

On Thu, Apr 19, 2012 at 12:01:48PM +0200, Sara Costa Romero wrote:
> Hi
> 
> I have some huge non-symmetric sparse systems I have to solve. Is there any
> LU factorization of sparse matrices or I have to consider them as dense? If
> there is not, can anyone give me an idea to minimize their fill factor
> preserving stability? It seems that if I preordering and stability do not
> match, one distrois the goods of the other.

Unfortunately PARI/GP does not provide efficient sparse linear algebra.
You shuld be better off using some other system like linbox or MUMPS.

Cheers,
Bill.