From 29cb20f6187aeed6da1f1a7efa3473ee46015ece Mon Sep 17 00:00:00 2001
From: Ian Jackson <ijackson@chiark.greenend.org.uk>
Date: Sun, 8 Apr 2018 15:46:07 +0100
Subject: [PATCH] simplex wip: use gsl_vector_get for X, for abandonment

This is too painful because it breaks POINT, which everyone expects to
be able to inspect with []

Signed-off-by: Ian Jackson <ijackson@chiark.greenend.org.uk>
---
 findcurve.c | 9 +++++----
 symbolic.py | 6 +++---
 2 files changed, 8 insertions(+), 7 deletions(-)

diff --git a/findcurve.c b/findcurve.c
index 055018f..1dd0941 100644
--- a/findcurve.c
+++ b/findcurve.c
@@ -26,14 +26,16 @@ static int NP;
 static double *INPUT; /* dyanmic array, on main's stack */
 static double PREP[NPREP];
 
-static void printcore(const double *X) {
+#define X(xi) gsl_vector_get(xg, (xi))
+
+static void printcore(gsl_vector *xg) {
   int i, j;
   DECLARE_F_G;
   CALCULATE_F_G;
   printf("[");
   for (i=0; i<NP; i++)
     for (j=0; j<3; j++)
-      printf(" %25.18g,", POINT(i)[j]);
+      printf(" %25.18g,", POINT(i,j));
   printf(" ]\n");
 }
 
@@ -44,8 +46,7 @@ static void prepare(double X[] /* startpoint */) {
 
 //#define DEBUG
 
-static double cb_Efunc(void *xp) {
-  const double *X = xp;
+static double cb_f(const gsl_vector *xg, void *params) {
   int P;
   DECLARE_F_G;
   CALCULATE_F_G;
diff --git a/symbolic.py b/symbolic.py
index 8680378..acc9657 100644
--- a/symbolic.py
+++ b/symbolic.py
@@ -237,11 +237,11 @@ def gen_point_references():
   cprint(' (PP) == 1    ? F :')
   cprint(' (PP) == NP-2 ? G :')
   cprint(' (PP) == NP-1 ? H :')
-  cprint(' &X[3*((PP)-2)]')
+  cprint(' &X(3*((PP)-2))')
   cprintraw(')')
 
-  cprintraw('#define EFlq X[ NX_DIRECT + 0 ]')
-  cprintraw('#define HGlq X[ NX_DIRECT + 1 ]')
+  cprintraw('#define EFlq X( NX_DIRECT + 0 )')
+  cprintraw('#define HGlq X( NX_DIRECT + 1 )')
   cprintraw('#define NX    ( NX_DIRECT + 2 )')
 
   for ai in range(0, len(abcd)):
-- 
2.30.2