diff --git a/include/viennals/lsCalculateVisibilities.hpp b/include/viennals/lsCalculateVisibilities.hpp
index a05d4301..8d6bfe09 100644
--- a/include/viennals/lsCalculateVisibilities.hpp
+++ b/include/viennals/lsCalculateVisibilities.hpp
@@ -15,7 +15,8 @@ template <class NumericType, int D> class CalculateVisibilities {
 
   CalculateVisibilities(
       const SmartPointer<Domain<NumericType, D>> &passedLevelSet,
-      int passedOpenBoundaryDirection, bool passedIsOpenBoundaryNegative)
+      int passedOpenBoundaryDirection = D - 1,
+      bool passedIsOpenBoundaryNegative = false)
       : levelSet(passedLevelSet),
         openBoundaryDirection(passedOpenBoundaryDirection),
         isOpenBoundaryNegative(passedIsOpenBoundaryNegative) {}