Fix mis-use of fabs in LAPACK

[alexxy/gromacs.git] / src / gromacs / linearalgebra / gmx_lapack / slagts.c

diff --git a/src/gromacs/linearalgebra/gmx_lapack/slagts.c b/src/gromacs/linearalgebra/gmx_lapack/slagts.c
index 364234ddc94fbf82dffd5a2e22c92c1ef90ed196..ff03ff732de3a3417653f8771ab5f92ab190577f 100644 (file)
--- a/src/gromacs/linearalgebra/gmx_lapack/slagts.c
+++ b/src/gromacs/linearalgebra/gmx_lapack/slagts.c
@@ -1,3 +1,4 @@
+#include <stdlib.h>
 #include <math.h>
 #include "types/simple.h"
 
@@ -31,7 +32,7 @@ F77_FUNC(slagts,SLAGTS)(int *job,
     --a;
 
     *info = 0;
-    if (fabs(*job) > 2 || *job == 0) {
+    if (abs(*job) > 2 || *job == 0) {
        *info = -1;
     } else if (*n < 0) {
        *info = -2;
@@ -77,7 +78,7 @@ F77_FUNC(slagts,SLAGTS)(int *job,
        }
     }
 
-    if (fabs(fabs(*job)-1.0)<GMX_FLOAT_MIN) {
+    if (1 == abs(*job)) {
        i__1 = *n;
        for (k = 2; k <= i__1; ++k) {
            if (in[k - 1] == 0) {