[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]
[Getfem-commits] r5275 - /trunk/getfem/src/getfem_mesh_fem.cc
From: |
Yves . Renard |
Subject: |
[Getfem-commits] r5275 - /trunk/getfem/src/getfem_mesh_fem.cc |
Date: |
Fri, 01 Apr 2016 11:25:43 -0000 |
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();
[Prev in Thread] |
Current Thread |
[Next in Thread] |
- [Getfem-commits] r5275 - /trunk/getfem/src/getfem_mesh_fem.cc,
Yves . Renard <=