#include <SparseDiagBlkMatrix.hpp>

Public Types
typedef SparseStorageType	StorageType

typedef SparseStorageType::DataType	DataType

typedef std::shared_ptr< SparseStorageType >	SparseStorageSharedPtr

typedef Array< OneD, DataType >	DataVectorType

typedef Array< OneD, const DataType >	ConstDataVectorType

typedef Array< OneD, SparseStorageSharedPtr >	SparseStorageSharedPtrVector

Public Member Functions
	NekSparseDiagBlkMatrix (const SparseStorageSharedPtrVector &sparseStoragePtrVector)

	NekSparseDiagBlkMatrix (const NekSparseDiagBlkMatrix &src)

	~NekSparseDiagBlkMatrix ()

IndexType	GetRows () const

IndexType	GetColumns () const

IndexType	GetNumNonZeroEntries ()

DataType	GetFillInRatio () const

IndexType	GetRows (int i) const
	number of rows at i-th submatrix More...

IndexType	GetColumns (int i) const
	number of columns at i-th submatrix More...

IndexType	GetNumberOfMatrixBlocks () const

IndexType	GetNumNonZeroEntries (int i) const

DataType	GetFillInRatio (int i) const

boost::call_traits< DataType >::const_reference	operator() (const IndexType row, const IndexType column) const

boost::call_traits< DataType >::const_reference	operator() (const IndexType block, const IndexType row, const IndexType column) const

void	Multiply (const DataVectorType &in, DataVectorType &out)

void	Multiply (const DataType in, DataType out)

void	MultiplySubMatrix (const IndexType blockNum, DataType in, DataType out)

size_t	GetMemoryFootprint ()

size_t	GetMemoryFootprint (IndexType i) const

unsigned long	GetMulCallsCounter () const

DataType	GetAvgRowDensity ()

DataType	GetAvgRowDensity (IndexType i) const

IndexType	GetBandwidth ()

IndexType	GetBandwidth (IndexType i)

COOMatTypeSharedPtr	GetCooStorage ()

COOMatTypeSharedPtr	GetCooStorage (IndexType i)

void	writeSparsityPatternTo (std::ostream &out, IndexType blockSize=64)

void	writeSubmatrixSparsityPatternTo (std::ostream &out, const IndexType subMatrixIdx, IndexType blockSize=64)

void	writeBlockSparsityPatternTo (std::ostream &out, const IndexType blk_row=0, const IndexType blk_col=0, IndexType blockSize=64)

void	writeSubmatrixBlockSparsityPatternTo (std::ostream &out, const IndexType subMatrixIdx, const IndexType blk_row=0, const IndexType blk_col=0, IndexType blockSize=64)
	Complementary routine to the previous. It generates exact non-zero pattern of a given block matrix entry to *this object. E.g., for a 6x6 matrix defined as A(i,i):=i, A(i,j):=0, block size = 2 and blk_row=blk_col=2 it produces 2x2 matrix with 5 and 6 along the diagonal. More...

Protected Attributes
IndexType	m_rows

IndexType	m_cols

IndexVector	m_rowoffset

unsigned long	m_mulCallsCounter

SparseStorageSharedPtrVector	m_submatrix

Detailed Description

template<typename SparseStorageType>
class Nektar::NekSparseDiagBlkMatrix< SparseStorageType >

Definition at line 61 of file SparseDiagBlkMatrix.hpp.

Member Typedef Documentation

◆ ConstDataVectorType

template<typename SparseStorageType>

typedef Array<OneD, const DataType> Nektar::NekSparseDiagBlkMatrix< SparseStorageType >::ConstDataVectorType

Definition at line 69 of file SparseDiagBlkMatrix.hpp.

◆ DataType

template<typename SparseStorageType>

typedef SparseStorageType::DataType Nektar::NekSparseDiagBlkMatrix< SparseStorageType >::DataType

Definition at line 66 of file SparseDiagBlkMatrix.hpp.

◆ DataVectorType

template<typename SparseStorageType>

typedef Array<OneD, DataType> Nektar::NekSparseDiagBlkMatrix< SparseStorageType >::DataVectorType

Definition at line 68 of file SparseDiagBlkMatrix.hpp.

◆ SparseStorageSharedPtr

template<typename SparseStorageType>

typedef std::shared_ptr<SparseStorageType> Nektar::NekSparseDiagBlkMatrix< SparseStorageType >::SparseStorageSharedPtr

Definition at line 67 of file SparseDiagBlkMatrix.hpp.

◆ SparseStorageSharedPtrVector

