Changeset e5ce1a in git

Timestamp:

Jul 20, 2022, 6:02:55 PM (22 months ago)

Author:

Hans Schoenemann <hannes@…>

Branches:

(u'spielwiese', 'fe61d9c35bf7c61f2b6cbf1b56e25e2f08d536cc')

Children:

b0600ca938adb77a8cb08f73867bd8009f2ca8b1

Parents:

9d65054e5b279af66d19c244318144b3989d4388

Message:

fix: spasm_kernel

File:

• Singular/dyn_modules/sispasm/sispasm.cc (modified) (2 diffs)

Singular/dyn_modules/sispasm/sispasm.cc

@@ -103 +103 @@
 
 spasm* sp_kernel(spasm* A, const ring R)
-{ /*from kernel.c*/
+{
   int n = A->n;
   int m = A->m;
   int*  p = (int*)spasm_malloc(n * sizeof(int));
   int * qinv = (int*)spasm_malloc(m * sizeof(int));
+#if 0
+  /*from kernel.c*/
   spasm_find_pivots(A, p, qinv);  /* this does some useless stuff, but
                                    * pushes zero rows to the bottom */
@@ -127 +129 @@
 
   spasm* K = spasm_kernel(A_t, qinv);
+  spasm_csr_free(A_t);
+#else
+  spasm_find_pivots(A, qinv, p);
+  spasm* K = spasm_kernel(A, qinv);
+#endif
   free(p);
   free(qinv);
-  spasm_csr_free(A_t);
   return K;
 }