stardis-solver

commit 89dd2f49aa916fadf28098473ecb5f8570bba466
parent eac0546a719853431772f8549c931d43dd1c6ca3
Author: Vincent Forest <vincent.forest@meso-star.com>
Date:   Fri, 19 Jan 2024 14:20:30 +0100

Correct the way numerical problems are handled during WoS

There was a problem in the ray tracing used to fix the diffusion
position at the interface when the inconsistency of the support was
detected. The origin of the ray traced was not the actual position, but
the invalid position.

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

1 file changed, 2 insertions(+), 2 deletions(-)
diff --git a/src/sdis_heat_path_conductive_wos_Xd.h b/src/sdis_heat_path_conductive_wos_Xd.h
@@ -233,7 +233,7 @@ XD(sample_next_position)
       float rt_dir[DIM] = {0};
       float rt_range[2] = {0, 0};
 
-      fX_set_dX(rt_pos, pos);
+      fX_set_dX(rt_pos, rwalk->vtx.P);
       fX_set_dX(rt_dir, dir);
       rt_range[0] = 0;
       rt_range[1] = (float)INF;
@@ -256,7 +256,7 @@ XD(sample_next_position)
         goto error;
       }
 
-      res = XD(setup_hit_rt)(scn, pos, dir, &hit, rwalk);
+      res = XD(setup_hit_rt)(scn, rwalk->vtx.P, dir, &hit, rwalk);
       if(res != RES_OK) goto error;
     }
   }