diff --git a/vendor/ddc/.github/workflows/pages.yml b/vendor/ddc/.github/workflows/pages.yml
index cd9d0bfb1..39ed4c6a9 100644
--- a/vendor/ddc/.github/workflows/pages.yml
+++ b/vendor/ddc/.github/workflows/pages.yml
@@ -31,6 +31,7 @@ jobs:
         echo ${{ secrets.GITHUB_TOKEN }} | docker login ghcr.io -u ${{ github.actor }} --password-stdin || true
         docker pull ghcr.io/cexa-project/ddc/doxygen || true
         DOCKER_BUILDKIT=1 docker build \
+          --build-arg BUILDKIT_INLINE_CACHE=1 \
           --cache-from ghcr.io/cexa-project/ddc/doxygen \
           -t ghcr.io/cexa-project/ddc/doxygen \
           -t ghcr.io/cexa-project/ddc/doxygen:${GITHUB_SHA:0:7} \
diff --git a/vendor/ddc/docs/CMakeLists.txt b/vendor/ddc/docs/CMakeLists.txt
index 08ceffd80..dc0b4c4b8 100644
--- a/vendor/ddc/docs/CMakeLists.txt
+++ b/vendor/ddc/docs/CMakeLists.txt
@@ -19,7 +19,8 @@ set(DOXYGEN_ENUM_VALUES_PER_LINE 1)
 set(DOXYGEN_EXAMPLE_PATH "${DDC_SOURCE_DIR}/examples")
 set(DOXYGEN_EXPAND_ONLY_PREDEF YES)
 set(DOXYGEN_EXTRACT_ALL YES)
-set(DOXYGEN_EXCLUDE_SYMBOLS "detail" "details")
+set(DOXYGEN_EXCLUDE_SYMBOLS "detail")
+set(DOXYGEN_EXCLUDE_PATTERNS "*/experimental/*;*/detail/*")
 set(DOXYGEN_EXTRACT_LOCAL_CLASSES YES)
 set(DOXYGEN_FULL_PATH_NAMES NO)
 set(DOXYGEN_GENERATE_TREEVIEW YES)
@@ -39,7 +40,7 @@ set(DOXYGEN_PREDEFINED "DOXYGEN_IGNORE")
 set(DOXYGEN_PROJECT_LOGO "${CMAKE_CURRENT_SOURCE_DIR}/_template/logo.png")
 set(DOXYGEN_PROJECT_NUMBER "${DDC_VERSION}")
 set(DOXYGEN_QT_AUTOBRIEF YES)
-set(DOXYGEN_RECURSIVE NO)
+set(DOXYGEN_RECURSIVE YES)
 set(DOXYGEN_SEARCHENGINE NO)
 set(DOXYGEN_SHOW_FILES NO)
 set(DOXYGEN_SHOW_INCLUDE_FILES NO)
diff --git a/vendor/ddc/include/ddc/chunk.hpp b/vendor/ddc/include/ddc/chunk.hpp
index 14b0cc812..a0541fc0d 100644
--- a/vendor/ddc/include/ddc/chunk.hpp
+++ b/vendor/ddc/include/ddc/chunk.hpp
@@ -11,6 +11,7 @@
 
 #include "ddc/chunk_common.hpp"
 #include "ddc/chunk_span.hpp"
+#include "ddc/chunk_traits.hpp"
 #include "ddc/kokkos_allocator.hpp"
 #include "ddc/parallel_deepcopy.hpp"
 
diff --git a/vendor/ddc/include/ddc/chunk_common.hpp b/vendor/ddc/include/ddc/chunk_common.hpp
index a60efd3d2..fa1414e33 100644
--- a/vendor/ddc/include/ddc/chunk_common.hpp
+++ b/vendor/ddc/include/ddc/chunk_common.hpp
@@ -13,48 +13,12 @@
 
 #include <Kokkos_Core.hpp>
 
