Nektar::NekMatrix< NekMatrix< DataType, InnerMatrixType >, BlockMatrixTag > Class Template Reference

#include <BlockMatrix.hpp>

Inheritance diagram for Nektar::NekMatrix< NekMatrix< DataType, InnerMatrixType >, BlockMatrixTag >:

Classes
class	iterator_base

Public Types
typedef NekMatrix< DataType, InnerMatrixType >	InnerType
typedef NekMatrix< InnerType, BlockMatrixTag >	ThisType
typedef InnerType::NumberType	NumberType
typedef ConstMatrix< NumberType >	BaseType
typedef InnerType::GetValueType	GetValueType
typedef InnerType::ConstGetValueType	ConstGetValueType
typedef iterator_base< ThisType >	iterator
typedef iterator_base< const ThisType >	const_iterator

Public Member Functions
	NekMatrix (MatrixStorage type=eFULL)
	NekMatrix (unsigned int numberOfBlockRows, unsigned int numberOfBlockColumns, unsigned int rowsPerBlock, unsigned int columnsPerBlock, MatrixStorage type=eFULL)
	NekMatrix (unsigned int numberOfBlockRows, unsigned int numberOfBlockColumns, const unsigned int *rowsPerBlock, const unsigned int *columnsPerBlock, MatrixStorage type=eFULL)
	NekMatrix (unsigned int numberOfBlockRows, unsigned int numberOfBlockColumns, const Array< OneD, const unsigned int > &rowsPerBlock, const Array< OneD, const unsigned int > &columnsPerBlock, MatrixStorage type=eFULL)
	NekMatrix (const Array< OneD, const unsigned int > &rowsPerBlock, const Array< OneD, const unsigned int > &columnsPerBlock, MatrixStorage type=eFULL)
	NekMatrix (const ThisType &rhs)
unsigned int	GetRequiredStorageSize () const
unsigned int	CalculateBlockIndex (unsigned int row, unsigned int column) const
const InnerType *	GetBlockPtr (unsigned int row, unsigned int column) const
std::shared_ptr< const InnerType >	GetBlock (unsigned int row, unsigned int column) const
std::shared_ptr< InnerType > &	GetBlock (unsigned int row, unsigned int column)
void	SetBlock (unsigned int row, unsigned int column, std::shared_ptr< InnerType > &m)
ConstGetValueType	operator() (unsigned int row, unsigned int col) const
unsigned int	GetStorageSize () const
unsigned int	GetNumberOfBlockRows () const
unsigned int	GetNumberOfBlockColumns () const
unsigned int	GetNumberOfRowsInBlockRow (unsigned int blockRow) const
unsigned int	GetNumberOfColumnsInBlockColumn (unsigned int blockCol) const
void	GetBlockSizes (Array< OneD, unsigned int > &rowSizes, Array< OneD, unsigned int > &colSizes) const
iterator	begin ()
iterator	end ()
const_iterator	begin () const
const_iterator	end () const
Public Member Functions inherited from Nektar::ConstMatrix< NekMatrix< DataType, InnerMatrixType >::NumberType >
virtual	~ConstMatrix ()
boost::call_traits< NekMatrix< DataType, InnerMatrixType >::NumberType >::value_type	operator() (unsigned int row, unsigned int column) const
unsigned int	GetStorageSize () const
MatrixStorage	GetType () const
MatrixStorage	GetStorageType () const
unsigned int	GetRows () const
unsigned int	GetTransposedRows (char transpose) const
unsigned int	GetColumns () const
unsigned int	GetTransposedColumns (char transpose) const
const unsigned int *	GetSize () const
void	Transpose ()
char	GetTransposeFlag () const

Static Public Member Functions
static ThisType	CreateWrapper (const ThisType &rhs)
static std::shared_ptr< ThisType >	CreateWrapper (const std::shared_ptr< ThisType > &rhs)
Static Public Member Functions inherited from Nektar::ConstMatrix< NekMatrix< DataType, InnerMatrixType >::NumberType >
static unsigned int	CalculateIndex (MatrixStorage type, unsigned int row, unsigned int col, unsigned int numRows, unsigned int numColumns, const char transpose='N', unsigned int numSubDiags=0, unsigned int numSuperDiags=0)
static unsigned int	GetRequiredStorageSize (MatrixStorage type, unsigned int rows, unsigned int columns, unsigned int subDiags=0, unsigned int superDiags=0)

