htrdr

Solving radiative transfer in heterogeneous media
git clone git://git.meso-star.fr/htrdr.git
Log | Files | Refs | README | LICENSE

commit e1bac4ab9dae78d16590fbf31201c25215843ba5
parent bb4e30ac937947a8f7c2ea1b097c8bd448118167
Author: Vincent Forest <vincent.forest@meso-star.com>
Date:   Tue, 24 Aug 2021 11:57:06 +0200

Fix an issue in combustion_compute_radiance_sw

At each scattering position, the range of the traced ray could be wrong.

Diffstat:
Msrc/combustion/htrdr_combustion_compute_radiance_sw.c | 3+--
1 file changed, 1 insertion(+), 2 deletions(-)

diff --git a/src/combustion/htrdr_combustion_compute_radiance_sw.c b/src/combustion/htrdr_combustion_compute_radiance_sw.c @@ -878,7 +878,6 @@ combustion_compute_radiance_sw d3_set(pos, pos_in); d3_set(dir, dir_in); - d2(range, 0, FLT_MAX); wlen = htrdr_combustion_laser_get_wavelength(cmd->laser); @@ -904,7 +903,7 @@ combustion_compute_radiance_sw /* Find the intersection with the combustion chamber geometry */ d3_set(rt_args.ray_org, pos); d3_set(rt_args.ray_dir, dir); - d2_set(rt_args.ray_range, range); + d2(rt_args.ray_range, 0, DBL_MAX); rt_args.hit_from = hit_prev; /* Avoid self intersection */ HTRDR(geometry_trace_ray(cmd->geom, &rt_args, &hit_curr));