DynRankView¶
Header file: <Kokkos_DynRankView.hpp>
Usage¶
DynRankView is a potential reference counted multidimensional array with
compile time layouts and memory space. Its semantics are similar to that of std::shared_ptr.
The DynRankView differs from the [[View|Kokkos::View]] in that its rank is not provided
with the DataType template parameter; it is determined dynamically based on the
number of extent arguments passed to the constructor. The rank has an upper bound of 7 dimensions.
Description¶
-
template<class DataType, class LayoutType, class MemorySpace, class MemoryTraits>
class DynRankView;¶ - Template Parameters:
DataType –
Defines the fundamental scalar type of the
DynRankViewAttention
This parameter is mandatory
The basic structure is
ScalarType. Examples:double: aDynRankViewofdouble, dimensions are passed as arguments to the constructor, the number of which determine the rank.
LayoutType –
Determines the mapping of indices into the underlying 1D memory storage
Important
This parameter is optional
Custom Layouts can be implemented, but Kokkos comes with some built-in ones:
LayoutRight: Strides increase from the right most to the left most dimension. The last dimension has a stride of one. This corresponds to how C multi dimensional arrays ([][][]) are laid out in memory.LayoutLeft: Strides increase from the left most to the right most dimension. The first dimension has a stride of one. This is the layout Fortran uses for its arrays.LayoutStride: Strides can be arbitrary for each dimension.
MemorySpace –
Controls the storage location
Important
This parameter is optional
If omitted, the default memory space of the default execution space is used (i.e.
Kokkos::DefaultExecutionSpace::memory_space)MemoryTraits –
Finer-grained control on the memory access
Important
This parameter is optional
Unmanaged: The DynRankView will not be reference counted. The allocation has to be provided to the constructor.Atomic: All accesses to the view will use atomic operations.RandomAccess: Hint that the view is used in a random access manner. If the view is alsoconst, this will trigger special load operations on GPUs (i.e. texture fetches).Restrict: There is no aliasing of the view by other data structures in the current scope.
Important
Template parameters other than
DataTypeare optional, but ordering is enforced. That means for example thatLayoutTypecan be omitted but if bothMemorySpaceandMemoryTraitsare specified,MemorySpacemust come beforeMemoryTraits.Public Static Variables
rank: Rank of the view (i.e. the dimensionality).rank_dynamic: Number of runtime determined dimensions.reference_type_is_lvalue_reference: Whether the reference type is a C++ lvalue reference.
Public Data Types Typedefs
-
type data_type¶
The
DataTypeof the DynRankView.
-
type const_data_type¶
Const version of
DataType, same asdata_typeif that is already const.
-
type non_const_data_type¶
Non-const version of
DataType, same asdata_typeif that is already non-const.
-
type scalar_array_type¶
If
DataTyperepresents some properly specialised array data type such as Sacado FAD types,scalar_array_typeis the underlying fundamental scalar type.
-
type const_scalar_array_type¶
The const version of
scalar_array_type, same asscalar_array_typeif that is already const
-
type non_const_scalar_array_type¶
The non-Const version of
scalar_array_type, same asscalar_array_typeif that is already non-const.
Public Scalar Typedefs
-
type value_type¶
The
data_typestripped of its array specifiers, i.e. the scalar type of the data the view is referencing (e.g. ifdata_typeisconst int*******,value_typeisconst int).
-
type const_value_type¶
Const version of
value_type.
-
type non_const_value_type¶
Non-const version of
value_type.
Public Spaces Typedefs
-
type execution_space¶
Execution space associated with the view, will be used for performing view initialization, and certain deep_copy operations.
-
type memory_space¶
Data storage location type.
-
type device_type¶
The compound type defined by
Device<execution_space,memory_space>.
-
type memory_traits¶
The memory traits of the view.
-
type host_mirror_space¶
Host accessible memory space used in
HostMirror.
Public View Typedefs
-
type non_const_type¶
This view type with all template parameters explicitly defined.
-
type const_type¶
This view type with all template parameters explicitly defined using a
constdata type.
-
type HostMirror¶
Compatible view type with the same
DataTypeandLayoutTypestored in host accessible memory space.
Public Data Handles Typedefs
-
type reference_type¶
Return type of the view access operators.
-
type pointer_type¶
Pointer to scalar type.
Other Public Typedefs
-
type array_layout¶
The layout of the
DynRankView.
-
type size_type¶
Index type associated with the memory space of this view.
-
type dimension¶
An integer array like type, able to represent the extents of the view.
-
type specialize¶
A specialization tag used for partial specialization of the mapping construct underlying a Kokkos
DynRankView.
Constructors
-
DynRankView()¶
Default constructor. No allocations are made, no reference counting happens. All extents are zero and its data pointer is
nullptrand its rank is set to 0.
-
DynRankView(const DynRankView<DT, Prop...> &rhs)¶
Copy constructor with compatible DynRankViews. Follows DynRankView assignment rules.
-
DynRankView(DynRankView &&rhs)¶
Move constructor.
-
DynRankView(const std::string &name, const IntType&... indices)¶
Requires:
array_layout::is_regular == trueStandard allocating constructor.
name: a user provided label, which is used for profiling and debugging purposes. Names are not required to be unique.indices: runtime dimensions of the view.
-
DynRankView(const std::string &name, const array_layout &layout)¶
Standard allocating constructor.
name: a user provided label, which is used for profiling and debugging purposes. Names are not required to be unique.layout: an instance of a layout class.
-
DynRankView(const AllocProperties &prop, const IntType&... indices)¶
Requires:
array_layout::is_regular == trueAllocating constructor with allocation properties. An allocation properties object is returned by the
view_allocfunction.indices: runtime dimensions of the view.
-
DynRankView(const AllocProperties &prop, const array_layout &layout)¶
Allocating constructor with allocation properties and a layout object.
layout: an instance of a layout class.
-
DynRankView(const pointer_type &ptr, const IntType&... indices)¶
Requires:
array_layout::is_regular == trueUnmanaged data wrapping constructor.
ptr: pointer to a user provided memory allocation. Must provide storage of sizeDynRankView::required_allocation_size(n0,...,nR).indices: runtime dimensions of the view.
-
DynRankView(const pointer_type &ptr, const array_layout &layout)¶
Unmanaged data wrapper constructor.
ptr: pointer to a user provided memory allocation. Must provide storage of sizeDynRankView::required_allocation_size(layout)(NEEDS TO BE IMPLEMENTED)layout: an instance of a layout class.
-
DynRankView(const ScratchSpace &space, const IntType&... indices)¶
Requires:
sizeof(IntType...)==rank_dynamic()andarray_layout::is_regular == trueConstructor which acquires memory from a Scratch Memory handle.
space: scratch memory handle. Typically returned fromteam_handlesinTeamPolicykernels.indices: runtime dimensions of the view.
-
DynRankView(const ScratchSpace &space, const array_layout &layout)¶
Constructor which acquires memory from a Scratch Memory handle.
space: scratch memory handle. Typically returned fromteam_handlesinTeamPolicykernels.layout: an instance of a layout class.
-
DynRankView(const DynRankView<DT, Prop...> &rhs, Args... args)¶
Subview constructor. See
subviewfunction for arguments.
Data Access Functions
-
reference_type operator()(const IntType&... indices) const¶
Returns a value of
reference_typewhich may or not be reference itself. The number of index arguments must match therankof the view. See notes onreference_typefor properties of the return type.
-
reference_type access(const IntType &i0 = 0, const IntType &i1 = 0, const IntType &i2 = 0, const IntType &i3 = 0, const IntType &i4 = 0, const IntType &i5 = 0, const IntType &i6 = 0) const¶
Returns a value of
reference_typewhich may or not be reference itself. The number of index arguments must be equal or larger than therankof the view. Index arguments beyondrankmust be0, which will be enforced ifKOKKOS_DEBUGis defined. See notes onreference_typefor properties of the return type.
Data Layout, Dimensions, Strides
-
constexpr array_layout layout() const¶
Returns the layout object. Can be used to to construct other views with the same dimensions.
-
template<class iType>
constexpr size_t extent(const iType &dim) const¶ Returns the extent of the specified dimension.
iTypemust be an integral type, anddimmust be smaller thanrank.
-
template<class iType>
constexpr int extent_int(const iType &dim) const¶ Returns the extent of the specified dimension as an
int.iTypemust be an integral type, anddimmust be smaller thanrank. Compared toextentthis function can be useful on architectures whereintoperations are more efficient thansize_t. It also may eliminate the need for type casts in applications which otherwise perform all index operations withint.
-
template<class iType>
constexpr size_t stride(const iType &dim) const¶ Returns the stride of the specified dimension.
iTypemust be an integral type, anddimmust be smaller thanrank. Example:a.stride(3) == (&a(i0,i1,i2,i3+1,i4)-&a(i0,i1,i2,i3,i4))
-
constexpr size_t stride_0() const¶
Return the stride of dimension 0.
-
constexpr size_t stride_1() const¶
Return the stride of dimension 1.
-
constexpr size_t stride_2() const¶
Return the stride of dimension 2.
-
constexpr size_t stride_3() const¶
Return the stride of dimension 3.
-
constexpr size_t stride_4() const¶
Return the stride of dimension 4.
-
constexpr size_t stride_5() const¶
Return the stride of dimension 5.
-
constexpr size_t stride_6() const¶
Return the stride of dimension 6.
-
constexpr size_t stride_7() const¶
Return the stride of dimension 7.
-
constexpr size_t span() const¶
Return the memory span in elements between the element with the lowest and the highest address. This can be larger than the product of extents due to padding, and or non-contiguous data layout as for example
LayoutStrideallows.
-
constexpr pointer_type data() const¶
Return the pointer to the underlying data allocation.
-
bool span_is_contiguous() const¶
Whether the span is contiguous (i.e. whether every memory location between in span belongs to the index space covered by the view).
-
static constexpr size_t required_allocation_size(size_t N0 = 0, size_t N1 = 0, size_t N2 = 0, size_t N3 = 0, size_t N4 = 0, size_t N5 = 0, size_t N6 = 0);¶
Returns the number of bytes necessary for an unmanaged view of the provided dimensions. This function is only valid if
array_layout::is_regular == true.
-
static constexpr size_t required_allocation_size(const array_layout &layout);¶
- Returns:
the number of bytes necessary for an unmanaged view of the provided layout.
Other Public Methods
-
int use_count() const¶
- Returns:
the current reference count of the underlying allocation.
-
const char *label() const;¶
- Returns:
the label of the
DynRankView.
-
constexpr unsigned rank() const¶
- Returns:
the dynamic rank of the
DynRankView.
-
constexpr bool is_allocated() const¶
- Returns:
true if the view points to a valid memory location. This function works for both managed and unmanaged views. With the unmanaged view, there is no guarantee that referenced address is valid, only that it is a non-null pointer.
Assignment Rules¶
Assignment rules cover the assignment operator as well as copy constructors. We aim at making all logically legal assignments possible, while intercepting illegal assignments if possible at compile time, otherwise at runtime. In the following, we use DstType and SrcType as the type of the destination view and source view respectively. dst_view and src_view refer to the runtime instances of the destination and source views, i.e.:
ScrType src_view(...);
DstType dst_view(src_view);
dst_view = src_view;
The following conditions must be met at and are evaluated at compile time:
DstType::rank == SrcType::rankDstType::non_const_value_typeis the same asSrcType::non_const_value_typeIf
std::is_const<SrcType::value_type>::value == truethanstd::is_const<DstType::value_type>::value == true.MemorySpaceAccess<DstType::memory_space,SrcType::memory_space>::assignable == true
Furthermore there are rules which must be met if DstType::array_layout is not the same as SrcType::array_layout. These rules only cover cases where both layouts are one of LayoutLeft , LayoutRight or LayoutStride
- If neither
DstType::array_layoutnorSrcType::array_layoutisLayoutStride: If
DstType::rank > 1thanDstType::array_layoutmust be the same asSrcType::array_layout.
- If neither
- If either
DstType::array_layoutorSrcType::array_layoutisLayoutStride For each dimension
kit must hold thatdst_view.extent(k) == src_view.extent(k)
- If either
Examples¶
#include<Kokkos_Core.hpp>
#include<cstdio>
int main(int argc, char* argv[]) {
Kokkos::initialize(argc,argv);
int N0 = atoi(argv[1]);
int N1 = atoi(argv[2]);
Kokkos::DynRankView<double> a("A",N0);
Kokkos::DynRankView<double> b("B",N1);
Kokkos::parallel_for("InitA", N0, KOKKOS_LAMBDA (const int& i) {
a(i) = i;
});
Kokkos::parallel_for("InitB", N1, KOKKOS_LAMBDA (const int& i) {
b(i) = i;
});
Kokkos::DynRankView<double,Kokkos::LayoutLeft> c("C",N0,N1);
{
Kokkos::DynRankView<const double> const_a(a);
Kokkos::DynRankView<const double> const_b(b);
Kokkos::parallel_for("SetC", Kokkos::MDRangePolicy<Kokkos::Rank<2,Kokkos::Iterate::Left>>({0,0},{N0,N1}),
KOKKOS_LAMBDA (const int& i0, const int& i1) {
c(i0,i1) = a(i0) * b(i1);
});
}
Kokkos::finalize();
}