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, B is an m by n matrix, A is a unit, or |
* where alpha is a scalar, B is an m by n matrix, 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. |
* |
* |
* Arguments |
* Arguments |
* ========== |
* ========== |
Line 49
|
Line 49
|
* |
* |
* 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 234
|
Line 234
|
END IF |
END IF |
ELSE |
ELSE |
* |
* |
* Form B := alpha*A'*B. |
* Form B := alpha*A**T*B. |
* |
* |
IF (UPPER) THEN |
IF (UPPER) THEN |
DO 110 J = 1,N |
DO 110 J = 1,N |
Line 300
|
Line 300
|
END IF |
END IF |
ELSE |
ELSE |
* |
* |
* Form B := alpha*B*A'. |
* Form B := alpha*B*A**T. |
* |
* |
IF (UPPER) THEN |
IF (UPPER) THEN |
DO 260 K = 1,N |
DO 260 K = 1,N |