Merge branch 'master' into blender2.8
[blender.git] / source / blender / blenlib / intern / math_vector_inline.c
index 4c40921edb6e2c39d0c39bc8b5a2e1ecf9612a06..189b94a6f1311abdaba75383b4c7357123832e0e 100644 (file)
@@ -753,6 +753,16 @@ MINLINE void cross_v3_v3v3(float r[3], const float a[3], const float b[3])
        r[2] = a[0] * b[1] - a[1] * b[0];
 }
 
+/* cross product suffers from severe precision loss when vectors are
+ * nearly parallel or opposite; doing the computation in double helps a lot */
+MINLINE void cross_v3_v3v3_hi_prec(float r[3], const float a[3], const float b[3])
+{
+       BLI_assert(r != a && r != b);
+       r[0] = (float)((double)a[1] * (double)b[2] - (double)a[2] * (double)b[1]);
+       r[1] = (float)((double)a[2] * (double)b[0] - (double)a[0] * (double)b[2]);
+       r[2] = (float)((double)a[0] * (double)b[1] - (double)a[1] * (double)b[0]);
+}
+
 /* Newell's Method */
 /* excuse this fairly specific function,
  * its used for polygon normals all over the place