stardis-solver

commit 996d5238dc4575cc58024b4f6c96f57e8927e7d7
parent 4873f81ef8f619e42c475de22254e9a51f0b1c50
Author: Christophe Coustet <christophe.coustet@meso-star.com>
Date:   Wed, 30 Sep 2020 11:53:45 +0200

Fix a test that had false expectations and should fail

Diffstat:
M
src/test_sdis_solve_boundary_flux.c
 | 
24
++++++------------------

1 file changed, 6 insertions(+), 18 deletions(-)
diff --git a/src/test_sdis_solve_boundary_flux.c b/src/test_sdis_solve_boundary_flux.c
@@ -447,29 +447,17 @@ main(int argc, char** argv)
   printf("Average values of the right side of the square = ");
   check_estimator(estimator, N, analyticT, analyticCF, analyticRF, analyticTF);
   OK(sdis_estimator_ref_put(estimator));
-  
-  /* Average temperature on the left side of the box */
+
+  /* Flux computation on Dirichlet boundaries is not available yet.
+   * Once available, the expected total flux is the same we expect on the right
+   * side (as the other sides are adiabatic). */
   prims[0] = 2;
   prims[1] = 3;
-
-  analyticT = Tb;
-  analyticCF = H * (analyticT - Tf);
-  analyticRF = Hrad * (analyticT - Trad);
-  analyticTF = analyticCF + analyticRF;
-
   bound_args.nprimitives = 2;
-  OK(SOLVE(box_scn, &bound_args, &estimator));
-  printf("Average values of the left side of the box = ");
-  check_estimator(estimator, N, analyticT, analyticCF, analyticRF, analyticTF);
-  OK(sdis_estimator_ref_put(estimator));
-
-  /* Average temperature on the left/right side of the square */
+  BA(SOLVE(box_scn, &bound_args, &estimator));
   prims[0] = 1;
   bound_args.nprimitives = 1;
-  OK(SOLVE(square_scn, &bound_args, &estimator));
-  printf("Average values of the left side of the square = ");
-  check_estimator(estimator, N, analyticT, analyticCF, analyticRF, analyticTF);
-  OK(sdis_estimator_ref_put(estimator));
+  BA(SOLVE(square_scn, &bound_args, &estimator));
   #undef SOLVE
 
   OK(sdis_scene_ref_put(box_scn));