Tomáš Oberhuber
--- a/src/TNL/Matrices/Dense.hpp
+++ b/src/TNL/Matrices/Dense.hpp
+          typename Device,
+          typename Index,
+          bool RowMajorOrder,
+          typename RealAllocator >
+void
+Dense< Real, Device, Index, RowMajorOrder, RealAllocator >::
+setElements( std::initializer_list< std::initializer_list< RealType > > data )
+{
+   IndexType rows = data.size();
+   IndexType columns = 0;
+   for( auto row : data )
+      columns = max( columns, row.size() );
+   this->setDimensions( rows, columns );
+   if( ! std::is_same< DeviceType, Devices::Host >::value )
+   {
+      Dense< RealType, Devices::Host, IndexType > hostDense( rows, columns );
--- a/src/TNL/Matrices/SparseMatrix.hpp
+++ b/src/TNL/Matrices/SparseMatrix.hpp
+          typename MatrixType,
+          template< typename, typename, typename > class Segments,
+          typename RealAllocator,
+          typename IndexAllocator >
+   template< typename MapIndex,
+             typename MapValue >
+void
+SparseMatrix< Real, Device, Index, MatrixType, Segments, RealAllocator, IndexAllocator >::
+setElements( const std::map< std::pair< MapIndex, MapIndex > , MapValue >& map )
+{
+   Containers::Vector< IndexType, Devices::Host, IndexType > rowsCapacities( this->getRows(), 0 );
+   for( auto element : map )
+      rowsCapacities[ element.first.first ]++;
+   if( !std::is_same< DeviceType, Devices::Host >::value )
+   {
+      SparseMatrix< Real, Devices::Host, Index, MatrixType, Segments > hostMatrix( this->getRows(), this->getColumns() );