Let ReferenceCell::get_nodal_type_quadrature() return a 'const' reference.

diff --git a/include/deal.II/grid/reference_cell.h b/include/deal.II/grid/reference_cell.h
index 12a6b8cd407840d1e241ebc7665eb4bf416515d2..63f216235e5dde91c8c4bf8097da077d5894c266 100644 (file)
--- a/include/deal.II/grid/reference_cell.h
+++ b/include/deal.II/grid/reference_cell.h
@@ -732,7 +732,7 @@ namespace ReferenceCell
    * @note The weights are not filled.
    */
   template <int dim>
-  Quadrature<dim> &
+  const Quadrature<dim> &
   get_nodal_type_quadrature(const Type &reference_cell);
 
   namespace internal

diff --git a/source/grid/reference_cell.cc b/source/grid/reference_cell.cc
index 9e473a05923de5f9ebc4483ecb659ba5e8aa4ce4..4643d48c7c489e6498a87de57ee2fadb4f523089 100644 (file)
--- a/source/grid/reference_cell.cc
+++ b/source/grid/reference_cell.cc
@@ -223,7 +223,7 @@ namespace ReferenceCell
   }
 
   template <int dim>
-  Quadrature<dim> &
+  const Quadrature<dim> &
   get_nodal_type_quadrature(const Type &reference_cell)
   {
     AssertDimension(dim, reference_cell.get_dimension());
@@ -237,22 +237,26 @@ namespace ReferenceCell
 
     if (reference_cell == get_hypercube(dim))
       {
-        static Quadrature<dim> quadrature = create_quadrature(reference_cell);
+        static const Quadrature<dim> quadrature =
+          create_quadrature(reference_cell);
         return quadrature;
       }
     else if (reference_cell == Type::Tri || reference_cell == Type::Tet)
       {
-        static Quadrature<dim> quadrature = create_quadrature(reference_cell);
+        static const Quadrature<dim> quadrature =
+          create_quadrature(reference_cell);
         return quadrature;
       }
     else if (reference_cell == Type::Pyramid)
       {
-        static Quadrature<dim> quadrature = create_quadrature(reference_cell);
+        static const Quadrature<dim> quadrature =
+          create_quadrature(reference_cell);
         return quadrature;
       }
     else if (reference_cell == Type::Wedge)
       {
-        static Quadrature<dim> quadrature = create_quadrature(reference_cell);
+        static const Quadrature<dim> quadrature =
+          create_quadrature(reference_cell);
         return quadrature;
       }
     else

diff --git a/source/grid/reference_cell.inst.in b/source/grid/reference_cell.inst.in
index 29eec7c46d7f7928d9bcdcf2cc5d2af6d914f839..e5b8045cdf02cc720b6d558e8408f7eb35ab4441 100644 (file)
--- a/source/grid/reference_cell.inst.in
+++ b/source/grid/reference_cell.inst.in
@@ -38,6 +38,6 @@ for (deal_II_dimension : DIMENSIONS)
   {
     template Quadrature<deal_II_dimension> get_gauss_type_quadrature(
       const Type &reference_cell, const unsigned n_points_1D);
-    template Quadrature<deal_II_dimension> &get_nodal_type_quadrature(
+    template const Quadrature<deal_II_dimension> &get_nodal_type_quadrature(
       const Type &reference_cell);
   }