Changeset 14914

Timestamp:

Sep 20, 2007, 10:39:28 AM (19 years ago)

Author:

eugene

Message:

report error in gaussj

File:

• trunk/Ohana/src/opihi/cmd.data/gaussj.c (modified) (3 diffs)

trunk/Ohana/src/opihi/cmd.data/gaussj.c

@@ -5 +5 @@
   float *m, *v;
   double **a, **b;
-  int i, j, N;
+  int i, j, N, status;
   Vector *B;
   Buffer *A;
@@ -33 +33 @@
     b[i][0] = v[i]; 
   }
-  dgaussj (a, N, b, 1);
+
+  status = dgaussj (a, N, b, 1);
 
   for (i = 0; i < N; i++) {
@@ -48 +49 @@
   free (a);
   free (b);
-  return (TRUE);
+
+  if (!status) {
+      gprint (GP_ERR, "failure in matrix solution\n");
+  }
+  return (status);
 
  usage: