dot-solver_stabilize_working_weights: Stabilize GLM Working Weights Before Linear Algebra

In lgspline: Lagrangian Multiplier Smoothing Splines for Smooth Function Estimation

Stabilize GLM Working Weights Before Linear Algebra

Description

Keeps iteratively reweighted least-squares weights finite and bounded away from zero so the weighted Gram and Woodbury updates stay numerically well posed.

Usage

.solver_stabilize_working_weights(
  W,
  min_weight = sqrt(.Machine$double.eps),
  max_weight = 10000
)

Arguments

W	Numeric working-weight vector.
min_weight	Lower bound applied after replacing non-finite values.
max_weight	Upper bound for extremely large finite weights.

Value

Numeric vector of stabilized working weights.

lgspline documentation built on May 8, 2026, 5:07 p.m.