Changeset 70a952 in git for Singular/iparith.cc

Timestamp:

May 29, 2012, 12:13:49 PM (11 years ago)

Author:

Hans Schoenemann <hannes@…>

Branches:

(u'spielwiese', '8e0ad00ce244dfd0756200662572aef8402f13d5')

Children:

9aac3e40e249db0a2411b4c4f36cdc9a5d791080

Parents:

5998134399af73455bd45b86419f1c2d580ea52d

Message:

fix: matrices in lusolve, ludecomp, luinverse must be constant (from master)

File:

• Singular/iparith.cc (modified) (4 diffs)

Singular/iparith.cc

@@ -4415 +4415 @@
      A list [P, L, U] is returned. */
   matrix mat = (const matrix)v->Data();
+  if (!idIsConstant((ideal)mat))
+  {
+    WerrorS("matrix must be constant");
+    return TRUE;
+  }
   matrix pMat;
   matrix lMat;
@@ -6816 +6821 @@
         return TRUE;
       }
+      if (!idIsConstant((ideal)aMat))
+      {
+        WerrorS("matrix must be constant");
+        return TRUE;
+      }
       invertible = luInverse(aMat, iMat);
     }
@@ -6836 +6846 @@
        return TRUE;
      }
+      if (!idIsConstant((ideal)pMat)
+      || (!idIsConstant((ideal)lMat))
+      || (!idIsConstant((ideal)uMat))
+      )
+      {
+        WerrorS("matricesx must be constant");
+        return TRUE;
+      }
      invertible = luInverseFromLUDecomp(pMat, lMat, uMat, iMat);
   }
@@ -6916 +6934 @@
     Werror("third matrix (%d x %d) and vector (%d x 1) do not fit",
            uMat->rows(), uMat->cols(), bVec->rows());
+    return TRUE;
+  }
+  if (!idIsConstant((ideal)pMat)
+  ||(!idIsConstant((ideal)lMat))
+  ||(!idIsConstant((ideal)uMat))
+  ||(!idIsConstant((ideal)bVec))
+  )
+  {
+    WerrorS("matrices must be constant");
     return TRUE;
   }