Changeset 8487b7 in git

Timestamp:

Jul 4, 2013, 11:11:20 AM (10 years ago)

Author:

Hans Schoenemann <hannes@…>

Branches:

(u'jengelh-datetime', 'ceac47cbc86fe4a15902392bdbb9bd2ae0ea02c6')(u'spielwiese', 'f875bbaccd0831e36aaed09ff6adeb3eb45aeb94')

Children:

90e49f7869b86be3e49012a3f3091d48b56c09ad

Parents:

358f9d7899436a4f0ac58f23dc33e07dce976d44

Message:

chg: simplified luDecomp

from master

File:

• kernel/linearAlgebra.cc (modified) (3 diffs)

kernel/linearAlgebra.cc

@@ -107 +107 @@
   int cc = aMat->cols();
   pMat = mpNew(rr, rr);
-  uMat = mpNew(rr, cc);
-  /* copy aMat into uMat: */
-  for (int r = 1; r <= rr; r++)
-    for (int c = 1; c <= cc; c++)
-      MATELEM(uMat, r, c) = p_Copy(aMat->m[c - 1 + (r - 1) * cc], R);
+  uMat = mp_Copy(aMat,R); /* copy aMat into uMat: */
 
   /* we use an int array to store all row permutations;
@@ -159 +155 @@
          is non-zero: */
       number pivotElement = pGetCoeff(MATELEM(uMat, r, r + cOffset));
-      poly p; number n;
+      poly p;
       for (int rGauss = r + 1; rGauss <= rr; rGauss++)
       {
@@ -165 +161 @@
         if (p != NULL)
         {
-          n = n_Div(pGetCoeff(p), pivotElement, R->cf);
+          number n = n_Div(pGetCoeff(p), pivotElement, R->cf);
           n_Normalize(n,R->cf);