Struct TensorChannelHessian 

pub struct TensorChannelHessian {
    pub h: Array3<f64>,
}

A FamilyChannelHessian backed directly by a pre-computed (n × K × K) tensor. Used by the default channel_hessian_at implementation and by tests.

This is the β-independent path: fill_subject reads from the frozen tensor without any recomputation.

Fields

h: Array3<f64>

Trait Implementations

impl FamilyChannelHessian for TensorChannelHessian

fn n_outputs(&self) -> usize

Number of output channels n_outputs (= K in the row Jacobian).

fn n_subjects(&self) -> usize

Number of subjects (rows).

fn fill_subject(&self, i: usize, out: &mut [f64])

Fill the n_outputs × n_outputs per-subject channel Hessian W_i into out (row-major, length n_outputs * n_outputs) for subject i. Negative eigenvalues must be clamped to zero (PSD projection) before or inside this call.

fn evaluate_full(&self) -> Array3<f64>

Materialise the full (n_subjects × n_outputs × n_outputs) tensor. Default implementation calls fill_subject for each row.

fn channel_hessian_at( &self, beta: &[f64], family_scalars: Option<&Arc<dyn Any + Send + Sync>>, ) -> Result<Arc<dyn FamilyChannelHessian>, String>

Return a refreshed W evaluated at beta using family_scalars when those scalars carry the per-row primary state at the current β.