From e7c21fcf1a902ec57ea629deb4cedcc38c2a7c00 Mon Sep 17 00:00:00 2001
From: Jonathan Youett <youett@mi.fu-berlin.de>
Date: Mon, 17 Oct 2011 12:46:47 +0000
Subject: [PATCH] Make the local bases of each tangent space only count as one
 global (block) dof. This way one gets a global block structure in the
 tangential stiffness matrix. Note that now the local assemblers have to know
 this special local structure.

[[Imported from SVN: r7938]]
---
 dune/gfe/globalgfetestfunctionbasis.hh | 11 +++++------
 1 file changed, 5 insertions(+), 6 deletions(-)

diff --git a/dune/gfe/globalgfetestfunctionbasis.hh b/dune/gfe/globalgfetestfunctionbasis.hh
index 0b495d9d..15bf079f 100644
--- a/dune/gfe/globalgfetestfunctionbasis.hh
+++ b/dune/gfe/globalgfetestfunctionbasis.hh
@@ -69,19 +69,18 @@ public:
             delete(lfe_);
     }
 
-    /** \brief Get the total number of global basis functions.
-     *
-     *  For each Lagrange point there are (dim TangentSpace) basis functions
+    /** \brief Get the total number of global (block) basis functions.
+     *  Only return the number of Lagrange points. For each Lagrange point there are (tangentDim) local shape functions
      */
     int size() const
     { 
-        return basis_.size()*tangentDim;
+        return basis_.size();
     }
 
-    /** \brief Get the global index of the basis function. */
+    /** \brief Get the global index of the (block) basis function. */
     int index(const Element& e, const int i) const
     {
-        return basis_.index(e,i/tangentDim)*tangentDim + (i%tangentDim);   
+        return basis_.index(e,i);   
     }
 
 
-- 
GitLab