fix photon sharing 0 output issue in negative patterns

fangq · Aug 7, 2023 · 9aaba97 · 9aaba97
1 parent 0800d0e
commit 9aaba97
Showing 1 changed file with 1 addition and 1 deletion.
diff --git a/src/mcx_core.cu b/src/mcx_core.cu
@@ -1986,7 +1986,7 @@ __global__ void mcx_main_loop(uint media[], OutputType field[], float genergy[],
 #endif
                         } else {
                             for (int i = 0; i < gcfg->srcnum; i++) {
-                                if (ppath[gcfg->w0offset + i] > 0.f) {
+                                if (fabs(ppath[gcfg->w0offset + i]) > 0.f) {
 #ifdef USE_DOUBLE
                                     atomicAdd(& field[(idx1dold + tshift * gcfg->dimlen.z)*gcfg->srcnum + i], weight * ppath[gcfg->w0offset + i]);
 #else