version 1.18, 2018/05/29 07:18:30
|
version 1.20, 2023/08/07 08:39:32
|
Line 134
|
Line 134
|
*> \author Univ. of Colorado Denver |
*> \author Univ. of Colorado Denver |
*> \author NAG Ltd. |
*> \author NAG Ltd. |
* |
* |
*> \date June 2016 |
|
* |
|
*> \ingroup complex16OTHERauxiliary |
*> \ingroup complex16OTHERauxiliary |
* |
* |
*> \par Further Details: |
*> \par Further Details: |
Line 169
|
Line 167
|
SUBROUTINE ZLATDF( IJOB, N, Z, LDZ, RHS, RDSUM, RDSCAL, IPIV, |
SUBROUTINE ZLATDF( IJOB, N, Z, LDZ, RHS, RDSUM, RDSCAL, IPIV, |
$ JPIV ) |
$ JPIV ) |
* |
* |
* -- LAPACK auxiliary routine (version 3.7.0) -- |
* -- LAPACK auxiliary routine -- |
* -- LAPACK is a software package provided by Univ. of Tennessee, -- |
* -- LAPACK is a software package provided by Univ. of Tennessee, -- |
* -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..-- |
* -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..-- |
* June 2016 |
|
* |
* |
* .. Scalar Arguments .. |
* .. Scalar Arguments .. |
INTEGER IJOB, LDZ, N |
INTEGER IJOB, LDZ, N |
Line 261
|
Line 258
|
* |
* |
* 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 ) |