Fix mis-use of fabs in LAPACK
authorMark Abraham <mark.j.abraham@gmail.com>
Wed, 20 Aug 2014 08:20:21 +0000 (10:20 +0200)
committerGerrit Code Review <gerrit@gerrit.gromacs.org>
Wed, 3 Sep 2014 20:44:07 +0000 (22:44 +0200)
clang-3.5 warns that fabs is not the right approach.

Change-Id: Ib7125fcad2f4fda99abf7b67f296f24d84cb34be
(cherry picked from commit d9c0e416e180cdd7e0a158d8d9b9e1fe59ef89b7)

src/gromacs/linearalgebra/gmx_lapack/dlagts.c
src/gromacs/linearalgebra/gmx_lapack/slagts.c

index bc342af697ae32f882e82fce933dc3c81a2504a5..3f784285884f9775378ff45d9ba53098deb0fb26 100644 (file)
@@ -1,3 +1,4 @@
+#include <stdlib.h>
 #include <math.h>
 #include "types/simple.h"
 
@@ -31,7 +32,7 @@ F77_FUNC(dlagts,DLAGTS)(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(dlagts,DLAGTS)(int *job,
        }
     }
 
-    if (fabs(fabs(*job)-1.0)<GMX_DOUBLE_MIN) {
+    if (1 == abs(*job)) {
        i__1 = *n;
        for (k = 2; k <= i__1; ++k) {
            if (in[k - 1] == 0) {
index 364234ddc94fbf82dffd5a2e22c92c1ef90ed196..ff03ff732de3a3417653f8771ab5f92ab190577f 100644 (file)
@@ -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) {