[Top][All Lists]
[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]
Re: [Help-gsl] integration w/ qagp, bessel_Jnu
From: |
Jorge Talamantes |
Subject: |
Re: [Help-gsl] integration w/ qagp, bessel_Jnu |
Date: |
Thu, 17 Nov 2005 12:21:58 -0800 |
>Split up the integral (or integrand) to compute the part near x=0
>using the asymptotic form of j(x,n) for small x to avoid underflow.
>Or disable the underflow error if it doesn't affect the final results.
>
>--
>Brian Gough
Hi. Thanks for your reply. I implemented a solution equivalent to your
second suggestion: it turns out that I can prevent gnu_sf_bessel_Jnu(x,n)
from collapsing with an underflow error if I ensure that (n-x)<~70. (n, x
> 0.) I did this by setting the integrand equal to zero if this condition
is not satisfied. Then I checked my integral for a few relevant values of
n and x. (I used matlab for cross-checking.) The integral did not change.
Thanks again.
Jorge
---
Jorge Talamantes, Physics and Geology, California State University,
Bakersfield, CA 93311, jtalamantes at csub dot edu, 661.654.2335