Protected Member Functions
boost::call_traits< NumberType >::value_type	v_GetValue (unsigned int row, unsigned int column) const override
unsigned int	v_GetStorageSize () const override
void	v_Transpose () override
Protected Member Functions inherited from Nektar::ConstMatrix< NekMatrix< DataType, InnerMatrixType >::NumberType >
	ConstMatrix (unsigned int rows, unsigned int columns, MatrixStorage policy=eFULL)
	ConstMatrix (const ConstMatrix< NekMatrix< DataType, InnerMatrixType >::NumberType > &rhs)
ConstMatrix< NekMatrix< DataType, InnerMatrixType >::NumberType > &	operator= (const ConstMatrix< NekMatrix< DataType, InnerMatrixType >::NumberType > &rhs)
void	Resize (unsigned int rows, unsigned int columns)
	Resets the rows and columns in the array. This method does not update the data storage to match the new row and column counts. More...
void	SetTransposeFlag (char newValue)
char	GetRawTransposeFlag () const
virtual boost::call_traits< NekMatrix< DataType, InnerMatrixType >::NumberType >::value_type	v_GetValue (unsigned int row, unsigned int column) const=0
virtual unsigned int	v_GetStorageSize () const=0
virtual void	v_Transpose ()
virtual char	v_GetTransposeFlag () const

Private Member Functions
void	Initialize (const unsigned int *rowsPerBlock, const unsigned int *columnsPerBlock)

Static Private Member Functions
static unsigned int	GetNumberOfElementsInBlock (unsigned int block, unsigned int totalBlocks, const Array< OneD, unsigned int > &sizes)

Private Attributes
Array< OneD, std::shared_ptr< InnerType > >	m_data
std::shared_ptr< InnerType >	m_nullBlockPtr
Array< OneD, unsigned int >	m_rowSizes
Array< OneD, unsigned int >	m_columnSizes
unsigned int	m_storageSize
unsigned int	m_numberOfBlockRows
unsigned int	m_numberOfBlockColumns

Static Private Attributes
static NumberType	m_zeroElement

Detailed Description

template<typename DataType, typename InnerMatrixType>
class Nektar::NekMatrix< NekMatrix< DataType, InnerMatrixType >, BlockMatrixTag >

Definition at line 49 of file BlockMatrix.hpp.

Member Typedef Documentation

BaseType

template<typename DataType , typename InnerMatrixType >

typedef ConstMatrix<NumberType> Nektar::NekMatrix< NekMatrix< DataType, InnerMatrixType >, BlockMatrixTag >::BaseType

Definition at line 57 of file BlockMatrix.hpp.

const_iterator

template<typename DataType , typename InnerMatrixType >

typedef iterator_base<const ThisType> Nektar::NekMatrix< NekMatrix< DataType, InnerMatrixType >, BlockMatrixTag >::const_iterator

Definition at line 125 of file BlockMatrix.hpp.

ConstGetValueType

template<typename DataType , typename InnerMatrixType >

typedef InnerType::ConstGetValueType Nektar::NekMatrix< NekMatrix< DataType, InnerMatrixType >, BlockMatrixTag >::ConstGetValueType

Definition at line 63 of file BlockMatrix.hpp.

GetValueType

template<typename DataType , typename InnerMatrixType >

typedef InnerType::GetValueType Nektar::NekMatrix< NekMatrix< DataType, InnerMatrixType >, BlockMatrixTag >::GetValueType

Definition at line 62 of file BlockMatrix.hpp.

InnerType

template<typename DataType , typename InnerMatrixType >

typedef NekMatrix<DataType, InnerMatrixType> Nektar::NekMatrix< NekMatrix< DataType, InnerMatrixType >, BlockMatrixTag >::InnerType

Definition at line 54 of file BlockMatrix.hpp.

iterator

template<typename DataType , typename InnerMatrixType >

typedef iterator_base<ThisType> Nektar::NekMatrix< NekMatrix< DataType, InnerMatrixType >, BlockMatrixTag >::iterator

Definition at line 124 of file BlockMatrix.hpp.

NumberType

template<typename DataType , typename InnerMatrixType >

typedef InnerType::NumberType Nektar::NekMatrix< NekMatrix< DataType, InnerMatrixType >, BlockMatrixTag >::NumberType

Definition at line 56 of file BlockMatrix.hpp.

ThisType

template<typename DataType , typename InnerMatrixType >

