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
This commit is contained in:
Benoit Jacob 2009-08-12 18:30:37 -04:00
parent 6bb6810c5e
commit 99802094e4

View File

@ -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);