|
|
|
@ -249,11 +249,6 @@ PTpwl(double arg, void *data) |
|
|
|
{ |
|
|
|
struct pwldata { int n; double vals[0]; } *thing = data; |
|
|
|
|
|
|
|
// fprintf(stderr, "%s called, data=%p, n=%d\n", __FUNCTION__, data, thing->n); |
|
|
|
// for(i=0; i<thing->n; i++) |
|
|
|
// fprintf(stderr, " %lf", thing->vals[i]); |
|
|
|
// fprintf(stderr, "\n"); |
|
|
|
|
|
|
|
double y; |
|
|
|
|
|
|
|
int k0 = 0; |
|
|
|
@ -271,8 +266,6 @@ PTpwl(double arg, void *data) |
|
|
|
(thing->vals[2*k1+1] - thing->vals[2*k0+1]) * |
|
|
|
(arg - thing->vals[2*k0]) / (thing->vals[2*k1] - thing->vals[2*k0]); |
|
|
|
|
|
|
|
// fprintf(stderr, "foo: interval %lf %lf %lf -> %lf\n", thing->vals[2*k0], arg, thing->vals[2*k1], y); |
|
|
|
|
|
|
|
return y; |
|
|
|
} |
|
|
|
|
|
|
|
@ -281,11 +274,6 @@ PTpwl_derivative(double arg, void *data) |
|
|
|
{ |
|
|
|
struct pwldata { int n; double vals[0]; } *thing = data; |
|
|
|
|
|
|
|
// fprintf(stderr, "%s called, data=%p, n=%d\n", __FUNCTION__, data, thing->n); |
|
|
|
// for(i=0; i<thing->n; i++) |
|
|
|
// fprintf(stderr, " %lf", thing->vals[i]); |
|
|
|
// fprintf(stderr, "\n"); |
|
|
|
|
|
|
|
double y; |
|
|
|
|
|
|
|
int k0 = 0; |
|
|
|
@ -303,7 +291,5 @@ PTpwl_derivative(double arg, void *data) |
|
|
|
(thing->vals[2*k1+1] - thing->vals[2*k0+1]) / |
|
|
|
(thing->vals[2*k1] - thing->vals[2*k0]); |
|
|
|
|
|
|
|
// fprintf(stderr, "bar: interval %lf %lf %lf -> %lf\n", thing->vals[2*k0], arg, thing->vals[2*k1], y); |
|
|
|
|
|
|
|
return y; |
|
|
|
} |