diff options
| author | james <james@8a072113-8704-0410-8d35-dd094bca7971> | 2012-06-29 11:47:00 +0000 |
|---|---|---|
| committer | james <james@8a072113-8704-0410-8d35-dd094bca7971> | 2012-06-29 11:47:00 +0000 |
| commit | 38febd8f674d1a3adf98334c422509d755b08a4b (patch) | |
| tree | 8192c262ac39f51edd72c64747804353301ae133 | |
| parent | fb89702b460bd02d29d453ecff8202f7180fc7da (diff) | |
Added fix to prevent dividing by T when T=0 provided by Alexander Kobotov of Intel.
| -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 ) / |