Signed-off-by: Ian Jackson <ijackson@chiark.greenend.org.uk>
return sqrt(s);
}
-static void cb_print(void *xp) {
- const double *x = xp;
- double F[N];
+static void printcore(const double *x) {
int i;
+ double F[N];
X_EXTRACT;
F_POPULATE;
- printf("\n[");
+ printf("[");
for (i=0; i<6; i++) printf(" %.18g,", x[i]);
for (i=0; i<6; i++) printf(" %.18g,", F[i]);
printf(" ]\n");
}
+static void __attribute__((unused)) cb_print(void *xp) {
+ const double *x = xp;
+ printf("\n");
+ printcore(x);
+}
+
static double scan1double(void) {
double v;
int r;
0, // cb_print,
0,0,0, sizeof(startpoint), siman_params);
- cb_print(startpoint);
+ printcore(startpoint);
printf("[]\n");
fflush(stdout);
hc.findcurve_result = l[0:6]
hc.threshold = l[0]**2
hc.total_dist = hc.threshold + l[1]**2
- vdbg().curve( hc.point_at_t )
+ #vdbg().curve( hc.point_at_t )
def point_at_t(hc, normalised_parameter):
dist = normalised_parameter * hc.total_dist