From: Wolfgang Bangerth Date: Sun, 7 Nov 2010 20:34:39 +0000 (+0000) Subject: Ensure that all values we work with are finite. X-Git-Tag: v8.0.0~5031 X-Git-Url: https://gitweb.dealii.org/cgi-bin/gitweb.cgi?a=commitdiff_plain;h=e281032f4c8c230eb56aee11ac9dc6ef77face52;p=dealii.git Ensure that all values we work with are finite. git-svn-id: https://svn.dealii.org/trunk@22626 0785d39b-7218-0410-832d-ea1e28bc413d --- diff --git a/deal.II/source/numerics/error_estimator.cc b/deal.II/source/numerics/error_estimator.cc index bb75944abc..e2efdb9795 100644 --- a/deal.II/source/numerics/error_estimator.cc +++ b/deal.II/source/numerics/error_estimator.cc @@ -385,7 +385,10 @@ namespace internal ExcInternalError()); for (unsigned int i=0; isecond.size(); ++i) - Assert (p->second[i] >= 0, ExcInternalError()); + { + Assert (numbers::is_finite(p->second[i]), ExcInternalError()); + Assert (p->second[i] >= 0, ExcInternalError()); + } face_integrals[p->first] = p->second; }