rpl/lapack/lapack/dlaqr5.f - diff

Return to dlaqr5.f CVS log

Up to [local] / rpl / lapack / lapack

version 1.8, 2010/12/21 13:53:31	version 1.9, 2011/07/22 07:38:07
Line 642	Line 642
CALL DLACPY( 'ALL', KNZ, JLEN, H( INCOL+1+J2, JCOL ),	CALL DLACPY( 'ALL', KNZ, JLEN, H( INCOL+1+J2, JCOL ),
$ LDH, WH( KZS+1, 1 ), LDWH )	$ LDH, WH( KZS+1, 1 ), LDWH )
*	*
* ==== Multiply by U21' ====	* ==== Multiply by U21**T ====
*	*
CALL DLASET( 'ALL', KZS, JLEN, ZERO, ZERO, WH, LDWH )	CALL DLASET( 'ALL', KZS, JLEN, ZERO, ZERO, WH, LDWH )
CALL DTRMM( 'L', 'U', 'C', 'N', KNZ, JLEN, ONE,	CALL DTRMM( 'L', 'U', 'C', 'N', KNZ, JLEN, ONE,
$ U( J2+1, 1+KZS ), LDU, WH( KZS+1, 1 ),	$ U( J2+1, 1+KZS ), LDU, WH( KZS+1, 1 ),
$ LDWH )	$ LDWH )
*	*
* ==== Multiply top of H by U11' ====	* ==== Multiply top of H by U11**T ====
*	*
CALL DGEMM( 'C', 'N', I2, JLEN, J2, ONE, U, LDU,	CALL DGEMM( 'C', 'N', I2, JLEN, J2, ONE, U, LDU,
$ H( INCOL+1, JCOL ), LDH, ONE, WH, LDWH )	$ H( INCOL+1, JCOL ), LDH, ONE, WH, LDWH )
Line 659	Line 659
CALL DLACPY( 'ALL', J2, JLEN, H( INCOL+1, JCOL ), LDH,	CALL DLACPY( 'ALL', J2, JLEN, H( INCOL+1, JCOL ), LDH,
$ WH( I2+1, 1 ), LDWH )	$ WH( I2+1, 1 ), LDWH )
*	*
* ==== Multiply by U21' ====	* ==== Multiply by U21**T ====
*	*
CALL DTRMM( 'L', 'L', 'C', 'N', J2, JLEN, ONE,	CALL DTRMM( 'L', 'L', 'C', 'N', J2, JLEN, ONE,
$ U( 1, I2+1 ), LDU, WH( I2+1, 1 ), LDWH )	$ U( 1, I2+1 ), LDU, WH( I2+1, 1 ), LDWH )