From: Martin Kronbichler <kronbichler@lnm.mw.tum.de>
Date: Tue, 14 May 2013 15:01:54 +0000 (+0000)
Subject: Avoid running into an assertion for very high order.
X-Git-Tag: v8.0.0~525
X-Git-Url: https://gitweb.dealii.org/cgi-bin/gitweb.cgi?a=commitdiff_plain;h=2cee4c74ac3d6dbb06f6305feb3d6ac5ecb96638;p=dealii.git

Avoid running into an assertion for very high order.

git-svn-id: https://svn.dealii.org/trunk@29511 0785d39b-7218-0410-832d-ea1e28bc413d
---

diff --git a/deal.II/source/fe/fe_q_base.cc b/deal.II/source/fe/fe_q_base.cc
index bb320f8a4c..331131daf3 100644
--- a/deal.II/source/fe/fe_q_base.cc
+++ b/deal.II/source/fe/fe_q_base.cc
@@ -185,14 +185,15 @@ struct FE_Q_Base<POLY,xdim,xspacedim>::Implementation
     const std::vector<unsigned int> &index_map_inverse =
       fe.poly_space.get_numbering_inverse();
     const std::vector<unsigned int> face_index_map =
-      FE_Q_Helper::face_lexicographic_to_hierarchic_numbering<2>(fe.degree);
-    Assert(fe.poly_space.compute_value(index_map_inverse[0],Point<2>()) == 1.,
+      FE_Q_Helper::face_lexicographic_to_hierarchic_numbering<dim>(fe.degree);
+    Assert(std::abs(fe.poly_space.compute_value(index_map_inverse[0],Point<dim>())
+                    - 1.) < 1e-14,
            ExcInternalError());
 
     for (unsigned int i=0; i<constraint_points.size(); ++i)
       for (unsigned int j=0; j<fe.degree+1; ++j)
         {
-          Point<2> p;
+          Point<dim> p;
           p[0] = constraint_points[i](0);
           fe.interface_constraints(i,face_index_map[j]) =
             fe.poly_space.compute_value(index_map_inverse[j], p);
@@ -304,8 +305,9 @@ struct FE_Q_Base<POLY,xdim,xspacedim>::Implementation
     const std::vector<unsigned int> &index_map_inverse =
       fe.poly_space.get_numbering_inverse();
     const std::vector<unsigned int> face_index_map =
-      FE_Q_Helper::face_lexicographic_to_hierarchic_numbering<3>(fe.degree);
-    Assert(fe.poly_space.compute_value(index_map_inverse[0],Point<3>()) == 1.,
+      FE_Q_Helper::face_lexicographic_to_hierarchic_numbering<dim>(fe.degree);
+    Assert(std::abs(fe.poly_space.compute_value(index_map_inverse[0],Point<dim>())
+                    - 1.) < 1e-14,
            ExcInternalError());
 
     fe.interface_constraints