Eigen::FullPivLU<RowMajorMatrixXf>luOfAtOmegaAReg(AtOmegaAReg_Eigen);// Calculate the full-pivoting LU decomposition of the regularized AtA. Note: We could also try FullPivHouseholderQR if our system is non-minimal (i.e. there are more constraints than unknowns).