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 ) |