Struct ConstrainedGaussianRemlForwardProblem 

pub struct ConstrainedGaussianRemlForwardProblem<'a> {
    pub x: ArrayView2<'a, f64>,
    pub y: ArrayView2<'a, f64>,
    pub penalty: ArrayView2<'a, f64>,
    pub weights: Option<ArrayView1<'a, f64>>,
    pub constraints: Option<&'a LinearInequalityConstraints>,
    pub init_lambda: Option<f64>,
}

Inputs to the single-penalty constrained Gaussian REML fit.

The public operation defines one criterion in both regimes: ordinary closed-form Gaussian REML when the unconstrained optimum is interior, and the exact affine-face restriction of that same criterion when constraints bind. There is deliberately no generic-GAM rho prior, shrinkage floor, or ALO correction hidden in this low-level primitive.

Fields

x: ArrayView2<'a, f64>

y: ArrayView2<'a, f64>

penalty: ArrayView2<'a, f64>

weights: Option<ArrayView1<'a, f64>>

constraints: Option<&'a LinearInequalityConstraints>

init_lambda: Option<f64>