rpl/lapack/blas/dsyrk.f - diff

Return to dsyrk.f CVS log

Up to [local] / rpl / lapack / blas

version 1.3, 2010/08/06 15:32:20	version 1.7, 2011/07/22 07:38:02
Line 13	Line 13
*	*
* DSYRK performs one of the symmetric rank k operations	* DSYRK performs one of the symmetric rank k operations
*	*
* C := alphaAA' + beta*C,	* C := alphaAA*T + betaC,
*	*
* or	* or
*	*
* C := alphaA'A + beta*C,	* C := alphaATA + 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 is an n by k matrix in the first case and a k by n matrix	* and A is an n by k matrix in the first case and a k by n matrix
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 := alphaAA' + beta*C.	* TRANS = 'N' or 'n' C := alphaAA*T + betaC.
*	*
* TRANS = 'T' or 't' C := alphaA'A + beta*C.	* TRANS = 'T' or 't' C := alphaATA + beta*C.
*	*
* TRANS = 'C' or 'c' C := alphaA'A + beta*C.	* TRANS = 'C' or 'c' C := alphaATA + beta*C.
*	*
* Unchanged on exit.	* Unchanged on exit.
*	*
Line 213	Line 213
*	*
IF (LSAME(TRANS,'N')) THEN	IF (LSAME(TRANS,'N')) THEN
*	*
* Form C := alphaAA' + beta*C.	* Form C := alphaAA*T + betaC.
*	*
IF (UPPER) THEN	IF (UPPER) THEN
DO 130 J = 1,N	DO 130 J = 1,N
Line 258	Line 258
END IF	END IF
ELSE	ELSE
*	*
* Form C := alphaA'A + beta*C.	* Form C := alphaATA + beta*C.
*	*
IF (UPPER) THEN	IF (UPPER) THEN
DO 210 J = 1,N	DO 210 J = 1,N