htrdr

commit fbcc6acaec441675842417f397999cf8f458d3ec
parent 7f857c562740aeec155d3023835c50dcc77f692a
Author: Vincent EYMET <vincent.eymet@meso-star.com>
Date:   Wed, 19 May 2021 10:45:45 +0200

Take into account null scattering properties field in htrdr-combustion

Diffstat:
M
src/combustion/htrdr_combustion_compute_radiance_sw.c
 | 
11
+++++++++--

1 file changed, 9 insertions(+), 2 deletions(-)
diff --git a/src/combustion/htrdr_combustion_compute_radiance_sw.c b/src/combustion/htrdr_combustion_compute_radiance_sw.c
@@ -354,12 +354,19 @@ sample_scattering_limited_hit_filter
    * current ray enters into the current voxel */
   traversal_dst = hit->distance[0];
 
-  /* First traversed leaf */
-  if(ctx->sampled_vox_collision_dst < 0) {
+  /* Tmax was not already computed */
+  if(ctx->Tmax < 0) {
     ctx->Tmax = (range[1] - range[0]) * ctx->ks_2hat;
     ctx->Ume = 1 - exp(-ctx->Tmax);
   }
 
+  /* No scattering into the whole traversed laser sheet */
+  if(ctx->Tmax == 0) {
+    pursue_traversal = 1;
+    return pursue_traversal;
+  }
+  ASSERT(ctx->Tmax > 0);
+
   for(;;) {
     atrstm_radcoefs_T radcoefs;
     double vox_dst;