Details | Last modification | View Log | RSS feed
| Rev | Author | Line No. | Line |
|---|---|---|---|
| 20 | pmbaty | 1 | #include "plane.h" |
| 2 | #include "CORE/STD/brmath.h" |
||
| 3 | #include "harness/trace.h" |
||
| 4 | #include <math.h> |
||
| 5 | |||
| 6 | // IDA: br_int_32 __cdecl BrPlaneEquation(br_vector4 *eqn, br_vector3 *v0, br_vector3 *v1, br_vector3 *v2) |
||
| 7 | br_int_32 BrPlaneEquation(br_vector4* eqn, br_vector3* v0, br_vector3* v1, br_vector3* v2) { |
||
| 8 | //br_vector3 a; // Pierre-Marie Baty -- unused variable |
||
| 9 | //br_vector3 b; // Pierre-Marie Baty -- unused variable |
||
| 10 | float ax; |
||
| 11 | float ay; |
||
| 12 | float az; |
||
| 13 | float bx; |
||
| 14 | float by; |
||
| 15 | float bz; |
||
| 16 | float nx; |
||
| 17 | float ny; |
||
| 18 | float nz; |
||
| 19 | float d; |
||
| 20 | float l; |
||
| 21 | |||
| 22 | ax = v1->v[0] - v0->v[0]; |
||
| 23 | ay = v1->v[1] - v0->v[1]; |
||
| 24 | az = v1->v[2] - v0->v[2]; |
||
| 25 | bx = v2->v[0] - v0->v[0]; |
||
| 26 | by = v2->v[1] - v0->v[1]; |
||
| 27 | bz = v2->v[2] - v0->v[2]; |
||
| 28 | nx = ay * bz - az * by; |
||
| 29 | ny = az * bx - ax * bz; |
||
| 30 | nz = ax * by - ay * bx; |
||
| 31 | l = nx * nx + ny * ny + nz * nz; |
||
| 32 | l = BrFloatSqrt(l); |
||
| 33 | if (isnan(l)) { |
||
| 34 | nx = 0.f; |
||
| 35 | ny = 0.f; |
||
| 36 | nz = 1.f; |
||
| 37 | } else { |
||
| 38 | nx = nx * (1.f / l); |
||
| 39 | ny = ny * (1.f / l); |
||
| 40 | nz = nz * (1.f / l); |
||
| 41 | } |
||
| 42 | d = nx * v0->v[0] + ny * v0->v[1] + nz * v0->v[2]; |
||
| 43 | eqn->v[0] = nx; |
||
| 44 | eqn->v[1] = ny; |
||
| 45 | eqn->v[2] = nz; |
||
| 46 | eqn->v[3] = d; |
||
| 47 | return !isnan(l); |
||
| 48 | } |