#include <MatrixBase.hpp>

Inheritance diagram for Nektar::ConstMatrix< DataType >:

Public Member Functions
virtual	~ConstMatrix ()

boost::call_traits< DataType > ::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 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
	ConstMatrix (unsigned int rows, unsigned int columns, MatrixStorage policy=eFULL)

	ConstMatrix (const ConstMatrix< DataType > &rhs)

ConstMatrix< DataType > &	operator= (const ConstMatrix< DataType > &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

Private Member Functions
virtual boost::call_traits < DataType >::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 Attributes
unsigned int	m_size [2]

char	m_transpose

MatrixStorage	m_storageType

Detailed Description

template<typename DataType>
singleton Nektar::ConstMatrix< DataType >

Definition at line 57 of file MatrixBase.hpp.

Constructor & Destructor Documentation

template<typename DataType >

Nektar::ConstMatrix< DataType >::~ConstMatrix ( )

virtual

Definition at line 38 of file MatrixBase.cpp.

38 {}

template<typename DataType >

Nektar::ConstMatrix< DataType >::ConstMatrix	(	unsigned int	rows,
		unsigned int	columns,
		MatrixStorage	policy = `eFULL`
	)

protected

Definition at line 226 of file MatrixBase.cpp.

References Nektar::ConstMatrix< DataType >::m_size.

                                                              :
         m_size(),
         m_transpose('N'),
         m_storageType(policy)
     {
         m_size[0] = rows;
         m_size[1] = columns;
     }

template<typename DataType>

Nektar::ConstMatrix< DataType >::ConstMatrix ( const ConstMatrix< DataType > & rhs )

protected

Definition at line 237 of file MatrixBase.cpp.

References Nektar::ConstMatrix< DataType >::m_size.

                                                                        :
         m_size(),
         m_transpose(rhs.m_transpose),
         m_storageType(rhs.m_storageType)
     {
         m_size[0] = rhs.m_size[0];
         m_size[1] = rhs.m_size[1];
     }

Member Function Documentation

template<typename DataType >

unsigned int Nektar::ConstMatrix< DataType >::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

Definition at line 120 of file MatrixBase.cpp.

References ASSERTL1, Nektar::BandedMatrixFuncs::CalculateIndex(), Nektar::FullMatrixFuncs::CalculateIndex(), Nektar::UpperTriangularMatrixFuncs::CalculateIndex(), Nektar::LowerTriangularMatrixFuncs::CalculateIndex(), Nektar::SymmetricMatrixFuncs::CalculateIndex(), Nektar::DiagonalMatrixFuncs::CalculateIndex(), Nektar::SymmetricBandedMatrixFuncs::CalculateIndex(), Nektar::eBANDED, Nektar::eDIAGONAL, ErrorUtil::efatal, Nektar::eFULL, Nektar::eLOWER_TRIANGULAR, Nektar::eLOWER_TRIANGULAR_BANDED, Nektar::ePOSITIVE_DEFINITE_SYMMETRIC, Nektar::ePOSITIVE_DEFINITE_SYMMETRIC_BANDED, Nektar::eSYMMETRIC, Nektar::eSYMMETRIC_BANDED, Nektar::eUPPER_TRIANGULAR, Nektar::eUPPER_TRIANGULAR_BANDED, and NEKERROR.

Referenced by Nektar::NekMatrix< DataType, StandardMatrixTag >::CalculateIndex().

     {
         if(transpose == 'T' )
         {
             std::swap(row, col);
         }
         switch(type)
         {
             case eFULL:
                 return FullMatrixFuncs::CalculateIndex(numRows, numColumns, row, col);
                 break;
             case eDIAGONAL:
                 return DiagonalMatrixFuncs::CalculateIndex(row, col);
                 break;
             case eUPPER_TRIANGULAR:
                 return UpperTriangularMatrixFuncs::CalculateIndex(row, col);
                 break;
             case eLOWER_TRIANGULAR:
                 return LowerTriangularMatrixFuncs::CalculateIndex(numColumns, row, col);
                 break;
             case eSYMMETRIC:
             case ePOSITIVE_DEFINITE_SYMMETRIC:
                 return SymmetricMatrixFuncs::CalculateIndex(row, col);
                 break;
             case eBANDED:
                 return BandedMatrixFuncs::CalculateIndex(numRows, numColumns,
                     row, col, numSubDiags, numSuperDiags);
                 break;
             case eSYMMETRIC_BANDED:
             case ePOSITIVE_DEFINITE_SYMMETRIC_BANDED:
                 {
                     ASSERTL1(numSubDiags==numSuperDiags,
                              std::string("Number of sub- and superdiagonals should ") +
                              std::string("be equal for a symmetric banded matrix"));
                     return SymmetricBandedMatrixFuncs::CalculateIndex(row, col,
                         numSuperDiags);
                 }
                 break;
             case eUPPER_TRIANGULAR_BANDED:
                 NEKERROR(ErrorUtil::efatal, "Not yet implemented.");
                 break;
             case eLOWER_TRIANGULAR_BANDED:
                 NEKERROR(ErrorUtil::efatal, "Not yet implemented.");
                 break;
 
             default:
                 NEKERROR(ErrorUtil::efatal, "Unhandled matrix type");
         }
 
         return std::numeric_limits<unsigned int>::max();
     }

template<typename DataType >

unsigned int Nektar::ConstMatrix< DataType >::GetColumns ( ) const

Definition at line 78 of file MatrixBase.cpp.

Referenced by Nektar::AddEqualNegatedLhs(), Nektar::MultiplyEqual(), Nektar::NegateInPlace(), Nektar::NekMatrix< DataType, StandardMatrixTag >::operator==(), Nektar::SubtractEqualNegatedLhs(), and Nektar::Transpose().

     {
         return GetTransposedColumns(GetTransposeFlag());
     }

template<typename DataType>

char Nektar::ConstMatrix< DataType >::GetRawTransposeFlag ( ) const

inlineprotected

Definition at line 115 of file MatrixBase.hpp.

             {
                 return m_transpose;
             }

template<typename DataType >

unsigned int Nektar::ConstMatrix< DataType >::GetRequiredStorageSize	(	MatrixStorage	type,
		unsigned int	rows,
		unsigned int	columns,
		unsigned int	subDiags = `0`,
		unsigned int	superDiags = `0`
	)

static

Definition at line 176 of file MatrixBase.cpp.

References ASSERTL1, Nektar::eBANDED, Nektar::eDIAGONAL, ErrorUtil::efatal, Nektar::eFULL, Nektar::eLOWER_TRIANGULAR, Nektar::eLOWER_TRIANGULAR_BANDED, Nektar::ePOSITIVE_DEFINITE_SYMMETRIC, Nektar::ePOSITIVE_DEFINITE_SYMMETRIC_BANDED, Nektar::eSYMMETRIC, Nektar::eSYMMETRIC_BANDED, Nektar::eUPPER_TRIANGULAR, Nektar::eUPPER_TRIANGULAR_BANDED, Nektar::BandedMatrixFuncs::GetRequiredStorageSize(), Nektar::FullMatrixFuncs::GetRequiredStorageSize(), Nektar::TriangularMatrixFuncs::GetRequiredStorageSize(), Nektar::DiagonalMatrixFuncs::GetRequiredStorageSize(), Nektar::TriangularBandedMatrixFuncs::GetRequiredStorageSize(), and NEKERROR.

Referenced by Nektar::NekMatrix< DataType, StandardMatrixTag >::GetRequiredStorageSize().

     {
         switch(type)
         {
             case eFULL:
                 return FullMatrixFuncs::GetRequiredStorageSize(rows, columns);
                 break;
             case eDIAGONAL:
                 return DiagonalMatrixFuncs::GetRequiredStorageSize(rows, columns);
                 break;
             case eUPPER_TRIANGULAR:
                 return UpperTriangularMatrixFuncs::GetRequiredStorageSize(rows, columns);
                 break;
             case eLOWER_TRIANGULAR:
                 return LowerTriangularMatrixFuncs::GetRequiredStorageSize(rows, columns);
                 break;
             case eSYMMETRIC:
             case ePOSITIVE_DEFINITE_SYMMETRIC:
                 return SymmetricMatrixFuncs::GetRequiredStorageSize(rows, columns);
                 break;
             case eBANDED:
                 return BandedMatrixFuncs::GetRequiredStorageSize(rows, columns,
                     subDiags, superDiags);
                 break;
             case eSYMMETRIC_BANDED:
             case ePOSITIVE_DEFINITE_SYMMETRIC_BANDED:
                 {
                     ASSERTL1(subDiags==superDiags,
                              std::string("Number of sub- and superdiagonals should ") +
                              std::string("be equal for a symmetric banded matrix"));
                     return SymmetricBandedMatrixFuncs::GetRequiredStorageSize(rows, columns,
                         superDiags);
                 }
                 break;
             case eUPPER_TRIANGULAR_BANDED:
                 NEKERROR(ErrorUtil::efatal, "Unhandled matrix type");
                 break;
             case eLOWER_TRIANGULAR_BANDED:
                 NEKERROR(ErrorUtil::efatal, "Unhandled matrix type");
                 break;
 
             default:
                 NEKERROR(ErrorUtil::efatal, "Unhandled matrix type");
         }
 
         return 0;
     }

template<typename DataType >

unsigned int Nektar::ConstMatrix< DataType >::GetRows ( ) const

Definition at line 59 of file MatrixBase.cpp.

Referenced by Nektar::AddEqualNegatedLhs(), Nektar::MultiplyEqual(), Nektar::NegateInPlace(), Nektar::NekMatrix< DataType, StandardMatrixTag >::operator==(), Nektar::SubtractEqualNegatedLhs(), and Nektar::Transpose().

     {
         return GetTransposedRows(GetTransposeFlag());
     }

template<typename DataType >

const unsigned int * Nektar::ConstMatrix< DataType >::GetSize ( ) const

Definition at line 97 of file MatrixBase.cpp.

97 { return m_size; }

Nektar::ConstMatrix::m_size

unsigned int m_size[2]

Definition: MatrixBase.hpp:126

template<typename DataType >

unsigned int Nektar::ConstMatrix< DataType >::GetStorageSize ( ) const

Definition at line 53 of file MatrixBase.cpp.

     {
         return v_GetStorageSize();
     }

template<typename DataType>

MatrixStorage Nektar::ConstMatrix< DataType >::GetStorageType ( ) const

inline

Definition at line 72 of file MatrixBase.hpp.

             {
                 return m_storageType;
             }

template<typename DataType >

unsigned int Nektar::ConstMatrix< DataType >::GetTransposedColumns ( char transpose ) const

Definition at line 84 of file MatrixBase.cpp.

     {
         if( transpose == 'N' )
         {
             return m_size[1];
         }
         else
         {
             return m_size[0];
         }
     }

template<typename DataType >

unsigned int Nektar::ConstMatrix< DataType >::GetTransposedRows ( char transpose ) const

Definition at line 65 of file MatrixBase.cpp.

     {
         if( transpose == 'N' )
         {
             return m_size[0];
         }
         else
         {
             return m_size[1];
         }
     }

template<typename DataType >

char Nektar::ConstMatrix< DataType >::GetTransposeFlag ( ) const

Definition at line 114 of file MatrixBase.cpp.

Referenced by Nektar::NekMatrix< DataType, StandardMatrixTag >::operator==().

     {
         return v_GetTransposeFlag();
     }

template<typename DataType>

MatrixStorage Nektar::ConstMatrix< DataType >::GetType ( ) const

inline

Definition at line 67 of file MatrixBase.hpp.

Referenced by Nektar::Transpose().

             {
                 return m_storageType;
             }

template<typename DataType >

boost::call_traits< DataType >::value_type Nektar::ConstMatrix< DataType >::operator()	(	unsigned int	row,
		unsigned int	column
	)		const

Definition at line 41 of file MatrixBase.cpp.

References ASSERTL2.

     {
         ASSERTL2(row < GetRows(), std::string("Row ") + boost::lexical_cast<std::string>(row) +
             std::string(" requested in a matrix with a maximum of ") + boost::lexical_cast<std::string>(GetRows()) +
             std::string(" rows"));
         ASSERTL2(column < GetColumns(), std::string("Column ") + boost::lexical_cast<std::string>(column) +
             std::string(" requested in a matrix with a maximum of ") + boost::lexical_cast<std::string>(GetColumns()) +
             std::string(" columns"));
         return v_GetValue(row, column);
     }

template<typename DataType>

ConstMatrix< DataType > & Nektar::ConstMatrix< DataType >::operator= ( const ConstMatrix< DataType > & rhs )

protected

Definition at line 247 of file MatrixBase.cpp.

References Nektar::ConstMatrix< DataType >::m_size, and Nektar::ConstMatrix< DataType >::m_transpose.

Referenced by Nektar::Matrix< DataType >::operator=().

     {
         m_size[0] = rhs.m_size[0];
         m_size[1] = rhs.m_size[1];
         m_transpose = rhs.m_transpose;
         m_storageType = m_storageType;
         return *this;
     }

template<typename DataType >

void Nektar::ConstMatrix< DataType >::Resize	(	unsigned int	rows,
		unsigned int	columns
	)

protected

Resets the rows and columns in the array. This method does not update the data storage to match the new row and column counts.

Definition at line 258 of file MatrixBase.cpp.

     {
         m_size[0] = rows;
         m_size[1] = columns;
     }

template<typename DataType >

void Nektar::ConstMatrix< DataType >::SetTransposeFlag ( char newValue )

protected

Definition at line 265 of file MatrixBase.cpp.

     {
         m_transpose = newValue;
     }

template<typename DataType >

void Nektar::ConstMatrix< DataType >::Transpose ( )

Definition at line 100 of file MatrixBase.cpp.

Referenced by Nektar::Transpose().

     {
         if( m_transpose == 'N' )
         {
             m_transpose = 'T';
         }
         else
         {
             m_transpose = 'N';
         }
         v_Transpose();
     }