diff --git a/src/geometryXVx/rhs/collisions_intra.cpp b/src/geometryXVx/rhs/collisions_intra.cpp index 5cb5eada1..726e13ac3 100644 --- a/src/geometryXVx/rhs/collisions_intra.cpp +++ b/src/geometryXVx/rhs/collisions_intra.cpp @@ -72,8 +72,12 @@ CollisionsIntra::CollisionsIntra(IDomainSpXVx const& mesh, double nustar0) int const ncells(ddc::select(mesh).size() - 1); if constexpr (uniform_edge_v) { double const step(ddc::step()); - ddc::init_discrete_space< - ddc::UniformPointSampling>(vx0 - step, vxN + step, ncells + 3); + ddc::init_discrete_space( + ddc::UniformPointSampling:: + init(ddc::Coordinate(vx0 - step), + ddc::Coordinate(vxN + step), + ddc::DiscreteVector>( + ncells + 3))); } else { int const npoints(ncells + 3); std::vector> breaks(npoints); @@ -87,8 +91,12 @@ CollisionsIntra::CollisionsIntra(IDomainSpXVx const& mesh, double nustar0) if constexpr (uniform_edge_v) { double const step(ddc::step()); - ddc::init_discrete_space< - ddc::UniformPointSampling>(vx0 - step, vxN + step, ncells + 2); + ddc::init_discrete_space( + ddc::UniformPointSampling:: + init(ddc::Coordinate(vx0 - step / 2), + ddc::Coordinate(vxN + step / 2), + ddc::DiscreteVector>( + ncells + 2))); } else { int const npoints(ncells + 2); std::vector> breaks(npoints); diff --git a/src/geometryXVx/rhs/collisions_intra.hpp b/src/geometryXVx/rhs/collisions_intra.hpp index f5327928e..2d2ebf19e 100644 --- a/src/geometryXVx/rhs/collisions_intra.hpp +++ b/src/geometryXVx/rhs/collisions_intra.hpp @@ -54,7 +54,7 @@ class CollisionsIntra : public IRightHandSide } static constexpr bool uniform_edge_v - = std::is_same_v>; + = std::is_same_v>; public: using ghosted_vx_point_sampling = std::conditional_t<