Template Numerical Library version\ main:4904c12
Loading...
Searching...
No Matches
TNL::Algorithms::Segments::BiEllpackBase< Device, Index, Organization, WarpSize > Class Template Reference

BiEllpackBase serves as a base class for TNL::Algorithms::Segments::BiEllpack and TNL::Algorithms::Segments::BiEllpackView. More...

#include <TNL/Algorithms/Segments/BiEllpackBase.h>

Inheritance diagram for TNL::Algorithms::Segments::BiEllpackBase< Device, Index, Organization, WarpSize >:
[legend]
Collaboration diagram for TNL::Algorithms::Segments::BiEllpackBase< Device, Index, Organization, WarpSize >:
[legend]

Public Types

using ConstOffsetsView = typename OffsetsView::ConstViewType
 The type for representing the constant vector view with offsets of the BiEllpack groups.
using DeviceType = Device
 The device where the segments are operating.
using IndexType = std::remove_const_t< Index >
 The type used for indexing of segments elements.
using OffsetsView = Containers::VectorView< Index, DeviceType, IndexType >
 The type for representing the vector view with offsets of the BiEllpack groups.
using SegmentViewType = BiEllpackSegmentView< IndexType, Organization, WarpSize >
 Accessor type for one particular segment.

Public Member Functions

__cuda_callable__ BiEllpackBase ()=default
 Default constructor with no parameters to create empty segments view.
__cuda_callable__ BiEllpackBase (BiEllpackBase &&) noexcept=default
 Move constructor.
__cuda_callable__ BiEllpackBase (const BiEllpackBase &)=default
 Copy constructor.
__cuda_callable__ BiEllpackBase (IndexType size, IndexType storageSize, OffsetsView segmentsPermutation, OffsetsView groupPointers)
 Constructor that initializes segments based all necessary data.
template<typename Function>
void forAllElements (Function &&function) const
template<typename Condition, typename Function>
void forAllElementsIf (Condition condition, Function function) const
template<typename Function>
void forAllSegments (Function &&function) const
template<typename Array, typename Function>
void forElements (const Array &segmentIndexes, Function function) const
template<typename Array, typename Function>
void forElements (const Array &segmentIndexes, Index begin, Index end, Function function) const
template<typename Function>
void forElements (IndexType begin, IndexType end, Function &&function) const
template<typename Condition, typename Function>
void forElementsIf (IndexType begin, IndexType end, Condition condition, Function function) const
template<typename Function>
void forSegments (IndexType begin, IndexType end, Function &&function) const
__cuda_callable__ IndexType getElementCount () const
 Returns the number of elements managed by all segments.
__cuda_callable__ IndexType getGlobalIndex (Index segmentIdx, Index localIdx) const
 Computes the global index of an element managed by the segments.
__cuda_callable__ OffsetsView getGroupPointersView ()
 Returns a view containing the pointers to the beginning of each segment.
__cuda_callable__ ConstOffsetsView getGroupPointersView () const
 Returns a view containing the pointers to the beginning of each segment.
__cuda_callable__ IndexType getSegmentCount () const
 Returns the number of segments.
__cuda_callable__ IndexType getSegmentsCount () const
 Returns the number of segments. Deprecated, use getSegmentCount().
__cuda_callable__ IndexType getSegmentSize (IndexType segmentIdx) const
 Returns the size of a particular segment denoted by segmentIdx.
__cuda_callable__ OffsetsView getSegmentsPermutationView ()
 Returns a view containing the permutation of segments.
__cuda_callable__ ConstOffsetsView getSegmentsPermutationView () const
 Returns a view containing the permutation of segments.
__cuda_callable__ SegmentViewType getSegmentView (IndexType segmentIdx) const
 Returns a segment view (i.e., a segment accessor) for the specified segment index.
__cuda_callable__ IndexType getSize () const
 Returns the number of elements managed by all segments.
__cuda_callable__ IndexType getStorageSize () const
 Returns number of elements that needs to be allocated by a container connected to this segments.
