bd_solver

modspectra.bd_solver(ell, xyz, z_sigma_lim, Hz, bd_max, el_constant1, el_constant2)

Function to solve for the ellipse equation to fit into form of ellipse_equation Chooses only to solve the equation numerically when necessary, avoiding the special cases. Funciton written in form to use with multiprocessing.pool and functools.partial

The defining characteristic of the Ellipse is set by the relation between the semi-major and semi-minor axis with the following equation from Liszt & Burton (1982): ad / bd = el_constant1 + el_constant2 * bd / bd_max

You can trick the “Ellipse” shape into a “Circle” by setting el_constant1 = 1., el_constant2 = 0.

Parameters:

ell: ‘int’

element number to iterate over

xyz: ‘ndarray with shape (3,N)’

xyz-coordinates

z_sigma_lim: ‘number’

sigma cuttoff to stop solving Ellipse equation for for z above a specified scale height threshold

Hz: ‘number’

Scale height along z axis

bd_max: ‘number’

Maximum semi-minor axis allowed within defined elliptical disk

el_constant1: ‘number’

First parameter for defining ellipse

el_constant2: ‘number’

second parameter for defining ellipse