Struct InnerSolutionBuilder 

pub struct InnerSolutionBuilder<'dp> { /* private fields */ }

Builder for InnerSolution that provides sensible defaults and auto-computes derived quantities (nullspace_dim).

Implementations

impl<'dp> InnerSolutionBuilder<'dp>

pub fn new( log_likelihood: f64, penalty_quadratic: f64, beta: Array1<f64>, n_observations: usize, hessian_op: Arc<dyn HessianOperator>, penalty_coords: Vec<PenaltyCoordinate>, penalty_logdet: PenaltyLogdetDerivs, dispersion: DispersionHandling, ) -> Self

Create a builder with the required core fields.

pub fn deriv_provider(self, p: Box<dyn HessianDerivativeProvider + 'dp>) -> Self

pub fn firth_term(self, term: Option<ExactJeffreysTerm>) -> Self

Install a pre-built Jeffreys/Firth term (Tier-A operator-backed via ExactJeffreysTerm::new, or the Tier-B value-only carrier via ExactJeffreysTerm::value_only).

pub fn hessian_logdet_correction(self, correction: f64) -> Self

pub fn penalty_subspace_trace( self, kernel: Option<Arc<PenaltySubspaceTrace>>, ) -> Self

Install the projected-logdet trace kernel that pairs with the hessian_logdet_correction on a rank-deficient penalty surface. See PenaltySubspaceTrace for the derivation and when it is required for gradient consistency.

pub fn rho_curvature_scale(self, scale: f64) -> Self

pub fn rho_prior(self, prior: RhoPrior) -> Self

pub fn nullspace_dim_override(self, dim: f64) -> Self

Override the auto-computed nullspace dimension.

By default, build() computes nullspace_dim as beta.len() - sum(penalty_coord.rank()). Use this when the caller has a different authoritative value (e.g. from stored per-penalty dims).

pub fn ext_coords(self, coords: Vec<HyperCoord>) -> Self

pub fn ext_coord_pair_fn( self, f: Box<dyn Fn(usize, usize) -> HyperCoordPair + Send + Sync>, ) -> Self

pub fn rho_ext_pair_fn( self, f: Box<dyn Fn(usize, usize) -> HyperCoordPair + Send + Sync>, ) -> Self

pub fn fixed_drift_deriv(self, f: FixedDriftDerivFn) -> Self

pub fn contracted_psi_second_order( self, f: Option<ContractedPsiSecondOrderFn>, ) -> Self

Install the direction-contracted second-order ψ hook (#740). When set, the outer-Hessian operator builder uses it instead of the K² per-pair base_h2 ψψ assembly. See ContractedPsiSecondOrderFn.

pub fn barrier_config(self, config: Option<BarrierConfig>) -> Self

pub fn kkt_residual(self, residual: Option<ProjectedKktResidual>) -> Self

pub fn active_constraints( self, block: Option<Arc<ActiveLinearConstraintBlock>>, ) -> Self

Stash the active linear-inequality constraint block carried alongside the inner solution. Used by PenaltySubspaceTrace::with_active_constraints at REML/LAML evaluation time to form the constraint-aware kernel K_T = K_S − K_S Aᵀ (A K_S Aᵀ)⁻¹ A K_S.

pub fn build(self) -> InnerSolution<'dp>

Build the InnerSolution, auto-computing nullspace_dim from penalty coordinates.