Hi everyone,
I have a program in which I would like to compute the interpolation
coefficients (cubic splines say...) and then I would like to have these
coefficients available globally to other functions.
So now I have done the following:
gsl_interp_accel *acc;
gsl_spline *spline;
before main
inside main I have done
gsl_interp_accel *acc
= gsl_interp_accel_alloc ();
gsl_spline *spline
= gsl_spline_alloc (gsl_interp_cspline,M);
gsl_spline_init (spline,tvalues.e,yvalues.e,M);
val = gsl_spline_eval (spline, 0.6, acc);
and val is correctly computed.
main then calls another function X (say) which in turn calls function Y
(say) a bunch of times. Inside Y I need to evaluate the same spline at a
bunch of different points. But if I do
val = gsl_spline_eval (spline, t, acc);
inside function Y the code terminates during runtime with a signal 11
termination. (It seems to compile fine...).
Any thoughts?
Thanks,
R Haynes
--
Dr. Ronald D. Haynes (PhD)
Assistant Professor
Mathematics and Statistics
Acadia University
Wolfville, N.S.
Using Opera's revolutionary e-mail client: http://www.opera.com/mail/
_______________________________________________
Help-gsl mailing list
[email protected]
http://lists.gnu.org/mailman/listinfo/help-gsl
