[Getfem-commits] r5275 - /trunk/getfem/src/getfem_mesh_fem.cc

[Yves . Renard]

Author: renard
Date: Fri Apr  1 13:25:43 2016
New Revision: 5275

URL: http://svn.gna.org/viewcvs/getfem?rev=5275&view=rev
Log:
minor improvement

Modified:
    trunk/getfem/src/getfem_mesh_fem.cc

Modified: trunk/getfem/src/getfem_mesh_fem.cc
URL: 
http://svn.gna.org/viewcvs/getfem/trunk/getfem/src/getfem_mesh_fem.cc?rev=5275&r1=5274&r2=5275&view=diff
==============================================================================
--- trunk/getfem/src/getfem_mesh_fem.cc (original)
+++ trunk/getfem/src/getfem_mesh_fem.cc Fri Apr  1 13:25:43 2016
@@ -360,17 +360,19 @@
         } else if (!dof_linkable(fd.pnd)) { // If the dof is not linkable
           itab[i] = nbdof;
           nbdof += Qdim / pf->target_dim();
-        } else {                            // for a standard linkable dof
+        } else {                            // For a standard linkable dof
           pgp->transform(linked_mesh().points_of_convex(cv), i, P);
          size_type idof = nbdof;
-         linked_mesh().neighbours_of_convex(cv, ftab[i], s);
+         if (ftab[i].size()) {
+           linked_mesh().neighbours_of_convex(cv, ftab[i], s);
          
-         for (size_type ncv : s) { // For each neighbour
-                                   // control if the dof already exists.
-           fd.ind_node = dof_nodes[ncv].search_node(P);
-           if (fd.ind_node != size_type(-1)) {
-             auto it = dof_sorts[ncv].find(fd);
-             if (it != dof_sorts[ncv].end()) { idof = it->second; break; }
+           for (size_type ncv : s) { // For each neighbour
+                                     // control if the dof already exists.
+             fd.ind_node = dof_nodes[ncv].search_node(P);
+             if (fd.ind_node != size_type(-1)) {
+               auto it = dof_sorts[ncv].find(fd);
+               if (it != dof_sorts[ncv].end()) { idof = it->second; break; }
+             }
            }
          }
          if (idof == nbdof) nbdof += Qdim / pf->target_dim();