Dear GSL, I have the following question. I need to pass thousands of times two arrays to a custom function in order to perform a spline interpolation. However allocating at each call in the function the dimension of the object and then freeing it via the commands
spline=gsl_spline_alloc(gsl_interp_cspline, BIN); gsl_spline_init(spline, array_x, array_y, BIN); \*commands /* gsl_spline_free(spline); is very time consuming. What I would like to do is allocating the object just once, i.e. when the function is called for the first time, and then proceed on up to the end of the loop. How can I do that ? Thank you in advance Best regards Ruben Farinelli