diff options
-rw-r--r-- | SRC/dlasq3.f | 4 | ||||
-rw-r--r-- | SRC/slasq3.f | 4 |
2 files changed, 4 insertions, 4 deletions
diff --git a/SRC/dlasq3.f b/SRC/dlasq3.f index d044b10a..731a1c5f 100644 --- a/SRC/dlasq3.f +++ b/SRC/dlasq3.f @@ -267,8 +267,8 @@ Z( NN-3 ) = Z( NN-7 ) Z( NN-7 ) = S END IF - IF( Z( NN-5 ).GT.Z( NN-3 )*TOL2 ) THEN - T = HALF*( ( Z( NN-7 )-Z( NN-3 ) )+Z( NN-5 ) ) + T = HALF*( ( Z( NN-7 )-Z( NN-3 ) )+Z( NN-5 ) ) + IF( Z( NN-5 ).GT.Z( NN-3 )*TOL2.AND.T.NE.ZERO ) THEN S = Z( NN-3 )*( Z( NN-5 ) / T ) IF( S.LE.T ) THEN S = Z( NN-3 )*( Z( NN-5 ) / diff --git a/SRC/slasq3.f b/SRC/slasq3.f index 3d158ee5..967e42fb 100644 --- a/SRC/slasq3.f +++ b/SRC/slasq3.f @@ -267,8 +267,8 @@ Z( NN-3 ) = Z( NN-7 ) Z( NN-7 ) = S END IF - IF( Z( NN-5 ).GT.Z( NN-3 )*TOL2 ) THEN - T = HALF*( ( Z( NN-7 )-Z( NN-3 ) )+Z( NN-5 ) ) + T = HALF*( ( Z( NN-7 )-Z( NN-3 ) )+Z( NN-5 ) ) + IF( Z( NN-5 ).GT.Z( NN-3 )*TOL2.AND.T.NE.ZERO ) THEN S = Z( NN-3 )*( Z( NN-5 ) / T ) IF( S.LE.T ) THEN S = Z( NN-3 )*( Z( NN-5 ) / |