Nektar::StdRegionsTests Namespace Reference

Functions
	BOOST_AUTO_TEST_CASE (TestLocCollapsedCoordsTet)

Function Documentation

BOOST_AUTO_TEST_CASE()

Nektar::StdRegionsTests::BOOST_AUTO_TEST_CASE

(

TestLocCollapsedCoordsTet

)

Definition at line 45 of file TestLocCollapsedCoords.cpp.

{
    unsigned int nq = 6;
    auto pkey0 =
        LibUtilities::PointsKey(nq, LibUtilities::eGaussLobattoLegendre);
    auto pkey1 =
        LibUtilities::PointsKey(nq - 1, LibUtilities::eGaussRadauMAlpha1Beta0);
    auto pkey2 =
        LibUtilities::PointsKey(nq - 1, LibUtilities::eGaussRadauMAlpha2Beta0);
 
    unsigned int nm = 4;
    auto bkey0 = LibUtilities::BasisKey(LibUtilities::eModified_A, nm, pkey0);
    auto bkey1 = LibUtilities::BasisKey(LibUtilities::eModified_B, nm, pkey1);
    auto bkey2 = LibUtilities::BasisKey(LibUtilities::eModified_C, nm, pkey2);
 
    // Create tet expansion
    auto stdTet = StdRegions::StdTetExp(bkey0, bkey1, bkey2);
 
    // Get some evenly distributed points
    auto pkey =
        LibUtilities::PointsKey(nq, LibUtilities::eNodalTetEvenlySpaced);
    auto evenlySpaced = LibUtilities::PointsManager()[pkey];
 
    int np = evenlySpaced->GetTotNumPoints();
    Array<OneD, NekDouble> xc(np), yc(np), zc(np);
    evenlySpaced->GetPoints(xc, yc, zc);
 
    double epsilon = 1.0e-8;
 
    // Loop over points and ensure that the xi -> eta -> xi mapping returns the
    // original coordinates in the tetrahedron.
    for (auto i = 0; i < np; ++i)
    {
        Array<OneD, NekDouble> xi(3), eta(3), tmp(3);
        xi[0] = xc[i];
        xi[1] = yc[i];
        xi[2] = zc[i];
 
        stdTet.LocCoordToLocCollapsed(xi, eta);
        stdTet.LocCollapsedToLocCoord(eta, tmp);
 
        BOOST_CHECK_CLOSE(tmp[0], xi[0], epsilon);
        BOOST_CHECK_CLOSE(tmp[1], xi[1], epsilon);
        BOOST_CHECK_CLOSE(tmp[2], xi[2], epsilon);
    }
}

References Nektar::LibUtilities::eGaussLobattoLegendre, Nektar::LibUtilities::eModified_A, Nektar::LibUtilities::eModified_B, Nektar::LibUtilities::eModified_C, Nektar::LibUtilities::eNodalTetEvenlySpaced, and Nektar::LibUtilities::PointsManager().