__cuda_callable__ IndexType getVirtualSegments () const
BiEllpackBaseoperator= (BiEllpackBase &&)=delete
 Move-assignment operator.
BiEllpackBaseoperator= (const BiEllpackBase &)=delete
 Copy-assignment operator.
void printStructure (std::ostream &str) const
 Prints the structure of the segments to the specified output stream.

Static Public Member Functions

static constexpr int getLogWarpSize ()
 Returns the log2 of the warp size used for the segments.
static constexpr ElementsOrganization getOrganization ()
 Returns the data layout for the biellpack format (it is always row-major order).
static std::string getSegmentsType ()
 Returns string with the segments type.
static std::string getSerializationType ()
 Returns string with the serialization type.
static constexpr int getWarpSize ()
 Returns the warp size used for the segments.
static constexpr bool havePadding ()
 This function denotes that the biellpack format uses padding elements.

Protected Member Functions

__cuda_callable__ void bind (IndexType size, IndexType storageSize, OffsetsView segmentsPermutation, OffsetsView groupPointers)
 Re-initializes the internal attributes of the base class.
__cuda_callable__ IndexType getVirtualSegments (IndexType segmentsCount) const

Protected Attributes

OffsetsView groupPointers
OffsetsView segmentsPermutation
IndexType size = 0
IndexType storageSize = 0

Detailed Description

template<typename Device, typename Index, ElementsOrganization Organization, int WarpSize>
class TNL::Algorithms::Segments::BiEllpackBase< Device, Index, Organization, WarpSize >

BiEllpackBase serves as a base class for TNL::Algorithms::Segments::BiEllpack and TNL::Algorithms::Segments::BiEllpackView.

Template Parameters
Deviceis type of device where the segments will be operating.
Indexis type for indexing of the elements managed by the segments.
Organizationis the organization of the elements in the segments—either row-major or column-major order.
WarpSizeis the warp size used for the segments.

Member Function Documentation

◆ bind()

template<typename Device, typename Index, ElementsOrganization Organization, int WarpSize>
__cuda_callable__ void TNL::Algorithms::Segments::BiEllpackBase< Device, Index, Organization, WarpSize >::bind ( IndexType size,
IndexType storageSize,
OffsetsView segmentsPermutation,
OffsetsView groupPointers )
protected

Re-initializes the internal attributes of the base class.

Note that this function is protected to ensure that the user cannot modify the base class of segments. For the same reason, in future code development we also need to make sure that all non-const functions in the base class return by value and not by reference.

◆ getGlobalIndex()

template<typename Device, typename Index, ElementsOrganization Organization, int WarpSize>
__cuda_callable__ IndexType TNL::Algorithms::Segments::BiEllpackBase< Device, Index, Organization, WarpSize >::getGlobalIndex ( Index segmentIdx,
Index localIdx ) const
nodiscard

Computes the global index of an element managed by the segments.

The global index serves as a reference to the element within its container.

Parameters
segmentIdxThe index of the segment containing the element.
localIdxThe local index of the element within the segment.
Returns
The global index of the element.

◆ getSegmentsType()

template<typename Device, typename Index, ElementsOrganization Organization, int WarpSize>
std::string TNL::Algorithms::Segments::BiEllpackBase< Device, Index, Organization, WarpSize >::getSegmentsType ( )
staticnodiscard

Returns string with the segments type.

