question

Author: walchko

src/squaternion.hpp

@@ -33,10 +33,12 @@ struct Quaternion {
         data_t m = sqrt(w * w + x * x + y * y + z * z);
 
         if (m < 1.0e-6) {
-            // throw QuaternionError();
             return 1;
         }
 
+        // use isinf(inv_m) instead?
+        // the above test will capture really large inv_m
+        // in addition to infinate
         data_t inv_m = 1.0 / m;
 
         this->w *= inv_m;