diff options
| author | Bastien Montagne <montagne29@wanadoo.fr> | 2018-07-20 13:13:58 +0300 |
|---|---|---|
| committer | Bastien Montagne <montagne29@wanadoo.fr> | 2018-07-20 13:13:58 +0300 |
| commit | ecfb74833d693a3cd7c49855369127cf4b4d168d (patch) | |
| tree | 7531c1d7672524feb976fa121fcfa8e78497306a /source/blender/blenlib/intern/math_vector_inline.c | |
| parent | 2de283615f081fa76fb05ac0bf81f831651f9e51 (diff) | |
| parent | 710a05954990c347932b9f540a1b220ea2c5c0da (diff) | |
Merge branch 'master' into blender2.8
Diffstat (limited to 'source/blender/blenlib/intern/math_vector_inline.c')
| -rw-r--r-- | source/blender/blenlib/intern/math_vector_inline.c | 10 |
1 files changed, 10 insertions, 0 deletions
diff --git a/source/blender/blenlib/intern/math_vector_inline.c b/source/blender/blenlib/intern/math_vector_inline.c index 4c40921edb6..189b94a6f13 100644 --- a/source/blender/blenlib/intern/math_vector_inline.c +++ b/source/blender/blenlib/intern/math_vector_inline.c @@ -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 |