From 30292ed3991238eeeb6a4700c4e6b723da3cc05f Mon Sep 17 00:00:00 2001
From: Daniel Arndt <daniel.arndt@iwr.uni-heidelberg.de>
Date: Fri, 15 Dec 2017 17:17:58 +0100
Subject: [PATCH] Fix fe_field_function_03

---
 include/deal.II/numerics/fe_field_function.templates.h | 6 ++++++
 tests/mpi/fe_field_function_03.cc                      | 4 ++++
 2 files changed, 10 insertions(+)

diff --git a/include/deal.II/numerics/fe_field_function.templates.h b/include/deal.II/numerics/fe_field_function.templates.h
index e26d997d1e..0071ed6ba5 100644
--- a/include/deal.II/numerics/fe_field_function.templates.h
+++ b/include/deal.II/numerics/fe_field_function.templates.h
@@ -283,6 +283,8 @@ namespace Functions
     // Now gather all the informations we need
     for (unsigned int i=0; i<n_cells; ++i)
       {
+        AssertThrow (!cells[i]->is_artificial(),
+                     VectorTools::ExcPointNotAvailableHere());
         fe_v.reinit(cells[i], i, 0);
         const unsigned int nq = qpoints[i].size();
         std::vector< Vector<typename VectorType::value_type> > vvalues (nq, Vector<typename VectorType::value_type>(this->n_components));
@@ -344,6 +346,8 @@ namespace Functions
     // Now gather all the informations we need
     for (unsigned int i=0; i<n_cells; ++i)
       {
+        AssertThrow (!cells[i]->is_artificial(),
+                     VectorTools::ExcPointNotAvailableHere());
         fe_v.reinit(cells[i], i, 0);
         const unsigned int nq = qpoints[i].size();
         std::vector< std::vector<Tensor<1,dim,typename VectorType::value_type> > >
@@ -409,6 +413,8 @@ namespace Functions
     // Now gather all the informations we need
     for (unsigned int i=0; i<n_cells; ++i)
       {
+        AssertThrow (!cells[i]->is_artificial(),
+                     VectorTools::ExcPointNotAvailableHere());
         fe_v.reinit(cells[i], i, 0);
         const unsigned int nq = qpoints[i].size();
         std::vector< Vector<typename VectorType::value_type> > vvalues (nq, Vector<typename VectorType::value_type>(this->n_components));
diff --git a/tests/mpi/fe_field_function_03.cc b/tests/mpi/fe_field_function_03.cc
index bf4bd61c39..1c8a5a0195 100644
--- a/tests/mpi/fe_field_function_03.cc
+++ b/tests/mpi/fe_field_function_03.cc
@@ -188,6 +188,10 @@ void test()
         {
           deallog << "  ExcPointNotAvailableHere" << std::endl;
         }
+      catch (std::exception &exc)
+        {
+          deallog << exc.what() << std::endl;
+        }
       catch (...)
         {
           deallog << "  Oh no! Some other error that we shouldn't get." << std::endl;
-- 
2.39.5