template<typename SparseStorageType>

typedef Array<OneD, SparseStorageSharedPtr> Nektar::NekSparseDiagBlkMatrix< SparseStorageType >::SparseStorageSharedPtrVector

Definition at line 70 of file SparseDiagBlkMatrix.hpp.

◆ StorageType

template<typename SparseStorageType>

typedef SparseStorageType Nektar::NekSparseDiagBlkMatrix< SparseStorageType >::StorageType

Definition at line 65 of file SparseDiagBlkMatrix.hpp.

Constructor & Destructor Documentation

◆ NekSparseDiagBlkMatrix() [1/2]

template<typename SparseStorageType >

Nektar::NekSparseDiagBlkMatrix< SparseStorageType >::NekSparseDiagBlkMatrix ( const SparseStorageSharedPtrVector & sparseStoragePtrVector )

Definition at line 53 of file SparseDiagBlkMatrix.cpp.

References Nektar::NekSparseDiagBlkMatrix< SparseStorageType >::m_cols, Nektar::NekSparseDiagBlkMatrix< SparseStorageType >::m_rowoffset, and Nektar::NekSparseDiagBlkMatrix< SparseStorageType >::m_rows.

                                                                     :
         m_rows(0),
         m_cols(0),
         m_rowoffset(sparseStoragePtrVector.num_elements()+1, 0.0),
         m_mulCallsCounter(0),
         m_submatrix(sparseStoragePtrVector.num_elements(),sparseStoragePtrVector)
     {
         for (int i = 0; i < sparseStoragePtrVector.num_elements(); i++)
         {
             const IndexType rows    = sparseStoragePtrVector[i]->GetRows();
             m_rows += rows;
             m_cols += sparseStoragePtrVector[i]->GetColumns();
             m_rowoffset[i+1] = m_rowoffset[i] + rows;
         }
     }

◆ NekSparseDiagBlkMatrix() [2/2]

template<typename SparseStorageType >

Nektar::NekSparseDiagBlkMatrix< SparseStorageType >::NekSparseDiagBlkMatrix ( const NekSparseDiagBlkMatrix< SparseStorageType > & src )

Definition at line 71 of file SparseDiagBlkMatrix.cpp.

                                                                                                       :
         m_rows(src.m_rows),
         m_cols(src.m_cols),
         m_rowoffset(src.m_rowoffset),
         m_mulCallsCounter(src.m_mulCallsCounter),
         m_submatrix(src.m_submatrix)
     {
     }

◆ ~NekSparseDiagBlkMatrix()

template<typename SparseStorageType >

Nektar::NekSparseDiagBlkMatrix< SparseStorageType >::~NekSparseDiagBlkMatrix ( )

Definition at line 81 of file SparseDiagBlkMatrix.cpp.

82 {

83 }

Member Function Documentation

◆ GetAvgRowDensity() [1/2]

template<typename SparseStorageType >

SparseStorageType::DataType Nektar::NekSparseDiagBlkMatrix< SparseStorageType >::GetAvgRowDensity ( )

Definition at line 249 of file SparseDiagBlkMatrix.cpp.

References Nektar::NekSparseDiagBlkMatrix< SparseStorageType >::GetNumNonZeroEntries(), Nektar::NekSparseDiagBlkMatrix< SparseStorageType >::GetRows(), and Nektar::NekSparseDiagBlkMatrix< SparseStorageType >::m_submatrix.

     {
         DataType avgRowDensity = 0.0;
         for (int i = 0; i < m_submatrix.num_elements(); i++)
         {
             avgRowDensity += (DataType) m_submatrix[i]->GetNumNonZeroEntries() /
                              (DataType) m_submatrix[i]->GetRows();
         }
         return avgRowDensity / m_submatrix.num_elements();
     }

◆ GetAvgRowDensity() [2/2]

template<typename SparseStorageType >

SparseStorageType::DataType Nektar::NekSparseDiagBlkMatrix< SparseStorageType >::GetAvgRowDensity ( IndexType i ) const

Definition at line 261 of file SparseDiagBlkMatrix.cpp.

References Nektar::NekSparseDiagBlkMatrix< SparseStorageType >::GetRows(), and Nektar::NekSparseDiagBlkMatrix< SparseStorageType >::m_submatrix.

     {
         return (DataType) m_submatrix[i]->GetNumNonZeroEntries() /
                (DataType) m_submatrix[i]->GetRows();
     }

◆ GetBandwidth() [1/2]

template<typename SparseStorageType >

IndexType Nektar::NekSparseDiagBlkMatrix< SparseStorageType >::GetBandwidth ( )

