From: Ian Jackson <ijackson@chiark.greenend.org.uk>
Date: Sun, 19 Nov 2017 20:14:06 +0000 (+0000)
Subject: findcurve: python output
X-Git-Url: http://www.chiark.greenend.org.uk/ucgi/~ian/git?a=commitdiff_plain;h=7bcf47dc67a95e6296c0cd71fad33de204bb94b8;p=moebius3.git

findcurve: python output

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

diff --git a/findcurve.c b/findcurve.c
index 169116a..9277db4 100644
--- 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");
   }
 }