From 2f5de047b6438de0f29052dcbc343343b291c9c7 Mon Sep 17 00:00:00 2001 From: langou Date: Mon, 31 Jan 2011 18:27:33 +0000 Subject: [PATCH] Change from Christof (Voemel) sent on Mon, 31 Jan 2011. ********************************************************************* change IF( USEDQD ) THEN * The initial SIGMA was to the outer end of the spectrum * the matrix is definite and we need not retreat. TAU = SPDIAM*EPS*N + TWO*PIVMIN ELSE into IF( USEDQD ) THEN * The initial SIGMA was to the outer end of the spectrum * the matrix is definite and we need not retreat. TAU = SPDIAM*EPS*N + TWO*PIVMIN TAU = MAX( TAU,TWO*EPS*ABS(SIGMA) ) ELSE that is, with an added line to prevent TAU from becoming zero TAU = MAX( TAU,TWO*EPS*ABS(SIGMA) ) ********************************************************************* Thanks Christof! --- SRC/dlarre.f | 1 + SRC/slarre.f | 1 + 2 files changed, 2 insertions(+) diff --git a/SRC/dlarre.f b/SRC/dlarre.f index f08043d..9e483f2 100644 --- a/SRC/dlarre.f +++ b/SRC/dlarre.f @@ -532,6 +532,7 @@ * The initial SIGMA was to the outer end of the spectrum * the matrix is definite and we need not retreat. TAU = SPDIAM*EPS*N + TWO*PIVMIN + TAU = MAX( TAU,TWO*EPS*ABS(SIGMA) ) ELSE IF(MB.GT.1) THEN CLWDTH = W(WEND) + WERR(WEND) - W(WBEGIN) - WERR(WBEGIN) diff --git a/SRC/slarre.f b/SRC/slarre.f index b2a439b..b3815b9 100644 --- a/SRC/slarre.f +++ b/SRC/slarre.f @@ -536,6 +536,7 @@ * The initial SIGMA was to the outer end of the spectrum * the matrix is definite and we need not retreat. TAU = SPDIAM*EPS*N + TWO*PIVMIN + TAU = MAX( TAU,TWO*EPS*ABS(SIGMA) ) ELSE IF(MB.GT.1) THEN CLWDTH = W(WEND) + WERR(WEND) - W(WBEGIN) - WERR(WBEGIN) -- 2.7.4