Nektar::MovementTests Namespace Reference

Functions
LibUtilities::SessionReaderSharedPtr	CreateSession ()
SpatialDomains::ZoneBaseShPtr	CreateZone (SpatialDomains::MovementType type, int zoneID, int domainID, LibUtilities::InterpreterSharedPtr interpreter)
	Produce dummy Zone objects, containing empty domain pointers. More...
SpatialDomains::InterfaceShPtr	CreateInterface (int interfaceID, std::vector< int > compositeIDs)
	Produce dummy Interface objects, containing empty domain pointers. More...
	BOOST_AUTO_TEST_CASE (TestAddGetZones)
	BOOST_AUTO_TEST_CASE (TestAddGetInterfaces)
	BOOST_AUTO_TEST_CASE (TestWriteMovement)

Variables
const std::string	angVelStr = "0.1*t"
const std::string	xEqnStr = "0.1*x - 0.1*t"
const std::string	yEqnStr = "0.1*y^2 - x"
const std::string	zEqnStr = "sqrt(t)"
const NekPoint< NekDouble >	origin = {1., 2., 3.}
const DNekVec	axis = {1., 2., 3.}
const std::vector< NekDouble >	velocity = {1., 1., 2.}

Function Documentation

BOOST_AUTO_TEST_CASE() [1/3]

Nektar::MovementTests::BOOST_AUTO_TEST_CASE

(

TestAddGetInterfaces

)

Definition at line 149 of file TestMovement.cpp.

{
    SpatialDomains::Movement m;
    SpatialDomains::InterfaceShPtr interface1 = CreateInterface(0, {0}),
                                   interface2 = CreateInterface(1, {1, 2, 3}),
                                   interface3 = CreateInterface(3, {4});
    m.AddInterface("north", interface1, interface2);
    m.AddInterface("east", interface2, interface3);
    m.AddInterface("south", interface3, interface1);
    SpatialDomains::InterfaceCollection interfaces = m.GetInterfaces();
    BOOST_TEST(interfaces.size() == 3);
    SpatialDomains::InterfacePairShPtr north = interfaces.at(
                                           std::make_pair(0, "north")),
                                       east = interfaces.at(
                                           std::make_pair(1, "east")),
                                       south = interfaces.at(
                                           std::make_pair(2, "south"));
    BOOST_TEST(north->m_leftInterface == interface1);
    BOOST_TEST(north->m_rightInterface == interface2);
    BOOST_TEST(east->m_leftInterface == interface2);
    BOOST_TEST(east->m_rightInterface == interface3);
    BOOST_TEST(south->m_leftInterface == interface3);
    BOOST_TEST(south->m_rightInterface == interface1);
}

References Nektar::SpatialDomains::Movement::AddInterface(), CreateInterface(), and Nektar::SpatialDomains::Movement::GetInterfaces().

BOOST_AUTO_TEST_CASE() [2/3]

Nektar::MovementTests::BOOST_AUTO_TEST_CASE

(

TestAddGetZones

)

Definition at line 128 of file TestMovement.cpp.

{
    LibUtilities::InterpreterSharedPtr interpreter =
        MemoryManager<LibUtilities::Interpreter>::AllocateSharedPtr();
    SpatialDomains::Movement m;
    SpatialDomains::ZoneBaseShPtr zone1 = CreateZone(
                                      SpatialDomains::MovementType::eFixed, 0,
                                      0, interpreter),
                                  zone2 = CreateZone(
                                      SpatialDomains::MovementType::eRotate, 2,
                                      2, interpreter);
    m.AddZone(zone1);
    m.AddZone(zone2);
    std::map<int, SpatialDomains::ZoneBaseShPtr> zones = m.GetZones();
    BOOST_TEST(zones.size() == 2);
    BOOST_TEST(zones.at(0).get() == zone1.get());
    BOOST_TEST(zones.at(2).get() == zone2.get());
    BOOST_TEST(zones.at(0) == zone1);
    BOOST_TEST(zones.at(2) == zone2);
}

References Nektar::SpatialDomains::Movement::AddZone(), Nektar::MemoryManager< DataType >::AllocateSharedPtr(), CreateZone(), Nektar::SpatialDomains::eFixed, Nektar::SpatialDomains::eRotate, and Nektar::SpatialDomains::Movement::GetZones().

