Eric Agol

> If you're interested in a Julia version, there has been some work on implementing at least some parts over in: https://github.com/JuliaAstro/Transits.jl Also https://github.com/rodluger/Limbdark.jl

And also join fit of direct-imaging and astrometric data: https://arxiv.org/abs/1910.10169

And, it looks like in their future work list they mention HMC.

I think 1. should suffice in most cases. 2. can often be described by a model for the variation of the orbital elements, which is more akin to photodynamics.

Thanks for developing this package, @montyvesselinov I am also interested in running this in parallel. I was not able to use the `addproc()` command. I did find the `setprocs()` command...

@langfzac This is still an issue which needs to be fixed. It causes a problem when `b=1-r,` which is encountered in some of the tests with Photodynamics.jl. The version without...

@langfzac Okay, I made changes 1, 2, 4, 9 & 10, and this fixed the issue with Limbdark.jl for `b=1-r` (for the case I tested with `r = 0.08740709353364863`). Now...

@langfzac Okay, so commit https://github.com/rodluger/Limbdark.jl/commit/9524317c1c5c04705fa4864d7ad46d21fbe5b638 seems to fix this. I also had to change the tolerance in the `M_n` test since the test wasn't passing. There are some other issues...

@rodluger What radius-ratio `r` do you use? I did find better than ppm precision for septic limb-darkening; you should be able to see this in the latest version of `occult_nonlinear.jl`....

@rodluger The idea is that given a basis of `\mu^j` terms in equation (59), and arbitrary `I(\mu)` function is a linear combination of these terms, and so the coefficients can...