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 207 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 226 of file MatrixFuncs.cpp.

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

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

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

CalculateIndex()

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

static

Definition at line 213 of file MatrixFuncs.cpp.

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