BOOST_AUTO_TEST_CASE() [3/3]

Nektar::MovementTests::BOOST_AUTO_TEST_CASE

(

TestWriteMovement

)

Definition at line 174 of file TestMovement.cpp.

{
    LibUtilities::InterpreterSharedPtr interpreter =
        MemoryManager<LibUtilities::Interpreter>::AllocateSharedPtr();
    SpatialDomains::Movement m;
    m.AddZone(
        CreateZone(SpatialDomains::MovementType::eFixed, 0, 10, interpreter));
    m.AddZone(CreateZone(SpatialDomains::MovementType::ePrescribe, 2, 12,
                         interpreter));
    m.AddZone(CreateZone(SpatialDomains::MovementType::eTranslate, 3, 13,
                         interpreter));
    m.AddZone(
        CreateZone(SpatialDomains::MovementType::eRotate, 1, 11, interpreter));
    m.AddInterface("north", CreateInterface(0, {0}),
                   CreateInterface(1, {1, 2, 3, 4}));
    m.AddInterface("south", CreateInterface(2, {5, 6, 7, 8, 9, 11}),
                   CreateInterface(3, {12}));
 
    SpatialDomains::InterfaceCollection interfaces     = m.GetInterfaces();
    std::map<int, SpatialDomains::ZoneBaseShPtr> zones = m.GetZones();
    TiXmlElement *nektar = new TiXmlElement("NEKTAR");
    m.WriteMovement(nektar);
    TiXmlElement *movement = nektar->FirstChildElement("MOVEMENT");
    BOOST_TEST(movement);
    TiXmlElement *xmlZones      = movement->FirstChildElement("ZONES"),
                 *xmlInterfaces = movement->FirstChildElement("INTERFACES");
    BOOST_TEST(xmlZones);
    BOOST_TEST(xmlInterfaces);
 
    int id, err;
    std::string attr;
    std::vector<NekDouble> vec;
 
    // Check the fixed zone
    TiXmlElement *zone = xmlZones->FirstChildElement();
    BOOST_TEST(zone->Value() == "F");
    err = zone->QueryIntAttribute("ID", &id);
    BOOST_TEST(err == TIXML_SUCCESS);
    BOOST_TEST(id == 0);
    err = zone->QueryStringAttribute("DOMAIN", &attr);
    BOOST_TEST(err == TIXML_SUCCESS);
    BOOST_TEST(attr == "D[10]");
 
    // Check the rotating zone
    zone = xmlZones->IterateChildren(zone)->ToElement();
    BOOST_TEST(zone->Value() == "R");
    err = zone->QueryIntAttribute("ID", &id);
    BOOST_TEST(err == TIXML_SUCCESS);
    BOOST_TEST(id == 1);
    err = zone->QueryStringAttribute("DOMAIN", &attr);
    BOOST_TEST(err == TIXML_SUCCESS);
    BOOST_TEST(attr == "D[11]");
    err = zone->QueryStringAttribute("ORIGIN", &attr);
    BOOST_TEST(err == TIXML_SUCCESS);
    ParseUtils::GenerateVector(attr, vec);
    BOOST_TEST(NekPoint<NekDouble>(vec.at(0), vec.at(1), vec.at(2)) == origin);
    err = zone->QueryStringAttribute("AXIS", &attr);
    BOOST_TEST(err == TIXML_SUCCESS);
    DNekVec xmlAxis(attr);
    BOOST_TEST(xmlAxis.GetDimension() == 3);
    for (int i = 0; i < 3; i++)
    {
        BOOST_TEST(xmlAxis[i] == axis[i]);
    }
    err = zone->QueryStringAttribute("ANGVEL", &attr);
    BOOST_TEST(err == TIXML_SUCCESS);
    BOOST_TEST(attr == angVelStr);
 
    // Check the prescribed zone
    zone = xmlZones->IterateChildren(zone)->ToElement();
    BOOST_TEST(zone->Value() == "P");
    err = zone->QueryIntAttribute("ID", &id);
    BOOST_TEST(err == TIXML_SUCCESS);
    BOOST_TEST(id == 2);
    err = zone->QueryStringAttribute("DOMAIN", &attr);
    BOOST_TEST(err == TIXML_SUCCESS);
    BOOST_TEST(attr == "D[12]");
    err = zone->QueryStringAttribute("XDEFORM", &attr);
    BOOST_TEST(err == TIXML_SUCCESS);
    BOOST_TEST(attr == xEqnStr);
    err = zone->QueryStringAttribute("YDEFORM", &attr);
    BOOST_TEST(err == TIXML_SUCCESS);
    BOOST_TEST(attr == yEqnStr);
    err = zone->QueryStringAttribute("ZDEFORM", &attr);
    BOOST_TEST(err == TIXML_SUCCESS);
    BOOST_TEST(attr == zEqnStr);
 
    // Check the translating zone
    zone = xmlZones->IterateChildren(zone)->ToElement();
    BOOST_TEST(zone->Value() == "T");
    err = zone->QueryIntAttribute("ID", &id);
    BOOST_TEST(err == TIXML_SUCCESS);
    BOOST_TEST(id == 3);
    err = zone->QueryStringAttribute("DOMAIN", &attr);
    BOOST_TEST(err == TIXML_SUCCESS);
    BOOST_TEST(attr == "D[13]");
    err = zone->QueryStringAttribute("VELOCITY", &attr);
    BOOST_TEST(err == TIXML_SUCCESS);
    vec.clear();
    ParseUtils::GenerateVector(attr, vec);
    BOOST_TEST(vec == velocity);
 
    BOOST_TEST(xmlZones->LastChild() == zone);
 
    TiXmlElement *intr;
 
    // Check the north interface
    TiXmlElement *interface = xmlInterfaces->FirstChildElement();
    BOOST_TEST(interface->Value() == "INTERFACE");
    err = interface->QueryStringAttribute("NAME", &attr);
    BOOST_TEST(err == TIXML_SUCCESS);
    BOOST_TEST(attr == "north");
    intr = interface->FirstChildElement();
    BOOST_TEST(intr->Value() == "L");
    err = intr->QueryIntAttribute("ID", &id);
    BOOST_TEST(err == TIXML_SUCCESS);
    BOOST_TEST(id == 0);
    err = intr->QueryStringAttribute("BOUNDARY", &attr);
    BOOST_TEST(err == TIXML_SUCCESS);
    BOOST_TEST(attr == "C[0]");
    intr = interface->IterateChildren(intr)->ToElement();
    BOOST_TEST(intr->Value() == "R");
    err = intr->QueryIntAttribute("ID", &id);
    BOOST_TEST(err == TIXML_SUCCESS);
    BOOST_TEST(id == 1);
    err = intr->QueryStringAttribute("BOUNDARY", &attr);
    BOOST_TEST(err == TIXML_SUCCESS);
    BOOST_TEST(attr == "C[1-4]");
 
    // Check the south interface
    interface = xmlInterfaces->IterateChildren(interface)->ToElement();
    BOOST_TEST(interface->Value() == "INTERFACE");
    err = interface->QueryStringAttribute("NAME", &attr);
    BOOST_TEST(err == TIXML_SUCCESS);
    BOOST_TEST(attr == "south");
    intr = interface->FirstChildElement();
    BOOST_TEST(intr->Value() == "L");
    err = intr->QueryIntAttribute("ID", &id);
    BOOST_TEST(err == TIXML_SUCCESS);
    BOOST_TEST(id == 2);
    err = intr->QueryStringAttribute("BOUNDARY", &attr);
    BOOST_TEST(err == TIXML_SUCCESS);
    BOOST_TEST(attr == "C[5-9,11]");
    intr = interface->IterateChildren(intr)->ToElement();
    BOOST_TEST(intr->Value() == "R");
    err = intr->QueryIntAttribute("ID", &id);
    BOOST_TEST(err == TIXML_SUCCESS);
    BOOST_TEST(id == 3);
    err = intr->QueryStringAttribute("BOUNDARY", &attr);
    BOOST_TEST(err == TIXML_SUCCESS);
    BOOST_TEST(attr == "C[12]");
 
    BOOST_TEST(xmlInterfaces->LastChild() == interface);
 
    delete nektar;
}

