From 49384f57d0fb83b1780ddeb976b9a2c575da1f35 Mon Sep 17 00:00:00 2001
From: Tobias Reiter <t.reiter1@live.de>
Date: Thu, 28 Nov 2024 15:03:06 +0100
Subject: [PATCH] Scale direction vector by grid delta

---
 include/viennals/lsCalculateVisibilities.hpp | 2 +-
 1 file changed, 1 insertion(+), 1 deletion(-)

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);