stardis-solver

commit 6e24ccd2685a5784bb6990d4689c307311c0fa5e
parent 911e0a4da0657bfaa511ad001449bfbf2a7c40f1
Author: Vincent Forest <vincent.forest@meso-star.com>
Date:   Tue, 26 Oct 2021 16:34:38 +0200

Fix the picard1 function

Diffstat:
M
src/sdis_heat_path_boundary_Xd_solid_fluid_picard1.h
 | 
10
++++++++--

1 file changed, 8 insertions(+), 2 deletions(-)
diff --git a/src/sdis_heat_path_boundary_Xd_solid_fluid_picard1.h b/src/sdis_heat_path_boundary_Xd_solid_fluid_picard1.h
@@ -173,6 +173,10 @@ XD(solid_fluid_boundary_picard1_path)
     ASSERT(fluid->type == SDIS_FLUID);
   }
 
+  /* Setup a fragment for the fluid side */
+  frag_fluid = *frag;
+  frag_fluid.side = fluid_side;
+
   /* Fetch the solid properties */
   lambda = solid_get_thermal_conductivity(solid, &rwalk->vtx);
   delta = solid_get_delta(solid, &rwalk->vtx);
@@ -262,7 +266,10 @@ XD(solid_fluid_boundary_picard1_path)
 
     /* Sample a radiative path and get the Tref at its end.
      * TODO handle the registration of the path geometry */
-    T_s = *T; rwalk_s = *rwalk;
+    T_s = *T;
+    rwalk_s = *rwalk;
+    rwalk_s.mdm = fluid;
+    rwalk_s.hit_side = fluid_side;
     res = XD(radiative_path)(scn, ctx, &rwalk_s, rng, &T_s);
     if(res != RES_OK) goto error;
 
@@ -295,7 +302,6 @@ XD(solid_fluid_boundary_picard1_path)
     /* Null-collision, looping at the beginning */
   }
 
-
 exit:
   return res;
 error: