diff options
| -rw-r--r-- | SRC/dbdsqr.f | 4 | ||||
| -rw-r--r-- | SRC/sbdsqr.f | 4 |
2 files changed, 4 insertions, 4 deletions
diff --git a/SRC/dbdsqr.f b/SRC/dbdsqr.f index 0e46e71c..d3630f34 100644 --- a/SRC/dbdsqr.f +++ b/SRC/dbdsqr.f @@ -411,12 +411,12 @@ 40 CONTINUE 50 CONTINUE SMINOA = SMINOA / SQRT( DBLE( N ) ) - THRESH = MAX( TOL*SMINOA, MAXITR*N*N*UNFL ) + THRESH = MAX( TOL*SMINOA, MAXITR*(N*(N*UNFL)) ) ELSE * * Absolute accuracy desired * - THRESH = MAX( ABS( TOL )*SMAX, MAXITR*N*N*UNFL ) + THRESH = MAX( ABS( TOL )*SMAX, MAXITR*(N*(N*UNFL)) ) END IF * * Prepare for main iteration loop for the singular values diff --git a/SRC/sbdsqr.f b/SRC/sbdsqr.f index 14560108..38fd1a05 100644 --- a/SRC/sbdsqr.f +++ b/SRC/sbdsqr.f @@ -410,12 +410,12 @@ 40 CONTINUE 50 CONTINUE SMINOA = SMINOA / SQRT( REAL( N ) ) - THRESH = MAX( TOL*SMINOA, MAXITR*N*N*UNFL ) + THRESH = MAX( TOL*SMINOA, MAXITR*(N*(N*UNFL)) ) ELSE * * Absolute accuracy desired * - THRESH = MAX( ABS( TOL )*SMAX, MAXITR*N*N*UNFL ) + THRESH = MAX( ABS( TOL )*SMAX, MAXITR*(N*(N*UNFL)) ) END IF * * Prepare for main iteration loop for the singular values |