version 1.8, 2011/11/21 20:43:04
|
version 1.14, 2016/08/27 15:34:39
|
Line 166
|
Line 166
|
*> \author Univ. of Colorado Denver |
*> \author Univ. of Colorado Denver |
*> \author NAG Ltd. |
*> \author NAG Ltd. |
* |
* |
*> \date November 2011 |
*> \date November 2015 |
* |
* |
*> \ingroup doubleOTHERcomputational |
*> \ingroup doubleOTHERcomputational |
* |
* |
Line 174
|
Line 174
|
SUBROUTINE DSTEIN( N, D, E, M, W, IBLOCK, ISPLIT, Z, LDZ, WORK, |
SUBROUTINE DSTEIN( N, D, E, M, W, IBLOCK, ISPLIT, Z, LDZ, WORK, |
$ IWORK, IFAIL, INFO ) |
$ IWORK, IFAIL, INFO ) |
* |
* |
* -- LAPACK computational routine (version 3.4.0) -- |
* -- LAPACK computational routine (version 3.6.0) -- |
* -- 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..-- |
* November 2011 |
* November 2015 |
* |
* |
* .. Scalar Arguments .. |
* .. Scalar Arguments .. |
INTEGER INFO, LDZ, M, N |
INTEGER INFO, LDZ, M, N |
Line 297
|
Line 297
|
BLKSIZ = BN - B1 + 1 |
BLKSIZ = BN - B1 + 1 |
IF( BLKSIZ.EQ.1 ) |
IF( BLKSIZ.EQ.1 ) |
$ GO TO 60 |
$ GO TO 60 |
GPIND = B1 |
GPIND = J1 |
* |
* |
* Compute reorthogonalization criterion and stopping criterion. |
* Compute reorthogonalization criterion and stopping criterion. |
* |
* |
Line 370
|
Line 370
|
* |
* |
* Normalize and scale the righthand side vector Pb. |
* Normalize and scale the righthand side vector Pb. |
* |
* |
|
JMAX = IDAMAX( BLKSIZ, WORK( INDRV1+1 ), 1 ) |
SCL = BLKSIZ*ONENRM*MAX( EPS, |
SCL = BLKSIZ*ONENRM*MAX( EPS, |
$ ABS( WORK( INDRV4+BLKSIZ ) ) ) / |
$ ABS( WORK( INDRV4+BLKSIZ ) ) ) / |
$ DASUM( BLKSIZ, WORK( INDRV1+1 ), 1 ) |
$ ABS( WORK( INDRV1+JMAX ) ) |
CALL DSCAL( BLKSIZ, SCL, WORK( INDRV1+1 ), 1 ) |
CALL DSCAL( BLKSIZ, SCL, WORK( INDRV1+1 ), 1 ) |
* |
* |
* Solve the system LU = Pb. |
* Solve the system LU = Pb. |