Definition at line 268 of file SparseDiagBlkMatrix.cpp.

References Nektar::NekSparseDiagBlkMatrix< SparseStorageType >::m_submatrix.

     {
         IndexType bandwidth = 0;
         for (int i = 0; i < m_submatrix.num_elements(); i++)
         {
             typename SparseStorageType::const_iterator entry = m_submatrix[i]->begin();
             for (; entry != m_submatrix[i]->end(); ++entry)
             {
                 bandwidth = (std::max)(static_cast<int>(bandwidth),
                     2*abs( static_cast<int>(entry->first.first - entry->first.second)+1));
             }
         }
         return bandwidth;
     }

◆ GetBandwidth() [2/2]

template<typename SparseStorageType >

IndexType Nektar::NekSparseDiagBlkMatrix< SparseStorageType >::GetBandwidth ( IndexType i )

Definition at line 284 of file SparseDiagBlkMatrix.cpp.

References Nektar::NekSparseDiagBlkMatrix< SparseStorageType >::m_submatrix.

     {
         IndexType bandwidth = 0;
         typename SparseStorageType::const_iterator entry = m_submatrix[i]->begin();
         for (; entry != m_submatrix[i]->end(); ++entry)
         {
             bandwidth = (std::max)(static_cast<int>(bandwidth),
                 2*abs( static_cast<int>(entry->first.first - entry->first.second)+1));
         }
         return bandwidth;
     }

◆ GetColumns() [1/2]

template<typename SparseStorageType >

IndexType Nektar::NekSparseDiagBlkMatrix< SparseStorageType >::GetColumns ( ) const

Definition at line 92 of file SparseDiagBlkMatrix.cpp.

References Nektar::NekSparseDiagBlkMatrix< SparseStorageType >::m_cols.

     {
         return m_cols;
     }

◆ GetColumns() [2/2]

template<typename SparseStorageType >

IndexType Nektar::NekSparseDiagBlkMatrix< SparseStorageType >::GetColumns ( int i ) const

number of columns at i-th submatrix

Definition at line 106 of file SparseDiagBlkMatrix.cpp.

References Nektar::NekSparseDiagBlkMatrix< SparseStorageType >::m_submatrix.

     {
         return m_submatrix[i]->GetColumns();
     }

◆ GetCooStorage() [1/2]

template<typename SparseStorageType >

COOMatTypeSharedPtr Nektar::NekSparseDiagBlkMatrix< SparseStorageType >::GetCooStorage ( )

Definition at line 311 of file SparseDiagBlkMatrix.cpp.

References Nektar::NekSparseDiagBlkMatrix< SparseStorageType >::m_submatrix, and Nektar::NekSparseDiagBlkMatrix< SparseStorageType >::writeSubmatrixSparsityPatternTo().

     {
         COOMatTypeSharedPtr coo (new COOMatType());
         IndexType row_offset = 0;
         IndexType col_offset = 0;
         for (IndexType i = 0; i < m_submatrix.num_elements(); i++)
         {
             typename SparseStorageType::const_iterator entry = m_submatrix[i]->begin();
             for (; entry != m_submatrix[i]->end(); entry++)
             {
                 IndexType loc_row = entry->first.first;
                 IndexType loc_col = entry->first.second;
                 (*coo)[std::make_pair(loc_row + row_offset, loc_col + col_offset) ] = entry->second;
             }
             row_offset += m_submatrix[i]->GetRows();
             col_offset += m_submatrix[i]->GetColumns();
         }
         return coo;
     }

◆ GetCooStorage() [2/2]

template<typename SparseStorageType >

COOMatTypeSharedPtr Nektar::NekSparseDiagBlkMatrix< SparseStorageType >::GetCooStorage ( IndexType i )

Definition at line 297 of file SparseDiagBlkMatrix.cpp.

References Nektar::NekSparseDiagBlkMatrix< SparseStorageType >::m_submatrix.

     {
         COOMatTypeSharedPtr coo (new COOMatType());
         typename SparseStorageType::const_iterator entry = m_submatrix[i]->begin();
         for (; entry != m_submatrix[i]->end(); entry++)
         {
             IndexType loc_row = entry->first.first;
             IndexType loc_col = entry->first.second;
             (*coo)[std::make_pair(loc_row, loc_col) ] = entry->second;
         }
         return coo;
     }

◆ GetFillInRatio() [1/2]

template<typename SparseStorageType >

NekSparseDiagBlkMatrix< SparseStorageType >::DataType Nektar::NekSparseDiagBlkMatrix< SparseStorageType >::GetFillInRatio ( ) const

