stardis-solver

commit d8bfb43cfc81d3e6d53dda2e299caa7e11dada0c
parent 4163547be0d54752613c42feac7eb23c0435fb27
Author: Vincent Forest <vincent.forest@meso-star.com>
Date:   Fri, 29 Mar 2024 16:01:25 +0100

Add geometry dump to unsteady 2D analytical test

Dump the super shape and sampled paths. The latter cheks path
registration.

Diffstat:
M
src/test_sdis_unsteady_analytic_profile_2d.c
 | 
46
+++++++++++++++++++++++++++++++++++++++++++++-

1 file changed, 45 insertions(+), 1 deletion(-)
diff --git a/src/test_sdis_unsteady_analytic_profile_2d.c b/src/test_sdis_unsteady_analytic_profile_2d.c
@@ -61,6 +61,32 @@ static double temperature(const double pos[2], const double time)
   return 280.0;
 }
 
+static void
+dump_heat_path
+  (FILE* fp,
+   struct sdis_scene* scn,
+   const enum sdis_diffusion_algorithm diff_algo,
+   const double pos[2],
+   const double time,
+   const size_t npaths)
+{
+  struct sdis_solve_probe_args args = SDIS_SOLVE_PROBE_ARGS_DEFAULT;
+  struct sdis_estimator* estimator = NULL;
+
+  args.nrealisations = npaths;
+  args.position[0] = pos[0];
+  args.position[1] = pos[1];
+  args.time_range[0] = time;
+  args.time_range[1] = time;
+  args.diff_algo = diff_algo;
+  args.register_paths = SDIS_HEAT_PATH_ALL;
+  OK(sdis_solve_probe(scn, &args, &estimator));
+
+  dump_heat_paths(fp, estimator);
+
+  OK(sdis_estimator_ref_put(estimator));
+}
+
 /*******************************************************************************
  * Solid, i.e. medium of the super shape
  ******************************************************************************/
@@ -270,7 +296,7 @@ check_probe
    const double time, /* [s] */
    const int green)
 {
-  struct sdis_solve_probe_args args  = SDIS_SOLVE_PROBE_ARGS_DEFAULT;
+  struct sdis_solve_probe_args args = SDIS_SOLVE_PROBE_ARGS_DEFAULT;
   struct sdis_mc T = SDIS_MC_NULL;
   struct sdis_estimator* estimator = NULL;
   double ref = 0;
@@ -314,6 +340,8 @@ main(int argc, char** argv)
   struct sdis_medium* dummy = NULL; /* Medium surrounding the solid */
   struct sdis_scene* scn = NULL;
 
+  /* Miscellaneous */
+  FILE* fp = NULL;
   struct super_shape sshape = SUPER_SHAPE_NULL;
   const double pos[3] = {0.2,0.3}; /* [m/fp_to_meter] */
   const double time = 5; /* [s] */
@@ -323,6 +351,12 @@ main(int argc, char** argv)
 
   sshape = create_super_shape();
 
+  /* Save the super shape geometry for debug and visualisation */
+  CHK(fp = fopen("super_shape_2d.obj", "w"));
+  dump_segments(fp, sshape.positions, super_shape_nvertices(&sshape),
+    sshape.indices, super_shape_nsegments(&sshape));
+  CHK(fclose(fp) == 0);
+
   solid = create_solid(sdis);
   dummy = create_dummy(sdis);
   interf = create_interface(sdis, solid, dummy);
@@ -332,6 +366,16 @@ main(int argc, char** argv)
   check_probe(scn, SDIS_DIFFUSION_WOS, pos, time, 0/*green*/);
   check_probe(scn, SDIS_DIFFUSION_WOS, pos, time, 1/*green*/);
 
+  /* Write 10 heat paths sampled by the delta sphere algorithm */
+  CHK(fp = fopen("delta_sphere_2d.vtk", "w"));
+  dump_heat_path(fp, scn, SDIS_DIFFUSION_DELTA_SPHERE, pos, time, 10);
+  CHK(fclose(fp) == 0);
+
+  /* Write 10 heat paths sampled by the WoS algorithm */
+  CHK(fp = fopen("wos_2d.vtk", "w"));
+  dump_heat_path(fp, scn, SDIS_DIFFUSION_WOS, pos, time, 10);
+  CHK(fclose(fp) == 0);
+
   release_super_shape(&sshape);
   OK(sdis_device_ref_put(sdis));
   OK(sdis_interface_ref_put(interf));