10#include <initializer_list>
18#include <Kokkos_Core.hpp>
20#include "ddc/coordinate.hpp"
21#include "ddc/discrete_domain.hpp"
22#include "ddc/discrete_element.hpp"
23#include "ddc/discrete_space.hpp"
24#include "ddc/discrete_vector.hpp"
30struct NonUniformPointSamplingBase
41 using continuous_dimension_type = CDim;
46 template <
class DDim,
class MemorySpace>
49 template <
class ODDim,
class OMemorySpace>
52 Kokkos::View<Coordinate<CDim>*, MemorySpace> m_points;
54 DiscreteElement<DDim> m_reference;
61 using discrete_element_type = DiscreteElement<DDim>;
68 Impl(std::initializer_list<Coordinate<CDim>>
const points)
69 :
Impl(points.begin(), points.end())
74 template <
class InputRange>
75 explicit Impl(InputRange
const& points) :
Impl(points.begin(), points.end())
80 template <
class InputIt>
81 Impl(InputIt
const points_begin, InputIt
const points_end)
82 : m_reference(create_reference_discrete_element<DDim>())
84 using view_type = Kokkos::View<Coordinate<CDim>*, MemorySpace>;
85 if (!std::is_sorted(points_begin, points_end)) {
86 throw std::runtime_error(
"Input points must be sorted");
89 std::vector<Coordinate<CDim>> host_points(points_begin, points_end);
90 m_points = view_type(
"NonUniformPointSampling::points", host_points.size());
91 Kokkos::deep_copy(m_points, view_type(host_points.data(), host_points.size()));
94 template <
class OriginMemorySpace>
95 explicit Impl(
Impl<DDim, OriginMemorySpace>
const& impl)
96 : m_points(Kokkos::create_mirror_view_and_copy(MemorySpace(), impl.m_points))
97 , m_reference(impl.m_reference)
113 return m_points.size();
124 discrete_element_type
const& icoord)
const noexcept
126 return m_points((icoord -
front()).value());
131
132
133
134
135 template <
class DDim,
class InputRange>
136 static std::tuple<
typename DDim::
template Impl<DDim, Kokkos::HostSpace>,
DiscreteDomain<DDim>>
137 init(InputRange
const& non_uniform_points)
139 assert(!non_uniform_points.empty());
141 typename DDim::
template Impl<DDim, Kokkos::HostSpace> disc(non_uniform_points);
143 return std::make_tuple(std::move(disc), std::move(domain));
147
148
149
150
151
152 template <
class DDim,
class InputRange>
154 typename DDim::
template Impl<DDim, Kokkos::HostSpace>,
160 InputRange
const& domain_r,
161 InputRange
const& pre_ghost_r,
162 InputRange
const& post_ghost_r)
164 assert(!domain_r.empty());
170 std::vector<
typename InputRange::value_type> full_domain;
172 std::copy(pre_ghost_r.begin(), pre_ghost_r.end(), std::back_inserter(full_domain));
173 std::copy(domain_r.begin(), domain_r.end(), std::back_inserter(full_domain));
174 std::copy(post_ghost_r.begin(), post_ghost_r.end(), std::back_inserter(full_domain));
176 typename DDim::
template Impl<DDim, Kokkos::HostSpace> disc(full_domain);
178 DiscreteDomain<DDim> ghosted_domain(disc.front(), n + n_ghosts_before + n_ghosts_after);
179 DiscreteDomain<DDim> pre_ghost = ghosted_domain.take_first(n_ghosts_before);
180 DiscreteDomain<DDim> main_domain = ghosted_domain.remove(n_ghosts_before, n_ghosts_after);
181 DiscreteDomain<DDim> post_ghost = ghosted_domain.take_last(n_ghosts_after);
182 return std::make_tuple(
184 std::move(main_domain),
185 std::move(ghosted_domain),
186 std::move(pre_ghost),
187 std::move(post_ghost));
193 :
public std::is_base_of<detail::NonUniformPointSamplingBase, DDim>::type
203 is_non_uniform_point_sampling_v<
typename DDimImpl::discrete_dimension_type>,
206std::ostream& operator<<(std::ostream& out, DDimImpl
const& mesh)
208 return out <<
"NonUniformPointSampling(" << mesh.size() <<
")";
211template <
class DDim, std::enable_if_t<is_non_uniform_point_sampling_v<DDim>,
int> = 0>
213 DiscreteElement<DDim>
const& c)
215 return discrete_space<DDim>().coordinate(c);
218template <
class DDim, std::enable_if_t<is_non_uniform_point_sampling_v<DDim>,
int> = 0>
220 DiscreteElement<DDim> i)
222 return coordinate(i) - coordinate(i - 1);
225template <
class DDim, std::enable_if_t<is_non_uniform_point_sampling_v<DDim>,
int> = 0>
227 DiscreteElement<DDim> i)
229 return coordinate(i + 1) - coordinate(i);
232template <
class DDim, std::enable_if_t<is_non_uniform_point_sampling_v<DDim>,
int> = 0>
236 return coordinate(d.front());
239template <
class DDim, std::enable_if_t<is_non_uniform_point_sampling_v<DDim>,
int> = 0>
243 return coordinate(d.back());
246template <
class DDim, std::enable_if_t<is_non_uniform_point_sampling_v<DDim>,
int> = 0>
250 return rmax(d) - rmin(d);
friend class DiscreteDomain
KOKKOS_FUNCTION constexpr bool operator!=(DiscreteVector< OTags... > const &rhs) const noexcept
The top-level namespace of DDC.
constexpr bool is_non_uniform_point_sampling_v
KOKKOS_FUNCTION Coordinate< typename DDim::continuous_dimension_type > rmax(DiscreteDomain< DDim > const &d)
KOKKOS_FUNCTION Coordinate< typename DDim::continuous_dimension_type > rlength(DiscreteDomain< DDim > const &d)
KOKKOS_FUNCTION Coordinate< typename DDim::continuous_dimension_type > distance_at_left(DiscreteElement< DDim > i)
KOKKOS_FUNCTION Coordinate< typename DDim::continuous_dimension_type > rmin(DiscreteDomain< DDim > const &d)
KOKKOS_FUNCTION Coordinate< typename DDim::continuous_dimension_type > coordinate(DiscreteElement< DDim > const &c)
KOKKOS_FUNCTION Coordinate< typename DDim::continuous_dimension_type > distance_at_right(DiscreteElement< DDim > i)