getfem-commits
[Top][All Lists]
Advanced

[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]

[Getfem-commits] r4943 - in /trunk/getfem/src: bgeot_node_tab.cc getfem_


From: Yves . Renard
Subject: [Getfem-commits] r4943 - in /trunk/getfem/src: bgeot_node_tab.cc getfem_mesh.cc
Date: Wed, 08 Apr 2015 16:02:35 -0000

Author: renard
Date: Wed Apr  8 18:02:35 2015
New Revision: 4943

URL: http://svn.gna.org/viewcvs/getfem?rev=4943&view=rev
Log:
fixed a small bug

Modified:
    trunk/getfem/src/bgeot_node_tab.cc
    trunk/getfem/src/getfem_mesh.cc

Modified: trunk/getfem/src/bgeot_node_tab.cc
URL: 
http://svn.gna.org/viewcvs/getfem/trunk/getfem/src/bgeot_node_tab.cc?rev=4943&r1=4942&r2=4943&view=diff
==============================================================================
--- trunk/getfem/src/bgeot_node_tab.cc  (original)
+++ trunk/getfem/src/bgeot_node_tab.cc  Wed Apr  8 18:02:35 2015
@@ -128,15 +128,18 @@
       for (size_type is = 0; is < sorters.size(); ++is) {
        if (existi) sorters[is].insert(j);
        if (existj) sorters[is].insert(i);
+        GMM_ASSERT3(sorters[is].size() == card(), "internal error");
       }
     }
   }
 
   void node_tab::sup_node(size_type i) {
+    resort();
     if (index().is_in(i)) {
       for (size_type is = 0; is < sorters.size(); ++is) {
        sorters[is].erase(i);
-       GMM_ASSERT3(sorters[is].size()+1 == card(), "internal error");
+        GMM_ASSERT3(sorters[is].size()+1 == card(), "Internal error");
+        // if (sorters[is].size()+1 != card()) { resort(); }
       }
       dal::dynamic_tas<base_node>::sup(i);
       

Modified: trunk/getfem/src/getfem_mesh.cc
URL: 
http://svn.gna.org/viewcvs/getfem/trunk/getfem/src/getfem_mesh.cc?rev=4943&r1=4942&r2=4943&view=diff
==============================================================================
--- trunk/getfem/src/getfem_mesh.cc     (original)
+++ trunk/getfem/src/getfem_mesh.cc     Wed Apr  8 18:02:35 2015
@@ -207,6 +207,7 @@
 #endif
 
   void mesh::optimize_structure() {
+    pts.resort();
     size_type i, j;
     j = nb_convex();
     for (i = 0; i < j; i++)




reply via email to

[Prev in Thread] Current Thread [Next in Thread]