version 1.4, 2010/08/07 13:22:09
|
version 1.7, 2011/07/22 07:38:02
|
Line 13
|
Line 13
|
* |
* |
* DSYR2K performs one of the symmetric rank 2k operations |
* DSYR2K performs one of the symmetric rank 2k operations |
* |
* |
* C := alpha*A*B' + alpha*B*A' + beta*C, |
* C := alpha*A*B**T + alpha*B*A**T + beta*C, |
* |
* |
* or |
* or |
* |
* |
* C := alpha*A'*B + alpha*B'*A + beta*C, |
* C := alpha*A**T*B + alpha*B**T*A + beta*C, |
* |
* |
* where alpha and beta are scalars, C is an n by n symmetric matrix |
* where alpha and beta are scalars, C is an n by n symmetric matrix |
* and A and B are n by k matrices in the first case and k by n |
* and A and B are n by k matrices in the first case and k by n |
Line 43
|
Line 43
|
* On entry, TRANS specifies the operation to be performed as |
* On entry, TRANS specifies the operation to be performed as |
* follows: |
* follows: |
* |
* |
* TRANS = 'N' or 'n' C := alpha*A*B' + alpha*B*A' + |
* TRANS = 'N' or 'n' C := alpha*A*B**T + alpha*B*A**T + |
* beta*C. |
* beta*C. |
* |
* |
* TRANS = 'T' or 't' C := alpha*A'*B + alpha*B'*A + |
* TRANS = 'T' or 't' C := alpha*A**T*B + alpha*B**T*A + |
* beta*C. |
* beta*C. |
* |
* |
* TRANS = 'C' or 'c' C := alpha*A'*B + alpha*B'*A + |
* TRANS = 'C' or 'c' C := alpha*A**T*B + alpha*B**T*A + |
* beta*C. |
* beta*C. |
* |
* |
* Unchanged on exit. |
* Unchanged on exit. |
Line 234
|
Line 234
|
* |
* |
IF (LSAME(TRANS,'N')) THEN |
IF (LSAME(TRANS,'N')) THEN |
* |
* |
* Form C := alpha*A*B' + alpha*B*A' + C. |
* Form C := alpha*A*B**T + alpha*B*A**T + C. |
* |
* |
IF (UPPER) THEN |
IF (UPPER) THEN |
DO 130 J = 1,N |
DO 130 J = 1,N |
Line 283
|
Line 283
|
END IF |
END IF |
ELSE |
ELSE |
* |
* |
* Form C := alpha*A'*B + alpha*B'*A + C. |
* Form C := alpha*A**T*B + alpha*B**T*A + C. |
* |
* |
IF (UPPER) THEN |
IF (UPPER) THEN |
DO 210 J = 1,N |
DO 210 J = 1,N |