Nektar::LowerTriangularMatrixFuncs Struct Reference

#include <MatrixFuncs.h>

Inheritance diagram for Nektar::LowerTriangularMatrixFuncs:

Static Public Member Functions
static unsigned int	CalculateIndex (unsigned int totalColumns, unsigned int curRow, unsigned int curColumn)
static std::tuple< unsigned int, unsigned int >	Advance (const unsigned int totalRows, const unsigned int totalColumns, const unsigned int curRow, const unsigned int curColumn, char transpose='N')
Static Public Member Functions inherited from Nektar::TriangularMatrixFuncs
static unsigned int	GetRequiredStorageSize (unsigned int rows, unsigned int columns)

Detailed Description

Definition at line 179 of file MatrixFuncs.h.

Member Function Documentation

Advance()

std::tuple< unsigned int, unsigned int > Nektar::LowerTriangularMatrixFuncs::Advance ( const unsigned int  totalRows, const unsigned int  totalColumns, const unsigned int  curRow, const unsigned int  curColumn, char  transpose = 'N'  )

static

Definition at line 220 of file MatrixFuncs.cpp.

References Nektar::UpperTriangularMatrixFuncs::Advance(), and ASSERTL1.

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

    {
        ASSERTL1(totalRows == totalColumns, "Triangular matrices must be square.");
        ASSERTL1(curRow < totalRows, "Attemping to iterate through an element on row " +
            std::to_string(curRow) + " of a (" +
            std::to_string(totalRows) + ", " +
            std::to_string(totalColumns) + " lower triangular matrix.");
        ASSERTL1(curColumn < totalColumns, "Attemping to iterate through an element on row " +
            std::to_string(curColumn) + " of a (" +
            std::to_string(totalRows) + ", " +
            std::to_string(totalColumns) + " lower triangular matrix.");
        ASSERTL1(curRow >= curColumn, "Attemping to iterate through element (" +
            std::to_string(curRow) + ", " +
            std::to_string(curColumn) + ") of a (" +
            std::to_string(totalRows) + ", " +
            std::to_string(totalColumns) + " lower triangular matrix.");

        if( transpose == 'T' )
        {
            return UpperTriangularMatrixFuncs::Advance(
                totalColumns, totalRows, curColumn, curRow);
        }

        unsigned int nextRow = curRow;
        unsigned int nextColumn = curColumn;

        if( nextRow < totalRows )
        {
            ++nextRow;
        }

        if( nextRow >= totalRows )
        {
            ++nextColumn;
            nextRow = nextColumn;
        }

        if( nextColumn >= totalColumns )
        {
            nextRow = std::numeric_limits<unsigned int>::max();
            nextColumn = std::numeric_limits<unsigned int>::max();
        }

        return std::tuple<unsigned int, unsigned int>(nextRow, nextColumn);
    }

CalculateIndex()

unsigned int Nektar::LowerTriangularMatrixFuncs::CalculateIndex ( unsigned int  totalColumns, unsigned int  curRow, unsigned int  curColumn  )

static

Definition at line 207 of file MatrixFuncs.cpp.

Referenced by Nektar::ConstMatrix< NekMatrix< DataType, InnerMatrixType >::NumberType >::CalculateIndex().

    {
        if( curRow >= curColumn )
        {
            return curRow + (2*totalColumns - curColumn - 1)*(curColumn)/2;
        }
        else
        {
            return std::numeric_limits<unsigned int>::max();
        }
    }