Diff for /rpl/lapack/blas/ztrsm.f between versions 1.6 and 1.7

version 1.6, 2010/12/21 13:51:27 version 1.7, 2011/07/22 07:38:03
Line 18 Line 18
 *  where alpha is a scalar, X and B are m by n matrices, A is a unit, or  *  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  *  non-unit,  upper or lower triangular matrix  and  op( A )  is one  of
 *  *
 *     op( A ) = A   or   op( A ) = A'   or   op( A ) = conjg( A' ).  *     op( A ) = A   or   op( A ) = A**T   or   op( A ) = A**H.
 *  *
 *  The matrix X is overwritten on B.  *  The matrix X is overwritten on B.
 *  *
Line 51 Line 51
 *  *
 *              TRANSA = 'N' or 'n'   op( A ) = A.  *              TRANSA = 'N' or 'n'   op( A ) = A.
 *  *
 *              TRANSA = 'T' or 't'   op( A ) = A'.  *              TRANSA = 'T' or 't'   op( A ) = A**T.
 *  *
 *              TRANSA = 'C' or 'c'   op( A ) = conjg( A' ).  *              TRANSA = 'C' or 'c'   op( A ) = A**H.
 *  *
 *           Unchanged on exit.  *           Unchanged on exit.
 *  *
Line 245 Line 245
               END IF                END IF
           ELSE            ELSE
 *  *
 *           Form  B := alpha*inv( A' )*B  *           Form  B := alpha*inv( A**T )*B
 *           or    B := alpha*inv( conjg( A' ) )*B.  *           or    B := alpha*inv( A**H )*B.
 *  *
               IF (UPPER) THEN                IF (UPPER) THEN
                   DO 140 J = 1,N                    DO 140 J = 1,N
Line 336 Line 336
               END IF                END IF
           ELSE            ELSE
 *  *
 *           Form  B := alpha*B*inv( A' )  *           Form  B := alpha*B*inv( A**T )
 *           or    B := alpha*B*inv( conjg( A' ) ).  *           or    B := alpha*B*inv( A**H ).
 *  *
               IF (UPPER) THEN                IF (UPPER) THEN
                   DO 330 K = N,1,-1                    DO 330 K = N,1,-1

Removed from v.1.6  
changed lines
  Added in v.1.7


CVSweb interface <joel.bertrand@systella.fr>