star-sp

commit 958cefa602645282b4ac22fcd1fc2f689a6c7f37
parent a1eecf0fa952a1af2f98d4464e7a2b217d7cce6b
Author: Christophe Coustet <christophe.coustet@meso-star.com>
Date:   Tue, 24 Sep 2019 16:58:29 +0200

Slight change in a test

Diffstat:
M
src/test_ssp_ran_triangle.h
 | 
15
++++++++-------

1 file changed, 8 insertions(+), 7 deletions(-)
diff --git a/src/test_ssp_ran_triangle.h b/src/test_ssp_ran_triangle.h
@@ -77,7 +77,7 @@ TEST()
   struct mem_allocator allocator;
   REAL samps[NSAMPS][3];
   REAL A[3], B[3], C[3];
-  REAL v0[3], v1[3], v2[3], v3[3], v4[3], v5[3];
+  REAL v0[3], v1[3], v2[3], m0[3], m1[3], m2[3];
   REAL plane[4];
   REAL pdf;
   size_t counter[2];
@@ -97,9 +97,9 @@ TEST()
   R3_SUB(v0, B, A);
   R3_SUB(v1, C, A);
   R3_SUB(v2, C, B);
-  R3_MINUS(v3, v0);
-  R3_MINUS(v4, v1);
-  R3_MINUS(v5, v2);
+  R3_MINUS(m0, v0);
+  R3_MINUS(m1, v1);
+  R3_MINUS(m2, v2);
   R3_CROSS(plane, v0, v1);
   plane[3] = -R3_DOT(plane, C);
 
@@ -115,14 +115,15 @@ TEST()
     dot = R3_DOT(R3_CROSS(tmp0, tmp0, v0), R3_CROSS(tmp1, v1, v0));
     CHK(sign(dot) == 1);
     R3_SUB(tmp0, samps[i], B);
-    dot = R3_DOT(R3_CROSS(tmp0, tmp0, v2), R3_CROSS(tmp1, v3, v2));
+    dot = R3_DOT(R3_CROSS(tmp0, tmp0, v2), R3_CROSS(tmp1, m0, v2));
     CHK(sign(dot) == 1);
     R3_SUB(tmp0, samps[i], C);
-    dot = R3_DOT(R3_CROSS(tmp0, tmp0, v4), R3_CROSS(tmp1, v5, v4));
+    dot = R3_DOT(R3_CROSS(tmp0, tmp0, m1), R3_CROSS(tmp1, m2, m1));
     CHK(sign(dot) == 1);
 
     area = R3_LEN(tmp1) * 0.5f;
-    CHK(EQ_EPS_R(1 / area, RAN_TRIANGLE_UNIFORM_PDF(A, B, C), (REAL)1.e-6) == 1);
+    CHK(EQ_EPS_R(pdf, RAN_TRIANGLE_UNIFORM_PDF(A, B, C), (REAL)1.e-8) == 1);
+    CHK(EQ_EPS_R(1 / area, pdf, (REAL)1.e-6) == 1);
   }
 
   nsteps = 10000;