References Nektar::SpatialDomains::Movement::AddInterface(), Nektar::SpatialDomains::Movement::AddZone(), Nektar::MemoryManager< DataType >::AllocateSharedPtr(), angVelStr, axis, CreateInterface(), CreateZone(), Nektar::SpatialDomains::eFixed, Nektar::SpatialDomains::ePrescribe, Nektar::SpatialDomains::eRotate, Nektar::SpatialDomains::eTranslate, Nektar::ParseUtils::GenerateVector(), Nektar::NekVector< DataType >::GetDimension(), Nektar::SpatialDomains::Movement::GetInterfaces(), Nektar::SpatialDomains::Movement::GetZones(), origin, velocity, Nektar::SpatialDomains::Movement::WriteMovement(), xEqnStr, yEqnStr, and zEqnStr.

CreateInterface()

SpatialDomains::InterfaceShPtr Nektar::MovementTests::CreateInterface	(	int	interfaceID,
		std::vector< int >	compositeIDs
	)

Produce dummy Interface objects, containing empty domain pointers.

Definition at line 114 of file TestMovement.cpp.

{
    SpatialDomains::CompositeMap edge;
    for (auto &id : compositeIDs)
    {
        edge[id] =
            MemoryManager<SpatialDomains::Composite>::AllocateSharedPtr();
    }
    return SpatialDomains::InterfaceShPtr(
        MemoryManager<SpatialDomains::Interface>::AllocateSharedPtr(interfaceID,
                                                                    edge));
}