typedef NekMatrix<InnerType, BlockMatrixTag> Nektar::NekMatrix< NekMatrix< DataType, InnerMatrixType >, BlockMatrixTag >::ThisType

Definition at line 55 of file BlockMatrix.hpp.

Constructor & Destructor Documentation

NekMatrix() [1/6]

template<typename DataType , typename InnerMatrixType >

Nektar::NekMatrix< NekMatrix< DataType, InnerMatrixType >, BlockMatrixTag >::NekMatrix ( MatrixStorage  type = eFULL )

explicit

Definition at line 42 of file BlockMatrix.cpp.

    : BaseType(0, 0, type), m_data(), m_rowSizes(), m_columnSizes(),
      m_storageSize(), m_numberOfBlockRows(0), m_numberOfBlockColumns(0)
{
}

NekMatrix() [2/6]

template<typename DataType , typename InnerMatrixType >

Nektar::NekMatrix< NekMatrix< DataType, InnerMatrixType >, BlockMatrixTag >::NekMatrix	(	unsigned int	numberOfBlockRows,
		unsigned int	numberOfBlockColumns,
		unsigned int	rowsPerBlock,
		unsigned int	columnsPerBlock,
		MatrixStorage	type = eFULL
	)

Definition at line 50 of file BlockMatrix.cpp.

    : BaseType(numberOfBlockRows * rowsPerBlock,
               numberOfBlockColumns * columnsPerBlock, type),
      m_data(), m_rowSizes(numberOfBlockRows),
      m_columnSizes(numberOfBlockColumns), m_storageSize(0),
      m_numberOfBlockRows(numberOfBlockRows),
      m_numberOfBlockColumns(numberOfBlockColumns)
{
    m_storageSize = GetRequiredStorageSize();
    m_data        = Array<OneD, std::shared_ptr<InnerType>>(
        m_storageSize, std::shared_ptr<InnerType>());
    for (unsigned int i = 1; i <= numberOfBlockRows; ++i)
    {
        m_rowSizes[i - 1] = i * rowsPerBlock - 1;
    }
 
    for (unsigned int i = 1; i <= numberOfBlockColumns; ++i)
    {
        m_columnSizes[i - 1] = i * columnsPerBlock - 1;
    }
}

NekMatrix() [3/6]

template<typename DataType , typename InnerMatrixType >

Nektar::NekMatrix< NekMatrix< DataType, InnerMatrixType >, BlockMatrixTag >::NekMatrix	(	unsigned int	numberOfBlockRows,
		unsigned int	numberOfBlockColumns,
		const unsigned int *	rowsPerBlock,
		const unsigned int *	columnsPerBlock,
		MatrixStorage	type = eFULL
	)

Definition at line 75 of file BlockMatrix.cpp.

    : BaseType(
          std::accumulate(rowsPerBlock, rowsPerBlock + numberOfBlockRows, 0),
          std::accumulate(columnsPerBlock,
                          columnsPerBlock + numberOfBlockColumns, 0),
          type),
      m_data(), m_rowSizes(numberOfBlockRows),
      m_columnSizes(numberOfBlockColumns), m_storageSize(0),
      m_numberOfBlockRows(numberOfBlockRows),
      m_numberOfBlockColumns(numberOfBlockColumns)
{
    m_storageSize = GetRequiredStorageSize();
    m_data        = Array<OneD, std::shared_ptr<InnerType>>(
        m_storageSize, std::shared_ptr<InnerType>());
    Initialize(rowsPerBlock, columnsPerBlock);
}

NekMatrix() [4/6]

template<typename DataType , typename InnerMatrixType >

Nektar::NekMatrix< NekMatrix< DataType, InnerMatrixType >, BlockMatrixTag >::NekMatrix	(	unsigned int	numberOfBlockRows,
		unsigned int	numberOfBlockColumns,
		const Array< OneD, const unsigned int > &	rowsPerBlock,
		const Array< OneD, const unsigned int > &	columnsPerBlock,
		MatrixStorage	type = eFULL
	)

