findcurve: python output

Signed-off-by: Ian Jackson <ijackson@chiark.greenend.org.uk>

diff --git a/findcurve.c b/findcurve.c
index 169116a826df4969cd9cfdab8c1ae9ad0c037829..9277db403c3b8eeb92ff56a83a2915db83aaba1e 100644 (file)
--- a/findcurve.c
+++ b/findcurve.c
@@ -18,7 +18,7 @@ static inline _Bool IS_SMALL(double v) {
 }
 
 static inline double sinc(double x) {
-  return gsl_sf_sinc(x / M_PI);  
+  return gsl_sf_sinc(x / M_PI);
 }
 
 static int cb_f(const gsl_vector *x, void *params, gsl_vector *f) {
@@ -89,6 +89,11 @@ int main(void) {
       r = gsl_multiroot_fdfsolver_iterate(solver);
       if (r) break;
 
+      printf("[");
+      for (i=0; i<6; i++) printf(" %.18g,", gsl_vector_get(solver->x, i));
+      for (i=0; i<6; i++) printf(" %.18g,", gsl_vector_get(solver->f, i));
+      printf(" ]\n");
+
       r = gsl_multiroot_test_residual(solver->f, epsilon);
       if (r != GSL_CONTINUE) break;
     }
@@ -97,5 +102,7 @@ int main(void) {
       fprintf(stderr,"ERROR %s\n",gsl_strerror(r));
       exit(-1);
     }
+
+    printf("None\n");
   }
 }