From: Wolfgang Bangerth <bangerth@math.tamu.edu>
Date: Thu, 15 Aug 2013 01:55:08 +0000 (+0000)
Subject: Patch by Juan Carlos Araujo Cabarcas: Relax an assertion.
X-Git-Tag: v8.1.0~1062
X-Git-Url: https://gitweb.dealii.org/cgi-bin/gitweb.cgi?a=commitdiff_plain;h=b18534590db90c6edeae730bd73b9ae59178cc10;p=dealii.git

Patch by Juan Carlos Araujo Cabarcas: Relax an assertion.

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

diff --git a/deal.II/doc/news/changes.h b/deal.II/doc/news/changes.h
index 75be54c8ea..89bdcf8984 100644
--- a/deal.II/doc/news/changes.h
+++ b/deal.II/doc/news/changes.h
@@ -56,6 +56,16 @@ inconvenience this causes.
 <h3>Specific improvements</h3>
 
 <ol>
+  <li>
+  Fixed: Creating objects of type MappingQ was previously only possible
+  for low order polynomials. For orders higher than around 6, one ran
+  into assertions that tested for internal consistency. These assertions
+  have now been appropriately relaxes for the growth of round-off errors
+  with growing polynomial degrees.
+  <br>
+  (Juan Carlos Araujo Cabarcas, Wolfgang Bangerth, 2013/08/14)
+  </li>
+
   <li>
   New: MappingQEulerian is now also instantiated for vector elements
   of type TrilinosWrappers::Vector as well as the MPI and block
diff --git a/deal.II/source/fe/mapping_q.cc b/deal.II/source/fe/mapping_q.cc
index 663d77482e..18cb2b1110 100644
--- a/deal.II/source/fe/mapping_q.cc
+++ b/deal.II/source/fe/mapping_q.cc
@@ -696,7 +696,7 @@ MappingQ<dim,spacedim>::set_laplace_on_quad_vector(Table<2,double> &loqvs) const
   // one. check this
   for (unsigned int unit_point=0; unit_point<loqvs.n_rows(); ++unit_point)
     Assert(std::fabs(std::accumulate(loqvs[unit_point].begin(),
-                                     loqvs[unit_point].end(),0.)-1)<1e-13,
+                                     loqvs[unit_point].end(),0.)-1)<1e-12*this->degree*this->degree,
            ExcInternalError());
 }