References Nektar::MemoryManager< DataType >::AllocateSharedPtr().

Referenced by BOOST_AUTO_TEST_CASE().

CreateSession()

LibUtilities::SessionReaderSharedPtr Nektar::MovementTests::CreateSession

(

)

Definition at line 46 of file TestMovement.cpp.

{
    char arg     = ' ';
    char *args[] = {&arg};
    return LibUtilities::SessionReader::CreateInstance(1, args);
}

References Nektar::LibUtilities::SessionReader::CreateInstance().

CreateZone()

SpatialDomains::ZoneBaseShPtr Nektar::MovementTests::CreateZone	(	SpatialDomains::MovementType	type,
		int	zoneID,
		int	domainID,
		LibUtilities::InterpreterSharedPtr	interpreter
	)

Produce dummy Zone objects, containing empty domain pointers.

Definition at line 60 of file TestMovement.cpp.

{
    SpatialDomains::CompositeMap domain;
    switch (type)
    {
        case SpatialDomains::MovementType::eFixed:
        {
            return SpatialDomains::ZoneFixedShPtr(
                MemoryManager<SpatialDomains::ZoneFixed>::AllocateSharedPtr(
                    zoneID, domainID, domain, 3));
        }
        break;
        case SpatialDomains::MovementType::eRotate:
        {
            LibUtilities::EquationSharedPtr angularVelEqn =
                std::make_shared<LibUtilities::Equation>(interpreter,
                                                         angVelStr);
            return SpatialDomains::ZoneRotateShPtr(
                MemoryManager<SpatialDomains::ZoneRotate>::AllocateSharedPtr(
                    zoneID, domainID, domain, 3, origin, axis, angularVelEqn));
        }
        break;
        case SpatialDomains::MovementType::eTranslate:
        {
            return SpatialDomains::ZoneTranslateShPtr(
                MemoryManager<SpatialDomains::ZoneTranslate>::AllocateSharedPtr(
                    zoneID, domainID, domain, 3, velocity));
        }
        break;
        case SpatialDomains::MovementType::ePrescribe:
        {
            LibUtilities::EquationSharedPtr
                xEqn = std::make_shared<LibUtilities::Equation>(interpreter,
                                                                xEqnStr),
                yEqn = std::make_shared<LibUtilities::Equation>(interpreter,
                                                                yEqnStr),
                zEqn = std::make_shared<LibUtilities::Equation>(interpreter,
                                                                zEqnStr);
            return SpatialDomains::ZonePrescribeShPtr(
                MemoryManager<SpatialDomains::ZonePrescribe>::AllocateSharedPtr(
                    zoneID, domainID, domain, 3, xEqn, yEqn, zEqn));
        }
        break;
        default:
        {
        }
        break;
    }
    return nullptr;
}

References angVelStr, axis, Nektar::SpatialDomains::eFixed, Nektar::SpatialDomains::ePrescribe, Nektar::SpatialDomains::eRotate, Nektar::SpatialDomains::eTranslate, origin, velocity, xEqnStr, yEqnStr, and zEqnStr.