Example
#include <iostream>
#include <functional>
#include <TNL/Algorithms/Segments/CSR.h>
#include <TNL/Algorithms/Segments/Ellpack.h>
#include <TNL/Devices/Host.h>
#include <TNL/Devices/Cuda.h>
template< typename Segments >
void
SegmentsExample()
{
/***
* Create segments and print the segments type.
*/
Segments segments;
std::cout << "The segments type is: " << segments.getSegmentsType() << '\n';
}
int
main( int argc, char* argv[] )
{
std::cout << "Example of CSR segments on host:\n";
SegmentsExample< TNL::Algorithms::Segments::CSR< TNL::Devices::Host, int > >();
std::cout << "Example of Ellpack segments on host:\n";
SegmentsExample< TNL::Algorithms::Segments::Ellpack< TNL::Devices::Host, int > >();
#ifdef __CUDACC__
std::cout << "Example of CSR segments on CUDA GPU:\n";
SegmentsExample< TNL::Algorithms::Segments::CSR< TNL::Devices::Cuda, int > >();
std::cout << "Example of Ellpack segments on CUDA GPU:\n";
SegmentsExample< TNL::Algorithms::Segments::Ellpack< TNL::Devices::Cuda, int > >();
#endif
return EXIT_SUCCESS;
}
Namespace for the segments data structures.
Definition _NamespaceDoxy.h:7
Output
Example of CSR segments on host:
The segments type is: CSR
Example of Ellpack segments on host:
The segments type is: Ellpack
Example of CSR segments on CUDA GPU:
The segments type is: CSR
Example of Ellpack segments on CUDA GPU:
The segments type is: Ellpack

◆ getSegmentView()

template<typename Device, typename Index, ElementsOrganization Organization, int WarpSize>
__cuda_callable__ SegmentViewType TNL::Algorithms::Segments::BiEllpackBase< Device, Index, Organization, WarpSize >::getSegmentView ( IndexType segmentIdx) const
nodiscard

Returns a segment view (i.e., a segment accessor) for the specified segment index.

Parameters
segmentIdxThe index of the requested segment.
Returns
The segment view of the specified segment.
Example
#include <iostream>
#include <TNL/Containers/Vector.h>
#include <TNL/Algorithms/Segments/CSR.h>
#include <TNL/Algorithms/Segments/Ellpack.h>
#include <TNL/Algorithms/SequentialFor.h>
#include <TNL/Devices/Host.h>
#include <TNL/Devices/Cuda.h>
template< typename Segments >
void
SegmentsExample()
{
using Device = typename Segments::DeviceType;
/***
* Create segments with given segments sizes.
*/
const int size( 5 );
Segments segments{ 1, 2, 3, 4, 5 };
auto view = segments.getView();
/***
* Print the elements mapping using segment view.
*/
std::cout << "Mapping of local indexes to global indexes:\n";
auto f = [ = ] __cuda_callable__( int segmentIdx )
{
printf( "Segment idx. %d: ", segmentIdx ); // printf works even in GPU kernels
auto segment = view.getSegmentView( segmentIdx );
for( auto element : segment )
printf( "%d -> %d \t", element.localIndex(), element.globalIndex() );
printf( "\n" );
};
}
int
main( int argc, char* argv[] )
{
std::cout << "Example of CSR segments on host:\n";
SegmentsExample< TNL::Algorithms::Segments::CSR< TNL::Devices::Host, int > >();
std::cout << "Example of Ellpack segments on host:\n";
SegmentsExample< TNL::Algorithms::Segments::Ellpack< TNL::Devices::Host, int > >();
#ifdef __CUDACC__
std::cout << "Example of CSR segments on CUDA GPU:\n";
SegmentsExample< TNL::Algorithms::Segments::CSR< TNL::Devices::Cuda, int > >();
std::cout << "Example of Ellpack segments on CUDA GPU:\n";
SegmentsExample< TNL::Algorithms::Segments::Ellpack< TNL::Devices::Cuda, int > >();
#endif
return EXIT_SUCCESS;
}
#define __cuda_callable__
Definition Macros.h:49
static void exec(Index start, Index end, Function f)
Static method for execution of the loop.
Definition SequentialFor.h:36
Output
Example of CSR segments on host:
Mapping of local indexes to global indexes:
Segment idx. 0: 0 -> 0
Segment idx. 1: 0 -> 1 1 -> 2
Segment idx. 2: 0 -> 3 1 -> 4 2 -> 5
Segment idx. 3: 0 -> 6 1 -> 7 2 -> 8 3 -> 9
Segment idx. 4: 0 -> 10 1 -> 11 2 -> 12 3 -> 13 4 -> 14
Example of Ellpack segments on host:
Mapping of local indexes to global indexes:
Segment idx. 0: 0 -> 0 1 -> 1 2 -> 2 3 -> 3 4 -> 4
Segment idx. 1: 0 -> 5 1 -> 6 2 -> 7 3 -> 8 4 -> 9
Segment idx. 2: 0 -> 10 1 -> 11 2 -> 12 3 -> 13 4 -> 14
Segment idx. 3: 0 -> 15 1 -> 16 2 -> 17 3 -> 18 4 -> 19
Segment idx. 4: 0 -> 20 1 -> 21 2 -> 22 3 -> 23 4 -> 24
Example of CSR segments on CUDA GPU:
Mapping of local indexes to global indexes:
Segment idx. 0: 0 -> 0
Segment idx. 1: 0 -> 1 1 -> 2
Segment idx. 2: 0 -> 3 1 -> 4 2 -> 5
Segment idx. 3: 0 -> 6 1 -> 7 2 -> 8 3 -> 9
Segment idx. 4: 0 -> 10 1 -> 11 2 -> 12 3 -> 13 4 -> 14
Example of Ellpack segments on CUDA GPU:
Mapping of local indexes to global indexes:
Segment idx. 0: 0 -> 0 1 -> 32 2 -> 64 3 -> 96 4 -> 128
Segment idx. 1: 0 -> 1 1 -> 33 2 -> 65 3 -> 97 4 -> 129
Segment idx. 2: 0 -> 2 1 -> 34 2 -> 66 3 -> 98 4 -> 130
Segment idx. 3: 0 -> 3 1 -> 35 2 -> 67 3 -> 99 4 -> 131
Segment idx. 4: 0 -> 4 1 -> 36 2 -> 68 3 -> 100 4 -> 132

