Re: [Help-gsl] solving linear system with large matrices

[Jordi Gutierrez Hermoso]

On 07/03/07, Jan Ehrhardt <address@hidden> wrote:
> I want to solve a linear system Ax=b, where A is symmetric, positive
> definite square matrix.  However, the size of the matrix is very
> large: approx. 65000x65000 or larger so an allocation of the whole
> matrix will fail.

> From chapter 13 of the GSL manual:

    The library provides simple linear algebra operations which
    operate directly on the gsl_vector and gsl_matrix objects. These
    are intended for use with "small" systems where simple algorithms
    are acceptable.

    Anyone interested in large systems will want to use the
    sophisticated routines found in LAPACK.

A convenient interface to LAPACK is found in GNU Octave, if you prefer
to avoid coding in Fortran, as I do. ;-)

HTH,
- Jordi G. H.