Add another safety check.

git-svn-id: https://svn.dealii.org/trunk@4909 0785d39b-7218-0410-832d-ea1e28bc413d

diff --git a/deal.II/deal.II/source/numerics/error_estimator.cc b/deal.II/deal.II/source/numerics/error_estimator.cc
index c9efeaa25d9c325c1cacff5ab1f6747a6a7b4734..7c6460dbf55333606576a0fefcc92e72ffdc36aa 100644 (file)
--- a/deal.II/deal.II/source/numerics/error_estimator.cc
+++ b/deal.II/deal.II/source/numerics/error_estimator.cc
@@ -921,6 +921,13 @@ integrate_over_regular_face (Data                       &data,
        for (unsigned int p=0; p<n_q_points; ++p)
          face_integral[n] += ::sqr(data.phi[n][p][component]) *
                              data.JxW_values[p];
+
+                                  // double check that the element
+                                  // already exists and that it was
+                                  // not already written to
+  Assert (data.face_integrals.find (face) != data.face_integrals.end(),
+         ExcInternalError());
+  Assert (data.face_integrals[face][0] < 0, ExcInternalError());
   
   data.face_integrals[face] = face_integral;
 };