Don't call LAPACK when LU-decomposing empty real/complex matrices

Otherwise, Macaulay2 crashes Fixes bug reported by Joel Louwsma in Zulip
Macaulay2 · May 15, 2024 · 2d9afc3 · 2d9afc3
1 parent 1b27bd8
commit 2d9afc3
Showing 1 changed file with 6 additions and 0 deletions.
diff --git a/M2/Macaulay2/e/dmat-lu-inplace.hpp b/M2/Macaulay2/e/dmat-lu-inplace.hpp
@@ -284,6 +284,9 @@ inline void DMatLUinPlace<M2::ARingRR>::computeLU()
   int info;
   int min = (rows <= cols) ? rows : cols;
 
+  if (min == 0)
+    return;
+
   // printf("entering DMatLUinPlace::computeLUNaive for RR\n");
 
   int* perm = newarray_atomic(int, min);
@@ -345,6 +348,9 @@ inline void DMatLUinPlace<M2::ARingCC>::computeLU()
   int info;
   int min = (rows <= cols) ? rows : cols;
 
+  if (min == 0)
+    return;
+
   // printf("entering DMatLUtemplate::computeLUNaive for RR\n");
 
   int* perm = newarray_atomic(int, min);