#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 59 of file SparseDiagBlkMatrix.hpp.

Member Typedef Documentation

◆ ConstDataVectorType

template<typename SparseStorageType >

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

Definition at line 66 of file SparseDiagBlkMatrix.hpp.

◆ DataType

template<typename SparseStorageType >

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

Definition at line 63 of file SparseDiagBlkMatrix.hpp.

◆ DataVectorType

template<typename SparseStorageType >

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

Definition at line 65 of file SparseDiagBlkMatrix.hpp.

◆ SparseStorageSharedPtr

template<typename SparseStorageType >

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

Definition at line 64 of file SparseDiagBlkMatrix.hpp.

◆ SparseStorageSharedPtrVector

template<typename SparseStorageType >

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

Definition at line 67 of file SparseDiagBlkMatrix.hpp.

◆ StorageType

template<typename SparseStorageType >

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

Definition at line 62 of file SparseDiagBlkMatrix.hpp.

Constructor & Destructor Documentation

◆ NekSparseDiagBlkMatrix() [1/2]

template<typename SparseStorageType >

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

Definition at line 51 of file SparseDiagBlkMatrix.cpp.

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

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

◆ NekSparseDiagBlkMatrix() [2/2]

template<typename SparseStorageType >

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

Definition at line 67 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 75 of file SparseDiagBlkMatrix.cpp.

76{

77}

Member Function Documentation

◆ GetAvgRowDensity() [1/2]

template<typename SparseStorageType >

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

Definition at line 251 of file SparseDiagBlkMatrix.cpp.

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

◆ GetAvgRowDensity() [2/2]

template<typename SparseStorageType >

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

Definition at line 264 of file SparseDiagBlkMatrix.cpp.

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

References Nektar::NekSparseDiagBlkMatrix< SparseStorageType >::GetNumNonZeroEntries().

◆ GetBandwidth() [1/2]

template<typename SparseStorageType >

IndexType Nektar::NekSparseDiagBlkMatrix< SparseStorageType >::GetBandwidth

Definition at line 271 of file SparseDiagBlkMatrix.cpp.

