diff --git a/SRC/dlaed2.f b/SRC/dlaed2.f
index c26f4f5ba..20e67a548 100644
--- a/SRC/dlaed2.f
+++ b/SRC/dlaed2.f
@@ -313,7 +313,7 @@
       IMAX = IDAMAX( N, Z, 1 )
       JMAX = IDAMAX( N, D, 1 )
       EPS = DLAMCH( 'Epsilon' )
-      TOL = EIGHT*EPS*MAX( ABS( D( JMAX ) ), ABS( Z( IMAX ) ) )
+      TOL = EIGHT*EIGHT*EPS*MAX( ABS( D( JMAX ) ), ABS( Z( IMAX ) ) )
 *
 *     If the rank-1 modifier is small enough, no more needs to be done
 *     except to reorganize Q so that its columns correspond with the
diff --git a/SRC/dlaed8.f b/SRC/dlaed8.f
index cfb97bf34..df6f0162a 100644
--- a/SRC/dlaed8.f
+++ b/SRC/dlaed8.f
@@ -355,7 +355,7 @@
       IMAX = IDAMAX( N, Z, 1 )
       JMAX = IDAMAX( N, D, 1 )
       EPS = DLAMCH( 'Epsilon' )
-      TOL = EIGHT*EPS*ABS( D( JMAX ) )
+      TOL = EIGHT*EIGHT*EPS*ABS( D( JMAX ) )
 *
 *     If the rank-1 modifier is small enough, no more needs to be done
 *     except to reorganize Q so that its columns correspond with the
diff --git a/SRC/slaed2.f b/SRC/slaed2.f
index 6acf91e19..7ebd202e1 100644
--- a/SRC/slaed2.f
+++ b/SRC/slaed2.f
@@ -313,7 +313,7 @@
       IMAX = ISAMAX( N, Z, 1 )
       JMAX = ISAMAX( N, D, 1 )
       EPS = SLAMCH( 'Epsilon' )
-      TOL = EIGHT*EPS*MAX( ABS( D( JMAX ) ), ABS( Z( IMAX ) ) )
+      TOL = EIGHT*EIGHT*EPS*MAX( ABS( D( JMAX ) ), ABS( Z( IMAX ) ) )
 *
 *     If the rank-1 modifier is small enough, no more needs to be done
 *     except to reorganize Q so that its columns correspond with the
diff --git a/SRC/slaed8.f b/SRC/slaed8.f
index bc6b3179a..b0d884b62 100644
--- a/SRC/slaed8.f
+++ b/SRC/slaed8.f
@@ -355,7 +355,7 @@
       IMAX = ISAMAX( N, Z, 1 )
       JMAX = ISAMAX( N, D, 1 )
       EPS = SLAMCH( 'Epsilon' )
-      TOL = EIGHT*EPS*ABS( D( JMAX ) )
+      TOL = EIGHT*EIGHT*EPS*ABS( D( JMAX ) )
 *
 *     If the rank-1 modifier is small enough, no more needs to be done
 *     except to reorganize Q so that its columns correspond with the