From d0f90cae9cecc77366d33b61bb517080bcd4c93a Mon Sep 17 00:00:00 2001 From: bangerth Date: Sun, 7 Nov 2010 20:34:39 +0000 Subject: [PATCH] Ensure that all values we work with are finite. git-svn-id: https://svn.dealii.org/trunk@22626 0785d39b-7218-0410-832d-ea1e28bc413d --- deal.II/source/numerics/error_estimator.cc | 5 ++++- 1 file changed, 4 insertions(+), 1 deletion(-) 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; } -- 2.39.5