search

fbartolic / caustic

Fitting single lens microlensing events with PyMC3
MIT License
1 stars 1 forks source link

Implement finite source effects #14

Open fbartolic opened 5 years ago

fbartolic commented 5 years ago

Finite source effects for single lens effects involve solving elliptical integrals which is costly. Current modeling codes use pre-computed tables for these elliptic integrals and interpolate a 1D function. If this approximation is fine it should be straightforward to implement in PyMC3, although instead of calling a table, I should fit a GP or a polynomial instead. Relevant references are here and here.

fbartolic commented 5 years ago

Implemented a simple Finite source model which relies on a pre-computed table for the elliptic integrals. Works well with single-band data, performance with multi-wavelength data is slow. Need to generalize this to other forms of limb darkening laws at some point.

fbartolic commented 4 years ago

Reopening this because I need to rethink parametrization options and limb darkening prescriptions.