https://github.com/dimforge/nalgebra/blob/main/src/linalg/svd.rs#L98 Better to use `T::EPSILON` rather than `T::RealField::default_epsilon()`?