Definition at line 96 of file BlockMatrix.cpp.

    : BaseType(std::accumulate(rowsPerBlock.data(),
                               rowsPerBlock.data() + numberOfBlockRows, 0),
               std::accumulate(columnsPerBlock.data(),
                               columnsPerBlock.data() + numberOfBlockColumns,
                               0),
               type),
      m_data(), m_rowSizes(numberOfBlockRows),
      m_columnSizes(numberOfBlockColumns), m_storageSize(0),
      m_numberOfBlockRows(numberOfBlockRows),
      m_numberOfBlockColumns(numberOfBlockColumns)
{
    m_storageSize = GetRequiredStorageSize();
    m_data        = Array<OneD, std::shared_ptr<InnerType>>(
        m_storageSize, std::shared_ptr<InnerType>());
    Initialize(rowsPerBlock.data(), columnsPerBlock.data());
}

NekMatrix() [5/6]

template<typename DataType , typename InnerMatrixType >

Nektar::NekMatrix< NekMatrix< DataType, InnerMatrixType >, BlockMatrixTag >::NekMatrix	(	const Array< OneD, const unsigned int > &	rowsPerBlock,
		const Array< OneD, const unsigned int > &	columnsPerBlock,
		MatrixStorage	type = eFULL
	)

Definition at line 118 of file BlockMatrix.cpp.

    : BaseType(
          std::accumulate(rowsPerBlock.begin(), rowsPerBlock.end(), 0),
          std::accumulate(columnsPerBlock.begin(), columnsPerBlock.end(), 0),
          type),
      m_data(), m_rowSizes(rowsPerBlock.size()),
      m_columnSizes(columnsPerBlock.size()), m_storageSize(0),
      m_numberOfBlockRows(rowsPerBlock.size()),
      m_numberOfBlockColumns(columnsPerBlock.size())
{
    m_storageSize = GetRequiredStorageSize();
    m_data        = Array<OneD, std::shared_ptr<InnerType>>(
        m_storageSize, std::shared_ptr<InnerType>());
    Initialize(rowsPerBlock.data(), columnsPerBlock.data());
}

NekMatrix() [6/6]

template<typename DataType , typename InnerMatrixType >

Nektar::NekMatrix< NekMatrix< DataType, InnerMatrixType >, BlockMatrixTag >::NekMatrix

(

const ThisType &

rhs

)

Member Function Documentation

begin() [1/2]

template<typename DataType , typename InnerMatrixType >

NekMatrix< NekMatrix< DataType, InnerMatrixType >, BlockMatrixTag >::iterator Nektar::NekMatrix< NekMatrix< DataType, InnerMatrixType >, BlockMatrixTag >::begin

Definition at line 417 of file BlockMatrix.cpp.

{
    return iterator(*this, 0, 0);
}

begin() [2/2]

template<typename DataType , typename InnerMatrixType >

NekMatrix< NekMatrix< DataType, InnerMatrixType >, BlockMatrixTag >::const_iterator Nektar::NekMatrix< NekMatrix< DataType, InnerMatrixType >, BlockMatrixTag >::begin

Definition at line 433 of file BlockMatrix.cpp.

{
    return const_iterator(*this, 0, 0);
}

CalculateBlockIndex()

template<typename DataType , typename InnerMatrixType >

unsigned int Nektar::NekMatrix< NekMatrix< DataType, InnerMatrixType >, BlockMatrixTag >::CalculateBlockIndex	(	unsigned int	row,
		unsigned int	column
	)		const

Definition at line 161 of file BlockMatrix.cpp.

{
    return BaseType::CalculateIndex(this->GetStorageType(), row, column,
                                    m_numberOfBlockRows, m_numberOfBlockColumns,
                                    this->GetTransposeFlag());
}

CreateWrapper() [1/2]

template<typename DataType , typename InnerMatrixType >

static std::shared_ptr< ThisType > Nektar::NekMatrix< NekMatrix< DataType, InnerMatrixType >, BlockMatrixTag >::CreateWrapper ( const std::shared_ptr< ThisType > &  rhs )

static

CreateWrapper() [2/2]

template<typename DataType , typename InnerMatrixType >

static ThisType Nektar::NekMatrix< NekMatrix< DataType, InnerMatrixType >, BlockMatrixTag >::CreateWrapper ( const ThisType &  rhs )

static

end() [1/2]

template<typename DataType , typename InnerMatrixType >

NekMatrix< NekMatrix< DataType, InnerMatrixType >, BlockMatrixTag >::iterator Nektar::NekMatrix< NekMatrix< DataType, InnerMatrixType >, BlockMatrixTag >::end

Definition at line 425 of file BlockMatrix.cpp.

{
    return iterator(*this);
}

end() [2/2]

