Revert initial points to surface

ViennaTools · Mar 26, 2024 · 4adeb1f · 4adeb1f
1 parent 7d85b25
commit 4adeb1f
Show file tree

Hide file tree

Showing 2 changed files with 5 additions and 16 deletions.
diff --git a/include/psDiffusion.hpp → include/viennaps/psDiffusion.hpp b/include/psDiffusion.hpp → include/viennaps/psDiffusion.hpp
diff --git a/include/viennaps/psMeanFreePath.hpp b/include/viennaps/psMeanFreePath.hpp
@@ -84,22 +84,11 @@ template <class NumericType, int D> class psMeanFreePath {
         auto particleSeed = rayInternal::tea<3>(idx, seed);
         rayRNG RngState(particleSeed);
 
-        rayTriple<NumericType> origin, direction;
-        int cellIdx = -1;
-
-        if constexpr (D == 3) {
-          auto pointIdx = pointDist(RngState);
-          direction = rayReflectionDiffuse<NumericType, 3>(
-              surfaceNormals[pointIdx], RngState);
-          cellIdx = getStartingCell(surfacePoints[pointIdx]);
-          origin = cellSet->getCellCenter(cellIdx);
-        } else {
-          auto cellIdx = cellDist(RngState);
-          while (
-              !psMaterialMap::isMaterial(materialIds->at(cellIdx), material)) {
-            cellIdx = cellDist(RngState);
-          }
-        }
+        auto pointIdx = pointDist(RngState);
+        auto direction = rayReflectionDiffuse<NumericType, 3>(
+            surfaceNormals[pointIdx], RngState);
+        int cellIdx = getStartingCell(surfacePoints[pointIdx]);
+        auto origin = cellSet->getCellCenter(cellIdx);
 
         NumericType distanceOffset = 0.;