stardis-solver

commit 198e24510e6bd88c0d77c017e329d2b2400f3a75
parent 83a4459b1ff420e862ea41901c953aeaa094ca9d
Author: Vincent Forest <vincent.forest@meso-star.com>
Date:   Thu, 11 Jan 2024 15:08:36 +0100

Correction to the calculation of incident diffuse external flux

A typo led to the corresponding weight being multiplied by the pdf of
the direction sampled towards the source, whereas it should be divided
by the pdf.

Diffstat:
M
src/sdis_heat_path_boundary_Xd_handle_external_net_flux.h
 | 
2
+-

1 file changed, 1 insertion(+), 1 deletion(-)
diff --git a/src/sdis_heat_path_boundary_Xd_handle_external_net_flux.h b/src/sdis_heat_path_boundary_Xd_handle_external_net_flux.h
@@ -333,7 +333,7 @@ XD(compute_incident_diffuse_flux)
       /* The source is above the surface */
       } else {
         const double Ld = XD(direct_contribution)(scn, &src_sample, pos, &hit);
-        L = Ld * cos_theta/PI * src_sample.pdf; /* [W/m^2] */
+        L = Ld * cos_theta / (PI * src_sample.pdf); /* [W/m^2] */
       }
     }
     incident_diffuse_flux += L;