Nektar::StdRegionsUnitTests Namespace Reference

Functions
void	CheckInvBwdTrans (const StdRegions::StdExpansionSharedPtr &exp)
	BOOST_AUTO_TEST_CASE (TestQuadInvBwdTransPseudoInverse)
	BOOST_AUTO_TEST_CASE (TestHexInvBwdTransPseudoInverse)

Function Documentation

BOOST_AUTO_TEST_CASE() [1/2]

Nektar::StdRegionsUnitTests::BOOST_AUTO_TEST_CASE

(

TestHexInvBwdTransPseudoInverse

)

Definition at line 99 of file TestInvBwdTrans.cpp.

{
    using namespace LibUtilities;
 
    const PointsKey pointsKey(5, eGaussLobattoLegendre);
    const BasisKey basisKey(eModified_A, 3, pointsKey);
 
    StdRegions::StdExpansionSharedPtr exp =
        MemoryManager<StdRegions::StdHexExp>::AllocateSharedPtr(
            basisKey, basisKey, basisKey);
    CheckInvBwdTrans(exp);
}

References Nektar::MemoryManager< DataType >::AllocateSharedPtr(), CheckInvBwdTrans(), Nektar::LibUtilities::eGaussLobattoLegendre, and Nektar::LibUtilities::eModified_A.

BOOST_AUTO_TEST_CASE() [2/2]

Nektar::StdRegionsUnitTests::BOOST_AUTO_TEST_CASE

(

TestQuadInvBwdTransPseudoInverse

)

Definition at line 86 of file TestInvBwdTrans.cpp.

{
    using namespace LibUtilities;
 
    const PointsKey pointsKey(6, eGaussLobattoLegendre);
    const BasisKey basisKey(eModified_A, 4, pointsKey);
 
    StdRegions::StdExpansionSharedPtr exp =
        MemoryManager<StdRegions::StdQuadExp>::AllocateSharedPtr(basisKey,
                                                                 basisKey);
    CheckInvBwdTrans(exp);
}

References Nektar::MemoryManager< DataType >::AllocateSharedPtr(), CheckInvBwdTrans(), Nektar::LibUtilities::eGaussLobattoLegendre, and Nektar::LibUtilities::eModified_A.

CheckInvBwdTrans()

void Nektar::StdRegionsUnitTests::CheckInvBwdTrans

(

const StdRegions::StdExpansionSharedPtr &

exp

)

Definition at line 44 of file TestInvBwdTrans.cpp.

{
    const StdRegions::StdMatrixKey bwdKey(StdRegions::eBwdTrans,
                                          exp->DetShapeType(), *exp);
    const StdRegions::StdMatrixKey invKey(StdRegions::eInvBwdTrans,
                                          exp->DetShapeType(), *exp);
    const StdRegions::StdMatrixKey invIProdKey(StdRegions::eInvIProductWRTBase,
                                               exp->DetShapeType(), *exp);
 
    DNekMatSharedPtr bwd      = exp->GetStdMatrix(bwdKey);
    DNekMatSharedPtr invBwd   = exp->GetStdMatrix(invKey);
    DNekMatSharedPtr invIProd = exp->GetStdMatrix(invIProdKey);
 
    BOOST_REQUIRE_EQUAL(bwd->GetRows(), exp->GetTotPoints());
    BOOST_REQUIRE_EQUAL(bwd->GetColumns(), exp->GetNcoeffs());
    BOOST_REQUIRE_EQUAL(invBwd->GetRows(), exp->GetNcoeffs());
    BOOST_REQUIRE_EQUAL(invBwd->GetColumns(), exp->GetTotPoints());
    BOOST_REQUIRE_EQUAL(invIProd->GetRows(), exp->GetTotPoints());
    BOOST_REQUIRE_EQUAL(invIProd->GetColumns(), exp->GetNcoeffs());
 
    const NekDouble epsilon = 1.0e-8;
    DNekMat identity        = (*invBwd) * (*bwd);
 
    for (unsigned int i = 0; i < identity.GetRows(); ++i)
    {
        for (unsigned int j = 0; j < identity.GetColumns(); ++j)
        {
            const NekDouble expected = i == j ? 1.0 : 0.0;
            BOOST_CHECK_SMALL(std::abs(*identity(i, j) - expected), epsilon);
        }
    }
 
    for (unsigned int i = 0; i < invBwd->GetRows(); ++i)
    {
        for (unsigned int j = 0; j < invBwd->GetColumns(); ++j)
        {
            BOOST_CHECK_SMALL(std::abs(*(*invIProd)(j, i) - *(*invBwd)(i, j)),
                              epsilon);
        }
    }
}

References Nektar::StdRegions::eBwdTrans, Nektar::StdRegions::eInvBwdTrans, and Nektar::StdRegions::eInvIProductWRTBase.

Referenced by BOOST_AUTO_TEST_CASE(), and BOOST_AUTO_TEST_CASE().