diff --git a/src/cosserat-continuum.cc b/src/cosserat-continuum.cc
index 45c0c5289a83428ad41cc7dd7d8c5005bd31122f..1bd4144a3c9a616571fffa131ce13bb40723b2ed 100644
--- a/src/cosserat-continuum.cc
+++ b/src/cosserat-continuum.cc
@@ -266,12 +266,10 @@ int main (int argc, char *argv[]) try
 
     for (auto&& vertex : vertices(gridView))
     {
-        bool isDirichlet;
-        pythonDirichletVertices.evaluate(vertex.geometry().corner(0), isDirichlet);
+        bool isDirichlet = pythonDirichletVertices(vertex.geometry().corner(0));
         dirichletVertices[indexSet.index(vertex)] = isDirichlet;
 
-        bool isNeumann;
-        pythonNeumannVertices.evaluate(vertex.geometry().corner(0), isNeumann);
+        bool isNeumann = pythonNeumannVertices(vertex.geometry().corner(0));
         neumannVertices[indexSet.index(vertex)] = isNeumann;
     }
 
diff --git a/src/gradient-flow.cc b/src/gradient-flow.cc
index e9c32fb13b43737564321e46db8845c226dfe9df..0ada4c171bd8e209fb95163e7a9d887552fc3a79 100644
--- a/src/gradient-flow.cc
+++ b/src/gradient-flow.cc
@@ -164,8 +164,7 @@ int main (int argc, char *argv[]) try
 
   for (auto&& vertex : vertices(grid->leafGridView()))
   {
-    bool isDirichlet;
-    pythonDirichletVertices.evaluate(vertex.geometry().corner(0), isDirichlet);
+    bool isDirichlet = pythonDirichletVertices(vertex.geometry().corner(0));
     dirichletVertices[indexSet.index(vertex)] = isDirichlet;
   }
 
diff --git a/src/harmonicmaps.cc b/src/harmonicmaps.cc
index 84a08d3fe277d95ef060979336a4a5f436df7de2..772fc06e6103f0b407841c6a42b496596f107bd5 100644
--- a/src/harmonicmaps.cc
+++ b/src/harmonicmaps.cc
@@ -224,9 +224,7 @@ int main (int argc, char *argv[])
 
     for (auto&& vertex : vertices(gridView))
     {
-      //bool isDirichlet;
       bool isDirichlet = pythonDirichletVertices(vertex.geometry().corner(0));
-      pythonDirichletVertices.evaluate(vertex.geometry().corner(0), isDirichlet);
       dirichletVertices[indexSet.index(vertex)] = isDirichlet;
     }