[Getfem-commits] r5288 - in /trunk/getfem/src: getfem/getfem_fem.h getfem_fem.cc

[Yves . Renard]

Author: renard
Date: Tue Apr  5 12:47:23 2016
New Revision: 5288

URL: http://svn.gna.org/viewcvs/getfem?rev=5288&view=rev
Log:
Avoid problem with global dofs

Modified:
    trunk/getfem/src/getfem/getfem_fem.h
    trunk/getfem/src/getfem_fem.cc

Modified: trunk/getfem/src/getfem/getfem_fem.h
URL: 
http://svn.gna.org/viewcvs/getfem/trunk/getfem/src/getfem/getfem_fem.h?rev=5288&r1=5287&r2=5288&view=diff
==============================================================================
--- trunk/getfem/src/getfem/getfem_fem.h        (original)
+++ trunk/getfem/src/getfem/getfem_fem.h        Tue Apr  5 12:47:23 2016
@@ -317,7 +317,7 @@
     const base_node &node_of_dof(size_type cv, size_type i) const
       { return (*(node_tab(cv)))[i];}
     virtual const std::vector<short_type> &
-      faces_of_dof(size_type /*cv*/, size_type i) const { return face_tab[i]; }
+    faces_of_dof(size_type /*cv*/, size_type i) const;
     bool is_on_real_element(void) const { return real_element_defined; }
     bool is_equivalent(void) const { return is_equiv; }
     bool need_G(void) const

Modified: trunk/getfem/src/getfem_fem.cc
URL: 
http://svn.gna.org/viewcvs/getfem/trunk/getfem/src/getfem_fem.cc?rev=5288&r1=5287&r2=5288&view=diff
==============================================================================
--- trunk/getfem/src/getfem_fem.cc      (original)
+++ trunk/getfem/src/getfem_fem.cc      Tue Apr  5 12:47:23 2016
@@ -221,7 +221,7 @@
   { c.hess_base_value(t, withM); }
 
   /* ******************************************************************** */
-  /*        Class for description of an interpolation dof.                    
*/
+  /*        Class for description of an interpolation dof.                */
   /* ******************************************************************** */
 
   enum ddl_type { LAGRANGE, NORMAL_DERIVATIVE, DERIVATIVE, MEAN_VALUE,
@@ -513,6 +513,12 @@
     return true;
   }
 
+
+  /* ******************************************************************** */
+  /*        Members methods of virtual_fem           .                    */
+  /* ******************************************************************** */
+
+
   void virtual_fem::add_node(const pdof_description &d, const base_node &pt,
                              const dal::bit_vector &faces) {
     short_type nb = cv_node.nb_points();
@@ -548,6 +554,12 @@
   void virtual_fem::unfreeze_cvs_node(void) {
     cv_node.structure() = cvs_node;
     pspt_valid = false;
+  }
+
+  const std::vector<short_type> &
+  virtual_fem::faces_of_dof(size_type /*cv*/, size_type i) const {
+    static const std::vector<short_type> no_faces;
+    return (i < face_tab.size()) ? face_tab[i] : no_faces;
   }
 
   void virtual_fem::copy(const virtual_fem &f) {