stardis

commit e139896e697204d0d6004f8d7fd1bc87588bee15
parent 9cb975b5a5a042d780bb7e43ac450111d966ce0c
Author: Christophe Coustet <christophe.coustet@meso-star.com>
Date:   Mon, 16 Sep 2019 18:43:51 +0200

Fix output for flux at boundary compute

Diffstat:
M
src/stardis-output.c
 | 
28
++++++++++++++++++++++++++--

1 file changed, 26 insertions(+), 2 deletions(-)
diff --git a/src/stardis-output.c b/src/stardis-output.c
@@ -694,12 +694,36 @@ print_single_MC_result
       result.E, /* Expected value */
       result.SE); /* Standard error */
     break;
-  case FLUX_BOUNDARY_COMPUTE:
-    printf("Flux at boundary %s at t=%g = %g +/- %g\n",
+  case FLUX_BOUNDARY_COMPUTE: {
+    enum sdis_estimator_type type;
+    sdis_estimator_get_type(estimator, &type);
+    if (res != RES_OK) goto error;
+    ASSERT(type == SDIS_ESTIMATOR_FLUX);
+
+    printf("Temperature at boundary %s at t=%g = %g +/- %g\n",
+      stardis->solve_name, stardis->probe[3],
+      result.E, /* Expected value */
+      result.SE); /* Standard error */
+    res = sdis_estimator_get_convective_flux(estimator, &result);
+    if (res != RES_OK) goto error;
+    printf("Convective flux at boundary %s at t=%g = %g +/- %g\n",
+      stardis->solve_name, stardis->probe[3],
+      result.E, /* Expected value */
+      result.SE); /* Standard error */
+    res = sdis_estimator_get_radiative_flux(estimator, &result);
+    if (res != RES_OK) goto error;
+    printf("Radiative flux at boundary %s at t=%g = %g +/- %g\n",
+      stardis->solve_name, stardis->probe[3],
+      result.E, /* Expected value */
+      result.SE); /* Standard error */
+    res = sdis_estimator_get_total_flux(estimator, &result);
+    if (res != RES_OK) goto error;
+    printf("Total flux Flux at boundary %s at t=%g = %g +/- %g\n",
       stardis->solve_name, stardis->probe[3],
       result.E, /* Expected value */
       result.SE); /* Standard error */
     break;
+  }
   default: FATAL("Invalid mode.");
   }
   printf("#failures: %zu/%zu\n", nfailures, stardis->N);