Made fe_q a public static member of MappingManifold.

diff --git a/include/deal.II/fe/mapping_manifold.h b/include/deal.II/fe/mapping_manifold.h
index 574a641aff907f4ae6e1cee9b0b2c9a9ab7bbb9f..0a6e712f8c3b347e7a16fb21407c2f5ba38f482d 100644 (file)
--- a/include/deal.II/fe/mapping_manifold.h
+++ b/include/deal.II/fe/mapping_manifold.h
@@ -304,11 +304,6 @@ public:
      */
     mutable std::vector<double> volume_elements;
 
-    /**
-     * A Q1 Finite element, to compute weights.
-     */
-    const FE_Q<dim> fe_q;
-
     /**
      * A pointer to the Manifold in use.
      *
@@ -371,12 +366,12 @@ public:
    * @}
    */
 
-protected:
-
   /**
    * An FE_Q object, used to compute weights for Manifold quadratures.
    */
-  const FE_Q<dim,spacedim> fe_q;
+  static const FE_Q<dim,spacedim> fe_q;
+
+protected:
 
   /**
    * Transforms the point @p p on the real cell to the corresponding point on

diff --git a/source/fe/mapping_manifold.cc b/source/fe/mapping_manifold.cc
index 8c9e2a1db99aa662cf57af2e4e3b46f7f0709ed6..0357c2e8b43cb3c4009562cf3c470b4d9aa8a9c2 100644 (file)
--- a/source/fe/mapping_manifold.cc
+++ b/source/fe/mapping_manifold.cc
@@ -41,9 +41,15 @@
 
 DEAL_II_NAMESPACE_OPEN
 
+
+template<int dim, int spacedim>
+const FE_Q<dim,spacedim>
+MappingManifold<dim,spacedim>::fe_q = FE_Q<dim,spacedim>(1);
+
+
+
 template<int dim, int spacedim>
-MappingManifold<dim,spacedim>::InternalData::InternalData () :
-  fe_q(1)
+MappingManifold<dim,spacedim>::InternalData::InternalData ()
 {}
 
 
@@ -176,18 +182,14 @@ initialize_face (const UpdateFlags      update_flags,
 
 
 template<int dim, int spacedim>
-MappingManifold<dim,spacedim>::MappingManifold () :
-  fe_q(1)
-  // support_point_weights_on_quad (compute_support_point_weights_on_quad<dim>(this->polynomial_degree)),
-  // support_point_weights_on_hex (compute_support_point_weights_on_hex<dim>(this->polynomial_degree)),
+MappingManifold<dim,spacedim>::MappingManifold ()
 {
 }
 
 
 
 template<int dim, int spacedim>
-MappingManifold<dim,spacedim>::MappingManifold (const MappingManifold<dim,spacedim> &mapping) :
-  fe_q(1)
+MappingManifold<dim,spacedim>::MappingManifold (const MappingManifold<dim,spacedim> &mapping)
 {}
 
 
@@ -452,7 +454,7 @@ namespace internal
 
                   // Get the weights to compute the np point in real space
                   for (unsigned int j=0; j<GeometryInfo<dim>::vertices_per_cell; ++j)
-                    weights[j] = data.fe_q.shape_value(j, np);
+                    weights[j] = MappingManifold<dim,spacedim>::fe_q.shape_value(j, np);
 
                   Point<spacedim> NP=data.manifold->
                                      get_new_point(Quadrature<spacedim>(data.vertices, weights));

diff --git a/tests/fe/mapping_q_manifold_01.cc b/tests/fe/mapping_q_manifold_01.cc
index 27667e66495928b877cb852d0330b64823170c74..a3336d13fc89052ba8df59df0011118517b03bac 100644 (file)
--- a/tests/fe/mapping_q_manifold_01.cc
+++ b/tests/fe/mapping_q_manifold_01.cc
@@ -13,7 +13,7 @@
 //
 // ---------------------------------------------------------------------
 
-// Test high order MappingQ on Manifolds.
+// Test high order MappingQ on a ChartManifold.
 
 #include "../tests.h"