template<typename DataType , typename InnerMatrixType >

NekMatrix< NekMatrix< DataType, InnerMatrixType >, BlockMatrixTag >::const_iterator Nektar::NekMatrix< NekMatrix< DataType, InnerMatrixType >, BlockMatrixTag >::end

Definition at line 441 of file BlockMatrix.cpp.

{
    return const_iterator(*this);
}

GetBlock() [1/2]

template<typename DataType , typename InnerMatrixType >

std::shared_ptr< typename NekMatrix< NekMatrix< DataType, InnerMatrixType >, BlockMatrixTag >::InnerType > & Nektar::NekMatrix< NekMatrix< DataType, InnerMatrixType >, BlockMatrixTag >::GetBlock	(	unsigned int	row,
		unsigned int	column
	)

Definition at line 231 of file BlockMatrix.cpp.

{
    ASSERTL2(
        this->GetTransposeFlag() == 'N' ? row < m_numberOfBlockRows
                                        : row < m_numberOfBlockColumns,
        std::string("Row ") + std::to_string(row) +
            std::string(" requested in a block matrix with a maximum of ") +
            std::to_string(m_numberOfBlockRows) + std::string(" rows"));
    ASSERTL2(
        this->GetTransposeFlag() == 'N' ? column < m_numberOfBlockColumns
                                        : column < m_numberOfBlockRows,
        std::string("Column ") + std::to_string(column) +
            std::string(" requested in a block matrix with a maximum of ") +
            std::to_string(m_numberOfBlockColumns) + std::string(" columns"));
    int x = CalculateBlockIndex(row, column);
    if (x == -1)
    {
        return m_nullBlockPtr;
    }
    else
    {
        return m_data[x];
    }
}

References ASSERTL2.

GetBlock() [2/2]

template<typename DataType , typename InnerMatrixType >

std::shared_ptr< const typename NekMatrix< NekMatrix< DataType, InnerMatrixType >, BlockMatrixTag >::InnerType > Nektar::NekMatrix< NekMatrix< DataType, InnerMatrixType >, BlockMatrixTag >::GetBlock	(	unsigned int	row,
		unsigned int	column
	)		const

Definition at line 202 of file BlockMatrix.cpp.

{
    ASSERTL2(
        this->GetTransposeFlag() == 'N' ? row < m_numberOfBlockRows
                                        : row < m_numberOfBlockColumns,
        std::string("Row ") + std::to_string(row) +
            std::string(" requested in a block matrix with a maximum of ") +
            std::to_string(m_numberOfBlockRows) + std::string(" rows"));
    ASSERTL2(
        this->GetTransposeFlag() == 'N' ? column < m_numberOfBlockColumns
                                        : column < m_numberOfBlockRows,
        std::string("Column ") + std::to_string(column) +
            std::string(" requested in a block matrix with a maximum of ") +
            std::to_string(m_numberOfBlockColumns) + std::string(" columns"));
    int x = CalculateBlockIndex(row, column);
    if (x < 0)
    {
        return std::shared_ptr<const InnerType>();
    }
    else
    {
        return m_data[x];
    }
}

References ASSERTL2.

GetBlockPtr()

template<typename DataType , typename InnerMatrixType >

const NekMatrix< NekMatrix< DataType, InnerMatrixType >, BlockMatrixTag >::InnerType * Nektar::NekMatrix< NekMatrix< DataType, InnerMatrixType >, BlockMatrixTag >::GetBlockPtr	(	unsigned int	row,
		unsigned int	column
	)		const

Definition at line 173 of file BlockMatrix.cpp.

{
    ASSERTL2(
        this->GetTransposeFlag() == 'N' ? row < m_numberOfBlockRows
                                        : row < m_numberOfBlockColumns,
        std::string("Row ") + std::to_string(row) +
            std::string(" requested in a block matrix with a maximum of ") +
            std::to_string(m_numberOfBlockRows) + std::string(" rows"));
    ASSERTL2(
        this->GetTransposeFlag() == 'N' ? column < m_numberOfBlockColumns
                                        : column < m_numberOfBlockColumns,
        std::string("Column ") + std::to_string(column) +
            std::string(" requested in a block matrix with a maximum of ") +
            std::to_string(m_numberOfBlockColumns) + std::string(" columns"));
    int x = CalculateBlockIndex(row, column);
    if (x == -1)
    {
        return nullptr;
    }
    else
    {
        return m_data[x].get();
    }
}

