[Top][All Lists]

[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]
## [Help-gsl] regular coulomb wave function array failing

**From**: |
Andrew Spott |

**Subject**: |
[Help-gsl] regular coulomb wave function array failing |

**Date**: |
Thu, 5 Apr 2012 17:43:51 -0600 |

I'm attempting to find the coulomb wave functions from l = 0 to l = 50 for the
points between x = 0 and x = 1000 for k's of ~0 (I get a bunch of underflow
values for k < .0099) to k = 1:
gsl_sf_coulomb_wave_F_array(0, lmax, 1. /k , k * j * dr , values, exponent)
where j goes from 0 to 10001, dr = .1, k goes from ~ .01 to 1 (in steps of .01
ish).
However, for certain regions in k * j * dr, I get loss of accuracy errors, and
I don't know how to make this function iterate for longer to make those errors
go away.
Any ideas (I've attached the error log out of my debug program which shows the
beginning and end of the loss of accuracy errors, and were they happen(in terms
of k * x = k * j * dr, and k).)
-Andrew

**
**`errors`

*Description:* Binary data

[Prev in Thread] |
**Current Thread** |
[Next in Thread] |

**[Help-gsl] regular coulomb wave function array failing**,
*Andrew Spott* **<=**