This is the mail archive of the
gsl-discuss@sources.redhat.com
mailing list for the GSL project.
THANKS: about QAWF integration -- Fourier integrals
- From: xie qiong <xieq at yahoo dot com>
- To: gsl-discuss at sources dot redhat dot com
- Date: Mon, 11 Mar 2002 11:21:43 -0800 (PST)
- Subject: THANKS: about QAWF integration -- Fourier integrals
Hi, Brian,
Thank you very much for your help.
Qiong
--- Brian Gough <bjg@network-theory.co.uk> wrote:
> xie qiong writes:
> > double myfunc(double x, void *params)
> > {
> > double a = *(double *)params;
>
> Hi,
> In this function add,
>
> if (x == 0.0) return 0.0;
>
> to avoid computing 0/0. The cutoff of 1e-30 is
> problematic with the
> current implementation based on the original
> QUADPACK code. This
> could be improved in a future version of GSL.
>
> >
> > return exp(- log(x) * log(x)/(2.0 * a * a)) /
> > (sqrt(2.0 * PI) * a * x);
> > }
>
> regards
> Brian
__________________________________________________
Do You Yahoo!?
Try FREE Yahoo! Mail - the world's greatest free email!
http://mail.yahoo.com/