stardis-solver

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

commit af831a26c8b2a3f72c082aab2cb4e726348fcf39
parent bf95ef4dc80dab2c44ca49cdfa6b156e9a0c242b
Author: Vincent Forest <vincent.forest@meso-star.com>
Date:   Fri,  1 Mar 2019 09:21:20 +0100

Small adjustment of the scene_get_medium_<2d|3d> functions

Diffstat:
Msrc/sdis_scene_Xd.h | 5+++--
1 file changed, 3 insertions(+), 2 deletions(-)

diff --git a/src/sdis_scene_Xd.h b/src/sdis_scene_Xd.h @@ -793,7 +793,7 @@ error: static INLINE res_T XD(scene_get_medium) (const struct sdis_scene* scn, - const double pos[2], + const double pos[DIM], struct get_medium_info* info, /* May be NULL */ struct sdis_medium** out_medium) { @@ -862,7 +862,8 @@ XD(scene_get_medium) fX(normalize)(N, hit.normal); cos_N_dir = fX(dot)(N, dir); - if(absf(cos_N_dir) > 1.e-1f) { /* Not roughly orthognonal */ + /* Not too close and not roughly orthognonal */ + if(hit.distance > 1.e-6 || absf(cos_N_dir) > 1.e-1f) { const struct sdis_interface* interf; interf = scene_get_interface(scn, hit.prim.prim_id); medium = interface_get_medium