version 1.18, 2018/05/29 07:18:30
|
version 1.19, 2020/05/21 21:46:09
|
Line 261
|
Line 261
|
* |
* |
* Solve for U- part, lockahead for RHS(N) = +-1. This is not done |
* Solve for U- part, lockahead for RHS(N) = +-1. This is not done |
* In BSOLVE and will hopefully give us a better estimate because |
* In BSOLVE and will hopefully give us a better estimate because |
* any ill-conditioning of the original matrix is transfered to U |
* any ill-conditioning of the original matrix is transferred to U |
* and not to L. U(N, N) is an approximation to sigma_min(LU). |
* and not to L. U(N, N) is an approximation to sigma_min(LU). |
* |
* |
CALL ZCOPY( N-1, RHS, 1, WORK, 1 ) |
CALL ZCOPY( N-1, RHS, 1, WORK, 1 ) |