Added method print to Segments for better use with output streams.

Documentation/Examples/Algorithms/Segments/SegmentsExample_CSR_constructor_1.cpp

@@ -35,7 +35,7 @@ void SegmentsExample()
     * Print the data managed by the segments.
     */
    auto fetch = [=] __cuda_callable__ ( int globalIdx ) -> double { return data_view[ globalIdx ]; };
-   printSegments( segments, fetch, std::cout );
+   std::cout << segments.print( fetch ) << std::endl;
 }
 
 int main( int argc, char* argv[] )

Documentation/Examples/Algorithms/Segments/SegmentsExample_CSR_constructor_2.cpp

@@ -34,7 +34,7 @@ void SegmentsExample()
     * Print the data managed by the segments.
     */
    auto fetch = [=] __cuda_callable__ ( int globalIdx ) -> double { return data_view[ globalIdx ]; };
-   printSegments( segments, fetch, std::cout );
+   std::cout << segments.print( fetch ) << std::endl;
 }
 
 int main( int argc, char* argv[] )

Documentation/Examples/Algorithms/Segments/SegmentsExample_CSR_forElements.cpp

@@ -33,7 +33,7 @@ void SegmentsExample()
     * Print the data managed by the segments.
     */
    auto fetch = [=] __cuda_callable__ ( int globalIdx ) -> double { return data_view[ globalIdx ]; };
-   printSegments( segments, fetch, std::cout );
+   std::cout << segments.print( fetch ) << std::endl;
 }
 
 int main( int argc, char* argv[] )

Documentation/Examples/Algorithms/Segments/SegmentsExample_CSR_forSegments.cpp

@@ -36,7 +36,7 @@ void SegmentsExample()
     * Print the data managed by the segments.
     */
    auto fetch = [=] __cuda_callable__ ( int globalIdx ) -> double { return data_view[ globalIdx ]; };
-   printSegments( segments, fetch, std::cout );
+   std::cout << segments.print( fetch ) << std::endl;
 }
 
 int main( int argc, char* argv[] )

Documentation/Examples/Algorithms/Segments/SegmentsExample_General.cpp

@@ -36,7 +36,7 @@ void SegmentsExample()
     * Print the data managed by the segments.
     */
    auto fetch = [=] __cuda_callable__ ( IndexType globalIdx ) -> double { return data_view[ globalIdx ]; };
-   printSegments( segments, fetch, std::cout );
+   std::cout << segments.print( fetch ) << std::endl;
 
    /***
     * Compute sums of elements in particular segments.

Documentation/Examples/Algorithms/Segments/SegmentsExample_forElements.cpp

@@ -34,7 +34,7 @@ void SegmentsExample()
    std::cout << "Data setup with no check ... " << std::endl;
    std::cout << "Array: " << data << std::endl;
    auto fetch = [=] __cuda_callable__ ( int globalIdx ) -> double { return data_view[ globalIdx ]; };
-   printSegments( segments, fetch, std::cout ) << std::endl;
+   std::cout << segments.print( fetch ) << std::endl;
 
    /***
     * Insert data into particular segments.
@@ -50,7 +50,7 @@ void SegmentsExample()
     */
    std::cout << "Data setup with check for padding elements..." << std::endl;
    std::cout << "Array: " << data << std::endl;
-   printSegments( segments, fetch, std::cout ) << std::endl;
+   std::cout << segments.print( fetch ) << std::endl;
 }
 
 int main( int argc, char* argv[] )

Documentation/Examples/Algorithms/Segments/SegmentsExample_forSegments-1.cpp

@@ -39,7 +39,7 @@ void SegmentsExample()
     * Print the data managed by the segments.
     */
    auto fetch = [=] __cuda_callable__ ( int globalIdx ) -> double { return data_view[ globalIdx ]; };
-   printSegments( segments, fetch, std::cout );
+   std::cout << segments.print( fetch ) << std::endl;
 }
 
 int main( int argc, char* argv[] )

