This is an automated email from the git hooks/post-receive script. renard pushed a commit to branch master in repository getfem.
The following commit(s) were added to refs/heads/master by this push: new 0f76b55 a small fix on neighbor transformation 0f76b55 is described below commit 0f76b552495b5112fab430da3b15b90ab1c7ffb3 Author: Yves Renard <yves.ren...@insa-lyon.fr> AuthorDate: Tue Oct 27 16:32:26 2020 +0100 a small fix on neighbor transformation --- src/getfem_generic_assembly_compile_and_exec.cc | 6 ++++-- 1 file changed, 4 insertions(+), 2 deletions(-) diff --git a/src/getfem_generic_assembly_compile_and_exec.cc b/src/getfem_generic_assembly_compile_and_exec.cc index 6eb3095..56ffbd1 100644 --- a/src/getfem_generic_assembly_compile_and_exec.cc +++ b/src/getfem_generic_assembly_compile_and_exec.cc @@ -4002,6 +4002,8 @@ namespace getfem { short_type f = ctx.face_num(); auto adj_face = m.adjacent_face(cv, f); if (adj_face.cv == size_type(-1)) { + GMM_WARNING2("Adjacent face not found, " + "probably an non-interior face"); inin.ctx.invalid_convex_num(); } else { gauss_pt_corresp gpc; @@ -4062,7 +4064,7 @@ namespace getfem { } } - if (inin.ctx.have_pgp()) { + if (inin.ctx.have_pgp() && inin.ctx.is_convex_num_valid()) { inin.ctx.set_ii(ipt); inin.pt_type = 1; inin.has_ctx = true; @@ -8266,7 +8268,7 @@ namespace getfem { } if (pspt != old_pspt) { first_gp = true; old_pspt = pspt; } if (pspt->size()) { - // iterations on Gauss points + // Iterations on Gauss points size_type first_ind = 0; if (v.f() != short_type(-1)) { gis.nbpt = pai->nb_points_on_face(v.f());