Nektar++
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Pages
Static Public Member Functions | List of all members
Nektar::LinearSystemSolver Struct Reference

#include <NekLinSys.hpp>

Static Public Member Functions

template<typename BVectorType , typename XVectorType >
static void Solve (const BVectorType &b, XVectorType &x, MatrixStorage m_matrixType, const Array< OneD, const int > &m_ipivot, unsigned int n, const Array< OneD, const double > &A, char m_transposeFlag, unsigned int m_numberOfSubDiagonals, unsigned int m_numberOfSuperDiagonals)
template<typename BVectorType , typename XVectorType >
static void SolveTranspose (const BVectorType &b, XVectorType &x, MatrixStorage m_matrixType, const Array< OneD, const int > &m_ipivot, unsigned int n, const Array< OneD, const double > &A, char m_transposeFlag, unsigned int m_numberOfSubDiagonals, unsigned int m_numberOfSuperDiagonals)

Detailed Description

Definition at line 72 of file NekLinSys.hpp.

Member Function Documentation

template<typename BVectorType , typename XVectorType >
static void Nektar::LinearSystemSolver::Solve ( const BVectorType &  b,
XVectorType &  x,
MatrixStorage  m_matrixType,
const Array< OneD, const int > &  m_ipivot,
unsigned int  n,
const Array< OneD, const double > &  A,
char  m_transposeFlag,
unsigned int  m_numberOfSubDiagonals,
unsigned int  m_numberOfSuperDiagonals 
)
inlinestatic

Definition at line 75 of file NekLinSys.hpp.

References ASSERTL0, 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 SolveTranspose().

{
switch(m_matrixType)
{
case eFULL:
{
x = b;
int info = 0;
Lapack::Dgetrs('N',n,1,A.get(),n,(int *)m_ipivot.get(),x.GetRawPtr(),n,info);
if( info < 0 )
{
std::string message = "ERROR: The " + boost::lexical_cast<std::string>(-info) + "th parameter had an illegal parameter for dgetrs";
ASSERTL0(false, message.c_str());
}
}
break;
case eDIAGONAL:
for(unsigned int i = 0; i < A.num_elements(); ++i)
{
x[i] = b[i]*A[i];
}
break;
case eUPPER_TRIANGULAR:
{
x = b;
int info = 0;
Lapack::Dtptrs('U', m_transposeFlag, 'N', n, 1, A.get(), x.GetRawPtr(), n, info);
if( info < 0 )
{
std::string message = "ERROR: The " + boost::lexical_cast<std::string>(-info) + "th parameter had an illegal parameter for dtrtrs";
ASSERTL0(false, message.c_str());
}
else if( info > 0 )
{
std::string message = "ERROR: The " + boost::lexical_cast<std::string>(-info) + "th diagonal element of A is 0 for dtrtrs";
ASSERTL0(false, message.c_str());
}
}
break;
case eLOWER_TRIANGULAR:
{
x = b;
int info = 0;
Lapack::Dtptrs('L', m_transposeFlag, 'N', n, 1, A.get(), x.GetRawPtr(), n, info);
if( info < 0 )
{
std::string message = "ERROR: The " + boost::lexical_cast<std::string>(-info) + "th parameter had an illegal parameter for dtrtrs";
ASSERTL0(false, message.c_str());
}
else if( info > 0 )
{
std::string message = "ERROR: The " + boost::lexical_cast<std::string>(-info) + "th diagonal element of A is 0 for dtrtrs";
ASSERTL0(false, message.c_str());
}
}
break;
case eSYMMETRIC:
{
x = b;
int info = 0;
Lapack::Dsptrs('U', n, 1, A.get(), m_ipivot.get(), x.GetRawPtr(), x.GetRows(), info);
if( info < 0 )
{
std::string message = "ERROR: The " + boost::lexical_cast<std::string>(-info) + "th parameter had an illegal parameter for dsptrs";
ASSERTL0(false, message.c_str());
}
}
break;
case ePOSITIVE_DEFINITE_SYMMETRIC:
{
x = b;
int info = 0;
Lapack::Dpptrs('U', n, 1, A.get(), x.GetRawPtr(), x.GetRows(), info);
if( info < 0 )
{
std::string message = "ERROR: The " + boost::lexical_cast<std::string>(-info) + "th parameter had an illegal parameter for dpptrs";
ASSERTL0(false, message.c_str());
}
}
break;
case eBANDED:
{
x = b;
int KL = m_numberOfSubDiagonals;
int KU = m_numberOfSuperDiagonals;
int info = 0;
Lapack::Dgbtrs(m_transposeFlag, n, KL, KU, 1, A.get(), 2*KL+KU+1, m_ipivot.get(), x.GetRawPtr(), n, info);
if( info < 0 )
{
std::string message = "ERROR: The " + boost::lexical_cast<std::string>(-info) + "th parameter had an illegal parameter for dgbtrs";
ASSERTL0(false, message.c_str());
}
}
break;
case ePOSITIVE_DEFINITE_SYMMETRIC_BANDED:
{
x = b;
int KU = m_numberOfSuperDiagonals;
int info = 0;
Lapack::Dpbtrs('U', n, KU, 1, A.get(), KU+1, x.GetRawPtr(), n, info);
if( info < 0 )
{
std::string message = "ERROR: The " + boost::lexical_cast<std::string>(-info) + "th parameter had an illegal parameter for dpbtrs";
ASSERTL0(false, message.c_str());
}
}
break;
case eSYMMETRIC_BANDED:
NEKERROR(ErrorUtil::efatal, "Unhandled matrix type");
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");
}
}
template<typename BVectorType , typename XVectorType >
static void Nektar::LinearSystemSolver::SolveTranspose ( const BVectorType &  b,
XVectorType &  x,
MatrixStorage  m_matrixType,
const Array< OneD, const int > &  m_ipivot,
unsigned int  n,
const Array< OneD, const double > &  A,
char  m_transposeFlag,
unsigned int  m_numberOfSubDiagonals,
unsigned int  m_numberOfSuperDiagonals 
)
inlinestatic

