From: bangerth Date: Thu, 15 Aug 2013 01:55:08 +0000 (+0000) Subject: Patch by Juan Carlos Araujo Cabarcas: Relax an assertion. X-Git-Url: https://gitweb.dealii.org/cgi-bin/gitweb.cgi?a=commitdiff_plain;h=57b7448602327c70e00cb69716547075d806eebb;p=dealii-svn.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.

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()); }