]> https://gitweb.dealii.org/ - dealii.git/commitdiff
Ensure that all values we work with are finite.
authorWolfgang Bangerth <bangerth@math.tamu.edu>
Sun, 7 Nov 2010 20:34:39 +0000 (20:34 +0000)
committerWolfgang Bangerth <bangerth@math.tamu.edu>
Sun, 7 Nov 2010 20:34:39 +0000 (20:34 +0000)
git-svn-id: https://svn.dealii.org/trunk@22626 0785d39b-7218-0410-832d-ea1e28bc413d

deal.II/source/numerics/error_estimator.cc

index bb75944abccaa326acc7f1edafdeb85d12bf0633..e2efdb979559fac1a91a7966a4202b806e2669df 100644 (file)
@@ -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;
        }

In the beginning the Universe was created. This has made a lot of people very angry and has been widely regarded as a bad move.

Douglas Adams


Typeset in Trocchi and Trocchi Bold Sans Serif.