stardis-solver

commit 402f34f44e40fd351ae14469f8834084d3b749b0
parent 2a9536eb711b57ca4a7cd6a0aec57ff2ea777555
Author: Vincent Forest <vincent.forest@meso-star.com>
Date:   Mon, 28 May 2018 09:45:34 +0200

Fix the 2D reinjection pattern

Diffstat:
M
src/sdis_solve_Xd.h
 | 
2
+-
M
src/test_sdis_volumic_power3_2d.c
 | 
5
+++--

2 files changed, 4 insertions(+), 3 deletions(-)
diff --git a/src/sdis_solve_Xd.h b/src/sdis_solve_Xd.h
@@ -525,7 +525,7 @@ XD(solid_solid_boundary_temperature)
     if(power != SDIS_VOLUMIC_POWER_NONE) {
       const double delta_in_meter = delta_boundary * fp_to_meter;
       const double lambda = solid_get_thermal_conductivity(mdm, &rwalk->vtx);
-      tmp = power * delta_in_meter * delta_in_meter / (2.0 * lambda);
+      tmp = power * delta_in_meter * delta_in_meter / (2.0 * DIM * lambda);
       T->value += tmp;
     }
 
diff --git a/src/test_sdis_volumic_power3_2d.c b/src/test_sdis_volumic_power3_2d.c
@@ -25,6 +25,8 @@
 #define T2 273.15
 #define H1 5.0
 #define H2 10.0
+#define MDb 1.0
+#define N 10000 /* #realisations */
 
 /*
  * The 2D scene is composed of 3 stacked solid slabs whose middle slab has a
@@ -155,7 +157,7 @@ solid_get_delta_boundary
   (const struct sdis_rwalk_vertex* vtx, struct sdis_data* data)
 {
   CHK(data != NULL && vtx != NULL);
-  return ((const struct solid*)sdis_data_cget(data))->delta * 2.1;
+  return ((const struct solid*)sdis_data_cget(data))->delta * MDb;
 }
 
 static double
@@ -246,7 +248,6 @@ main(int argc, char** argv)
   struct sdis_interface* interfaces[10/*#segment*/];
   struct sdis_mc T = SDIS_MC_NULL;
   double pos[2];
-  const size_t N = 10000;
   size_t i;
   (void)argc, (void)argv;