From e281032f4c8c230eb56aee11ac9dc6ef77face52 Mon Sep 17 00:00:00 2001
From: Wolfgang Bangerth <bangerth@math.tamu.edu>
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; i<p->second.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