Definition at line 142 of file SparseDiagBlkMatrix.cpp.

References Nektar::NekSparseDiagBlkMatrix< SparseStorageType >::m_submatrix.

     {
         IndexType stored = 0;
         IndexType nnz = 0;
         for (int i = 0; i < m_submatrix.num_elements(); i++)
         {
             stored += m_submatrix[i]->GetNumStoredDoubles();
             nnz    += m_submatrix[i]->GetNumNonZeroEntries();
         }
         return (DataType)stored/(DataType)nnz;
     }

◆ GetFillInRatio() [2/2]

template<typename SparseStorageType >

NekSparseDiagBlkMatrix< SparseStorageType >::DataType Nektar::NekSparseDiagBlkMatrix< SparseStorageType >::GetFillInRatio ( int i ) const

Definition at line 136 of file SparseDiagBlkMatrix.cpp.

References Nektar::NekSparseDiagBlkMatrix< SparseStorageType >::m_submatrix.

     {
         return m_submatrix[i]->GetFillInRatio();
     }

◆ GetMemoryFootprint() [1/2]

template<typename SparseStorageType >

size_t Nektar::NekSparseDiagBlkMatrix< SparseStorageType >::GetMemoryFootprint ( )

Definition at line 220 of file SparseDiagBlkMatrix.cpp.

References Nektar::NekSparseDiagBlkMatrix< SparseStorageType >::m_rowoffset, and Nektar::NekSparseDiagBlkMatrix< SparseStorageType >::m_submatrix.

     {
         size_t bytes =
                        sizeof(IndexType)*2 +      // sizes
                        sizeof(unsigned long)  +   // mulCallsCounter
                        sizeof(IndexVector) +
                        sizeof(IndexType)*m_rowoffset.capacity() + 
                        sizeof(SparseStorageSharedPtrVector) + 
                        sizeof(SparseStorageSharedPtr)*m_submatrix.capacity();
         for (int i = 0; i < m_submatrix.num_elements(); i++)
         {
             bytes += m_submatrix[i]->GetMemoryUsage();
         }
         return bytes;
     }

◆ GetMemoryFootprint() [2/2]

template<typename SparseStorageType >

size_t Nektar::NekSparseDiagBlkMatrix< SparseStorageType >::GetMemoryFootprint ( IndexType i ) const

Definition at line 237 of file SparseDiagBlkMatrix.cpp.

References Nektar::NekSparseDiagBlkMatrix< SparseStorageType >::m_submatrix.

     {
         return m_submatrix[i]->GetMemoryUsage();
     }

◆ GetMulCallsCounter()

template<typename SparseStorageType >

unsigned long Nektar::NekSparseDiagBlkMatrix< SparseStorageType >::GetMulCallsCounter ( ) const

Definition at line 243 of file SparseDiagBlkMatrix.cpp.

References Nektar::NekSparseDiagBlkMatrix< SparseStorageType >::m_mulCallsCounter.

     {
         return m_mulCallsCounter;
     }

◆ GetNumberOfMatrixBlocks()

template<typename SparseStorageType >

IndexType Nektar::NekSparseDiagBlkMatrix< SparseStorageType >::GetNumberOfMatrixBlocks ( ) const

Definition at line 112 of file SparseDiagBlkMatrix.cpp.

References Nektar::NekSparseDiagBlkMatrix< SparseStorageType >::m_submatrix.

     {
         return m_submatrix.num_elements();
     }

◆ GetNumNonZeroEntries() [1/2]

template<typename SparseStorageType >

IndexType Nektar::NekSparseDiagBlkMatrix< SparseStorageType >::GetNumNonZeroEntries ( )

Definition at line 118 of file SparseDiagBlkMatrix.cpp.

References Nektar::NekSparseDiagBlkMatrix< SparseStorageType >::m_submatrix.

Referenced by Nektar::NekSparseDiagBlkMatrix< SparseStorageType >::GetAvgRowDensity().

     {
         IndexType nnz = 0;
         for (int i = 0; i < m_submatrix.num_elements(); i++)
         {
             nnz += m_submatrix[i]->GetNumNonZeroEntries();
         }
         return nnz;
     }

◆ GetNumNonZeroEntries() [2/2]

template<typename SparseStorageType >

IndexType Nektar::NekSparseDiagBlkMatrix< SparseStorageType >::GetNumNonZeroEntries ( int i ) const

Definition at line 130 of file SparseDiagBlkMatrix.cpp.

