curveopt: symbolic: wip, before go back to conditional in C

[moebius3.git] / findcurve.c

diff --git a/findcurve.c b/findcurve.c
index f6df6e6fd28f15ba856b2ee5d0f598f393464d40..0fa5c10a7691f085d7d1d53b6e2302a94a12e2a2 100644 (file)
--- a/findcurve.c
+++ b/findcurve.c
@@ -26,6 +26,15 @@ static double target[N];
 
 static double cb_Efunc(void *xp) {
   const double *x = xp;
+
+  e = 0;
+  for (i=0; i<N-3; i++) {
+    // A is point i, B i+1, C i+2, etc.
+    if (i == 0) {
+    } else if (i==N-2) {
+    } else {
+      E_CALCULATE_ED;
+
   double F[N], e;
   int i;
   X_EXTRACT;
@@ -40,6 +49,11 @@ static double cb_Efunc(void *xp) {
   return e;
 }
 
+static int cb_fdf(const gsl_vector *x, void *params,
+                 gsl_vector *f, gsl_matrix *J) {
+  
+}
+
 static void cb_step(const gsl_rng *rng, void *xp, double step_size) {
   double *x = xp;
   int i;