Definition at line 210 of file NekLinSys.hpp.

References ASSERTL0, 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, NEKERROR, and Solve().

{
switch(m_matrixType)
{
case eFULL:
{
x = b;
int info = 0;
Lapack::Dgetrs('T',n,1,A.get(),n,(int *)m_ipivot.get(),x.GetRawPtr(), n,info);
if( info < 0 )
{
std::string message = "ERROR: The " + boost::lexical_cast<std::string>(-info) + "th parameter had an illegal parameter for dgetrs";
ASSERTL0(false, message.c_str());
}
}
break;
case eDIAGONAL:
Solve(b, x, m_matrixType, m_ipivot, n, A, m_transposeFlag, m_numberOfSubDiagonals, m_numberOfSuperDiagonals);
break;
case eUPPER_TRIANGULAR:
{
char trans = m_transposeFlag;
if( trans == 'N' )
{
trans = 'T';
}
else
{
trans = 'N';
}
x = b;
int info = 0;
Lapack::Dtptrs('U', trans, 'N', n, 1, A.get(), x.GetRawPtr(), n, info);
if( info < 0 )
{
std::string message = "ERROR: The " + boost::lexical_cast<std::string>(-info) + "th parameter had an illegal parameter for dtrtrs";
ASSERTL0(false, message.c_str());
}
else if( info > 0 )
{
std::string message = "ERROR: The " + boost::lexical_cast<std::string>(-info) + "th diagonal element of A is 0 for dtrtrs";
ASSERTL0(false, message.c_str());
}
}
break;
case eLOWER_TRIANGULAR:
{
char trans = m_transposeFlag;
if( trans == 'N' )
{
trans = 'T';
}
else
{
trans = 'N';
}
x = b;
int info = 0;
Lapack::Dtptrs('L', trans, 'N', n, 1, A.get(), x.GetRawPtr(), n, info);
if( info < 0 )
{
std::string message = "ERROR: The " + boost::lexical_cast<std::string>(-info) + "th parameter had an illegal parameter for dtrtrs";
ASSERTL0(false, message.c_str());
}
else if( info > 0 )
{
std::string message = "ERROR: The " + boost::lexical_cast<std::string>(-info) + "th diagonal element of A is 0 for dtrtrs";
ASSERTL0(false, message.c_str());
}
}
break;
case eSYMMETRIC:
case ePOSITIVE_DEFINITE_SYMMETRIC:
case ePOSITIVE_DEFINITE_SYMMETRIC_BANDED:
Solve(b, x, m_matrixType, m_ipivot, n, A, m_transposeFlag, m_numberOfSubDiagonals, m_numberOfSuperDiagonals);
break;
case eBANDED:
{
x = b;
int KL = m_numberOfSubDiagonals;
int KU = m_numberOfSuperDiagonals;
int info = 0;
Lapack::Dgbtrs(m_transposeFlag, n, KL, KU, 1, A.get(), 2*KL+KU+1, m_ipivot.get(), x.GetRawPtr(), n, info);
if( info < 0 )
{
std::string message = "ERROR: The " + boost::lexical_cast<std::string>(-info) + "th parameter had an illegal parameter for dgbtrs";
ASSERTL0(false, message.c_str());
}
}
break;
case eSYMMETRIC_BANDED:
NEKERROR(ErrorUtil::efatal, "Unhandled matrix type");
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");
}
}
Generated on Sun Mar 15 2015 00:11:48 for Nektar++ by   doxygen 1.8.1.2