diff options
| author |  Benoit Jacob <jacob.benoit.1@gmail.com> | 2009-08-12 18:30:37 -0400 |
|---|---|---|
| committer | Benoit Jacob <jacob.benoit.1@gmail.com> | 2009-08-12 18:30:37 -0400 |
| commit | 99802094e471bb64cc3c3e3d0b2391ba2a1e75a2 (patch) | |
| tree | 2d0a587044f096130ce2aadbfbb1d84d1cffbac0 /Eigen/src/Jacobi | |
| parent | 6bb6810c5e6b0b44310791a3cc19ff0d2bf42835 (diff) | |
do without an empirical homemade formula that i wasn't comfortable about...
turns out it's not needed anymore and removing it seems to only increase the precision
Diffstat (limited to 'Eigen/src/Jacobi')
| -rw-r--r-- | Eigen/src/Jacobi/Jacobi.h | 2 |
1 files changed, 1 insertions, 1 deletions
diff --git a/Eigen/src/Jacobi/Jacobi.h b/Eigen/src/Jacobi/Jacobi.h index 48195c1ec..785dd26ab 100644 --- a/Eigen/src/Jacobi/Jacobi.h +++ b/Eigen/src/Jacobi/Jacobi.h @@ -50,7 +50,7 @@ void MatrixBase<Derived>::applyJacobiOnTheRight(int p, int q, Scalar c, Scalar s template<typename Scalar> bool ei_makeJacobi(Scalar x, Scalar y, Scalar z, Scalar *c, Scalar *s) { - if(ei_abs(y) <= ei_abs(z-x) * 0.5 * machine_epsilon<Scalar>()) + if(y == 0) { *c = Scalar(1); *s = Scalar(0); |