Re: [Help-gsl] gsl_linalg_LU_decomp Function

[Jordi Gutiérrez Hermoso]

On 07/08/07, Matthew Boulton <address@hidden> wrote:
>    gsl_permutation *indx = gsl_permutation_calloc( ncomp_total + 1 ); /* Has

Change this line to

     gsl_permutation *indx = gsl_permutation_calloc( mat_size+1 );

but it looks like you're trying to do an LU permutation of a
non-square matrix. Unfortunately, the current GSL implementation does
not allow this.

- Jordi G. H.