{
    IndexType bandwidth = 0;
    for (int i = 0; i < m_submatrix.size(); 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;
}

References tinysimd::abs().

◆ GetBandwidth() [2/2]

template<typename SparseStorageType >

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

Definition at line 291 of file SparseDiagBlkMatrix.cpp.

{
    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;
}

References tinysimd::abs().

◆ GetColumns() [1/2]

template<typename SparseStorageType >

IndexType Nektar::NekSparseDiagBlkMatrix< SparseStorageType >::GetColumns

Definition at line 86 of file SparseDiagBlkMatrix.cpp.

{
    return m_cols;
}

Referenced by Nektar::operator<<().

◆ GetColumns() [2/2]

template<typename SparseStorageType >

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

number of columns at i-th submatrix

Definition at line 100 of file SparseDiagBlkMatrix.cpp.

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

◆ GetCooStorage() [1/2]

template<typename SparseStorageType >

COOMatTypeSharedPtr Nektar::NekSparseDiagBlkMatrix< SparseStorageType >::GetCooStorage

Definition at line 321 of file SparseDiagBlkMatrix.cpp.

{
    COOMatTypeSharedPtr coo(new COOMatType());
    IndexType row_offset = 0;
    IndexType col_offset = 0;
    for (IndexType i = 0; i < m_submatrix.size(); 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 306 of file SparseDiagBlkMatrix.cpp.

{
    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

Definition at line 140 of file SparseDiagBlkMatrix.cpp.

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

References Nektar::NekSparseDiagBlkMatrix< SparseStorageType >::GetNumNonZeroEntries().

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

◆ GetFillInRatio() [2/2]

template<typename SparseStorageType >

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

Definition at line 133 of file SparseDiagBlkMatrix.cpp.

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

References Nektar::NekSparseDiagBlkMatrix< SparseStorageType >::GetFillInRatio().

◆ GetMemoryFootprint() [1/2]

template<typename SparseStorageType >

size_t Nektar::NekSparseDiagBlkMatrix< SparseStorageType >::GetMemoryFootprint

Definition at line 220 of file SparseDiagBlkMatrix.cpp.

{
    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.size(); 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 236 of file SparseDiagBlkMatrix.cpp.

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

◆ GetMulCallsCounter()

template<typename SparseStorageType >

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

Definition at line 243 of file SparseDiagBlkMatrix.cpp.

{
    return m_mulCallsCounter;
}

◆ GetNumberOfMatrixBlocks()

template<typename SparseStorageType >

IndexType Nektar::NekSparseDiagBlkMatrix< SparseStorageType >::GetNumberOfMatrixBlocks

Definition at line 106 of file SparseDiagBlkMatrix.cpp.

{
    return m_submatrix.size();
}

◆ GetNumNonZeroEntries() [1/2]

template<typename SparseStorageType >

IndexType Nektar::NekSparseDiagBlkMatrix< SparseStorageType >::GetNumNonZeroEntries

Definition at line 113 of file SparseDiagBlkMatrix.cpp.

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

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

◆ GetNumNonZeroEntries() [2/2]

template<typename SparseStorageType >

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

Definition at line 125 of file SparseDiagBlkMatrix.cpp.

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

◆ GetRows() [1/2]

template<typename SparseStorageType >

IndexType Nektar::NekSparseDiagBlkMatrix< SparseStorageType >::GetRows

Definition at line 80 of file SparseDiagBlkMatrix.cpp.

{
    return m_rows;
}

Referenced by Nektar::operator<<().

◆ GetRows() [2/2]

template<typename SparseStorageType >

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

number of rows at i-th submatrix

Definition at line 93 of file SparseDiagBlkMatrix.cpp.

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

◆ Multiply() [1/2]

template<typename SparseStorageType >

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

Definition at line 200 of file SparseDiagBlkMatrix.cpp.

{
    for (int i = 0; i < m_submatrix.size(); ++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 DataVectorType &	in,
		DataVectorType &	out
	)

Definition at line 189 of file SparseDiagBlkMatrix.cpp.

{
    for (int i = 0; i < m_submatrix.size(); ++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.

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

◆ operator()() [1/2]

template<typename SparseStorageType >

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

Definition at line 155 of file SparseDiagBlkMatrix.cpp.

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

◆ operator()() [2/2]

template<typename SparseStorageType >

boost::call_traits< typenameSparseStorageType::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.

{
    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);
}

◆ 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 482 of file SparseDiagBlkMatrix.cpp.

{
    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.size(); 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 388 of file SparseDiagBlkMatrix.cpp.

{
    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.size(); 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 437 of file SparseDiagBlkMatrix.cpp.

{
    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 351 of file SparseDiagBlkMatrix.cpp.

{
    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 124 of file SparseDiagBlkMatrix.hpp.

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

◆ m_mulCallsCounter

template<typename SparseStorageType >

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

protected

Definition at line 126 of file SparseDiagBlkMatrix.hpp.

◆ m_rowoffset

template<typename SparseStorageType >

IndexVector Nektar::NekSparseDiagBlkMatrix< SparseStorageType >::m_rowoffset

protected

Definition at line 125 of file SparseDiagBlkMatrix.hpp.

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

◆ m_rows

template<typename SparseStorageType >

IndexType Nektar::NekSparseDiagBlkMatrix< SparseStorageType >::m_rows

protected

Definition at line 123 of file SparseDiagBlkMatrix.hpp.

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

◆ m_submatrix

template<typename SparseStorageType >

SparseStorageSharedPtrVector Nektar::NekSparseDiagBlkMatrix< SparseStorageType >::m_submatrix

protected

Definition at line 127 of file SparseDiagBlkMatrix.hpp.

Public Types

Public Member Functions

Protected Attributes

Detailed Description

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