.

DrTimothyAldenDavis · Feb 20, 2024 · 1ed91eb · 1ed91eb
1 parent fa147f6
commit 1ed91eb
Showing 1 changed file with 3 additions and 0 deletions.
diff --git a/SPEX/SPEX_Cholesky/Source/spex_cholesky_up_triangular_solve.c b/SPEX/SPEX_Cholesky/Source/spex_cholesky_up_triangular_solve.c
@@ -283,6 +283,9 @@ SPEX_info spex_cholesky_up_triangular_solve
         {
             // x[k] = x[k] / rho[j-1]
 HERE
+// THIS FAILS with SPEX_PANIC: rhos[j] is zero,
+// then SPEX segfaults
+fprintf (stderr, "------ exact divide by rhos[j]: j = %" PRId64" of n %" PRId64 "\n", j, n) ;
             SPEX_MPZ_DIVEXACT(x->x.mpz[k],x->x.mpz[k],
                                             rhos->x.mpz[j-1]);
 HERE