|
|
@@ -804,7 +804,7 @@ |
|
|
CALL DLARTGP( B22BULGE, B22E(I-1), RWORK(IU2SN+I-1), |
|
|
CALL DLARTGP( B22BULGE, B22E(I-1), RWORK(IU2SN+I-1), |
|
|
$ RWORK(IU2CS+I-1), R ) |
|
|
$ RWORK(IU2CS+I-1), R ) |
|
|
ELSE IF( NU .LT. MU ) THEN |
|
|
ELSE IF( NU .LT. MU ) THEN |
|
|
CALL DLARTGS( B21E(I), B21E(I+1), NU, RWORK(IU2CS+I-1), |
|
|
|
|
|
|
|
|
CALL DLARTGS( B21E(I), B21D(I+1), NU, RWORK(IU2CS+I-1), |
|
|
$ RWORK(IU2SN+I-1) ) |
|
|
$ RWORK(IU2SN+I-1) ) |
|
|
ELSE |
|
|
ELSE |
|
|
CALL DLARTGS( B22D(I), B22E(I), MU, RWORK(IU2CS+I-1), |
|
|
CALL DLARTGS( B22D(I), B22E(I), MU, RWORK(IU2CS+I-1), |
|
|
|