version 1.8, 2010/12/21 13:53:51
|
version 1.9, 2011/07/22 07:38:17
|
Line 635
|
Line 635
|
CALL ZLACPY( 'ALL', KNZ, JLEN, H( INCOL+1+J2, JCOL ), |
CALL ZLACPY( '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**H ==== |
* |
* |
CALL ZLASET( 'ALL', KZS, JLEN, ZERO, ZERO, WH, LDWH ) |
CALL ZLASET( 'ALL', KZS, JLEN, ZERO, ZERO, WH, LDWH ) |
CALL ZTRMM( 'L', 'U', 'C', 'N', KNZ, JLEN, ONE, |
CALL ZTRMM( '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**H ==== |
* |
* |
CALL ZGEMM( 'C', 'N', I2, JLEN, J2, ONE, U, LDU, |
CALL ZGEMM( '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 652
|
Line 652
|
CALL ZLACPY( 'ALL', J2, JLEN, H( INCOL+1, JCOL ), LDH, |
CALL ZLACPY( 'ALL', J2, JLEN, H( INCOL+1, JCOL ), LDH, |
$ WH( I2+1, 1 ), LDWH ) |
$ WH( I2+1, 1 ), LDWH ) |
* |
* |
* ==== Multiply by U21' ==== |
* ==== Multiply by U21**H ==== |
* |
* |
CALL ZTRMM( 'L', 'L', 'C', 'N', J2, JLEN, ONE, |
CALL ZTRMM( '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 ) |