References Nektar::NekSparseDiagBlkMatrix< SparseStorageType >::m_submatrix.

     {
         return m_submatrix[i]->GetNumNonZeroEntries();
     }

◆ GetRows() [1/2]

template<typename SparseStorageType >

IndexType Nektar::NekSparseDiagBlkMatrix< SparseStorageType >::GetRows ( ) const

Definition at line 86 of file SparseDiagBlkMatrix.cpp.

References Nektar::NekSparseDiagBlkMatrix< SparseStorageType >::m_rows.

Referenced by Nektar::NekSparseDiagBlkMatrix< SparseStorageType >::GetAvgRowDensity(), Nektar::NekSparseDiagBlkMatrix< SparseStorageType >::operator()(), Nektar::NekSparseDiagBlkMatrix< SparseStorageType >::writeBlockSparsityPatternTo(), Nektar::NekSparseDiagBlkMatrix< SparseStorageType >::writeSparsityPatternTo(), and Nektar::NekSparseDiagBlkMatrix< SparseStorageType >::writeSubmatrixBlockSparsityPatternTo().

     {
         return m_rows;
     }

◆ GetRows() [2/2]

template<typename SparseStorageType >

IndexType Nektar::NekSparseDiagBlkMatrix< SparseStorageType >::GetRows ( int i ) const

number of rows at i-th submatrix

Definition at line 99 of file SparseDiagBlkMatrix.cpp.

References Nektar::NekSparseDiagBlkMatrix< SparseStorageType >::m_submatrix.

     {
         return m_submatrix[i]->GetRows();
     }

◆ Multiply() [1/2]

template<typename SparseStorageType >

void Nektar::NekSparseDiagBlkMatrix< SparseStorageType >::Multiply	(	const DataVectorType &	in,
		DataVectorType &	out
	)

Definition at line 187 of file SparseDiagBlkMatrix.cpp.

References Nektar::NekSparseDiagBlkMatrix< SparseStorageType >::m_mulCallsCounter, Nektar::NekSparseDiagBlkMatrix< SparseStorageType >::m_rowoffset, and Nektar::NekSparseDiagBlkMatrix< SparseStorageType >::m_submatrix.

     {
         for (int i = 0; i < m_submatrix.num_elements(); ++i)
         {
             m_submatrix[i]->Multiply(&in[m_rowoffset[i]], &out[m_rowoffset[i]]);
         }
         m_mulCallsCounter++;
     }

◆ Multiply() [2/2]

template<typename SparseStorageType >

void Nektar::NekSparseDiagBlkMatrix< SparseStorageType >::Multiply	(	const DataType *	in,
		DataType *	out
	)

Definition at line 199 of file SparseDiagBlkMatrix.cpp.

References Nektar::NekSparseDiagBlkMatrix< SparseStorageType >::m_mulCallsCounter, Nektar::NekSparseDiagBlkMatrix< SparseStorageType >::m_rowoffset, Nektar::NekSparseDiagBlkMatrix< SparseStorageType >::m_submatrix, and Nektar::NekSparseDiagBlkMatrix< SparseStorageType >::MultiplySubMatrix().

     {
         for (int i = 0; i < m_submatrix.num_elements(); ++i)
         {
             m_submatrix[i]->Multiply(&in[m_rowoffset[i]], &out[m_rowoffset[i]]);
         }
         m_mulCallsCounter++;
     }

◆ MultiplySubMatrix()

template<typename SparseStorageType >

void Nektar::NekSparseDiagBlkMatrix< SparseStorageType >::MultiplySubMatrix	(	const IndexType	blockNum,
		DataType *	in,
		DataType *	out
	)

Definition at line 211 of file SparseDiagBlkMatrix.cpp.

References Nektar::NekSparseDiagBlkMatrix< SparseStorageType >::m_rowoffset, and Nektar::NekSparseDiagBlkMatrix< SparseStorageType >::m_submatrix.

Referenced by Nektar::NekSparseDiagBlkMatrix< SparseStorageType >::Multiply().

     {
         m_submatrix[blockNum]->Multiply(in + m_rowoffset[blockNum], out + m_rowoffset[blockNum]);
         // m_mulCallsCounter++;
     }

◆ operator()() [1/2]

template<typename SparseStorageType >

boost::call_traits< typename SparseStorageType::DataType >::const_reference Nektar::NekSparseDiagBlkMatrix< SparseStorageType >::operator()	(	const IndexType	row,
		const IndexType	column
	)		const

Todo:: double check, might be a bug when local_col > local column

Definition at line 164 of file SparseDiagBlkMatrix.cpp.