+#include "ddc/chunk_traits.hpp"
 #include "ddc/detail/macros.hpp"
 #include "ddc/discrete_domain.hpp"
 
 namespace ddc {
 
-template <class T>
-inline constexpr bool enable_borrowed_chunk = false;
-
-template <class T>
-inline constexpr bool enable_chunk = false;
-
-template <class T>
-inline constexpr bool is_chunk_v = enable_chunk<std::remove_const_t<std::remove_reference_t<T>>>;
-
-template <class T>
-inline constexpr bool is_borrowed_chunk_v
-        = is_chunk_v<
-                  T> && (std::is_lvalue_reference_v<T> || enable_borrowed_chunk<std::remove_cv_t<std::remove_reference_t<T>>>);
-
-template <class T>
-struct chunk_traits
-{
-    static_assert(is_chunk_v<T>);
-    using value_type
-            = std::remove_cv_t<std::remove_pointer_t<decltype(std::declval<T>().data_handle())>>;
-    using pointer_type = decltype(std::declval<T>().data_handle());
-    using reference_type = decltype(*std::declval<T>().data_handle());
-};
-
-template <class T>
-using chunk_value_t = typename chunk_traits<T>::value_type;
-
-template <class T>
-using chunk_pointer_t = typename chunk_traits<T>::pointer_type;
-
-template <class T>
-using chunk_reference_t = typename chunk_traits<T>::reference_type;
-
-template <class T>
-inline constexpr bool is_writable_chunk_v
-        = !std::is_const_v<std::remove_pointer_t<chunk_pointer_t<T>>>;
-
 /** Access the domain (or subdomain) of a view
  * @param[in]  chunk the view whose domain to access
  * @return the domain of view in the queried dimensions
@@ -351,5 +315,4 @@ class ChunkCommon<ElementType, DiscreteDomain<DDims...>, LayoutStridedPolicy>
     }
 };
 
-
 } // namespace ddc
diff --git a/vendor/ddc/include/ddc/chunk_traits.hpp b/vendor/ddc/include/ddc/chunk_traits.hpp
new file mode 100644
index 000000000..ee1f7f114
--- /dev/null
+++ b/vendor/ddc/include/ddc/chunk_traits.hpp
@@ -0,0 +1,49 @@
+// Copyright (C) The DDC development team, see COPYRIGHT.md file
+//
+// SPDX-License-Identifier: MIT
+
+#pragma once
+
+#include <type_traits>
+#include <utility>
+
+namespace ddc {
+
+template <class T>
+inline constexpr bool enable_borrowed_chunk = false;
+
+template <class T>
+inline constexpr bool enable_chunk = false;
+
+template <class T>
+inline constexpr bool is_chunk_v = enable_chunk<std::remove_const_t<std::remove_reference_t<T>>>;
+
+template <class T>
+inline constexpr bool is_borrowed_chunk_v
+        = is_chunk_v<
+                  T> && (std::is_lvalue_reference_v<T> || enable_borrowed_chunk<std::remove_cv_t<std::remove_reference_t<T>>>);
+
+template <class T>
+struct chunk_traits
+{
+    static_assert(is_chunk_v<T>);
+    using value_type
+            = std::remove_cv_t<std::remove_pointer_t<decltype(std::declval<T>().data_handle())>>;
+    using pointer_type = decltype(std::declval<T>().data_handle());
+    using reference_type = decltype(*std::declval<T>().data_handle());
+};
+
+template <class T>
+using chunk_value_t = typename chunk_traits<T>::value_type;
+
+template <class T>
+using chunk_pointer_t = typename chunk_traits<T>::pointer_type;
+
+template <class T>
+using chunk_reference_t = typename chunk_traits<T>::reference_type;
+
+template <class T>
+inline constexpr bool is_writable_chunk_v
+        = !std::is_const_v<std::remove_pointer_t<chunk_pointer_t<T>>>;
+
+} // namespace ddc
diff --git a/vendor/ddc/include/ddc/coordinate.hpp b/vendor/ddc/include/ddc/coordinate.hpp
index 1182bba3a..68fdf69eb 100644
--- a/vendor/ddc/include/ddc/coordinate.hpp
+++ b/vendor/ddc/include/ddc/coordinate.hpp
@@ -4,7 +4,10 @@
 
 #pragma once
 
+#include <type_traits>
+
 #include "ddc/detail/tagged_vector.hpp"
+#include "ddc/discrete_element.hpp"
 #include "ddc/real_type.hpp"
 
 namespace ddc {
@@ -22,9 +25,24 @@ using CoordinateElement = Real;
 template <class... CDims>
 using Coordinate = detail::TaggedVector<CoordinateElement, CDims...>;
 
+template <class... DDim, std::enable_if_t<(sizeof...(DDim) > 1), int> = 0>
+KOKKOS_FUNCTION Coordinate<typename DDim::continuous_dimension_type...> coordinate(
+        DiscreteElement<DDim...> const& c)
+{
+    return Coordinate<typename DDim::continuous_dimension_type...>(
+            coordinate(DiscreteElement<DDim>(c))...);
+}
+
+// Gives access to the type of the coordinates of a discrete element
+// Example usage : "using Coords = coordinate_of_t<DElem>;"
 template <class T>
-struct coordinate_of;
+struct coordinate_of
+{
+    static_assert(is_discrete_element_v<T>, "Parameter T must be of type DiscreteElement");
+    using type = decltype(coordinate(std::declval<T>()));
+};
 
+/// Helper type of \ref ddc::coordinate_of
 template <class T>
 using coordinate_of_t = typename coordinate_of<T>::type;
 
diff --git a/vendor/ddc/include/ddc/coordinate_md.hpp b/vendor/ddc/include/ddc/coordinate_md.hpp
deleted file mode 100644
index 18244074e..000000000
--- a/vendor/ddc/include/ddc/coordinate_md.hpp
+++ /dev/null
@@ -1,21 +0,0 @@
-// Copyright (C) The DDC development team, see COPYRIGHT.md file
-//
-// SPDX-License-Identifier: MIT
-
-#pragma once
-
-#include <type_traits>
-
-#include "ddc/coordinate.hpp"
-#include "ddc/discrete_element.hpp"
-
-namespace ddc {
-
-template <class... DDim, std::enable_if_t<(sizeof...(DDim) > 1), int> = 0>
-KOKKOS_FUNCTION Coordinate<typename DDim::continuous_dimension_type...> coordinate(
-        DiscreteElement<DDim...> const& c)
-{
-    return Coordinate<typename DDim::continuous_dimension_type...>(coordinate(select<DDim>(c))...);
-}
-
-} // namespace ddc
diff --git a/vendor/ddc/include/ddc/ddc.hpp b/vendor/ddc/include/ddc/ddc.hpp
index c0668bbdd..dae35b18a 100644
--- a/vendor/ddc/include/ddc/ddc.hpp
+++ b/vendor/ddc/include/ddc/ddc.hpp
@@ -15,10 +15,10 @@
 #include "ddc/aligned_allocator.hpp"
 #include "ddc/chunk.hpp"
 #include "ddc/chunk_span.hpp"
+#include "ddc/chunk_traits.hpp"
 #include "ddc/kokkos_allocator.hpp"
 
 // Discretizations
-#include "ddc/coordinate_md.hpp"
 #include "ddc/discrete_domain.hpp"
 #include "ddc/discrete_element.hpp"
 #include "ddc/discrete_space.hpp"
diff --git a/vendor/ddc/include/ddc/dual_discretization.hpp b/vendor/ddc/include/ddc/detail/dual_discretization.hpp
similarity index 97%
rename from vendor/ddc/include/ddc/dual_discretization.hpp
rename to vendor/ddc/include/ddc/detail/dual_discretization.hpp
index c2e955c52..ae8e091c8 100644
--- a/vendor/ddc/include/ddc/dual_discretization.hpp
+++ b/vendor/ddc/include/ddc/detail/dual_discretization.hpp
@@ -16,7 +16,7 @@
 #include <hip/hip_runtime.h>
 #endif
 
-namespace ddc {
+namespace ddc::detail {
 
 template <class DDim>
 class DualDiscretization
@@ -80,4 +80,4 @@ class DualDiscretization
     }
 };
 
-} // namespace ddc
+} // namespace ddc::detail
diff --git a/vendor/ddc/include/ddc/detail/tagged_vector.hpp b/vendor/ddc/include/ddc/detail/tagged_vector.hpp
index 9100bbfb6..a6a49c3bb 100644
--- a/vendor/ddc/include/ddc/detail/tagged_vector.hpp
+++ b/vendor/ddc/include/ddc/detail/tagged_vector.hpp
@@ -21,17 +21,17 @@ template <class, class...>
 class TaggedVector;
 
 template <class T>
-struct IsTaggedVector : std::false_type
+struct is_tagged_vector : std::false_type
 {
 };
 
 template <class ElementType, class... Tags>
-struct IsTaggedVector<TaggedVector<ElementType, Tags...>> : std::true_type
+struct is_tagged_vector<TaggedVector<ElementType, Tags...>> : std::true_type
 {
 };
 
 template <class T>
-inline constexpr bool is_tagged_vector_v = IsTaggedVector<T>::value;
+inline constexpr bool is_tagged_vector_v = is_tagged_vector<T>::value;
 
 template <class ElementType, class... Tags>
 struct ToTypeSeq<TaggedVector<ElementType, Tags...>>
@@ -224,12 +224,12 @@ KOKKOS_FUNCTION constexpr auto const& take(
 
 
 template <class T>
-class ConversionOperators
+class TaggedVectorConversionOperators
 {
 };
 
 template <class ElementType, class Tag>
-class ConversionOperators<TaggedVector<ElementType, Tag>>
+class TaggedVectorConversionOperators<TaggedVector<ElementType, Tag>>
 {
 public:
     KOKKOS_FUNCTION constexpr operator ElementType const &() const noexcept
@@ -244,9 +244,9 @@ class ConversionOperators<TaggedVector<ElementType, Tag>>
 };
 
 template <class ElementType, class... Tags>
-class TaggedVector : public ConversionOperators<TaggedVector<ElementType, Tags...>>
+class TaggedVector : public TaggedVectorConversionOperators<TaggedVector<ElementType, Tags...>>
 {
-    friend class ConversionOperators<TaggedVector<ElementType, Tags...>>;
+    friend class TaggedVectorConversionOperators<TaggedVector<ElementType, Tags...>>;
 
     using tags_seq = detail::TypeSeq<Tags...>;
 
diff --git a/vendor/ddc/include/ddc/discrete_domain.hpp b/vendor/ddc/include/ddc/discrete_domain.hpp
index e1ab56746..0f927eb38 100644
--- a/vendor/ddc/include/ddc/discrete_domain.hpp
+++ b/vendor/ddc/include/ddc/discrete_domain.hpp
@@ -23,17 +23,17 @@ template <class... DDims>
 class DiscreteDomain;
 
 template <class T>
-struct IsDiscreteDomain : std::false_type
+struct is_discrete_domain : std::false_type
 {
 };
 
 template <class... Tags>
-struct IsDiscreteDomain<DiscreteDomain<Tags...>> : std::true_type
+struct is_discrete_domain<DiscreteDomain<Tags...>> : std::true_type
 {
 };
 
 template <class T>
-inline constexpr bool is_discrete_domain_v = IsDiscreteDomain<T>::value;
+inline constexpr bool is_discrete_domain_v = is_discrete_domain<T>::value;
 
 
 namespace detail {
diff --git a/vendor/ddc/include/ddc/discrete_element.hpp b/vendor/ddc/include/ddc/discrete_element.hpp
index e9d00c23f..15b963d6f 100644
--- a/vendor/ddc/include/ddc/discrete_element.hpp
+++ b/vendor/ddc/include/ddc/discrete_element.hpp
@@ -10,7 +10,6 @@
 #include <type_traits>
 #include <utility>
 
-#include "ddc/coordinate.hpp"
 #include "ddc/detail/macros.hpp"
 #include "ddc/detail/type_seq.hpp"
 #include "ddc/discrete_vector.hpp"
@@ -21,17 +20,17 @@ template <class...>
 class DiscreteElement;
 
 template <class T>
-struct IsDiscreteElement : std::false_type
+struct is_discrete_element : std::false_type
 {
 };
 
 template <class... Tags>
-struct IsDiscreteElement<DiscreteElement<Tags...>> : std::true_type
+struct is_discrete_element<DiscreteElement<Tags...>> : std::true_type
 {
 };
 
 template <class T>
-inline constexpr bool is_discrete_element_v = IsDiscreteElement<T>::value;
+inline constexpr bool is_discrete_element_v = is_discrete_element<T>::value;
 
 
 namespace detail {
@@ -443,13 +442,4 @@ KOKKOS_FUNCTION constexpr DiscreteVector<Tags...> operator-(
     return DiscreteVector<Tags...>((uid<Tags>(lhs) - uid<Tags>(rhs))...);
 }
 
-// Gives access to the type of the coordinates of a discrete element
-// Example usage : "using Coords = coordinate_of_t<DElem>;"
-template <class... DDims>
-struct coordinate_of<ddc::DiscreteElement<DDims...>>
-{
-    // maybe a static_assert on DDims ?
-    using type = Coordinate<typename DDims::continuous_dimension_type...>;
-};
-
 } // namespace ddc
diff --git a/vendor/ddc/include/ddc/discrete_space.hpp b/vendor/ddc/include/ddc/discrete_space.hpp
index f207d3461..1ddf03830 100644
--- a/vendor/ddc/include/ddc/discrete_space.hpp
+++ b/vendor/ddc/include/ddc/discrete_space.hpp
@@ -17,7 +17,9 @@
 
 #include <Kokkos_Core.hpp>
 
+#include "ddc/detail/dual_discretization.hpp"
 #include "ddc/detail/macros.hpp"
+
 #if defined(__CUDACC__)
 #include <sstream>
 
@@ -29,7 +31,6 @@
 #include <hip/hip_runtime.h>
 #endif
 
-#include "ddc/dual_discretization.hpp"
 
 namespace ddc {
 
diff --git a/vendor/ddc/include/ddc/discrete_vector.hpp b/vendor/ddc/include/ddc/discrete_vector.hpp
index aac322b5c..446f98f2c 100644
--- a/vendor/ddc/include/ddc/discrete_vector.hpp
+++ b/vendor/ddc/include/ddc/discrete_vector.hpp
@@ -21,17 +21,17 @@ template <class...>
 class DiscreteVector;
 
 template <class T>
-struct IsDiscreteVector : std::false_type
+struct is_discrete_vector : std::false_type
 {
 };
 
 template <class... Tags>
-struct IsDiscreteVector<DiscreteVector<Tags...>> : std::true_type
+struct is_discrete_vector<DiscreteVector<Tags...>> : std::true_type
 {
 };
 
 template <class T>
-inline constexpr bool is_discrete_vector_v = IsDiscreteVector<T>::value;
+inline constexpr bool is_discrete_vector_v = is_discrete_vector<T>::value;
 
 
 namespace detail {
@@ -205,13 +205,15 @@ KOKKOS_FUNCTION constexpr auto const& take(HeadDVect const& head, TailDVects con
     DDC_IF_NVCC_THEN_POP
 }
 
+namespace detail {
+
 template <class T>
-class ConversionOperators
+class DiscreteVectorConversionOperators
 {
 };
 
 template <class Tag>
-class ConversionOperators<DiscreteVector<Tag>>
+class DiscreteVectorConversionOperators<DiscreteVector<Tag>>
 {
 public:
     KOKKOS_FUNCTION constexpr operator DiscreteVectorElement const &() const noexcept
@@ -225,8 +227,6 @@ class ConversionOperators<DiscreteVector<Tag>>
     }
 };
 
-namespace detail {
-
 /// Returns a reference to the underlying `std::array`
 template <class... Tags>
 KOKKOS_FUNCTION constexpr std::array<DiscreteVectorElement, sizeof...(Tags)>& array(
@@ -250,9 +250,9 @@ KOKKOS_FUNCTION constexpr std::array<DiscreteVectorElement, sizeof...(Tags)> con
  * Each is tagged by its associated dimensions.
  */
 template <class... Tags>
-class DiscreteVector : public ConversionOperators<DiscreteVector<Tags...>>
+class DiscreteVector : public detail::DiscreteVectorConversionOperators<DiscreteVector<Tags...>>
 {
-    friend class ConversionOperators<DiscreteVector<Tags...>>;
+    friend class detail::DiscreteVectorConversionOperators<DiscreteVector<Tags...>>;
 
     friend KOKKOS_FUNCTION constexpr std::array<DiscreteVectorElement, sizeof...(Tags)>& detail::
             array<Tags...>(DiscreteVector<Tags...>& v) noexcept;
diff --git a/vendor/ddc/include/ddc/kernels/splines.hpp b/vendor/ddc/include/ddc/kernels/splines.hpp
index 4ad1848c8..fc67f83fd 100644
--- a/vendor/ddc/include/ddc/kernels/splines.hpp
+++ b/vendor/ddc/include/ddc/kernels/splines.hpp
@@ -4,7 +4,6 @@
 
 #pragma once
 
-#include "splines/bspline.hpp"
 #include "splines/bsplines_non_uniform.hpp"
 #include "splines/bsplines_uniform.hpp"
 #include "splines/constant_extrapolation_rule.hpp"
diff --git a/vendor/ddc/include/ddc/kernels/splines/bspline.hpp b/vendor/ddc/include/ddc/kernels/splines/bspline.hpp
deleted file mode 100644
index 83f90c70c..000000000
--- a/vendor/ddc/include/ddc/kernels/splines/bspline.hpp
+++ /dev/null
@@ -1,11 +0,0 @@
-// Copyright (C) The DDC development team, see COPYRIGHT.md file
-//
-// SPDX-License-Identifier: MIT
-
-#pragma once
-namespace ddc {
-template <class Tag>
-struct BSpline
-{
-};
-} // namespace ddc
diff --git a/vendor/ddc/include/ddc/kernels/splines/bsplines_non_uniform.hpp b/vendor/ddc/include/ddc/kernels/splines/bsplines_non_uniform.hpp
index f5d5dd1c8..8160de039 100644
--- a/vendor/ddc/include/ddc/kernels/splines/bsplines_non_uniform.hpp
+++ b/vendor/ddc/include/ddc/kernels/splines/bsplines_non_uniform.hpp
@@ -11,7 +11,6 @@
 
 #include <ddc/ddc.hpp>
 
-#include "bspline.hpp"
 #include "view.hpp"
 
 namespace ddc {
@@ -47,8 +46,6 @@ class NonUniformBSplines
 public:
     using tag_type = Tag;
 
-    using continuous_dimension_type = BSpline<Tag>;
-
 
     using discrete_dimension_type = NonUniformBSplines;
 
diff --git a/vendor/ddc/include/ddc/kernels/splines/bsplines_uniform.hpp b/vendor/ddc/include/ddc/kernels/splines/bsplines_uniform.hpp
index ecdd0494f..26a781f68 100644
--- a/vendor/ddc/include/ddc/kernels/splines/bsplines_uniform.hpp
+++ b/vendor/ddc/include/ddc/kernels/splines/bsplines_uniform.hpp
@@ -10,7 +10,6 @@
 
 #include <ddc/ddc.hpp>
 
-#include "bspline.hpp"
 #include "math_tools.hpp"
 #include "view.hpp"
 
@@ -46,8 +45,6 @@ class UniformBSplines
 public:
     using tag_type = Tag;
 
-    using continuous_dimension_type = BSpline<Tag>;
-
 
     using discrete_dimension_type = UniformBSplines;
 
diff --git a/vendor/ddc/include/ddc/kernels/splines/spline_evaluator_2d.hpp b/vendor/ddc/include/ddc/kernels/splines/spline_evaluator_2d.hpp
index 792bab49a..33fc2d4ef 100644
--- a/vendor/ddc/include/ddc/kernels/splines/spline_evaluator_2d.hpp
+++ b/vendor/ddc/include/ddc/kernels/splines/spline_evaluator_2d.hpp
@@ -651,47 +651,47 @@ class SplineEvaluator2D
      */
     template <class Layout, class... CoordsDims>
     KOKKOS_INLINE_FUNCTION double eval(
-            ddc::Coordinate<CoordsDims...> const& coord_eval,
+            ddc::Coordinate<CoordsDims...> coord_eval,
             ddc::ChunkSpan<double const, bsplines_domain_type, Layout, memory_space> const
                     spline_coef) const
     {
-        ddc::Coordinate<typename interpolation_mesh_type1::continuous_dimension_type>
-                coord_eval_interpolation1(coord_eval);
-        ddc::Coordinate<typename interpolation_mesh_type2::continuous_dimension_type>
-                coord_eval_interpolation2(coord_eval);
+        using Dim1 = typename interpolation_mesh_type1::continuous_dimension_type;
+        using Dim2 = typename interpolation_mesh_type2::continuous_dimension_type;
         if constexpr (bsplines_type1::is_periodic()) {
-            if (coord_eval_interpolation1 < ddc::discrete_space<bsplines_type1>().rmin()
-                || coord_eval_interpolation1 > ddc::discrete_space<bsplines_type1>().rmax()) {
-                coord_eval_interpolation1
+            if (ddc::get<Dim1>(coord_eval) < ddc::discrete_space<bsplines_type1>().rmin()
+                || ddc::get<Dim1>(coord_eval) > ddc::discrete_space<bsplines_type1>().rmax()) {
+                ddc::get<Dim1>(coord_eval)
                         -= Kokkos::floor(
-                                   (coord_eval_interpolation1
+                                   (ddc::get<Dim1>(coord_eval)
                                     - ddc::discrete_space<bsplines_type1>().rmin())
                                    / ddc::discrete_space<bsplines_type1>().length())
                            * ddc::discrete_space<bsplines_type1>().length();
             }
-        } else {
-            if (coord_eval_interpolation1 < ddc::discrete_space<bsplines_type1>().rmin()) {
-                return m_left_extrap_rule_1(coord_eval, spline_coef);
-            }
-            if (coord_eval_interpolation1 > ddc::discrete_space<bsplines_type1>().rmax()) {
-                return m_right_extrap_rule_1(coord_eval, spline_coef);
-            }
         }
         if constexpr (bsplines_type2::is_periodic()) {
-            if (coord_eval_interpolation2 < ddc::discrete_space<bsplines_type2>().rmin()
-                || coord_eval_interpolation2 > ddc::discrete_space<bsplines_type2>().rmax()) {
-                coord_eval_interpolation2
+            if (ddc::get<Dim2>(coord_eval) < ddc::discrete_space<bsplines_type2>().rmin()
+                || ddc::get<Dim2>(coord_eval) > ddc::discrete_space<bsplines_type2>().rmax()) {
+                ddc::get<Dim2>(coord_eval)
                         -= Kokkos::floor(
-                                   (coord_eval_interpolation2
+                                   (ddc::get<Dim2>(coord_eval)
                                     - ddc::discrete_space<bsplines_type2>().rmin())
                                    / ddc::discrete_space<bsplines_type2>().length())
                            * ddc::discrete_space<bsplines_type2>().length();
             }
-        } else {
-            if (coord_eval_interpolation2 < ddc::discrete_space<bsplines_type2>().rmin()) {
+        }
+        if constexpr (!bsplines_type1::is_periodic()) {
+            if (ddc::get<Dim1>(coord_eval) < ddc::discrete_space<bsplines_type1>().rmin()) {
+                return m_left_extrap_rule_1(coord_eval, spline_coef);
+            }
+            if (ddc::get<Dim1>(coord_eval) > ddc::discrete_space<bsplines_type1>().rmax()) {
+                return m_right_extrap_rule_1(coord_eval, spline_coef);
+            }
+        }
+        if constexpr (!bsplines_type2::is_periodic()) {
+            if (ddc::get<Dim2>(coord_eval) < ddc::discrete_space<bsplines_type2>().rmin()) {
                 return m_left_extrap_rule_2(coord_eval, spline_coef);
             }
-            if (coord_eval_interpolation2 > ddc::discrete_space<bsplines_type2>().rmax()) {
+            if (ddc::get<Dim2>(coord_eval) > ddc::discrete_space<bsplines_type2>().rmax()) {
                 return m_right_extrap_rule_2(coord_eval, spline_coef);
             }
         }
@@ -699,28 +699,22 @@ class SplineEvaluator2D
                 ddc::Coordinate<
                         typename interpolation_mesh_type1::continuous_dimension_type,
                         typename interpolation_mesh_type2::continuous_dimension_type>(
-                        coord_eval_interpolation1,
-                        coord_eval_interpolation2),
+                        ddc::get<Dim1>(coord_eval),
+                        ddc::get<Dim2>(coord_eval)),
                 spline_coef);
     }
 
     /**
      * @brief Evaluate the function or its derivative at the coordinate given.
      *
-     * @param[in] coord_eval1
-     * 			The coordinate on the first dimension where we want to evaluate.
-     * @param[in] coord_eval2
-     * 			The coordinate on the second dimension where we want to evaluate.
+     * @param[in] coord_eval
+     * 			The coordinate where we want to evaluate.
      * @param[in] splne_coef
      * 			The B-splines coefficients of the function we want to evaluate.
-     * @param[out] vals1
-     * 			A ChunkSpan with the not-null values of each function of the spline in the first dimension.
-     * @param[out] vals2
-     * 			A ChunkSpan with the not-null values of each function of the spline in the second dimension.
-     * @param[in] eval_type_1
+     * @tparam EvalType1
      * 			A flag indicating if we evaluate the function or its derivative in the first dimension.
      * 			The type of this object is either `eval_type` or `eval_deriv_type`.
-     * @param[in] eval_type_2
+     * @tparam EvalType2
      * 			A flag indicating if we evaluate the function or its derivative in the second dimension.
      *          The type of this object is either `eval_type` or `eval_deriv_type`.
      */
diff --git a/vendor/ddc/include/ddc/parallel_deepcopy.hpp b/vendor/ddc/include/ddc/parallel_deepcopy.hpp
index 127b20df5..b395e9153 100644
--- a/vendor/ddc/include/ddc/parallel_deepcopy.hpp
+++ b/vendor/ddc/include/ddc/parallel_deepcopy.hpp
@@ -8,7 +8,7 @@
 
 #include <Kokkos_Core.hpp>
 
-#include "ddc/chunk_common.hpp"
+#include "ddc/chunk_traits.hpp"
 
 namespace ddc {
 
diff --git a/vendor/ddc/include/ddc/parallel_fill.hpp b/vendor/ddc/include/ddc/parallel_fill.hpp
index a0b1b050e..52020e5ea 100644
--- a/vendor/ddc/include/ddc/parallel_fill.hpp
+++ b/vendor/ddc/include/ddc/parallel_fill.hpp
@@ -8,7 +8,7 @@
 
 #include <Kokkos_Core.hpp>
 
-#include "ddc/chunk_common.hpp"
+#include "ddc/chunk_traits.hpp"
 
 namespace ddc {
 
diff --git a/vendor/ddc/include/ddc/pdi.hpp b/vendor/ddc/include/ddc/pdi.hpp
index 39058a56e..1f64ca80b 100644
--- a/vendor/ddc/include/ddc/pdi.hpp
+++ b/vendor/ddc/include/ddc/pdi.hpp
@@ -12,7 +12,7 @@
 
 #include <pdi.h>
 
-#include "ddc/chunk_common.hpp"
+#include "ddc/chunk_traits.hpp"
 
 namespace ddc {
 
diff --git a/vendor/ddc/tests/splines/extrapolation_rule.cpp b/vendor/ddc/tests/splines/extrapolation_rule.cpp
index 9c355fe99..864b51f76 100644
--- a/vendor/ddc/tests/splines/extrapolation_rule.cpp
+++ b/vendor/ddc/tests/splines/extrapolation_rule.cpp
@@ -333,13 +333,13 @@ static void ExtrapolationRuleSplineTest()
             coords_eval.domain(),
             KOKKOS_LAMBDA(Index<IDim<X, I1, I2>...> const e) {
                 coords_eval(e) = ddc::coordinate(e);
-                // Set coords_eval outside of the domain
+                // Set coords_eval outside of the domain (+1 to ensure left bound is outside domain)
                 ddc::get<I1>(coords_eval(e))
                         = xN<I1>() + (ddc::select<I1>(ddc::coordinate(e)) - x0<I1>()) + 1;
-#if defined(BC_GREVILLE)
+                // Set coords_eval outside of the domain (this point should be found on the grid in
+                // the periodic case)
                 ddc::get<I2>(coords_eval(e))
-                        = xN<I2>() + (ddc::select<I2>(ddc::coordinate(e)) - x0<I2>()) + 1;
-#endif
+                        = 2 * xN<I2>() + (ddc::select<I2>(ddc::coordinate(e)) - x0<I2>());
             });
 
 
@@ -370,10 +370,16 @@ static void ExtrapolationRuleSplineTest()
                         discrete_element_type e_batch(e);
                 double tmp;
                 if (ddc::select<I2>(coords_eval(e)) > xN<I2>()) {
+#if defined(BC_PERIODIC)
+                    tmp = vals(ddc::DiscreteElement<IDim<X, I1, I2>...>(
+                            vals.template domain<IDim<I1, I1, I2>>().back(),
+                            e_without_interest));
+#else
                     tmp = vals(ddc::DiscreteElement<IDim<X, I1, I2>...>(
                             vals.template domain<IDim<I1, I1, I2>>().back(),
                             vals.template domain<IDim<I2, I1, I2>>().back(),
                             e_batch));
+#endif
                 } else {
                     tmp = vals(ddc::DiscreteElement<IDim<X, I1, I2>...>(
                             vals.template domain<IDim<I1, I1, I2>>().back(),