Struct ReparamResult 

pub struct ReparamResult {
    pub s_transformed: Array2<f64>,
    pub log_det: f64,
    pub det1: Array1<f64>,
    pub qs: Array2<f64>,
    pub canonical_transformed: Vec<CanonicalPenalty>,
    pub e_transformed: Array2<f64>,
    pub u_truncated: Array2<f64>,
    pub penalty_shrinkage_ridge: f64,
}

Result of the stable reparameterization algorithm from Wood (2011) Appendix B

Fields

s_transformed: Array2<f64>

Penalty matrix in TRANSFORMED coefficient coordinates.

This must be compatible with beta_transformed and X_transformed = X * Qs.

log_det: f64

Log-determinant of the penalty matrix (stable computation)

det1: Array1<f64>

First derivatives of log-determinant w.r.t. log-smoothing parameters

qs: Array2<f64>

Orthogonal transformation matrix Qs

canonical_transformed: Vec<CanonicalPenalty>

Canonical penalties in the TRANSFORMED coordinate frame. The single source of truth for penalty roots in the transformed frame. Downstream consumers use these for block-local PenaltyCoordinate construction, TK correction, and ext-coord paths.

e_transformed: Array2<f64>

Lambda-dependent penalty square root in TRANSFORMED coordinates (rank x p matrix). This is used for applying the actual penalty in the least squares solve.

u_truncated: Array2<f64>

Truncated eigenvectors (p × m where m = p - structural_rank).

Coordinate frame note:

• This matrix is stored in the TRANSFORMED coefficient frame (post-Qs), i.e. it is compatible with canonical_transformed, beta_transformed, and transformed Hessians without additional coordinate mapping.

These vectors span the structural null space used by positive-part log-determinant conventions.

penalty_shrinkage_ridge: f64

The rho-independent shrinkage ridge magnitude that was added to each eigenvalue of the penalized block. Zero means no shrinkage was applied.

Trait Implementations

impl Clone for ReparamResult

fn clone(&self) -> ReparamResult

Returns a duplicate of the value.

fn clone_from(&mut self, source: &Self)

Performs copy-assignment from source.