diff --git a/include/viennals/lsCalculateVisibilities.hpp b/include/viennals/lsCalculateVisibilities.hpp
index a1872f12..c1a7ee1e 100644
--- a/include/viennals/lsCalculateVisibilities.hpp
+++ b/include/viennals/lsCalculateVisibilities.hpp
@@ -52,7 +52,7 @@ template <class NumericType, int D> class CalculateVisibilities {
     //****************************
 
     // Invert the vector
-    auto dir = Normalize(Inv(direction));
+    auto dir = Normalize(Inv(direction)) * domiain.getGrid().getGridDelta();
 
     auto numDefinedPoints = domain.getNumberOfPoints();
     std::vector<NumericType> visibilities(numDefinedPoints);