Nektar::SymmetricMatrixFuncs Struct Reference

#include <MatrixFuncs.h>

Inheritance diagram for Nektar::SymmetricMatrixFuncs:

Collaboration diagram for Nektar::SymmetricMatrixFuncs:

Static Public Member Functions
static unsigned int	CalculateIndex (unsigned int curRow, unsigned int curColumn)
template<typename DataType >
static void	Invert (unsigned int rows, unsigned int columns, Array< OneD, DataType > &data)
static boost::tuples::tuple < unsigned int, unsigned int >	Advance (const unsigned int totalRows, const unsigned int totalColumns, const unsigned int curRow, const unsigned int curColumn)

Additional Inherited Members
Static Private Member Functions inherited from Nektar::TriangularMatrixFuncs
static unsigned int	GetRequiredStorageSize (unsigned int rows, unsigned int columns)

Detailed Description

Definition at line 199 of file MatrixFuncs.h.

Member Function Documentation

boost::tuples::tuple< unsigned int, unsigned int > Nektar::SymmetricMatrixFuncs::Advance ( const unsigned int  totalRows, const unsigned int  totalColumns, const unsigned int  curRow, const unsigned int  curColumn  )

static

Definition at line 296 of file MatrixFuncs.cpp.

References ASSERTL1.

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

    {
        ASSERTL1(totalRows == totalColumns, "Symmetric matrices must be square.");
        ASSERTL1(curRow < totalRows, "Attemping to iterate through an element on row " +
            boost::lexical_cast<std::string>(curRow) + " of a (" +
            boost::lexical_cast<std::string>(totalRows) + ", " +
            boost::lexical_cast<std::string>(totalColumns) + " symmetric matrix.");
        ASSERTL1(curColumn < totalColumns, "Attemping to iterate through an element on row " +
            boost::lexical_cast<std::string>(curColumn) + " of a (" +
            boost::lexical_cast<std::string>(totalRows) + ", " +
            boost::lexical_cast<std::string>(totalColumns) + " symmetric matrix.");

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

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

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

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

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

unsigned int Nektar::SymmetricMatrixFuncs::CalculateIndex ( unsigned int  curRow, unsigned int  curColumn  )

static

Definition at line 283 of file MatrixFuncs.cpp.

Referenced by Nektar::ConstMatrix< DataType >::CalculateIndex().

    {
        if( curRow <= curColumn )
        {
            return curRow + curColumn*(curColumn+1)/2;
        }
        else
        {
            return curColumn + curRow*(curRow + 1)/2;
        }
    }

template<typename DataType >

static void Nektar::SymmetricMatrixFuncs::Invert ( unsigned int  rows, unsigned int  columns, Array< OneD, DataType > &  data  )

inlinestatic

Definition at line 206 of file MatrixFuncs.h.

References ASSERTL0, and Nektar::Array< OneD, DataType >::get().

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

        {
#ifdef NEKTAR_USING_BLAS
            ASSERTL0(rows==columns, "Only square matrices can be inverted.");

            int n = columns;
            int info = 0;
            Array<OneD, int> ipivot(n);
            Array<OneD, DataType> work(n);

            Lapack::Dsptrf('U', n, data.get(), ipivot.get(), info);

            if( info < 0 )
            {
                std::string message = "ERROR: The " + boost::lexical_cast<std::string>(-info) + "th parameter had an illegal parameter for dsptrf";
                ASSERTL0(false, message.c_str());
            }
            else if( info > 0 )
            {
                std::string message = "ERROR: Element u_" + boost::lexical_cast<std::string>(info) +   boost::lexical_cast<std::string>(info) + " is 0 from dsptrf";
                ASSERTL0(false, message.c_str());
            }

            Lapack::Dsptri('U', n, data.get(), ipivot.get(),
                           work.get(), info);

            if( info < 0 )
            {
                std::string message = "ERROR: The " + boost::lexical_cast<std::string>(-info) + "th parameter had an illegal parameter for dsptri";
                ASSERTL0(false, message.c_str());
            }
            else if( info > 0 )
            {
                std::string message = "ERROR: Element u_" + boost::lexical_cast<std::string>(info) +   boost::lexical_cast<std::string>(info) + " is 0 from dsptri";
                ASSERTL0(false, message.c_str());
            }

#else
            // error Full matrix inversion not supported without blas.
            BOOST_STATIC_ASSERT(sizeof(DataType) == 0);
#endif
        }