Changed 1e-20 to 1e-16. 1e-20 was cutting it too close on 32-bit system, resulting in divide by zero later
parent
5b717c801f
commit
e1c0ad8335
|
@ -151,9 +151,9 @@ Vector Unit3::localCoordinates(const Unit3& y) const {
|
||||||
double dot = p.dot(q);
|
double dot = p.dot(q);
|
||||||
|
|
||||||
// Check for special cases
|
// Check for special cases
|
||||||
if (std::abs(dot - 1.0) < 1e-20)
|
if (std::abs(dot - 1.0) < 1e-16)
|
||||||
return (Vector(2) << 0, 0);
|
return (Vector(2) << 0, 0);
|
||||||
else if (std::abs(dot + 1.0) < 1e-20)
|
else if (std::abs(dot + 1.0) < 1e-16)
|
||||||
return (Vector(2) << M_PI, 0);
|
return (Vector(2) << M_PI, 0);
|
||||||
else {
|
else {
|
||||||
// no special case
|
// no special case
|
||||||
|
|
Loading…
Reference in New Issue