From aa901f0f8cb9f926f01e04ebb38382831d4a491b Mon Sep 17 00:00:00 2001
From: Emanuele De Rubeis <e.derubeis@ira.inaf.it>
Date: Fri, 17 May 2024 11:56:11 +0000
Subject: [PATCH] Update w-stacking.cu
---
w-stacking.cu | 10 +++++-----
1 file changed, 5 insertions(+), 5 deletions(-)
diff --git a/w-stacking.cu b/w-stacking.cu
index f52889a..9aa0d66 100755
--- a/w-stacking.cu
+++ b/w-stacking.cu
@@ -159,7 +159,7 @@ __global__ void convolve_g(
for (k = kmin; k <= kmax; k++)
{
- double v_dist = (double)k - pos_v;
+ double v_dist = (double)k+0.5 - pos_v;
int increaseprecision = 5;
for (j = jmin; j <= jmax; j++)
@@ -472,13 +472,13 @@ void wstack(
for (k = kmin; k <= kmax; k++)
{
- //double v_dist = (double)k+0.5 - pos_v;
- double v_dist = (double)k - pos_v;
+ double v_dist = (double)k+0.5 - pos_v;
+ //double v_dist = (double)k - pos_v;
for (j = jmin; j <= jmax; j++)
{
- //double u_dist = (double)j+0.5 - pos_u;
- double u_dist = (double)j - pos_u;
+ double u_dist = (double)j+0.5 - pos_u;
+ //double u_dist = (double)j - pos_u;
uint iKer = 2 * (j + k*grid_size_x + grid_w*grid_size_x*grid_size_y);
int jKer = (int)(increaseprecision * (fabs(u_dist+(double)KernelLen)));
int kKer = (int)(increaseprecision * (fabs(v_dist+(double)KernelLen)));
--
GitLab