stardis-solver

Solve coupled heat transfers
git clone git://git.meso-star.fr/stardis-solver.git
Log | Files | Refs | README | LICENSE

commit 97612c490aa06adaf8bf7c58016233c7a8dacbe4
parent 34f84f757a0f562f8561da1d5525e429864b9abd
Author: Christophe Coustet <christophe.coustet@meso-star.com>
Date:   Tue, 23 Oct 2018 13:56:51 +0200

Spacing

Diffstat:
Msrc/sdis_solve_Xd.h | 12++++++------
Msrc/test_sdis_convection.c | 4++--
Msrc/test_sdis_convection_non_uniform.c | 6+++---
3 files changed, 11 insertions(+), 11 deletions(-)

diff --git a/src/sdis_solve_Xd.h b/src/sdis_solve_Xd.h @@ -241,7 +241,7 @@ sample_time { ASSERT(time_range && time_range[0] >= 0 && time_range[1] >= time_range[0] && rng); - if (time_range[0] == time_range[1]) return time_range[0]; + if(time_range[0] == time_range[1]) return time_range[0]; return ssp_rng_uniform_double(rng, time_range[0], time_range[1]); } #define SAMPLE_TIME_DEFINED @@ -1701,18 +1701,18 @@ XD(interface_prebuild_fragment) #if SDIS_SOLVE_DIMENSION == 2 #define SET_PARAM(Dest, Src) (Dest).u = (Src); - st = (float) uv[0]; + st = (float)uv[0]; #else #define SET_PARAM(Dest, Src) f2_set((Dest).uv, (Src)); f2_set_d2(st, uv); #endif res = sXd(scene_view_get_primitive(scn->sXd(view), iprim, &prim)); - if (res != RES_OK) return res; + if(res != RES_OK) return res; res = sXd(primitive_get_attrib(&prim, SXD_POSITION, st, &attr)); - if (res != RES_OK) return res; + if(res != RES_OK) return res; dX_set_fX(vtx.P, attr.value); res = sXd(primitive_get_attrib(&prim, SXD_GEOMETRY_NORMAL, st, &attr)); - if (res != RES_OK) return res; + if(res != RES_OK) return res; fX(set)(hit.normal, attr.value); hit.distance = 0; @@ -2171,7 +2171,7 @@ XD(solve_boundary_flux) ++N; } } - if (res != RES_OK) goto error; + if(res != RES_OK) goto error; setup_estimator(estimator, nrealisations, N, weight_t, sqr_weight_t); setup_estimator_flux(estimator, FLUX_CONVECTIVE__, weight_fc, sqr_weight_fc); diff --git a/src/test_sdis_convection.c b/src/test_sdis_convection.c @@ -260,7 +260,7 @@ main(int argc, char** argv) Tinf = (H*(T0 + T1 + T2 + T3 + T4 + T5)) / (6 * H); printf("Temperature of the box at (%g %g %g)\n", SPLIT3(pos)); FOR_EACH(i, 0, 5) { - double time = i ? (double) i / nu : INF; + double time = i ? (double)i / nu : INF; double time_range[2]; time_range[0] = time_range[1] = time; ref = Tf_0 * exp(-nu * time) + Tinf * (1 - exp(-nu * time)); @@ -283,7 +283,7 @@ main(int argc, char** argv) Tinf = (H * (T0 + T1 + T2 + T3)) / (4 * H); printf("Temperature of the square at (%g %g)\n", SPLIT2(pos)); FOR_EACH(i, 0, 5) { - double time = i ? (double) i / nu : INF; + double time = i ? (double)i / nu : INF; double time_range[2]; time_range[0] = time_range[1] = time; ref = Tf_0 * exp(-nu * time) + Tinf * (1 - exp(-nu * time)); diff --git a/src/test_sdis_convection_non_uniform.c b/src/test_sdis_convection_non_uniform.c @@ -275,7 +275,7 @@ main(int argc, char** argv) / (HC0 + HC1 + HC2 + HC3 + HC4 + HC5); printf("Temperature of the box at (%g %g %g)\n", SPLIT3(pos)); FOR_EACH(i, 0, 5) { - double time = i ? (double) i / nu : INF; + double time = i ? (double)i / nu : INF; double time_range[2]; time_range[0] = time_range[1] = time; ref = Tf_0 * exp(-nu * time) + Tinf * (1 - exp(-nu * time)); @@ -298,7 +298,7 @@ main(int argc, char** argv) Tinf = (HC0 * T0 + HC1 * T1 + HC2 * T2 + HC3 * T3) / (HC0 + HC1 + HC2 + HC3); printf("Temperature of the square at (%g %g)\n", SPLIT2(pos)); FOR_EACH(i, 0, 5) { - double time = i ? (double) i / nu : INF; + double time = i ? (double)i / nu : INF; double time_range[2]; time_range[0] = time_range[1] = time; ref = Tf_0 * exp(-nu * time) + Tinf * (1 - exp(-nu * time)); @@ -310,7 +310,7 @@ main(int argc, char** argv) OK(sdis_estimator_get_temperature(estimator, &T)); OK(sdis_estimator_ref_put(estimator)); printf(" t=%g : %g ~ %g +/- %g\n", time, ref, T.E, T.SE); - if (nfails) + if(nfails) printf("#failures = %lu/%lu\n", (unsigned long)nfails,(unsigned long)N); CHK(eq_eps(T.E, ref, T.SE * 3)); }