[Top][All Lists]

[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]

Re: [Help-gsl] gsl_blas_dtrsm

From: Brian Gough
Subject: Re: [Help-gsl] gsl_blas_dtrsm
Date: Tue, 30 May 2006 12:18:31 +0100

Remi Vieux writes:
 > Does the matrix A has to be triangular to use this function? What
 > exactly is TransA for? Same question for the parameter Diag: what do you
 > mean by "the diagonal of A is used or taken as unity"?


The BLAS tr functions only work with triangular matrices.

I suggest taking an LU decomposition of A and solving in place
on the columns of B. e.g.


  for (j = 0; j < B->size2 ;j++) {
    gsl_vector_view x = gsl_matrix_column (B, j);
    gsl_linalg_LU_svx (A, p, &x.vector);

See the manual for details of the functions involved.

Brian Gough

Network Theory Ltd,
Publishing the GSL Manual -

reply via email to

[Prev in Thread] Current Thread [Next in Thread]