stardis-solver

commit 98097d811e42e6e1e0cbe21a06fe9c2baddd1822
parent 6647d8884babcb1f9e41e3959e7250527d57f0c9
Author: Vincent Forest <vincent.forest@meso-star.com>
Date:   Wed, 30 May 2018 10:56:01 +0200

Minor adjustments on unit tests

Diffstat:
M
src/test_sdis_conducto_radiative_2d.c
 | 
2
+-
M
src/test_sdis_solve_probe2.c
 | 
2
+-
M
src/test_sdis_solve_probe_boundary.c
 | 
2
+-
M
src/test_sdis_volumic_power.c
 | 
2
+-
M
src/test_sdis_volumic_power2_2d.c
 | 
14
+++++++-------

5 files changed, 11 insertions(+), 11 deletions(-)
diff --git a/src/test_sdis_conducto_radiative_2d.c b/src/test_sdis_conducto_radiative_2d.c
@@ -417,7 +417,7 @@ main(int argc, char** argv)
     printf("Temperature at (%g, %g)  = %g ~ %g +/- %g\n",
       SPLIT2(pos), ref, T.E, T.SE);
 
-    CHK(eq_eps(T.E, ref, 2*T.SE) == 1);
+    CHK(eq_eps(T.E, ref, 3*T.SE) == 1);
 
     CHK(sdis_estimator_ref_put(estimator) == RES_OK);
   }
diff --git a/src/test_sdis_solve_probe2.c b/src/test_sdis_solve_probe2.c
@@ -265,7 +265,7 @@ main(int argc, char** argv)
 
   /* Check the results */
   CHK(nfails + nreals == N);
-  CHK(eq_eps(T.E, ref, T.SE));
+  CHK(eq_eps(T.E, ref, 3*T.SE));
 
   /* Release data */
   CHK(sdis_estimator_ref_put(estimator) == RES_OK);
diff --git a/src/test_sdis_solve_probe_boundary.c b/src/test_sdis_solve_probe_boundary.c
@@ -366,7 +366,7 @@ main(int argc, char** argv)
   printf("Boundary temperature of the square at (%g %g) = %g ~ %g +/- %g\n",
     SPLIT2(pos), ref, T.E, T.SE);
   printf("#failures = %lu/%lu\n", (unsigned long)nfails, (unsigned long)N);
-  CHK(eq_eps(T.E, ref, T.SE));
+  CHK(eq_eps(T.E, ref, 3*T.SE));
   #undef SOLVE
 
   CHK(sdis_scene_ref_put(box_scn) == RES_OK);
diff --git a/src/test_sdis_volumic_power.c b/src/test_sdis_volumic_power.c
@@ -324,7 +324,7 @@ main(int argc, char** argv)
   printf("Temperature of the box at (%g %g %g) = %g ~ %g +/- %g\n",
     SPLIT3(pos), ref, T.E, T.SE);
   printf("#failures = %lu/%lu\n", (unsigned long)nfails, (unsigned long)N);
-  CHK(eq_eps(T.E, ref, T.SE*2));
+  CHK(eq_eps(T.E, ref, 3*T.SE));
 
   /* Solve in 2D */
   CHK(sdis_solve_probe(square_scn, N, pos, INF, 1.0, 0, 0, &estimator) == RES_OK);
diff --git a/src/test_sdis_volumic_power2_2d.c b/src/test_sdis_volumic_power2_2d.c
@@ -17,27 +17,27 @@
 #include "test_sdis_utils.h"
 #include <rsys/math.h>
 
-#define N 40000 /* #realisations */
+#define N 10000 /* #realisations */
 #define Pw 10000 /* Volumic power */
 #define MDb sqrt(2.0) /* Multiplier applied to delta to define delta_boundary */
 #define NONE -1
 
 /* H delta T: expected 286.83 C */
-/*#define Tboundary1 NONE
-#define Tboundary2 NONE*/
+#define Tboundary1 NONE
+#define Tboundary2 NONE
 /* Db = 2.1*D:           339.091 +/- 2.734920; #failures: 0
  * Db = 0.5*D:           286.182 +/- 2.202150; #failures: 0
  * Db = 0.5*D 100K real: 283.29  +/- 0.688329; #failures: 1
  * Db = 2.1*D wo Pw:     69.12   +/- 0.461998; #failures: 0
  * Db = 0.5*D wo Pw:     69.764  +/- 0.459326; #failures: 0 */
-/*#define DELTA 0.01*/
+#define DELTA 0.005
 
 /* Dirichlets: expected 246.93 C */
-#define Tboundary1 373.15
-#define Tboundary2 273.15
+/*#define Tboundary1 373.15*/
+/*#define Tboundary2 273.15*/
 /* Db = 2.1*D: 288.487 +/- 2.17576; #failures: 0
  * Db = 0.5*D: 244.382 +/- 1.71547; #failures: 1 */
-#define DELTA 0.01
+/*#define DELTA 0.01*/
 /* Db = 2.1*D: 273.532 +/- 1.98965; #failures: 0
  * Db = 0.5*D: 243.219 +/- 1.70725: #failures: 1 */
 /*#define DELTA 0.005*/