Documentation/Examples/Algorithms/Segments/SegmentsExample_forSegments-2.cpp

@@ -33,7 +33,7 @@ void SegmentsExample()
     */
    std::cout << "Values of elements after intial setup: " << std::endl;
    auto fetch = [=] __cuda_callable__ ( int globalIdx ) -> double { return data_view[ globalIdx ]; };
-   printSegments( segments, fetch, std::cout );
+   std::cout << segments.print( fetch );
 
    /***
     * Divide elements in each segment by a sum of all elements in the segment
@@ -55,7 +55,7 @@ void SegmentsExample()
     * Print the data managed by the segments.
     */
    std::cout << "Value of elements after dividing by sum in each segment:" << std::endl;
-   printSegments( segments, fetch, std::cout );
+   std::cout << segments.print( fetch ) << std::endl;
 }
 
 int main( int argc, char* argv[] )

Documentation/Examples/Algorithms/Segments/SegmentsExample_reduceSegments.cpp

@@ -36,7 +36,7 @@ void SegmentsExample()
     */
    std::cout << "Values of elements after intial setup: " << std::endl;
    auto fetch = [=] __cuda_callable__ ( int globalIdx ) -> double { return data_view[ globalIdx ]; };
-   printSegments( segments, fetch, std::cout );
+   std::cout << segments.print( fetch ) << std::endl;
 
    /***
     * Compute sums of elements in each segment.

Documentation/Examples/Algorithms/Segments/SegmentsPrintingExample-2.cpp

@@ -32,7 +32,7 @@ void SegmentsExample()
     */
    auto data_view = data.getView();
    auto fetch = [=] __cuda_callable__ ( int globalIdx ) -> double { return data_view[ globalIdx ]; };
-   printSegments( segments, fetch, std::cout );
+   std::cout << segments.print( fetch ) << std::endl;
    std::cout << std::endl;
 }
 

src/TNL/Algorithms/Segments/BiEllpack.h

@@ -135,6 +135,9 @@ class BiEllpack
 
       void load(File &file);
 
+      template< typename Fetch >
+      SegmentsPrinter< BiEllpack, Fetch > print( Fetch&& fetch ) const;
+
       void printStructure(std::ostream &str) const;
 
       // TODO: nvcc needs this public because of lambda function used inside

src/TNL/Algorithms/Segments/BiEllpack.hpp

@@ -581,6 +581,19 @@ load( File& file )
         >> this->groupPointers;
 }
 
+template< typename Device,
+          typename Index,
+          typename IndexAllocator,
+          ElementsOrganization Organization,
+          int WarpSize >
+      template< typename Fetch >
+auto
+BiEllpack< Device, Index, IndexAllocator, Organization, WarpSize >::
+print( Fetch&& fetch ) const -> SegmentsPrinter< BiEllpack, Fetch >
+{
+   return SegmentsPrinter< BiEllpack, Fetch >( *this, fetch );
+}
+
 template< typename Device,
           typename Index,
           typename IndexAllocator,

src/TNL/Algorithms/Segments/BiEllpackView.h

@@ -138,6 +138,9 @@ class BiEllpackView
 
       void load( File& file );
 
+      template< typename Fetch >
+      SegmentsPrinter< BiEllpackView, Fetch > print( Fetch&& fetch ) const;
+
       void printStructure( std::ostream& str ) const;
 
    protected:

src/TNL/Algorithms/Segments/BiEllpackView.hpp

@@ -477,6 +477,18 @@ save( File& file ) const
         << this->groupPointers;
 }
 
+template< typename Device,
+          typename Index,
+          ElementsOrganization Organization,
+          int WarpSize >
+      template< typename Fetch >
+auto
+BiEllpackView< Device, Index, Organization, WarpSize >::
+print( Fetch&& fetch ) const -> SegmentsPrinter< BiEllpackView, Fetch >
+{
+   return SegmentsPrinter< BiEllpackView, Fetch >( *this, fetch );
+}
+
 template< typename Device,
           typename Index,
           ElementsOrganization Organization,

