stardis-solver

commit 927016507d8c894ba7c1fe1512dcbcc833702090
parent 6bdba7dfd02d1c0daabf16d3b2dfd34ad22dd5ff
Author: Vincent Forest <vincent.forest@meso-star.com>
Date:   Wed, 30 May 2018 16:09:03 +0200

Small update of the volumic_power2_2d test

Diffstat:
M
src/test_sdis_volumic_power2_2d.c
 | 
7
+++++--

1 file changed, 5 insertions(+), 2 deletions(-)
diff --git a/src/test_sdis_volumic_power2_2d.c b/src/test_sdis_volumic_power2_2d.c
@@ -21,15 +21,17 @@
 #define Pw 10000 /* Volumic power */
 #define NONE -1
 
-/* H delta T. Expected temperature 286.83 C */
+/* H delta T */
 #define Tboundary1 NONE
 #define Tboundary2 NONE
 #define DELTA 0.01
+#define Tref 286.83 /* In celcius. Computed by EDF with Syrthes */
 
-/* Dirichlets. Expected temperature 246.93 C */
+/* Dirichlets */
 /*#define Tboundary1 373.15*/
 /*#define Tboundary2 273.15*/
 /*#define DELTA 0.01*/
+/*#define Tref 246.93*/ /* In celcius. Computed by EDF with Syrthes */
 
 /*
  *           _\  T1
@@ -378,6 +380,7 @@ main(int argc, char** argv)
   printf("Temperature at (%g %g) = %g +/- %g\n", SPLIT2(pos), T.E-273.15, T.SE);
   printf("#realisations: %lu; #failures: %lu\n",
     (unsigned long)nreals, (unsigned long)nfails);
+  CHK(eq_eps(T.E-273.15, Tref, T.SE*3));
   CHK(sdis_estimator_ref_put(estimator) == RES_OK);
 
   CHK(sdis_scene_ref_put(scn) == RES_OK);