debug test reflection

Author: JLBegin

pytissueoptics/rayscattering/opencl/src/propagation.c

@@ -64,7 +64,9 @@ void roulette(float weightThreshold, __global Photon *photons, __global uint *se
 }
 
 void reflect(FresnelIntersection *fresnelIntersection, __global Photon *photons, uint photonID){
+    printf("Initial direction: %f, %f, %f\n", photons[photonID].direction.x, photons[photonID].direction.y, photons[photonID].direction.z);
     rotateAround(&photons[photonID].direction, &fresnelIntersection->incidencePlane, fresnelIntersection->angleDeflection);
+    printf("Reflected direction: %f, %f, %f\n", photons[photonID].direction.x, photons[photonID].direction.y, photons[photonID].direction.z);
 }
 
 void refract(FresnelIntersection *fresnelIntersection, __global Photon *photons, uint photonID){
@@ -104,6 +106,11 @@ float reflectOrRefract(Intersection *intersection, __global Photon *photons, __c
     FresnelIntersection fresnelIntersection = computeFresnelIntersection(photons[photonID].direction, intersection,
                                                                          materials, surfaces, seeds, gid);
 
+    printf("FresnelIntersection.angleDeflection: %f\n", fresnelIntersection.angleDeflection);
+    printf("FresnelIntersection.isReflected: %d\n", fresnelIntersection.isReflected);
+    printf("FresnelIntersection.nextMaterialID: %d\n", fresnelIntersection.nextMaterialID);
+    printf("Intersection.isSmooth: %d\n", intersection->isSmooth);
+
     if (fresnelIntersection.isReflected) {
         if (intersection->isSmooth) {
             // Prevent reflection from crossing the raw surface.

pytissueoptics/rayscattering/opencl/src/vectorOperators.c

@@ -43,6 +43,9 @@ void rotateAroundAxisGlobal(__global float3 *mainVector, __global float3 *axisVe
     }
 
 void rotateAround(__global float3 *mainVector, float3 *axisVector, float theta){
+    printf("rotateAround: %f, %f, %f\n", axisVector->x, axisVector->y, axisVector->z);
+    printf("rotateAround mainVector: %f, %f, %f\n", mainVector->x, mainVector->y, mainVector->z);
+    printf("rotateAround theta: %f\n", theta);
 //    normalizeVectorLocal(axisVector);
     float sint = sin(theta);
     float cost = cos(theta);

pytissueoptics/rayscattering/tests/opencl/src/testCLPhoton.py

@@ -47,9 +47,6 @@ def setUp(self):
         sourcePath = os.path.join(OPENCL_SOURCE_DIR, "propagation.c")
         self.program = CLProgram(sourcePath)
 
-    def tearDown(self):
-        self.program.release()
-
     def testWhenMoveBy_shouldMovePhotonByTheGivenDistanceTowardsItsDirection(self):
         photonResult = self._photonFunc("moveBy", 10)
         self.assertEqual(self.INITIAL_POSITION + self.INITIAL_DIRECTION * 10, photonResult.position)
@@ -420,6 +417,7 @@ def _photonFunc(self, funcName: str, *args) -> PhotonResult:
                                 materialID=0, solidID=self.INITIAL_SOLID_ID, weight=self.INITIAL_WEIGHT)
         npArgs = [np.float32(arg) if isinstance(arg, (float, int)) else arg for arg in args]
         npArgs = [cl.cltypes.make_float3(*arg.array) if isinstance(arg, Vector) else arg for arg in npArgs]
+        print("npArgs:", npArgs)
         self.program.launchKernel(kernelName=funcName + "Kernel", N=1,
                                   arguments=npArgs + [photonBuffer, np.int32(0)])
         photonResult = self._getPhotonResult(photonBuffer)

tox.ini

@@ -8,4 +8,4 @@ setenv =
     PTO_N_WORK_UNITS = 128
     PTO_MAX_MEMORY_MB = 1024
 commands =
-    python -m unittest pytissueoptics/rayscattering/tests/opencl/src/testCLPhoton.py -v
+    python -m unittest -v pytissueoptics/rayscattering/tests/opencl/src/testCLPhoton.py -k testWhenReflectOrRefractWithReflectingIntersection_shouldReflectPhoton