From 71f288c5d00262c60e27766f2e7d90cf6dc54030 Mon Sep 17 00:00:00 2001
From: Ian Jackson <ijackson@chiark.greenend.org.uk>
Date: Sat, 7 Apr 2018 23:41:34 +0100
Subject: [PATCH] curveopt: wip

Signed-off-by: Ian Jackson <ijackson@chiark.greenend.org.uk>
---
 curveopt.py |  1 +
 findcurve.c | 20 +++++++++++++++++---
 2 files changed, 18 insertions(+), 3 deletions(-)

diff --git a/curveopt.py b/curveopt.py
index b0e725c..5ebe166 100644
--- a/curveopt.py
+++ b/curveopt.py
@@ -24,6 +24,7 @@ class OptimisedCurve():
     db = DiscreteBezier(cp, nt, bezier_constructor=BezierSegment)
 
     fc_input = map(db.point_at_it, range(0, nt+1))
+    dbg(repr(fc_input))
 
     for end in (False,True):
       ei = nt if end else 0
diff --git a/findcurve.c b/findcurve.c
index a8dfc46..0b13a21 100644
--- a/findcurve.c
+++ b/findcurve.c
@@ -42,23 +42,37 @@ static void prepare(double X[] /* startpoint */) {
   PREPARE;
 }
 
+#define DEBUG
+
 static double cb_Efunc(void *xp) {
   const double *X = xp;
   DECLARE_F_G;
   CALCULATE_F_G;
 
-  printf(" Efunc ");
-  printcore(X);
+#ifdef DEBUG
+  int i,j;
+  printf(" Efunc\n");
+  for (j=0; j<3; j++) {
+    for (i=0; i<NP; i++)
+      printf(" %7.4f", POINT(i)[j]);
+    printf("\n");
+  }
+  printf("             ");
+#endif
 
   double e = 0;
   int P;
   for (P=0; P<NP-3; P++) {
     double P_cost;
     CALCULATE_COST;
-    printf(" %.18g", P_cost);
+#ifdef DEBUG
+    printf(" %7.4f", P_cost);
+#endif
     e += P_cost;
   }
+#ifdef DEBUG
   printf("\n");
+#endif
   return e;
 }
 
-- 
2.30.2