debugging in cpoly

PDLPorters · Feb 13, 2024 · b9c81b5 · b9c81b5
1 parent cbdf65b
commit b9c81b5
Showing 1 changed file with 9 additions and 0 deletions.
diff --git a/Basic/Math/cpoly.c b/Basic/Math/cpoly.c
@@ -258,6 +258,7 @@ char *cpoly(double opr[], double opi[], int degree,
   }
 
   while (!failreason) {
+fprintf(stderr, "NOFAIL yet\n");
 
     /* Start the algorithm for one zero */
     if (nn < 3) {
@@ -318,6 +319,14 @@ char *cpoly(double opr[], double opi[], int degree,
   if (pc) free(pc);
   /* The zerofinder has failed on two major passes
      Return empty handed */
+if (failreason) {
+  fprintf(stderr, "FAIL data (got to %d):\n", nn);
+  int i;
+  for (i=0; i<degree+1; i++) {
+    fprintf(stderr, "  %20g%+20gi\n", opr[i], opi[i]);
+  }
+  fflush(stderr);
+}
   return failreason;
 }