[Top][All Lists]
[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]
[Bug-gsl] Chebyshev series derivatives
From: |
Jerry Gagelman |
Subject: |
[Bug-gsl] Chebyshev series derivatives |
Date: |
Wed, 3 Mar 2010 10:03:27 +0100 |
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] Chebyshev series derivatives,
Jerry Gagelman <=