References Nektar::NekSparseDiagBlkMatrix< SparseStorageType >::GetRows(), and Nektar::NekSparseDiagBlkMatrix< SparseStorageType >::m_submatrix.

     {
         IndexType i = 0;
         static DataType defaultReturnValue = 0;
 
         signed int local_row = glob_row;
         signed int local_col = glob_column;
 
         while ((local_row >= m_submatrix[i]->GetRows()) &&
                (local_col >= 0))
         {
             local_row -= m_submatrix[i]->GetRows();
             local_col -= m_submatrix[i]->GetColumns();
             i++;
         }
         /// \todo double check, might be a bug when local_col > local column
         if (local_col < 0) return defaultReturnValue;
 
         return m_submatrix[i]->GetValue(local_row, local_col);
     }

◆ operator()() [2/2]

template<typename SparseStorageType >

boost::call_traits< typename SparseStorageType::DataType >::const_reference Nektar::NekSparseDiagBlkMatrix< SparseStorageType >::operator()	(	const IndexType	block,
		const IndexType	row,
		const IndexType	column
	)		const

Definition at line 156 of file SparseDiagBlkMatrix.cpp.

References Nektar::NekSparseDiagBlkMatrix< SparseStorageType >::m_submatrix.

     {
         return m_submatrix[block]->GetValue(loc_row, loc_column);
     }

◆ writeBlockSparsityPatternTo()

template<typename SparseStorageType >

void Nektar::NekSparseDiagBlkMatrix< SparseStorageType >::writeBlockSparsityPatternTo	(	std::ostream &	out,
		const IndexType	blk_row = `0`,
		const IndexType	blk_col = `0`,
		IndexType	blockSize = `64`
	)

Definition at line 466 of file SparseDiagBlkMatrix.cpp.

References Nektar::NekSparseDiagBlkMatrix< SparseStorageType >::GetRows(), and Nektar::NekSparseDiagBlkMatrix< SparseStorageType >::m_submatrix.

     {
         blockSize = (std::min)(blockSize, GetRows());
         std::vector< std::vector<int> > grid (blockSize);
         for (int row = 0; row < blockSize; row++)
         {
             grid[row].resize(blockSize,0.0);
         }
 
         IndexType row_offset = 0;
         IndexType col_offset = 0;
         for (int i = 0; i < m_submatrix.num_elements(); i++)
         {
             typename SparseStorageType::const_iterator entry =
                                         m_submatrix[i]->begin();
             typename SparseStorageType::const_iterator stop  =
                                         m_submatrix[i]->end();
             for (; entry != stop; ++entry)
             {
                 const IndexType row = entry->first.first  + row_offset;
                 const IndexType col = entry->first.second + col_offset;
 
                 if (blk_row != row / blockSize ) continue;
                 if (blk_col != col / blockSize ) continue;
                 grid[row % blockSize][col % blockSize]++;
             }
             row_offset += m_submatrix[i]->GetRows();
             col_offset += m_submatrix[i]->GetColumns();
         }
 
         for (int row = 0; row < blockSize; row++)
         {
             for (int col = 0; col < blockSize; col++)
             {
                 out << grid[row][col] << " ";
             }
             out << std::endl;
         }
     }

◆ writeSparsityPatternTo()

template<typename SparseStorageType >

void Nektar::NekSparseDiagBlkMatrix< SparseStorageType >::writeSparsityPatternTo	(	std::ostream &	out,
		IndexType	blockSize = `64`
	)

Definition at line 379 of file SparseDiagBlkMatrix.cpp.

References Nektar::NekSparseDiagBlkMatrix< SparseStorageType >::GetRows(), Nektar::NekSparseDiagBlkMatrix< SparseStorageType >::m_submatrix, and Nektar::NekSparseDiagBlkMatrix< SparseStorageType >::writeSubmatrixBlockSparsityPatternTo().

