search

rodluger / Limbdark.jl

Analytical transit light curves for limb darkened stars
MIT License
14 stars 4 forks source link

Possible error in cel_bulirsch for p<0 #27

Closed ericagol closed 6 years ago

ericagol commented 6 years ago

The numerical integration is not agreeing with the analytic formula in test_cel_vec.jl. This may be due to singularity at \phi = \sin^{-1}\sqrt{1/(1-p)}.

ericagol commented 6 years ago

Indeed; the singularity was causing problems. When I do a logarithmic spacing in phi on either side, I get good agreement between the analytic cel integral and a numeric integral.