References ASSERTL2.

GetBlockSizes()

template<typename DataType , typename InnerMatrixType >

void Nektar::NekMatrix< NekMatrix< DataType, InnerMatrixType >, BlockMatrixTag >::GetBlockSizes	(	Array< OneD, unsigned int > &	rowSizes,
		Array< OneD, unsigned int > &	colSizes
	)		const

Definition at line 398 of file BlockMatrix.cpp.

{
    if (this->GetTransposeFlag() == 'T')
    {
        rowSizes = m_columnSizes;
        colSizes = m_rowSizes;
    }
    else
    {
        rowSizes = m_rowSizes;
        colSizes = m_columnSizes;
    }
}

GetNumberOfBlockColumns()

template<typename DataType , typename InnerMatrixType >

unsigned int Nektar::NekMatrix< NekMatrix< DataType, InnerMatrixType >, BlockMatrixTag >::GetNumberOfBlockColumns

Definition at line 353 of file BlockMatrix.cpp.

{
    if (this->GetTransposeFlag() == 'N')
    {
        return m_numberOfBlockColumns;
    }
    else
    {
        return m_numberOfBlockRows;
    }
}

GetNumberOfBlockRows()

template<typename DataType , typename InnerMatrixType >

unsigned int Nektar::NekMatrix< NekMatrix< DataType, InnerMatrixType >, BlockMatrixTag >::GetNumberOfBlockRows

Definition at line 339 of file BlockMatrix.cpp.

{
    if (this->GetTransposeFlag() == 'N')
    {
        return m_numberOfBlockRows;
    }
    else
    {
        return m_numberOfBlockColumns;
    }
}

GetNumberOfColumnsInBlockColumn()

template<typename DataType , typename InnerMatrixType >

unsigned int Nektar::NekMatrix< NekMatrix< DataType, InnerMatrixType >, BlockMatrixTag >::GetNumberOfColumnsInBlockColumn

(

unsigned int

blockCol

)

const

Definition at line 382 of file BlockMatrix.cpp.

{
    if (this->GetTransposeFlag() == 'T')
    {
        return GetNumberOfElementsInBlock(blockCol, m_numberOfBlockRows,
                                          m_rowSizes);
    }
    else
    {
        return GetNumberOfElementsInBlock(blockCol, m_numberOfBlockColumns,
                                          m_columnSizes);
    }
}

GetNumberOfElementsInBlock()

template<typename DataType , typename InnerMatrixType >

unsigned int Nektar::NekMatrix< NekMatrix< DataType, InnerMatrixType >, BlockMatrixTag >::GetNumberOfElementsInBlock ( unsigned int  block, unsigned int  totalBlocks, const Array< OneD, unsigned int > &  sizes  )

staticprivate

Definition at line 465 of file BlockMatrix.cpp.

{
    ASSERTL2(block < totalBlocks,
             std::string("Block Element ") + std::to_string(block) +
                 std::string(" requested in a matrix with a maximum of ") +
                 std::to_string(totalBlocks) + std::string(" blocks."));
    if (block == 0)
    {
        return sizes[block] + 1;
    }
    else
    {
        return sizes[block] - sizes[block - 1];
    }
}

References ASSERTL2.

GetNumberOfRowsInBlockRow()

template<typename DataType , typename InnerMatrixType >

unsigned int Nektar::NekMatrix< NekMatrix< DataType, InnerMatrixType >, BlockMatrixTag >::GetNumberOfRowsInBlockRow

(

unsigned int

blockRow

)

const

Definition at line 366 of file BlockMatrix.cpp.

{
    if (this->GetTransposeFlag() == 'N')
    {
        return GetNumberOfElementsInBlock(blockRow, m_numberOfBlockRows,
                                          m_rowSizes);
    }
    else
    {
        return GetNumberOfElementsInBlock(blockRow, m_numberOfBlockColumns,
                                          m_columnSizes);
    }
}

GetRequiredStorageSize()

template<typename DataType , typename InnerMatrixType >

unsigned int Nektar::NekMatrix< NekMatrix< DataType, InnerMatrixType >, BlockMatrixTag >::GetRequiredStorageSize

Definition at line 152 of file BlockMatrix.cpp.

{
    return BaseType::GetRequiredStorageSize(this->GetStorageType(),
                                            this->GetNumberOfBlockRows(),
                                            this->GetNumberOfBlockColumns());
}

