]> https://gitweb.dealii.org/ - dealii-svn.git/commitdiff
Patch by Juan Carlos Araujo Cabarcas: Relax an assertion.
authorbangerth <bangerth@0785d39b-7218-0410-832d-ea1e28bc413d>
Thu, 15 Aug 2013 01:55:08 +0000 (01:55 +0000)
committerbangerth <bangerth@0785d39b-7218-0410-832d-ea1e28bc413d>
Thu, 15 Aug 2013 01:55:08 +0000 (01:55 +0000)
git-svn-id: https://svn.dealii.org/trunk@30312 0785d39b-7218-0410-832d-ea1e28bc413d

deal.II/doc/news/changes.h
deal.II/source/fe/mapping_q.cc

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

In the beginning the Universe was created. This has made a lot of people very angry and has been widely regarded as a bad move.

Douglas Adams


Typeset in Trocchi and Trocchi Bold Sans Serif.