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 75be54c8ea08acdcbcf158f5dbcffe26e9d8070a..89bdcf89845ffd226d82f70daa79dc2eef431e5a 100644 (file)
--- 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 663d77482e3418f5a37b1f4ba7c2d2b01bf26f5d..18cb2b111045fd59d847b1ffcd54f81ffe349603 100644 (file)
--- 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());
 }