stardis-solver

commit 312238434d4392b2dbba78bf75ce48f6cff805d6
parent f7d4355e899579b82cef12cfda227281ea90ce1d
Author: Vincent Forest <vincent.forest@meso-star.com>
Date:   Thu,  4 Nov 2021 11:02:09 +0100

Fix the check of Tref in picard1

In picard1 when the path goes to the infinity we asserted that the
returned value was the radiative ambient temperature. It was effectively
the case excepted with a volumic power or a fixed flux: in this case the
Monte-Carlo weight stores alse the volumic power and/or the fixed flux
also.

In picard1, when the path goes to infinity, we asserted that the value
returned was the radiative ambient temperature. This was indeed the case
except with a volumic power or an imposed flow rate: in this case the
Monte-Carlo weight also stores the volumic power and/or the flux.

Diffstat:
M
src/sdis_heat_path_boundary_Xd_solid_fluid_picard1.h
 | 
1
-

1 file changed, 0 insertions(+), 1 deletion(-)
diff --git a/src/sdis_heat_path_boundary_Xd_solid_fluid_picard1.h b/src/sdis_heat_path_boundary_Xd_solid_fluid_picard1.h
@@ -78,7 +78,6 @@ XD(rwalk_get_Tref)
      * fetches the ambient radiative temperature. We do not use the limit
      * conditions as the reference temperature to make the sampled paths
      * independant of them. */
-    ASSERT(T->value == scn->trad.temperature);
     Tref = scn->trad.reference;
   } else {
     struct sdis_interface_fragment frag;