stardis-solver

commit 9bb7e5d12f5677ca474e603b13b63d5e21671aef
parent 89d1b58ce74c17d4a669b148d4a66d889cd62c75
Author: Vincent Forest <vincent.forest@meso-star.com>
Date:   Fri,  1 Mar 2019 11:18:16 +0100

Tune the test_sdis_volumic_power test

Reduce the delta to match the analytic solution.

Diffstat:
M
src/test_sdis_volumic_power.c
 | 
10
+++++-----

1 file changed, 5 insertions(+), 5 deletions(-)
diff --git a/src/test_sdis_volumic_power.c b/src/test_sdis_volumic_power.c
@@ -50,7 +50,7 @@
 #define T0 320
 #define LAMBDA 0.1
 #define P0 10
-#define DELTA 1.0/20.0
+#define DELTA 1.0/40.0
 
 /*******************************************************************************
  * Media
@@ -181,12 +181,12 @@ solve(struct sdis_scene* scn, const double pos[])
 
   switch(dim) {
     case SDIS_SCENE_2D:
-      printf("Temperature at (%g %g) = %g ~ %g +/- %g\n",
-        SPLIT2(pos), ref, T.E, T.SE);
+      printf("Temperature at (%g %g) = %g ~ %g +/- %g [%g, %g]\n",
+        SPLIT2(pos), ref, T.E, T.SE, T.E-3*T.SE, T.E+3*T.SE);
       break;
     case SDIS_SCENE_3D:
-      printf("Temperature at (%g %g %g) = %g ~ %g +/- %g\n",
-        SPLIT3(pos), ref, T.E, T.SE);
+      printf("Temperature at (%g %g %g) = %g ~ %g +/- %g [%g, %g]\n",
+        SPLIT3(pos), ref, T.E, T.SE, T.E -3*T.SE, T.E + 3*T.SE);
       break;
     default: FATAL("Unreachable code.\n"); break;
   }