version 1.6, 2010/12/21 13:51:25
|
version 1.7, 2011/07/22 07:38:02
|
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'. |
* op( A ) = A or op( A ) = A**T. |
* |
* |
* 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 ) = A'. |
* TRANSA = 'C' or 'c' op( A ) = A**T. |
* |
* |
* Unchanged on exit. |
* Unchanged on exit. |
* |
* |
Line 243
|
Line 243
|
END IF |
END IF |
ELSE |
ELSE |
* |
* |
* Form B := alpha*inv( A' )*B. |
* Form B := alpha*inv( A**T )*B. |
* |
* |
IF (UPPER) THEN |
IF (UPPER) THEN |
DO 130 J = 1,N |
DO 130 J = 1,N |
Line 319
|
Line 319
|
END IF |
END IF |
ELSE |
ELSE |
* |
* |
* Form B := alpha*B*inv( A' ). |
* Form B := alpha*B*inv( A**T ). |
* |
* |
IF (UPPER) THEN |
IF (UPPER) THEN |
DO 310 K = N,1,-1 |
DO 310 K = N,1,-1 |