version 1.19, 2018/05/29 07:18:01
|
version 1.20, 2020/05/21 21:46:01
|
Line 85
|
Line 85
|
*> RHS is DOUBLE PRECISION array, dimension (N) |
*> RHS is DOUBLE PRECISION array, dimension (N) |
*> On entry, RHS contains contributions from other subsystems. |
*> On entry, RHS contains contributions from other subsystems. |
*> On exit, RHS contains the solution of the subsystem with |
*> On exit, RHS contains the solution of the subsystem with |
*> entries acoording to the value of IJOB (see above). |
*> entries according to the value of IJOB (see above). |
*> \endverbatim |
*> \endverbatim |
*> |
*> |
*> \param[in,out] RDSUM |
*> \param[in,out] RDSUM |
Line 260
|
Line 260
|
* |
* |
* Solve for U-part, look-ahead for RHS(N) = +-1. This is not done |
* Solve for U-part, look-ahead 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 DCOPY( N-1, RHS, 1, XP, 1 ) |
CALL DCOPY( N-1, RHS, 1, XP, 1 ) |