From: Ian Jackson <ijackson@chiark.greenend.org.uk>
Date: Sun, 8 Apr 2018 14:47:57 +0000 (+0100)
Subject: Revert "simplex wip: use gsl_vector_get for X, for abandonment"
X-Git-Url: https://www.chiark.greenend.org.uk/ucgi/~ian/git?a=commitdiff_plain;h=937b1b9ee1006655de9fe02de812b81ce836fb65;p=moebius3.git

Revert "simplex wip: use gsl_vector_get for X, for abandonment"

This reverts commit 29cb20f6187aeed6da1f1a7efa3473ee46015ece.
---

diff --git a/findcurve.c b/findcurve.c
index 1dd0941..055018f 100644
--- a/findcurve.c
+++ b/findcurve.c
@@ -26,16 +26,14 @@ static int NP;
 static double *INPUT; /* dyanmic array, on main's stack */
 static double PREP[NPREP];
 
-#define X(xi) gsl_vector_get(xg, (xi))
-
-static void printcore(gsl_vector *xg) {
+static void printcore(const double *X) {
   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");
 }
 
@@ -46,7 +44,8 @@ static void prepare(double X[] /* startpoint */) {
 
 //#define DEBUG
 
-static double cb_f(const gsl_vector *xg, void *params) {
+static double cb_Efunc(void *xp) {
+  const double *X = xp;
   int P;
   DECLARE_F_G;
   CALCULATE_F_G;
diff --git a/symbolic.py b/symbolic.py
index acc9657..8680378 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)):