Struct DmDt

pub struct DmDt<T>
where
    T: Float,
{
    pub dt_grid: Grid<T>,
    pub dm_grid: Grid<T>,
}

dm–dt map plotter

Fields

dt_grid: Grid<T>

dm_grid: Grid<T>

Implementations

impl<T> DmDt<T>

where T: Float,

pub fn from_grids<Gdt, Gdm>(dt_grid: Gdt, dm_grid: Gdm) -> Self

where Gdt: Into<Grid<T>>, Gdm: Into<Grid<T>>,

Create new DmDt

pub fn from_lgdt_dm_limits( min_lgdt: T, max_lgdt: T, lgdt_size: usize, max_abs_dm: T, dm_size: usize, ) -> Self

Create new DmDt with logarithmic dt grid and linear dm grid

dt grid will have borders [10^min_lgdt, 10^max_lgdt), dm grid will have borders [-max_abs_dm, max_abs_dm)

pub fn shape(&self) -> (usize, usize)

N dt by N dm

pub fn points(&self, t: &[T], m: &[T]) -> Array2<u64>

Represents each pair of (t, m) points as a unity value in dm-dt map

t must be an ascending slice

pub fn gausses<Erf>(&self, t: &[T], m: &[T], err2: &[T]) -> Array2<T>

where T: ErfFloat, Erf: ErrorFunction<T>,

Represents each pair of (t, m, err2) points as a Gaussian distribution in dm-dt map

t must be an ascending slice.

Each observation is assumed to happen at time moment t_i and have Gaussian distribution of its magnitude N(m_i, err2_i). Each pair of observations (t_1, m_1, err2_1), (t_2, m_2, err2_2) is represented by 1-D Gaussian in the dm-dt space having constant dt and dm ~ N(m2-m1, err2_1 + err2_2). This distribution is integrated over each cell using Erf struct implementing ErrorFunction.

pub fn dt_points(&self, t: &[T]) -> Array1<u64>

Count dt in the each dt grid cell

pub fn cond_prob<Erf>(&self, t: &[T], m: &[T], err2: &[T]) -> Array2<T>

where T: ErfFloat, Erf: ErrorFunction<T>,

Conditional probability p(m2-m1|t2-t1)

Technically it is optimized version of DmDt::gausses() normalized by DmDt::dt_points but with better performance. Mathematically it represents the distribution of conditional probability p(m2-m1|t2-t1), see Soraisam et al. 2020 for details.

Trait Implementations

impl<T> Clone for DmDt<T>

where T: Float + Clone,

fn clone(&self) -> DmDt<T>

Returns a duplicate of the value.

fn clone_from(&mut self, source: &Self)

Performs copy-assignment from source.

impl<T> Debug for DmDt<T>

where T: Float + Debug,

fn fmt(&self, f: &mut Formatter<'_>) -> Result

Formats the value using the given formatter.

impl<'de, T> Deserialize<'de> for DmDt<T>

where T: Float + Deserialize<'de>,

fn deserialize<__D>(__deserializer: __D) -> Result<Self, __D::Error>

where __D: Deserializer<'de>,

Deserialize this value from the given Serde deserializer.

impl<T> Serialize for DmDt<T>

where T: Float + Serialize,

fn serialize<__S>(&self, __serializer: __S) -> Result<__S::Ok, __S::Error>

where __S: Serializer,

Serialize this value into the given Serde serializer.