[Help-gsl] multiple calls to gsl_interp

[Ronald Haynes]

Hi, if I create interp objects

gsl_interp_accel *acc;
gsl_spline *spline;

is there any reason that I can't call use them multiple times within a code, 
for example


  acc = gsl_interp_accel_alloc ();
                      spline = gsl_spline_alloc (gsl_interp_cspline,M)          
               gsl_spline_init (spline,tvalues.e,yvalues.e,M);


and then later within the same function

  acc = gsl_interp_accel_alloc ();
                      spline = gsl_spline_alloc (gsl_interp_cspline,N)          
               gsl_spline_init (spline,tvalues.e,yvalues.e,N);

where N is usually not the same as M (could be smaller or larger).  The 
dimension of tvalues,yvalues is >=max(M,N)

Just wondering if there is any problem with memory issues.

Thanks,
R Haynes



Dr. R. D. Haynes
Assistant Professor
Mathematics and Statistics
Acadia University
Wolfville, N.S. B4P2R6
Canada