center() to return the FEM expansion coefficients of the mean field, …

…removed pointwise centering routine
fdaPDE · Jan 23, 2024 · aa697d2 · aa697d2
1 parent adf30a0
commit aa697d2
Show file tree

Hide file tree

Showing 3 changed files with 1 addition and 180,011 deletions.
diff --git a/fdaPDE/models/functional/center.h b/fdaPDE/models/functional/center.h
@@ -38,14 +38,7 @@ CenterReturnType center(const DMatrix<double>& X, SmootherType_&& smoother, Cali
     smoother.init();
     smoother.solve();
     // compute mean matrix and return
-    DMatrix<double> mean = smoother.fitted().replicate(1, X.rows()).transpose();
-    return {X - mean, mean};
-}
-
-// pointwise centering of the data matrix X
-CenterReturnType center(const DMatrix<double>& X) {
-    DMatrix<double> mean = (X.colwise().sum().transpose() / X.rows()).replicate(1, X.rows()).transpose();
-    return {X - mean, mean};
+    return {X - smoother.fitted().replicate(1, X.rows()).transpose(), smoother.f()};
 }
 
 }   // namespace models