Struct FitSensitivity 

pub struct FitSensitivity<'a> { /* private fields */ }

The one sensitivity operator built at the optimum. See module docs.

Implementations

impl<'a> FitSensitivity<'a>

pub fn from_faer_cholesky(factor: &'a FaerCholeskyFactor, dim: usize) -> Self

pub fn from_projected( basis: &'a Array2<f64>, reduced_inverse: &'a Array2<f64>, ) -> Self

pub fn dim(&self) -> usize

Coefficient dimension p the operator acts on.

pub fn apply(&self, rhs: &Array1<f64>) -> Array1<f64>

H⁻¹ · rhs (pseudo-inverse action for the projected variant).

pub fn apply_multi(&self, rhs: ArrayView2<'_, f64>) -> Array2<f64>

H⁻¹ · RHS for a (p × m) block of right-hand sides — the batched form every multi-channel consumer should use (outer ρ-pair solves, ALO’s H⁻¹Xᵀ leverage block) so the factor is traversed once per block instead of once per column.

pub fn mode_response(&self, dg_dt: ArrayView2<'_, f64>) -> Option<Array2<f64>>

The IFT mode response ∂β̂/∂t = −H⁻¹ · ∂g/∂t for a (p × m) block of score perturbations — THE object of #935.

Returns None if any solved entry is non-finite (the factored curvature was unusable for this channel); callers must not substitute an approximation, matching the contract of the deleted ift_dbeta_drho_from_solver.

pub fn mode_response_coned( &self, hessian: ArrayView2<'_, f64>, dg_dt: ArrayView2<'_, f64>, col_supports: &[Range<usize>], ) -> Option<Array2<f64>>

Cone-of-influence mode response (#779), the lazy/local form of Self::mode_response. Each perturbation column ∂g/∂t_a is structurally supported only within col_supports[a], so its response −H⁻¹ ∂g/∂t_a is exactly zero outside the coupling component of hessian containing that support. Columns whose support is empty (a structurally inactive channel) are skipped with no solve; the active columns are solved as ONE batched block through Self::apply_multi — strictly better than the per-column BLAS-2 loop this replaces — and each result confined to its cone. On a fully coupled hessian every cone is the whole space and the result equals Self::mode_response bit-for-bit.

hessian must be the same curvature this operator inverts; a dimension mismatch (or any non-finite solved entry) returns None rather than silently substituting an approximation.