From: Martin Kronbichler <kronbichler@lnm.mw.tum.de>
Date: Mon, 11 Jan 2016 11:42:47 +0000 (+0100)
Subject: Add instantiation
X-Git-Tag: v8.4.0-rc2~76^2~8
X-Git-Url: https://gitweb.dealii.org/cgi-bin/gitweb.cgi?a=commitdiff_plain;h=dcd13d34391c17c8321c435b03ed5fcded8fca08;p=dealii.git

Add instantiation
---

diff --git a/include/deal.II/matrix_free/shape_info.templates.h b/include/deal.II/matrix_free/shape_info.templates.h
index 8b3aa6fa81..6b6a2d63a5 100644
--- a/include/deal.II/matrix_free/shape_info.templates.h
+++ b/include/deal.II/matrix_free/shape_info.templates.h
@@ -314,25 +314,25 @@ namespace internal
         for (unsigned int q=0; q<stride; ++q)
           {
             shape_val_evenodd[i*stride+q] =
-              0.5 * (shape_values[i*n_q_points_1d+q] +
-                     shape_values[i*n_q_points_1d+n_q_points_1d-1-q]);
+              Number(0.5) * (shape_values[i*n_q_points_1d+q] +
+                             shape_values[i*n_q_points_1d+n_q_points_1d-1-q]);
             shape_val_evenodd[(fe_degree-i)*stride+q] =
-              0.5 * (shape_values[i*n_q_points_1d+q] -
-                     shape_values[i*n_q_points_1d+n_q_points_1d-1-q]);
+              Number(0.5) * (shape_values[i*n_q_points_1d+q] -
+                             shape_values[i*n_q_points_1d+n_q_points_1d-1-q]);
 
             shape_gra_evenodd[i*stride+q] =
-              0.5 * (shape_gradients[i*n_q_points_1d+q] +
-                     shape_gradients[i*n_q_points_1d+n_q_points_1d-1-q]);
+              Number(0.5) * (shape_gradients[i*n_q_points_1d+q] +
+                             shape_gradients[i*n_q_points_1d+n_q_points_1d-1-q]);
             shape_gra_evenodd[(fe_degree-i)*stride+q] =
-              0.5 * (shape_gradients[i*n_q_points_1d+q] -
-                     shape_gradients[i*n_q_points_1d+n_q_points_1d-1-q]);
+              Number(0.5) * (shape_gradients[i*n_q_points_1d+q] -
+                             shape_gradients[i*n_q_points_1d+n_q_points_1d-1-q]);
 
             shape_hes_evenodd[i*stride+q] =
-              0.5 * (shape_hessians[i*n_q_points_1d+q] +
-                     shape_hessians[i*n_q_points_1d+n_q_points_1d-1-q]);
+              Number(0.5) * (shape_hessians[i*n_q_points_1d+q] +
+                             shape_hessians[i*n_q_points_1d+n_q_points_1d-1-q]);
             shape_hes_evenodd[(fe_degree-i)*stride+q] =
-              0.5 * (shape_hessians[i*n_q_points_1d+q] -
-                     shape_hessians[i*n_q_points_1d+n_q_points_1d-1-q]);
+              Number(0.5) * (shape_hessians[i*n_q_points_1d+q] -
+                             shape_hessians[i*n_q_points_1d+n_q_points_1d-1-q]);
           }
       if (fe_degree % 2 == 0)
         for (unsigned int q=0; q<stride; ++q)
diff --git a/source/matrix_free/matrix_free.cc b/source/matrix_free/matrix_free.cc
index 14a9177eca..4fab36c845 100644
--- a/source/matrix_free/matrix_free.cc
+++ b/source/matrix_free/matrix_free.cc
@@ -26,5 +26,6 @@ DEAL_II_NAMESPACE_OPEN
 
 template struct internal::MatrixFreeFunctions::ShapeInfo<double>;
 template struct internal::MatrixFreeFunctions::ShapeInfo<float>;
+template struct internal::MatrixFreeFunctions::ShapeInfo<long double>;
 
 DEAL_II_NAMESPACE_CLOSE
diff --git a/source/matrix_free/matrix_free.inst.in b/source/matrix_free/matrix_free.inst.in
index c4d2888d0a..47aeeaa5ba 100644
--- a/source/matrix_free/matrix_free.inst.in
+++ b/source/matrix_free/matrix_free.inst.in
@@ -39,6 +39,8 @@ for (deal_II_dimension : DIMENSIONS)
   template void internal::MatrixFreeFunctions::ShapeInfo<float>::reinit
   <deal_II_dimension>(const Quadrature<1> &, const FiniteElement
                       <deal_II_dimension,deal_II_dimension> &, const unsigned int);
+  template void internal::MatrixFreeFunctions::ShapeInfo<long double>::reinit
+  <deal_II_dimension>(const Quadrature<1> &, const FiniteElement
+                      <deal_II_dimension,deal_II_dimension> &, const unsigned int);
 #endif
 }
-