Referenced by Nektar::NekSparseDiagBlkMatrix< SparseStorageType >::writeSubmatrixSparsityPatternTo().

     {
         const int matRows = GetRows();
         const int gridRows = matRows / blockSize + (matRows % blockSize > 0);
         const int gridCols = gridRows;
 
         std::vector< std::vector<int> > grid (gridRows);
         for (int row = 0; row < gridRows; row++)
         {
             grid[row].resize(gridCols,0.0);
         }
 
         IndexType row_offset = 0;
         IndexType col_offset = 0;
         for (int i = 0; i < m_submatrix.num_elements(); i++)
         {
             typename SparseStorageType::const_iterator entry = m_submatrix[i]->begin();
             typename SparseStorageType::const_iterator stop  = m_submatrix[i]->end();
             for (; entry != stop; ++entry)
             {
                 const IndexType row = entry->first.first  + row_offset;
                 const IndexType col = entry->first.second + col_offset;
                 const int gridRow = row / blockSize;
                 const int gridCol = col / blockSize;
                 grid[gridRow][gridCol]++;
             }
             row_offset += m_submatrix[i]->GetRows();
             col_offset += m_submatrix[i]->GetColumns();
         }
 
         for (int row = 0; row < gridRows; row++)
         {
             for (int col = 0; col < gridCols; col++)
             {
                 out << grid[row][col] << " ";
             }
             out << std::endl;
         }
     }

◆ writeSubmatrixBlockSparsityPatternTo()

template<typename SparseStorageType >

void Nektar::NekSparseDiagBlkMatrix< SparseStorageType >::writeSubmatrixBlockSparsityPatternTo	(	std::ostream &	out,
		const IndexType	subMatrixIdx,
		const IndexType	blk_row = `0`,
		const IndexType	blk_col = `0`,
		IndexType	blockSize = `64`
	)

Complementary routine to the previous. It generates exact non-zero pattern of a given block matrix entry to *this object. E.g., for a 6x6 matrix defined as A(i,i):=i, A(i,j):=0, block size = 2 and blk_row=blk_col=2 it produces 2x2 matrix with 5 and 6 along the diagonal.

Definition at line 428 of file SparseDiagBlkMatrix.cpp.

References Nektar::NekSparseDiagBlkMatrix< SparseStorageType >::GetRows(), and Nektar::NekSparseDiagBlkMatrix< SparseStorageType >::m_submatrix.

Referenced by Nektar::NekSparseDiagBlkMatrix< SparseStorageType >::writeSparsityPatternTo().

     {
         blockSize = (std::min)(blockSize, m_submatrix[subMatrixIdx]->GetRows());
         std::vector< std::vector<int> > grid (blockSize);
         for (int row = 0; row < blockSize; row++)
         {
             grid[row].resize(blockSize,0.0);
         }
 
         typename SparseStorageType::const_iterator entry =
                                     m_submatrix[subMatrixIdx]->begin();
         typename SparseStorageType::const_iterator stop  =
                                     m_submatrix[subMatrixIdx]->end();
         for (; entry != stop; ++entry)
         {
             const IndexType row = entry->first.first;
             const IndexType col = entry->first.second;
 
             if (blk_row != row / blockSize ) continue;
             if (blk_col != col / blockSize ) continue;
             grid[row % blockSize][col % blockSize]++;
         }
 
         for (int row = 0; row < blockSize; row++)
         {
             for (int col = 0; col < blockSize; col++)
             {
                 out << grid[row][col] << " ";
             }
             out << std::endl;
         }
     }

◆ writeSubmatrixSparsityPatternTo()

template<typename SparseStorageType >

void Nektar::NekSparseDiagBlkMatrix< SparseStorageType >::writeSubmatrixSparsityPatternTo	(	std::ostream &	out,
		const IndexType	subMatrixIdx,
		IndexType	blockSize = `64`
	)

Definition at line 340 of file SparseDiagBlkMatrix.cpp.

References Nektar::NekSparseDiagBlkMatrix< SparseStorageType >::m_submatrix, and Nektar::NekSparseDiagBlkMatrix< SparseStorageType >::writeSparsityPatternTo().

Referenced by Nektar::NekSparseDiagBlkMatrix< SparseStorageType >::GetCooStorage().

     {
         const int matRows = m_submatrix[subMatrixIdx]->GetRows();
         const int gridRows = matRows / blockSize + (matRows % blockSize > 0);
         const int gridCols = gridRows;
 
         std::vector< std::vector<int> > grid (gridRows);
         for (int row = 0; row < gridRows; row++)
         {
             grid[row].resize(gridCols,0.0);
         }
 
         typename SparseStorageType::const_iterator entry =
                                     m_submatrix[subMatrixIdx]->begin();
         typename SparseStorageType::const_iterator stop  =
                                     m_submatrix[subMatrixIdx]->end();
         for (; entry != stop; ++entry)
         {
             const IndexType row = entry->first.first;
             const IndexType col = entry->first.second;
             const int gridRow = row / blockSize;
             const int gridCol = col / blockSize;
             grid[gridRow][gridCol]++;
         }
 
         for (int row = 0; row < gridRows; row++)
         {
             for (int col = 0; col < gridCols; col++)
             {
                 out << grid[row][col] << " ";
             }
             out << std::endl;
         }
     }

