fix nan in detected photon data when using hyperboloid src, see https…

…://groups.google.com/g/mcx-users/c/vyfHbzEO-0M/m/JzlpEZ3TBwAJ
fangq · Oct 5, 2022 · 63ffc1e · 63ffc1e
1 parent 10330ef
commit 63ffc1e
Showing 1 changed file with 1 addition and 1 deletion.
diff --git a/src/mcx_core.cu b/src/mcx_core.cu
@@ -1295,7 +1295,7 @@ __device__ inline int launchnewphoton(MCXpos* p, MCXdir* v, Stokes* s, MCXtime*
                     float phi = TWO_PI * rand_uniform01(t);
                     sincosf(phi, &sphi, &cphi);
 
-                    float r = sqrtf(-0.5f * logf(rand_uniform01(t))) * gcfg->srcparam1.x;
+                    float r = sqrtf(0.5f * rand_next_scatlen(t)) * gcfg->srcparam1.x;
 
                     /** parameter to generate photon path from coordinates at focus (depends on focal distance and rayleigh range) */
                     float tt = -gcfg->srcparam1.y / gcfg->srcparam1.z;