13#include <Kokkos_Core.hpp>
15#include "ddc/coordinate.hpp"
16#include "ddc/discrete_domain.hpp"
17#include "ddc/discrete_element.hpp"
18#include "ddc/discrete_space.hpp"
19#include "ddc/discrete_vector.hpp"
20#include "ddc/real_type.hpp"
25
30 using continuous_dimension_type = CDim;
32 using continuous_element_type = Coordinate<CDim>;
44 template <
class MemorySpace>
47 template <
class OMemorySpace>
51 continuous_element_type m_origin {0};
62 template <
class OriginMemorySpace>
63 explicit Impl(
Impl<OriginMemorySpace>
const& impl)
64 : m_origin(impl.m_origin)
72
73
74
75
76 Impl(continuous_element_type origin, Real step) : m_origin(origin), m_step(step)
92 return discrete_element_type {0};
103 coordinate(discrete_element_type
const& icoord)
const noexcept
105 return m_origin + continuous_element_type(icoord.uid()) * m_step;
110
111
112
113
114
115
116 static std::tuple<
Impl<Kokkos::HostSpace>, discrete_domain_type>
init(
117 continuous_element_type a,
118 continuous_element_type b,
119 discrete_vector_type n)
123 Impl<Kokkos::HostSpace> disc(a, continuous_element_type {(b - a) / (n - 1)});
124 discrete_domain_type domain {disc.front(), n};
125 return std::make_tuple(std::move(disc), std::move(domain));
129
130
131
132
133
134
135
136
138 Impl<Kokkos::HostSpace>,
139 discrete_domain_type,
140 discrete_domain_type,
141 discrete_domain_type,
142 discrete_domain_type>
144 continuous_element_type a,
145 continuous_element_type b,
146 discrete_vector_type n,
147 discrete_vector_type n_ghosts_before,
148 discrete_vector_type n_ghosts_after)
150 using discrete_domain_type = discrete_domain_type;
153 Real discretization_step {(b - a) / (n - 1)};
154 Impl<Kokkos::HostSpace>
155 disc(a - n_ghosts_before.value() * discretization_step, discretization_step);
156 discrete_domain_type ghosted_domain
157 = discrete_domain_type(disc.front(), n + n_ghosts_before + n_ghosts_after);
158 discrete_domain_type pre_ghost
159 = discrete_domain_type(ghosted_domain.front(), n_ghosts_before);
160 discrete_domain_type main_domain
161 = discrete_domain_type(ghosted_domain.front() + n_ghosts_before, n);
162 discrete_domain_type post_ghost
163 = discrete_domain_type(main_domain.back() + 1, n_ghosts_after);
164 return std::make_tuple(
166 std::move(main_domain),
167 std::move(ghosted_domain),
168 std::move(pre_ghost),
169 std::move(post_ghost));
173
174
175
176
177
178
179
181 Impl<Kokkos::HostSpace>,
182 discrete_domain_type,
183 discrete_domain_type,
184 discrete_domain_type,
185 discrete_domain_type>
187 continuous_element_type a,
188 continuous_element_type b,
189 discrete_vector_type n,
190 discrete_vector_type n_ghosts)
192 return init_ghosted(a, b, n, n_ghosts, n_ghosts);
213 is_uniform_sampling_v<
typename DDimImpl::discrete_dimension_type>,
215std::ostream& operator<<(std::ostream& out, DDimImpl
const& mesh)
217 return out <<
"UniformPointSampling( origin=" << mesh.origin() <<
", step=" << mesh.step()
233 return discrete_space<DDim>().origin();
241 return discrete_space<DDim>().front();
248 return discrete_space<DDim>().step();
266 return coordinate(d.front());
272 return coordinate(d.back());
278 return rmax(d) - rmin(d);
286template <
class... DDims>
288 : std::conditional_t<(is_uniform_sampling_v<DDims> && ...), std::true_type, std::false_type>
friend class DiscreteDomain
Definition: discrete_domain.hpp:53
Definition: aligned_allocator.hpp:11
KOKKOS_FUNCTION std::enable_if_t< is_uniform_sampling_v< DDim >, typename DDim::discrete_element_type > front() noexcept
Lower bound index of the mesh.
Definition: uniform_point_sampling.hpp:239
KOKKOS_FUNCTION std::enable_if_t< is_uniform_sampling_v< DDim >, typename DDim::continuous_element_type > origin() noexcept
Lower bound index of the mesh.
Definition: uniform_point_sampling.hpp:231
constexpr bool is_uniform_sampling_v
Definition: uniform_point_sampling.hpp:207
constexpr bool is_uniform_domain_v
Definition: uniform_point_sampling.hpp:293
KOKKOS_FUNCTION Coordinate< CDim > distance_at_left(DiscreteElement< UniformPointSampling< CDim > >)
Definition: uniform_point_sampling.hpp:252
KOKKOS_FUNCTION std::enable_if_t< is_uniform_sampling_v< DDim >, Real > step() noexcept
Spacing step of the mesh.
Definition: uniform_point_sampling.hpp:246
KOKKOS_FUNCTION Coordinate< CDim > coordinate(DiscreteElement< UniformPointSampling< CDim > > const &c)
Definition: uniform_point_sampling.hpp:222
KOKKOS_FUNCTION Coordinate< CDim > rmin(DiscreteDomain< UniformPointSampling< CDim > > const &d)
Definition: uniform_point_sampling.hpp:264
KOKKOS_FUNCTION Coordinate< CDim > distance_at_right(DiscreteElement< UniformPointSampling< CDim > >)
Definition: uniform_point_sampling.hpp:258
KOKKOS_FUNCTION Coordinate< CDim > rlength(DiscreteDomain< UniformPointSampling< CDim > > const &d)
Definition: uniform_point_sampling.hpp:276
KOKKOS_FUNCTION Coordinate< CDim > rmax(DiscreteDomain< UniformPointSampling< CDim > > const &d)
Definition: uniform_point_sampling.hpp:270
Definition: uniform_point_sampling.hpp:283