[Getfem-commits] [getfem-commits] branch master updated: a small fix on neighbor transformation

[Yves Renard]

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.Renard@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());