harmonica.jax

harmonica.jax.harmonica_transit_quad_ld(times, t0, period, a, inc, ecc=0., omega=0., u1=0., u2=0., r=jnp.array([0.1]))

Harmonica transits with jax – quadratic limb darkening.

Parameters:

• times (ndarray) – 1D array of model evaluation times [days].

• t0 (float) – Time of transit [days].

• period (float) – Orbital period [days].

• a (float) – Semi-major axis [stellar radii].

• inc (float) – Orbital inclination [radians].

• ecc (float, optional) – Eccentricity [], 0 <= ecc < 1. Default=0.

• omega (float, optional) – Argument of periastron [radians]. Default=0.

• u1 (float.) – Quadratic limb-darkening coefficient.

• u2 (float.) – Quadratic limb-darkening coefficient.

• r (ndarray) –

  Transmission string coefficients. 1D array of N Fourier coefficients that specify the planet radius as a function of angle in the sky-plane. The length of r must be odd, and the final two coefficients must not both be zero.

  \[r_{\rm{p}}(\theta) = \sum_{n=0}^N a_n \cos{(n \theta)} + \sum_{n=1}^N b_n \sin{(n \theta)}\]

  The input array is given as r=[a_0, a_1, b_1, a_2, b_2,..].

Returns:

Normalised transit light curve fluxes [].

Return type:

array

harmonica.jax.harmonica_transit_nonlinear_ld(times, t0, period, a, inc, ecc=0., omega=0., u1=0., u2=0., u3=0., u4=0., r=jnp.array([0.1]))

Harmonica transits with jax – non-linear limb darkening.

Parameters:

• times (ndarray) – 1D array of model evaluation times [days].

• t0 (float) – Time of transit [days].

• period (float) – Orbital period [days].

• a (float) – Semi-major axis [stellar radii].

• inc (float) – Orbital inclination [radians].

• ecc (float, optional) – Eccentricity [], 0 <= ecc < 1. Default=0.

• omega (float, optional) – Argument of periastron [radians]. Default=0.

• u1 (float.) – Non-linear limb-darkening coefficient.

• u2 (float.) – Non-linear limb-darkening coefficient.

• u3 (float.) – Non-linear limb-darkening coefficient.

• u4 (float.) – Non-linear limb-darkening coefficient.

• r (ndarray) –

  Transmission string coefficients. 1D array of N Fourier coefficients that specify the planet radius as a function of angle in the sky-plane. The length of r must be odd, and the final two coefficients must not both be zero.

  \[r_{\rm{p}}(\theta) = \sum_{n=0}^N a_n \cos{(n \theta)} + \sum_{n=1}^N b_n \sin{(n \theta)}\]

  The input array is given as r=[a_0, a_1, b_1, a_2, b_2,..].

Returns:

Normalised transit light curve fluxes [].

Return type:

array