Member Data Documentation

◆ m_cols

template<typename SparseStorageType>

IndexType Nektar::NekSparseDiagBlkMatrix< SparseStorageType >::m_cols

protected

Definition at line 126 of file SparseDiagBlkMatrix.hpp.

Referenced by Nektar::NekSparseDiagBlkMatrix< SparseStorageType >::GetColumns(), and Nektar::NekSparseDiagBlkMatrix< SparseStorageType >::NekSparseDiagBlkMatrix().

◆ m_mulCallsCounter

template<typename SparseStorageType>

unsigned long Nektar::NekSparseDiagBlkMatrix< SparseStorageType >::m_mulCallsCounter

protected

Definition at line 128 of file SparseDiagBlkMatrix.hpp.

Referenced by Nektar::NekSparseDiagBlkMatrix< SparseStorageType >::GetMulCallsCounter(), and Nektar::NekSparseDiagBlkMatrix< SparseStorageType >::Multiply().

◆ m_rowoffset

template<typename SparseStorageType>

IndexVector Nektar::NekSparseDiagBlkMatrix< SparseStorageType >::m_rowoffset

protected

Definition at line 127 of file SparseDiagBlkMatrix.hpp.

Referenced by Nektar::NekSparseDiagBlkMatrix< SparseStorageType >::GetMemoryFootprint(), Nektar::NekSparseDiagBlkMatrix< SparseStorageType >::Multiply(), Nektar::NekSparseDiagBlkMatrix< SparseStorageType >::MultiplySubMatrix(), and Nektar::NekSparseDiagBlkMatrix< SparseStorageType >::NekSparseDiagBlkMatrix().

◆ m_rows

template<typename SparseStorageType>

IndexType Nektar::NekSparseDiagBlkMatrix< SparseStorageType >::m_rows

protected

Definition at line 125 of file SparseDiagBlkMatrix.hpp.

Referenced by Nektar::NekSparseDiagBlkMatrix< SparseStorageType >::GetRows(), and Nektar::NekSparseDiagBlkMatrix< SparseStorageType >::NekSparseDiagBlkMatrix().

◆ m_submatrix

template<typename SparseStorageType>

SparseStorageSharedPtrVector Nektar::NekSparseDiagBlkMatrix< SparseStorageType >::m_submatrix

protected

Definition at line 129 of file SparseDiagBlkMatrix.hpp.

Public Types

Public Member Functions

Protected Attributes

Detailed Description

template<typename SparseStorageType> class Nektar::NekSparseDiagBlkMatrix< SparseStorageType >

Member Typedef Documentation

◆ ConstDataVectorType

◆ DataType

◆ DataVectorType

◆ SparseStorageSharedPtr

◆ SparseStorageSharedPtrVector

◆ StorageType

Constructor & Destructor Documentation

◆ NekSparseDiagBlkMatrix() [1/2]

◆ NekSparseDiagBlkMatrix() [2/2]

◆ ~NekSparseDiagBlkMatrix()

Member Function Documentation

◆ GetAvgRowDensity() [1/2]

◆ GetAvgRowDensity() [2/2]

◆ GetBandwidth() [1/2]

◆ GetBandwidth() [2/2]

◆ GetColumns() [1/2]

◆ GetColumns() [2/2]

◆ GetCooStorage() [1/2]

◆ GetCooStorage() [2/2]

◆ GetFillInRatio() [1/2]

◆ GetFillInRatio() [2/2]

◆ GetMemoryFootprint() [1/2]

◆ GetMemoryFootprint() [2/2]

◆ GetMulCallsCounter()

◆ GetNumberOfMatrixBlocks()

◆ GetNumNonZeroEntries() [1/2]

◆ GetNumNonZeroEntries() [2/2]

◆ GetRows() [1/2]

◆ GetRows() [2/2]

◆ Multiply() [1/2]

◆ Multiply() [2/2]

◆ MultiplySubMatrix()

◆ operator()() [1/2]

◆ operator()() [2/2]

◆ writeBlockSparsityPatternTo()

◆ writeSparsityPatternTo()

◆ writeSubmatrixBlockSparsityPatternTo()

◆ writeSubmatrixSparsityPatternTo()

Member Data Documentation

◆ m_cols

◆ m_mulCallsCounter

◆ m_rowoffset

◆ m_rows

◆ m_submatrix

template<typename SparseStorageType>
class Nektar::NekSparseDiagBlkMatrix< SparseStorageType >