GetStorageSize()

template<typename DataType , typename InnerMatrixType >

unsigned int Nektar::NekMatrix< NekMatrix< DataType, InnerMatrixType >, BlockMatrixTag >::GetStorageSize

Definition at line 332 of file BlockMatrix.cpp.

{
    return m_storageSize;
}

Initialize()

template<typename DataType , typename InnerMatrixType >

void Nektar::NekMatrix< NekMatrix< DataType, InnerMatrixType >, BlockMatrixTag >::Initialize ( const unsigned int *  rowsPerBlock, const unsigned int *  columnsPerBlock  )

private

Definition at line 486 of file BlockMatrix.cpp.

{
    m_storageSize = this->GetRows() * this->GetColumns();
    if (this->GetRows() > 0)
    {
        m_rowSizes[0] = rowsPerBlock[0] - 1;
        for (unsigned int i = 1; i < m_numberOfBlockRows; ++i)
        {
            m_rowSizes[i] = rowsPerBlock[i] + m_rowSizes[i - 1];
        }
    }
    if (this->GetColumns() > 0)
    {
        m_columnSizes[0] = columnsPerBlock[0] - 1;
        for (unsigned int i = 1; i < m_numberOfBlockColumns; ++i)
        {
            m_columnSizes[i] = columnsPerBlock[i] + m_columnSizes[i - 1];
        }
    }
}

operator()()

template<typename DataType , typename InnerMatrixType >

NekMatrix< NekMatrix< DataType, InnerMatrixType >, BlockMatrixTag >::ConstGetValueType Nektar::NekMatrix< NekMatrix< DataType, InnerMatrixType >, BlockMatrixTag >::operator()	(	unsigned int	row,
		unsigned int	col
	)		const

Definition at line 279 of file BlockMatrix.cpp.

{
    ASSERTL2(row < this->GetRows(),
             std::string("Row ") + std::to_string(row) +
                 std::string(" requested in a matrix with a maximum of ") +
                 std::to_string(this->GetRows()) + std::string(" rows"));
    ASSERTL2(col < this->GetColumns(),
             std::string("Column ") + std::to_string(col) +
                 std::string(" requested in a matrix with a maximum of ") +
                 std::to_string(this->GetColumns()) + std::string(" columns"));
 
    const Array<OneD, unsigned int> *rowSizes    = &m_rowSizes;
    const Array<OneD, unsigned int> *columnSizes = &m_columnSizes;
 
    if (this->GetTransposeFlag() == 'T')
    {
        std::swap(rowSizes, columnSizes);
    }
 
    unsigned int blockRow =
        std::lower_bound(rowSizes->begin(), rowSizes->end(), row) -
        rowSizes->begin();
    unsigned int blockColumn =
        std::lower_bound(columnSizes->begin(), columnSizes->end(), col) -
        columnSizes->begin();
    const std::shared_ptr<const InnerType> block =
        GetBlock(blockRow, blockColumn);
 
    unsigned int actualRow = row;
    if (blockRow > 0)
    {
        actualRow = row - ((*rowSizes)[blockRow - 1]) - 1;
    }
 
    unsigned int actualCol = col;
    if (blockColumn > 0)
    {
        actualCol = col - ((*columnSizes)[blockColumn - 1]) - 1;
    }
 
    if (block)
    {
        return (*block)(actualRow, actualCol);
    }
    else
    {
        return m_zeroElement;
    }
}

References ASSERTL2.

SetBlock()

template<typename DataType , typename InnerMatrixType >

void Nektar::NekMatrix< NekMatrix< DataType, InnerMatrixType >, BlockMatrixTag >::SetBlock	(	unsigned int	row,
		unsigned int	column,
		std::shared_ptr< InnerType > &	m
	)

Definition at line 258 of file BlockMatrix.cpp.

{
    ASSERTL2(
        this->GetTransposeFlag() == 'N' ? row < m_numberOfBlockRows
                                        : row < m_numberOfBlockColumns,
        std::string("Row ") + std::to_string(row) +
            std::string(" requested in a block matrix with a maximum of ") +
            std::to_string(m_numberOfBlockRows) + std::string(" rows"));
    ASSERTL2(
        this->GetTransposeFlag() == 'N' ? column < m_numberOfBlockColumns
                                        : column < m_numberOfBlockRows,
        std::string("Column ") + std::to_string(column) +
            std::string(" requested in a block matrix with a maximum of ") +
            std::to_string(m_numberOfBlockColumns) + std::string(" columns"));
    m_data[CalculateBlockIndex(row, column)] = InnerType::CreateWrapper(m);
}

