fix maximum time gate rounding bug

fangq · Jul 25, 2019 · a3b7714 · a3b7714
1 parent eb109e0
commit a3b7714
Show file tree

Hide file tree

Showing 2 changed files with 2 additions and 2 deletions.
diff --git a/src/mmc_core.cl b/src/mmc_core.cl
@@ -489,7 +489,7 @@ float branchless_badouel_raytet(ray *r, __constant MCXParam *gcfg,__constant int
 	    r->Lmove=((r->isend) ? r->Lmove : Lmin);
 	    r->pout=r->p0+FL3(Lmin)*r->vec;
 
-	    if((int)((r->photontimer+r->Lmove*(prop.n*R_C0)-gcfg->tstart)*gcfg->Rtstep)>(int)((gcfg->tend-gcfg->tstart)*gcfg->Rtstep)){ /*exit time window*/
+	    if((int)((r->photontimer+r->Lmove*(prop.n*R_C0)-gcfg->tstart)*gcfg->Rtstep)>gcfg->maxgate-1){ /*exit time window*/
 	       r->faceid=-2;
 	       r->pout.x=MMC_UNDEFINED;
 	       r->Lmove=(gcfg->tend-r->photontimer)/(prop.n*R_C0)-1e-4f;

diff --git a/src/tettracing.c b/src/tettracing.c
@@ -969,7 +969,7 @@ float branchless_badouel_raytet(ray *r, raytracer *tracer, mcconfig *cfg, visito
 	    T = _mm_add_ps(T, S);
 	    _mm_store_ps(&(r->pout.x),T);
 
-	    if((int)((r->photontimer+r->Lmove*rc-cfg->tstart)*visit->rtstep)>=(int)((cfg->tend-cfg->tstart)*visit->rtstep)){ /*exit time window*/
+	    if((int)((r->photontimer+r->Lmove*rc-cfg->tstart)*visit->rtstep)>cfg->maxgate-1){ /*exit time window*/
 	       r->faceid=-2;
 	       r->pout.x=MMC_UNDEFINED;
 	       r->Lmove=(cfg->tend-r->photontimer)/(prop->n*R_C0)-1e-4f;