Workaround: Add compatibility instantiations

diff --git a/source/lac/affine_constraints.cc b/source/lac/affine_constraints.cc
index aef54fb03aae10fb76009256a53af2c3a3185c87..39c51d7dbcc0e27782ca27c185d463ca10169802 100644 (file)
--- a/source/lac/affine_constraints.cc
+++ b/source/lac/affine_constraints.cc
@@ -21,17 +21,6 @@ DEAL_II_NAMESPACE_OPEN
 
 #include "affine_constraints.inst"
 
-/*
- * FIXME: This mixed variant is needed by the multigrid - it currently only
- * supports AffineConstraints<double>. If someone has the patience to
- * templateify the 10 custom wrapper classes around AffineConstraints,
- * please go ahead, fix the multigrid classes and remove this
- * instantiation:
- */
-template void
-dealii::AffineConstraints<double>::condense<dealii::BlockVector<float>>(
-  dealii::BlockVector<float> &) const;
-
 /*
  * Note: You probably do not want to add your custom instantiation down
  * here but use affine_constraints.inst.in instead. We use the following

diff --git a/source/lac/affine_constraints.inst.in b/source/lac/affine_constraints.inst.in
index 5d47c6796dc2213c84053a3f9aa207450730ac2c..ad54c969d00bb11f439b46f359775bd7f5eae0a4 100644 (file)
--- a/source/lac/affine_constraints.inst.in
+++ b/source/lac/affine_constraints.inst.in
@@ -336,3 +336,16 @@ for (V : EXTERNAL_PARALLEL_VECTORS)
 
     template void AffineConstraints<V::value_type>::set_zero<V>(V &) const;
   }
+
+
+//
+// FIXME: These mixed variants are needed for multigrid and matrix free.
+//
+
+for (T : DEAL_II_VEC_TEMPLATES)
+  {
+    template void dealii::AffineConstraints<double>::condense<T<float>>(
+      T<float> &) const;
+    template void dealii::AffineConstraints<double>::distribute<T<float>>(
+      T<float> &) const;
+  }