src/TNL/Algorithms/Segments/CSR.h

@@ -486,6 +486,27 @@ class CSR
        */
       void load( File& file );
 
+      /**
+       * \brief Return simple proxy object for insertion to output stream.
+       *
+       * The proxy object serves for wrapping segments with lambda function mediating access to data managed by the segments.
+       *
+       * \tparam Fetch is type of lambda function for data access.
+       * \param fetch is an instance of lambda function for data access. It is supposed to be defined as
+       *
+       * ```
+       * auto fetch = [=] __cuda_callable__ ( IndexType globalIdx ) -> ValueType { return data_view[ globalIdx ]; };
+       * ```
+       * \return Proxy object for insertion to output stream.
+       *
+       * \par Example
+       * \include Algorithms/Segments/SegmentsPrintingExample-2.cpp
+       * \par Output
+       * \include SegmentsPrintingExample-2.out
+       */
+      template< typename Fetch >
+      SegmentsPrinter< CSR, Fetch > print( Fetch&& fetch ) const;
+
    protected:
 
       OffsetsContainer offsets;

src/TNL/Algorithms/Segments/CSR.hpp

@@ -368,6 +368,18 @@ load( File& file )
    this->kernel.init( this->offsets );
 }
 
+template< typename Device,
+          typename Index,
+          typename Kernel,
+          typename IndexAllocator >
+      template< typename Fetch >
+auto
+CSR< Device, Index, Kernel, IndexAllocator >::
+print( Fetch&& fetch ) const -> SegmentsPrinter< CSR, Fetch >
+{
+   return SegmentsPrinter< CSR, Fetch >( *this, fetch );
+}
+
       } // namespace Segments
    }  // namespace Algorithms
 } // namespace TNL

src/TNL/Algorithms/Segments/CSRView.h

@@ -139,6 +139,9 @@ class CSRView
 
       void load( File& file );
 
+      template< typename Fetch >
+      SegmentsPrinter< CSRView, Fetch > print( Fetch&& fetch ) const;
+
    protected:
 
       OffsetsView offsets;

src/TNL/Algorithms/Segments/CSRView.hpp

@@ -318,6 +318,18 @@ load( File& file )
    this->kernel.init( this->offsets );
 }
 
+template< typename Device,
+          typename Index,
+          typename Kernel >
+      template< typename Fetch >
+auto
+CSRView< Device, Index, Kernel >::
+print( Fetch&& fetch ) const -> SegmentsPrinter< CSRView, Fetch >
+{
+   return SegmentsPrinter< CSRView, Fetch >( *this, fetch );
+}
+
+
       } // namespace Segments
    }  // namespace Containers
 } // namespace TNL

src/TNL/Algorithms/Segments/ChunkedEllpack.h

@@ -131,6 +131,9 @@ class ChunkedEllpack
 
       void load( File& file );
 
+      template< typename Fetch >
+      SegmentsPrinter< ChunkedEllpack, Fetch > print( Fetch&& fetch ) const;
+
       void printStructure( std::ostream& str ); // TODO const;
 
    protected:

src/TNL/Algorithms/Segments/ChunkedEllpack.hpp

@@ -534,6 +534,18 @@ load( File& file )
    file.load( &this->numberOfSlices );
 }
 
+template< typename Device,
+          typename Index,
+          typename IndexAllocator,
+          ElementsOrganization Organization >
+      template< typename Fetch >
+auto
+ChunkedEllpack< Device, Index, IndexAllocator, Organization >::
+print( Fetch&& fetch ) const -> SegmentsPrinter< ChunkedEllpack, Fetch >
+{
+   return SegmentsPrinter< ChunkedEllpack, Fetch >( *this, fetch );
+}
+
 template< typename Device,
           typename Index,
           typename IndexAllocator,