[Help-gsl] gsl interpolation problem

[Ronald Haynes]

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/