stardis-solver

commit aee9f6c7e4b84062e127cf93125839450eb59cb5
parent df0ff90c75cb418873c3303c70d6e15a9525d8bf
Author: Vincent Forest <vincent.forest@meso-star.com>
Date:   Fri, 13 Apr 2018 11:35:07 +0200

Fix the fluid temperature

The observation time was always set to 0

Diffstat:
M
src/sdis_solve_Xd.h
 | 
2
+-

1 file changed, 1 insertion(+), 1 deletion(-)
diff --git a/src/sdis_solve_Xd.h b/src/sdis_solve_Xd.h
@@ -442,7 +442,7 @@ XD(fluid_temperature)
    * FIXME we assume that hc is constant for the whole enclosure */
   mu = hc / (rho * cp) * enc->S_over_V;
   tau = ssp_ran_exp(rng, mu);
-  rwalk->vtx.time = MMIN(rwalk->vtx.time - tau, 0);
+  rwalk->vtx.time = MMAX(rwalk->vtx.time - tau, 0);
 
   /* Check the initial condition */
   tmp = fluid_get_temperature(rwalk->mdm, &rwalk->vtx);