pub fn fit( knots: &Knots, points: &DataPoints, params: &Parameters, penalization: Option<Penalization> ) -> Result<MatD, FitError>