version 1.3, 2010/08/06 15:32:22
|
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 |