Struct MeasureJetAnisotropyJets 

pub struct MeasureJetAnisotropyJets {
    pub q: Array2<f64>,
    pub indices: Vec<LIndex>,
    pub d_first: Vec<Array2<f64>>,
    pub d_second: Vec<Array2<f64>>,
}

The anisotropic energy together with its exact first and second jets with respect to the lower-triangular Cholesky factor entries of L.

indices[a] names the (row, col) of the a-th active lower-triangular entry (column-major over the lower triangle: for each column j, rows j..d). d_first[a] = ∂Q/∂L_{indices[a]}, and d_second[(a, b)] (stored for the full pair grid, symmetric in a, b) is ∂²Q/∂L_{indices[a]}∂L_{indices[b]}.

Fields

q: Array2<f64>

The det-normalized anisotropic energy form (m×m, symmetric PSD).

indices: Vec<LIndex>

Active lower-triangular L-entry indices, in the derivative order.

d_first: Vec<Array2<f64>>

First derivatives ∂Q/∂L_a, one m×m form per active index.

d_second: Vec<Array2<f64>>

Second derivatives ∂²Q/∂L_a∂L_b, indexed by a*n + b over the n = indices.len() active entries (full symmetric grid).

Implementations

impl MeasureJetAnisotropyJets

pub fn n_active(&self) -> usize

Number of active lower-triangular derivative channels.

pub fn second(&self, a: usize, b: usize) -> &Array2<f64>

Borrow the second-derivative form ∂²Q/∂L_a∂L_b.