From: Martin Kronbichler <kronbichler@lnm.mw.tum.de>
Date: Thu, 28 Apr 2016 21:09:57 +0000 (+0200)
Subject: Fix test of exception in FE_Q at degree 0
X-Git-Tag: v8.5.0-rc1~1068^2
X-Git-Url: https://gitweb.dealii.org/cgi-bin/gitweb.cgi?a=commitdiff_plain;h=refs%2Fpull%2F2563%2Fhead;p=dealii.git

Fix test of exception in FE_Q at degree 0
---

diff --git a/source/fe/fe_q.cc b/source/fe/fe_q.cc
index 343d5f5839..b1ad9218e5 100644
--- a/source/fe/fe_q.cc
+++ b/source/fe/fe_q.cc
@@ -25,7 +25,6 @@ DEAL_II_NAMESPACE_OPEN
 
 namespace
 {
-  template <int dim, int spacedim>
   std::vector<Point<1> >
   get_QGaussLobatto_points (const unsigned int degree)
   {
@@ -33,7 +32,7 @@ namespace
       return QGaussLobatto<1>(degree+1).get_points();
     else
       {
-        typedef FE_Q_Base<TensorProductPolynomials<dim>, dim, spacedim> FEQ;
+        typedef FE_Q_Base<TensorProductPolynomials<1>, 1, 1> FEQ;
         AssertThrow(false, typename FEQ::ExcFEQCannotHaveDegree0());
       }
     return std::vector<Point<1> >();
@@ -46,13 +45,13 @@ template <int dim, int spacedim>
 FE_Q<dim,spacedim>::FE_Q (const unsigned int degree)
   :
   FE_Q_Base<TensorProductPolynomials<dim>, dim, spacedim>
-  (TensorProductPolynomials<dim>(Polynomials::generate_complete_Lagrange_basis(get_QGaussLobatto_points<dim,spacedim>(degree))),
+  (TensorProductPolynomials<dim>(Polynomials::generate_complete_Lagrange_basis(get_QGaussLobatto_points(degree))),
    FiniteElementData<dim>(this->get_dpo_vector(degree),
                           1, degree,
                           FiniteElementData<dim>::H1),
    std::vector<bool> (1, false))
 {
-  this->initialize(get_QGaussLobatto_points<dim,spacedim>(degree));
+  this->initialize(get_QGaussLobatto_points(degree));
 }
 
 
diff --git a/source/fe/fe_q_base.cc b/source/fe/fe_q_base.cc
index dcb28803f2..f8f9e64194 100644
--- a/source/fe/fe_q_base.cc
+++ b/source/fe/fe_q_base.cc
@@ -1120,7 +1120,8 @@ template <class PolynomialType, int dim, int spacedim>
 std::vector<unsigned int>
 FE_Q_Base<PolynomialType,dim,spacedim>::get_dpo_vector(const unsigned int degree)
 {
-  AssertThrow(degree>0, ExcFEQCannotHaveDegree0());
+  typedef FE_Q_Base<PolynomialType,dim,spacedim> FEQ;
+  AssertThrow(degree>0, typename FEQ::ExcFEQCannotHaveDegree0());
   std::vector<unsigned int> dpo(dim+1, 1U);
   for (unsigned int i=1; i<dpo.size(); ++i)
     dpo[i]=dpo[i-1]*(degree-1);
diff --git a/tests/fe/fe_q_01.output b/tests/fe/fe_q_01.output
index 158d3d7463..e6fb08a7c4 100644
--- a/tests/fe/fe_q_01.output
+++ b/tests/fe/fe_q_01.output
@@ -1,2 +1,2 @@
 
-DEAL::ExcFEQCannotHaveDegree0()
+DEAL::typename FEQ::ExcFEQCannotHaveDegree0()