From b18534590db90c6edeae730bd73b9ae59178cc10 Mon Sep 17 00:00:00 2001 From: Wolfgang Bangerth Date: Thu, 15 Aug 2013 01:55:08 +0000 Subject: [PATCH] Patch by Juan Carlos Araujo Cabarcas: Relax an assertion. git-svn-id: https://svn.dealii.org/trunk@30312 0785d39b-7218-0410-832d-ea1e28bc413d --- deal.II/doc/news/changes.h | 10 ++++++++++ deal.II/source/fe/mapping_q.cc | 2 +- 2 files changed, 11 insertions(+), 1 deletion(-) 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.

Specific improvements

    +
  1. + 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. +
    + (Juan Carlos Araujo Cabarcas, Wolfgang Bangerth, 2013/08/14) +
  2. +
  3. 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::set_laplace_on_quad_vector(Table<2,double> &loqvs) const // one. check this for (unsigned int unit_point=0; unit_pointdegree*this->degree, ExcInternalError()); } -- 2.39.5