◆ getSerializationType()

template<typename Device, typename Index, ElementsOrganization Organization, int WarpSize>
std::string TNL::Algorithms::Segments::BiEllpackBase< Device, Index, Organization, WarpSize >::getSerializationType ( )
staticnodiscard

Returns string with the serialization type.

Example
#include <iostream>
#include <functional>
#include <TNL/Algorithms/Segments/CSR.h>
#include <TNL/Algorithms/Segments/Ellpack.h>
#include <TNL/Devices/Host.h>
#include <TNL/Devices/Cuda.h>
template< typename Segments >
void
SegmentsExample()
{
/***
* Create segments and print the serialization type.
*/
Segments segments;
std::cout << "The serialization type is: " << segments.getSerializationType() << '\n';
}
int
main( int argc, char* argv[] )
{
std::cout << "Example of CSR segments on host:\n";
SegmentsExample< TNL::Algorithms::Segments::CSR< TNL::Devices::Host, int > >();
std::cout << "Example of Ellpack segments on host:\n";
SegmentsExample< TNL::Algorithms::Segments::Ellpack< TNL::Devices::Host, int > >();
#ifdef __CUDACC__
std::cout << "Example of CSR segments on CUDA GPU:\n";
SegmentsExample< TNL::Algorithms::Segments::CSR< TNL::Devices::Cuda, int > >();
std::cout << "Example of Ellpack segments on CUDA GPU:\n";
SegmentsExample< TNL::Algorithms::Segments::Ellpack< TNL::Devices::Cuda, int > >();
#endif
return EXIT_SUCCESS;
}
Output
Example of CSR segments on host:
The serialization type is: CSR< int >
Example of Ellpack segments on host:
The serialization type is: Ellpack< int, RowMajorOrder, 32 >
Example of CSR segments on CUDA GPU:
The serialization type is: CSR< int >
Example of Ellpack segments on CUDA GPU:
The serialization type is: Ellpack< int, ColumnMajorOrder, 32 >

The documentation for this class was generated from the following file: