help-gsl
[Top][All Lists]
Advanced

[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

Attachment: errors
Description: Binary data


reply via email to

[Prev in Thread] Current Thread [Next in Thread]