Dear bug list,
I have found that the function gsl_cheb_calc_deriv() issues a segmentation
fault when both the Chebyshev series and the series for its derivative are
allocated for order N=1. No other value of N seems to effect the same error.
There are two different installations of the GSL that I use at home versus
at work: GSL 1.11 on Linux OpenSUSE 11.0, and GSL 1.12 on OS X 10.6
(installed via macports). The test code that I have included below produces
the same error on both systems.
Yours,
Jerry Gagelman
/* gsl_cheb_bug.c */
#include <stdlib.h>
#include <gsl/gsl_chebyshev.h>
int
main(int argc, char** argv)
{
int N = 1;
double c[2] = {0, 1}; /* initialize series as T_N */
gsl_cheb_series S = {c, N, -1, 1};
gsl_cheb_series* DS = gsl_cheb_alloc(N);
gsl_cheb_calc_deriv(DS, &S);
gsl_cheb_free(DS);
return 0;
}
_______________________________________________
Bug-gsl mailing list
[email protected]
http://lists.gnu.org/mailman/listinfo/bug-gsl
