LAPACK/ScaLAPACK Development

Hallo,
for me it seems to be too complicated to solve X*A=B, like A'*X'=B' with DTRTRS. Moreover A is triangular and with B are parts of one matrix. How to solve this "standardized" problem efficient?

Hi,
Did you try to call the BLAS routine: DTRSM directly?

Code: Select all

*  DTRSM  solves one of the matrix equations
*
*     op( A )*X = alpha*B,   or   X*op( A ) = alpha*B,
*
*  where alpha is a scalar, X and B are m by n matrices, A is a unit, or
*  non-unit,  upper or lower triangular matrix  and  op( A )  is one  of
*
*     op( A ) = A   or   op( A ) = A**T.
*
*  The matrix X is overwritten on B.


Julie

See: viewtopic.php?f=2&t=2675