add two Asserts in fe_enriched_step-36 unit test

diff --git a/tests/fe/fe_enriched_step-36.cc b/tests/fe/fe_enriched_step-36.cc
index a22183a967221f1fffa6a53451ecde3ef94bb723..22702ae43e9ad2c5836e031455f6f10d913e3414 100644 (file)
--- a/tests/fe/fe_enriched_step-36.cc
+++ b/tests/fe/fe_enriched_step-36.cc
@@ -77,6 +77,8 @@ template<int dim>
 double PotentialFunction<dim>::value(const dealii::Point<dim> &p,
                                      const unsigned int ) const
 {
+  Assert (p.square() > 0.,
+          ExcDivideByZero());
   return -1.0 / std::sqrt(p.square());
 }
 
@@ -114,6 +116,8 @@ public:
   {
     Tensor<1,dim> dist = p-origin;
     const double r = dist.norm();
+    Assert (r > 0.,
+            ExcDivideByZero());
     dist/=r;
     return -Z*std::exp(-Z*r)*dist;
   }