diff --git a/include/ddc/kernels/fft.hpp b/include/ddc/kernels/fft.hpp
index 74a01083d..884581968 100644
--- a/include/ddc/kernels/fft.hpp
+++ b/include/ddc/kernels/fft.hpp
@@ -332,29 +332,7 @@ int N(ddc::DiscreteDomain<DDimX...> x_mesh)
     static_assert(
             (is_uniform_point_sampling_v<DDimX> && ...),
             "DDimX dimensions should derive from UniformPointSampling");
-    return ddc::get<DDim>(x_mesh.extents());
-}
-
-template <typename DDim, typename... DDimX>
-double a(ddc::DiscreteDomain<DDimX...> x_mesh)
-{
-    static_assert(
-            (is_uniform_point_sampling_v<DDimX> && ...),
-            "DDimX dimensions should derive from UniformPointSampling");
-    return ((2 * N<DDim>(x_mesh) - 1) * coordinate(ddc::select<DDim>(x_mesh).front())
-            - coordinate(ddc::select<DDim>(x_mesh).back()))
-           / 2 / (N<DDim>(x_mesh) - 1);
-}
-
-template <typename DDim, typename... DDimX>
-double b(ddc::DiscreteDomain<DDimX...> x_mesh)
-{
-    static_assert(
-            (is_uniform_point_sampling_v<DDimX> && ...),
-            "DDimX dimensions should derive from UniformPointSampling");
-    return ((2 * N<DDim>(x_mesh) - 1) * coordinate(ddc::select<DDim>(x_mesh).back())
-            - coordinate(ddc::select<DDim>(x_mesh).front()))
-           / 2 / (N<DDim>(x_mesh) - 1);
+    return static_cast<int>(x_mesh.template extent<DDim>());
 }
 
 // core
@@ -595,7 +573,10 @@ typename DDimFx::template Impl<DDimFx, Kokkos::HostSpace> init_fourier_space(
             ddc::Coordinate<typename DDimFx::continuous_dimension_type>(0),
             ddc::Coordinate<typename DDimFx::continuous_dimension_type>(
                     2 * (ddc::detail::fft::N<DDimX>(x_mesh) - 1)
-                    / (ddc::detail::fft::b<DDimX>(x_mesh) - ddc::detail::fft::a<DDimX>(x_mesh))
+                    * (ddc::detail::fft::N<DDimX>(x_mesh) - 1)
+                    / static_cast<double>(
+                            ddc::detail::fft::N<DDimX>(x_mesh)
+                            * (ddc::coordinate(x_mesh.back()) - ddc::coordinate(x_mesh.front())))
                     * Kokkos::numbers::pi),
             ddc::DiscreteVector<DDimFx>(ddc::detail::fft::N<DDimX>(x_mesh)),
             ddc::DiscreteVector<DDimFx>(ddc::detail::fft::N<DDimX>(x_mesh)));
diff --git a/include/ddc/periodic_sampling.hpp b/include/ddc/periodic_sampling.hpp
index 5b73b65ec..b0882a3be 100644
--- a/include/ddc/periodic_sampling.hpp
+++ b/include/ddc/periodic_sampling.hpp
@@ -127,8 +127,8 @@ class PeriodicSampling : detail::PeriodicSamplingBase
         {
             return m_origin
                    + continuous_element_type(
-                             static_cast<int>((icoord.uid() + m_n_period / 2 - 1) % m_n_period)
-                             - static_cast<int>(m_n_period / 2 - 1))
+                             static_cast<int>((icoord.uid() + m_n_period / 2) % m_n_period)
+                             - static_cast<int>(m_n_period / 2))
                              * m_step;
         }
     };
diff --git a/tests/fft.cpp b/tests/fft.cpp
index f2cfd4727..183b4b28b 100644
--- a/tests/fft.cpp
+++ b/tests/fft.cpp
@@ -29,6 +29,34 @@ struct DFDim : ddc::PeriodicSampling<Kx>
 {
 };
 
+template <typename X>
+static void test_fourier_mesh(std::size_t Nx)
+{
+    double const a = -10;
+    double const b = 10;
+
+    DDom<DDim<X>> const x_mesh(ddc::init_discrete_space<DDim<X>>(DDim<X>::template init<DDim<X>>(
+            ddc::Coordinate<X>(a + (b - a) / Nx / 2),
+            ddc::Coordinate<X>(b - (b - a) / Nx / 2),
+            DVect<DDim<X>>(Nx))));
+    ddc::init_discrete_space<DFDim<ddc::Fourier<X>>>(
+            ddc::init_fourier_space<DFDim<ddc::Fourier<X>>>(ddc::DiscreteDomain<DDim<X>>(x_mesh)));
+    DDom<DFDim<ddc::Fourier<X>>> const k_mesh
+            = ddc::FourierMesh<DFDim<ddc::Fourier<X>>>(x_mesh, true);
+
+    double const epsilon = 1e-14;
+    ddc::DiscreteElement<DFDim<ddc::Fourier<X>>> const k_front = k_mesh.front();
+    for (ddc::DiscreteElement<DFDim<ddc::Fourier<X>>> k : k_mesh) {
+        double const ka = -Kokkos::numbers::pi * Nx / (b - a);
+        double const kb = Kokkos::numbers::pi * Nx / (b - a);
+        double const k_pred = 2 * (k - k_front) * Kokkos::numbers::pi / (b - a);
+        EXPECT_NEAR(
+                ddc::coordinate(k),
+                k_pred - (kb - ka) * Kokkos::floor((k_pred - ka) / (kb - ka)),
+                epsilon);
+    }
+}
+
 // TODO:
 // - FFT multidim but according to a subset of dimensions
 template <typename ExecSpace, typename MemorySpace, typename Tin, typename Tout, typename... X>
@@ -199,6 +227,16 @@ struct RDimX;
 struct RDimY;
 struct RDimZ;
 
+TEST(FourierMesh, Even)
+{
+    test_fourier_mesh<RDimX>(16);
+}
+
+TEST(FourierMesh, Odd)
+{
+    test_fourier_mesh<RDimX>(17);
+}
+
 #if fftw_serial_AVAIL
 TEST(FFTNorm, OFF)
 {