Referenced by BOOST_AUTO_TEST_CASE().

Variable Documentation

angVelStr

const std::string Nektar::MovementTests::angVelStr = "0.1*t"

Definition at line 53 of file TestMovement.cpp.

Referenced by BOOST_AUTO_TEST_CASE(), and CreateZone().

axis

const DNekVec Nektar::MovementTests::axis = {1., 2., 3.}

Definition at line 56 of file TestMovement.cpp.

Referenced by BOOST_AUTO_TEST_CASE(), CreateZone(), Nektar::SpatialDomains::Movement::ReadZones(), and ZoneRotate_Init().

origin

const NekPoint<NekDouble> Nektar::MovementTests::origin = {1., 2., 3.}

Definition at line 55 of file TestMovement.cpp.

Referenced by BOOST_AUTO_TEST_CASE(), CreateZone(), Nektar::Collections::Helmholtz_IterPerExp::operator()(), Nektar::SpatialDomains::Movement::ReadZones(), and ZoneRotate_Init().

velocity

const std::vector<NekDouble> Nektar::MovementTests::velocity = {1., 1., 2.}

Definition at line 57 of file TestMovement.cpp.

Referenced by BOOST_AUTO_TEST_CASE(), Nektar::SolverUtils::FilterAeroForces::CalculateForces(), CreateZone(), Nektar::IncNavierStokes::EvaluateAdvectionTerms(), Nektar::MMFAdvection::EvaluateAdvectionVelocity(), Nektar::VariableConverter::GetAbsoluteVelocity(), Nektar::NonlinearPeregrine::GetFluxVector(), Nektar::NonlinearSWE::GetFluxVector(), Nektar::CompressibleFlowSystem::GetFluxVector(), Nektar::CompressibleFlowSystem::GetFluxVectorDeAlias(), Nektar::SolverUtils::FluidInterface::GetVelocity(), Nektar::VariableConverter::GetVelocityVector(), Nektar::LinearSWE::GetVelocityVector(), Nektar::NonlinearPeregrine::GetVelocityVector(), Nektar::NonlinearSWE::GetVelocityVector(), Nektar::MultiRegions::ExpList::LinearAdvectionReactionSolve(), main(), OUTPUT(), Nektar::SpatialDomains::Movement::ReadZones(), Nektar::AdjointAdvection::v_Advect(), Nektar::AlternateSkewAdvection::v_Advect(), Nektar::LinearisedAdvection::v_Advect(), Nektar::NavierStokesAdvection::v_Advect(), Nektar::SkewSymmetricAdvection::v_Advect(), Nektar::CompressibleFlowSystem::v_ExtraFldOutput(), Nektar::NavierStokesCFE::v_ExtraFldOutput(), Nektar::AcousticSystem::v_GetMaxStdVelocity(), Nektar::CompressibleFlowSystem::v_GetMaxStdVelocity(), Nektar::CompressibleFlowSystem::v_GetVelocity(), Nektar::IncNavierStokes::v_GetVelocity(), Nektar::SolverUtils::FileSolution::v_GetVelocity(), Nektar::CoupledLinearNS::v_InitObject(), Nektar::MultiRegions::ContField::v_LinearAdvectionReactionSolve(), Nektar::FieldUtils::ProcessWSS::v_Process(), Nektar::VCSMapping::v_SetUpPressureForcing(), Nektar::VCSImplicit::v_SetUpViscousForcing(), and ZoneTranslate_Init().

xEqnStr

const std::string Nektar::MovementTests::xEqnStr = "0.1*x - 0.1*t"

Definition at line 53 of file TestMovement.cpp.

Referenced by BOOST_AUTO_TEST_CASE(), and CreateZone().

yEqnStr

const std::string Nektar::MovementTests::yEqnStr = "0.1*y^2 - x"

Definition at line 54 of file TestMovement.cpp.

Referenced by BOOST_AUTO_TEST_CASE(), and CreateZone().

zEqnStr

const std::string Nektar::MovementTests::zEqnStr = "sqrt(t)"

Definition at line 54 of file TestMovement.cpp.

Referenced by BOOST_AUTO_TEST_CASE(), and CreateZone().