stardis-solver

commit e7d8bc730d00a99b034147f939c45894322455c2
parent ec2005eed6500a0909747efc97b600a69367c540
Author: Vincent Forest <vincent.forest@meso-star.com>
Date:   Thu,  2 Dec 2021 16:22:54 +0100

Update the checks on a solid/fluid boundary with a fixed flux

A fixed flux != 0 is not supported when the Picard order is greater > 1

Diffstat:
M
src/sdis_heat_path_boundary_Xd_fixed_flux.h
 | 
4
++--

1 file changed, 2 insertions(+), 2 deletions(-)
diff --git a/src/sdis_heat_path_boundary_Xd_fixed_flux.h b/src/sdis_heat_path_boundary_Xd_fixed_flux.h
@@ -44,7 +44,7 @@ XD(solid_boundary_with_flux_path)
   struct XD(reinjection_step) reinject_step = XD(REINJECTION_STEP_NULL);
 
   /* Reinjection arguments */
-  struct XD(solid_reinjection_args) solid_reinject_args = 
+  struct XD(solid_reinjection_args) solid_reinject_args =
     XD(SOLID_REINJECTION_ARGS_NULL);
 
   /* Data attached to the boundary */
@@ -69,7 +69,7 @@ XD(solid_boundary_with_flux_path)
    * radiative temperature is linearized, i.e. when the picard order is equal
    * to 1 */
   picard_order = get_picard_order(ctx);
-  if(picard_order > 1 && phi > 0) {
+  if(picard_order > 1 && phi != 0) {
     log_err(scn->dev,
       "%s: invalid flux '%g' W/m^2. Could not manage a flux != 0 when the "
       "picard order is not equal to 1; Picard order is currently set to %lu.\n",