References ASSERTL2.

v_GetStorageSize()

template<typename DataType , typename InnerMatrixType >

unsigned int Nektar::NekMatrix< NekMatrix< DataType, InnerMatrixType >, BlockMatrixTag >::v_GetStorageSize

overrideprotectedvirtual

Implements Nektar::ConstMatrix< NekMatrix< DataType, InnerMatrixType >::NumberType >.

Definition at line 520 of file BlockMatrix.cpp.

{
    return this->GetStorageSize();
}

v_GetValue()

template<typename DataType , typename InnerMatrixType >

boost::call_traits< typenameNekMatrix< NekMatrix< DataType, InnerMatrixType >, BlockMatrixTag >::NumberType >::value_type Nektar::NekMatrix< NekMatrix< DataType, InnerMatrixType >, BlockMatrixTag >::v_GetValue ( unsigned int  row, unsigned int  column  ) const

overrideprotectedvirtual

Implements Nektar::ConstMatrix< NekMatrix< DataType, InnerMatrixType >::NumberType >.

Definition at line 512 of file BlockMatrix.cpp.

{
    return (*this)(row, column);
}

v_Transpose()

template<typename DataType , typename InnerMatrixType >

void Nektar::NekMatrix< NekMatrix< DataType, InnerMatrixType >, BlockMatrixTag >::v_Transpose

overrideprotectedvirtual

Reimplemented from Nektar::ConstMatrix< NekMatrix< DataType, InnerMatrixType >::NumberType >.

Definition at line 527 of file BlockMatrix.cpp.

{
    for (auto &ptr : m_data)
    {
        if (ptr.get() != nullptr)
        {
            ptr->Transpose();
        }
    }
}

Member Data Documentation

m_columnSizes

template<typename DataType , typename InnerMatrixType >

Array<OneD, unsigned int> Nektar::NekMatrix< NekMatrix< DataType, InnerMatrixType >, BlockMatrixTag >::m_columnSizes

private

Definition at line 222 of file BlockMatrix.hpp.

m_data

template<typename DataType , typename InnerMatrixType >

Array<OneD, std::shared_ptr<InnerType> > Nektar::NekMatrix< NekMatrix< DataType, InnerMatrixType >, BlockMatrixTag >::m_data

private

Definition at line 219 of file BlockMatrix.hpp.

m_nullBlockPtr

template<typename DataType , typename InnerMatrixType >

std::shared_ptr<InnerType> Nektar::NekMatrix< NekMatrix< DataType, InnerMatrixType >, BlockMatrixTag >::m_nullBlockPtr

private

Definition at line 220 of file BlockMatrix.hpp.

m_numberOfBlockColumns

template<typename DataType , typename InnerMatrixType >

unsigned int Nektar::NekMatrix< NekMatrix< DataType, InnerMatrixType >, BlockMatrixTag >::m_numberOfBlockColumns

private

Definition at line 225 of file BlockMatrix.hpp.

m_numberOfBlockRows

template<typename DataType , typename InnerMatrixType >

unsigned int Nektar::NekMatrix< NekMatrix< DataType, InnerMatrixType >, BlockMatrixTag >::m_numberOfBlockRows

private

Definition at line 224 of file BlockMatrix.hpp.

m_rowSizes

template<typename DataType , typename InnerMatrixType >

Array<OneD, unsigned int> Nektar::NekMatrix< NekMatrix< DataType, InnerMatrixType >, BlockMatrixTag >::m_rowSizes

private

Definition at line 221 of file BlockMatrix.hpp.

m_storageSize

template<typename DataType , typename InnerMatrixType >

unsigned int Nektar::NekMatrix< NekMatrix< DataType, InnerMatrixType >, BlockMatrixTag >::m_storageSize

private

Definition at line 223 of file BlockMatrix.hpp.

m_zeroElement

template<typename DataType , typename InnerMatrixType >

NekMatrix< NekMatrix< DataType, InnerMatrixType >, BlockMatrixTag >::NumberType Nektar::NekMatrix< NekMatrix< DataType, InnerMatrixType >, BlockMatrixTag >::m_zeroElement

staticprivate

Definition at line 226 of file BlockMatrix.hpp.