[Bug-gsl] Chebyshev series derivatives

[Jerry Gagelman]

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;
}