meepmeep.numba3d.solve3d_d

meepmeep.numba3d.solve3d_d(te, p, a, i, e, w, lan: float = 0.0) → tuple[ndarray[tuple[Any, ...], dtype[_ScalarT]], ndarray[tuple[Any, ...], dtype[_ScalarT]]]

Calculate Taylor expansion coefficients and their parameter derivatives for the 3D position around a given expansion-point time relative to the transit centre.

Parameters:

tefloat

Expansion-point time: the time of the Taylor-series expansion [days], measured relative to the transit centre (time of inferior conjunction). te=0 expands at the transit centre.

pfloat

Orbital period [days].

afloat

Semi-major axis of the orbit [R_star].

ifloat

Inclination of the orbit [rad].

efloat

Eccentricity of the orbit.

wfloat

Argument of periastron [rad].

lanfloat, optional

Longitude of the ascending node [rad]. A constant counterclockwise rotation of the sky-plane (x, y) coordinates about the line of sight; the line-of-sight (z) coordinate is unaffected. Defaults to 0.0.

Returns:

cfndarray (3, 5)

Position Taylor coefficients (identical to solve3d output).

dcfndarray (7, 3, 5)

Parameter derivative coefficients. dcf[k] = d(cf)/d(theta_k) for theta = (tc, p, a, i, e, w, lan). Row 0 is the derivative with respect to the transit-centre time tc (dM/dtc = -n); row 6 is the derivative with respect to the longitude of the ascending node.