From: Guido Kanschat <dr.guido.kanschat@gmail.com>
Date: Wed, 21 Aug 2013 20:22:22 +0000 (+0000)
Subject: add missing initialize for constraints
X-Git-Tag: v8.1.0~1019
X-Git-Url: https://gitweb.dealii.org/cgi-bin/gitweb.cgi?a=commitdiff_plain;h=bc9c94e83ddd4b24f603591e99b5d61e8223452a;p=dealii.git

add missing initialize for constraints

git-svn-id: https://svn.dealii.org/trunk@30386 0785d39b-7218-0410-832d-ea1e28bc413d
---

diff --git a/deal.II/include/deal.II/meshworker/simple.h b/deal.II/include/deal.II/meshworker/simple.h
index 0a5c70c6cb..dc9d1264ad 100644
--- a/deal.II/include/deal.II/meshworker/simple.h
+++ b/deal.II/include/deal.II/meshworker/simple.h
@@ -458,6 +458,15 @@ namespace MeshWorker
        */
       void initialize(MATRIX &m, VECTOR &rhs);
 
+      /**
+       * Initialize the constraints. After this function has been
+       * called with a valid ConstraintMatrix, the function
+       * ConstraintMatrix::distribute_local_to_global() will be used
+       * by assemble() to distribute the cell and face matrices into a
+       * global sparse matrix.
+       */
+      void initialize(const ConstraintMatrix &constraints);
+
       /**
        * Initialize the local data
        * in the
@@ -1189,7 +1198,15 @@ namespace MeshWorker
       ResidualSimple<VECTOR>::initialize(data);
     }
 
+    template <class MATRIX, class VECTOR>
+    inline void
+    MatrixSimple<MATRIX>::initialize(const ConstraintMatrix &c)
+    {
+      MatrixSimple<MATRIX>::initialize(c);
+      ResidualSimple<VECTOR>::initialize(c);
+    }
 
+    
     template <class MATRIX, class VECTOR>
     template <class DOFINFO>
     inline void