Nektar::Utilities::InputMCF Class Reference

#include <InputMCF.h>

Inheritance diagram for Nektar::Utilities::InputMCF:

Collaboration diagram for Nektar::Utilities::InputMCF:

Public Member Functions
	InputMCF (NekMeshUtils::MeshSharedPtr m)
	Set up InputCAD object. More...
virtual	~InputMCF ()
virtual void	Process ()
void	ParseFile (std::string nm)
Public Member Functions inherited from Nektar::NekMeshUtils::InputModule
NEKMESHUTILS_EXPORT	InputModule (MeshSharedPtr p_m)
NEKMESHUTILS_EXPORT void	OpenStream ()
	Open a file for input. More...
Public Member Functions inherited from Nektar::NekMeshUtils::Module
NEKMESHUTILS_EXPORT	Module (MeshSharedPtr p_m)
NEKMESHUTILS_EXPORT void	RegisterConfig (std::string key, std::string value)
	Register a configuration option with a module. More...
NEKMESHUTILS_EXPORT void	PrintConfig ()
	Print out all configuration options for a module. More...
NEKMESHUTILS_EXPORT void	SetDefaults ()
	Sets default configuration options for those which have not been set. More...
NEKMESHUTILS_EXPORT MeshSharedPtr	GetMesh ()
virtual NEKMESHUTILS_EXPORT void	ProcessVertices ()
	Extract element vertices. More...
virtual NEKMESHUTILS_EXPORT void	ProcessEdges (bool ReprocessEdges=true)
	Extract element edges. More...
virtual NEKMESHUTILS_EXPORT void	ProcessFaces (bool ReprocessFaces=true)
	Extract element faces. More...
virtual NEKMESHUTILS_EXPORT void	ProcessElements ()
	Generate element IDs. More...
virtual NEKMESHUTILS_EXPORT void	ProcessComposites ()
	Generate composites. More...
virtual NEKMESHUTILS_EXPORT void	ClearElementLinks ()

Static Public Member Functions
static NekMeshUtils::ModuleSharedPtr	create (NekMeshUtils::MeshSharedPtr m)
	Creates an instance of this class. More...

Static Public Attributes
static NekMeshUtils::ModuleKey	className
	ModuleKey for class. More...

Private Attributes
std::string	m_minDelta
std::string	m_maxDelta
std::string	m_eps
std::string	m_cadfile
std::string	m_order
std::string	m_blsurfs
std::string	m_blthick
std::string	m_blprog
std::string	m_bllayers
std::string	m_refinement
std::string	m_nacadomain
bool	m_makeBL
bool	m_surfopti
bool	m_varopti
bool	m_refine
bool	m_woct
bool	m_2D
bool	m_splitBL
bool	m_naca
bool	m_adjust
bool	m_adjustall
bool	m_manifold

Additional Inherited Members
Protected Member Functions inherited from Nektar::NekMeshUtils::InputModule
NEKMESHUTILS_EXPORT void	PrintSummary ()
	Print summary of elements. More...
Protected Member Functions inherited from Nektar::NekMeshUtils::Module
NEKMESHUTILS_EXPORT void	ReorderPrisms (PerMap &perFaces)
	Reorder node IDs so that prisms and tetrahedra are aligned correctly. More...
NEKMESHUTILS_EXPORT void	PrismLines (int prism, PerMap &perFaces, std::set< int > &prismsDone, std::vector< ElementSharedPtr > &line)
Protected Attributes inherited from Nektar::NekMeshUtils::InputModule
io::filtering_istream	m_mshFile
	Input stream. More...
std::ifstream	m_mshFileStream
	Input stream. More...
Protected Attributes inherited from Nektar::NekMeshUtils::Module
MeshSharedPtr	m_mesh
	Mesh object. More...
std::map< std::string, ConfigOption >	m_config
	List of configuration values. More...

Detailed Description

Definition at line 46 of file InputMCF.h.

Constructor & Destructor Documentation

Nektar::Utilities::InputMCF::InputMCF

(

NekMeshUtils::MeshSharedPtr

m

)

Set up InputCAD object.

Definition at line 59 of file InputMCF.cpp.

                                  : InputModule(m)
{
}

Nektar::Utilities::InputMCF::~InputMCF ( )

virtual

Definition at line 63 of file InputMCF.cpp.

{
}

Member Function Documentation

static NekMeshUtils::ModuleSharedPtr Nektar::Utilities::InputMCF::create ( NekMeshUtils::MeshSharedPtr  m )

inlinestatic

Creates an instance of this class.

Definition at line 54 of file InputMCF.h.

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

    {
        return MemoryManager<InputMCF>::AllocateSharedPtr(m);
    }

void Nektar::Utilities::InputMCF::ParseFile

(

std::string

nm

)

Definition at line 67 of file InputMCF.cpp.

References ASSERTL0, Nektar::LibUtilities::SessionReader::CreateInstance(), Nektar::iterator, L, m_2D, m_bllayers, m_blprog, m_blsurfs, m_blthick, m_cadfile, m_eps, m_makeBL, m_manifold, m_maxDelta, m_minDelta, m_naca, m_nacadomain, m_order, m_refine, m_refinement, m_splitBL, m_surfopti, m_varopti, m_woct, and class_topology::P.

Referenced by Process().

{
    vector<string> filename;
    filename.push_back(nm);
    LibUtilities::SessionReaderSharedPtr pSession =
        LibUtilities::SessionReader::CreateInstance(0, NULL, filename);

    ASSERTL0(pSession->DefinesElement("NEKTAR/MESHING"), "no meshing tag");
    ASSERTL0(pSession->DefinesElement("NEKTAR/MESHING/INFORMATION"),
             "no information tag");
    ASSERTL0(pSession->DefinesElement("NEKTAR/MESHING/PARAMETERS"),
             "no parameters tag");

    TiXmlElement *mcf = pSession->GetElement("NEKTAR/MESHING");

    TiXmlElement *info = mcf->FirstChildElement("INFORMATION");
    TiXmlElement *I    = info->FirstChildElement("I");
    map<string, string> information;
    while (I)
    {
        string tmp1, tmp2;
        I->QueryStringAttribute("PROPERTY", &tmp1);
        I->QueryStringAttribute("VALUE", &tmp2);
        information[tmp1] = tmp2;
        I                 = I->NextSiblingElement("I");
    }

    TiXmlElement *param = mcf->FirstChildElement("PARAMETERS");
    TiXmlElement *P     = param->FirstChildElement("P");
    map<string, string> parameters;
    while (P)
    {
        string tmp1, tmp2;
        P->QueryStringAttribute("PARAM", &tmp1);
        P->QueryStringAttribute("VALUE", &tmp2);
        parameters[tmp1] = tmp2;
        P                = P->NextSiblingElement("P");
    }

    set<string> boolparameters;

    if (pSession->DefinesElement("NEKTAR/MESHING/BOOLPARAMETERS"))
    {
        TiXmlElement *bparam = mcf->FirstChildElement("BOOLPARAMETERS");
        TiXmlElement *BP     = bparam->FirstChildElement("P");

        while (BP)
        {
            string tmp;
            BP->QueryStringAttribute("VALUE", &tmp);
            boolparameters.insert(tmp);
            BP = BP->NextSiblingElement("P");
        }
    }

    set<string> refinement;
    if (pSession->DefinesElement("NEKTAR/MESHING/REFINEMENT"))
    {
        TiXmlElement *refine = mcf->FirstChildElement("REFINEMENT");
        TiXmlElement *L      = refine->FirstChildElement("LINE");

        while (L)
        {
            stringstream ss;
            TiXmlElement *T = L->FirstChildElement("X1");
            ss << T->GetText() << ",";
            T = L->FirstChildElement("Y1");
            ss << T->GetText() << ",";
            T = L->FirstChildElement("Z1");
            ss << T->GetText() << ",";
            T = L->FirstChildElement("X2");
            ss << T->GetText() << ",";
            T = L->FirstChildElement("Y2");
            ss << T->GetText() << ",";
            T = L->FirstChildElement("Z2");
            ss << T->GetText() << ",";
            T = L->FirstChildElement("R");
            ss << T->GetText() << ",";
            T = L->FirstChildElement("D");
            ss << T->GetText();

            refinement.insert(ss.str());

            L = L->NextSiblingElement("LINE");
        }
    }

    map<string, string>::iterator it;

    it = information.find("CADFile");
    ASSERTL0(it != information.end(), "no cadfile defined");
    m_cadfile = it->second;

    it = information.find("MeshType");
    ASSERTL0(it != information.end(), "no meshtype defined");
    m_makeBL = it->second == "3DBndLayer";
    m_2D     = it->second == "2D";
    m_manifold = it->second == "Manifold";
    if (it->second == "2DBndLayer")
    {
        m_makeBL = true;
        m_2D     = true;
    }
    if (!m_makeBL && !m_2D && !m_manifold)
    {
        ASSERTL0(it->second == "3D", "unsure on MeshType")
    }

    it = parameters.find("MinDelta");
    ASSERTL0(it != parameters.end(), "no mindelta defined");
    m_minDelta = it->second;

    it = parameters.find("MaxDelta");
    ASSERTL0(it != parameters.end(), "no maxdelta defined");
    m_maxDelta = it->second;

    it = parameters.find("EPS");
    ASSERTL0(it != parameters.end(), "no eps defined");
    m_eps = it->second;

    it = parameters.find("Order");
    ASSERTL0(it != parameters.end(), "no order defined");
    m_order = it->second;

    if (m_makeBL)
    {
        it = parameters.find("BndLayerSurfaces");
        ASSERTL0(it != parameters.end(), "no BndLayersurfs defined");
        m_blsurfs = it->second;

        it = parameters.find("BndLayerThickness");
        ASSERTL0(it != parameters.end(), "no BndLayerthick defined");
        m_blthick = it->second;

        it        = parameters.find("BndLayerLayers");
        m_splitBL = it != parameters.end();
        if (m_splitBL)
        {
            m_bllayers = it->second;
            it         = parameters.find("BndLayerProgression");
            m_blprog   = it != parameters.end() ? it->second : "2.0";
        }
    }

    m_naca = false;
    if (m_2D && m_cadfile.find('.') == std::string::npos)
    {
        m_naca = true;

        stringstream ss;
        it = parameters.find("Xmin");
        ASSERTL0(it != parameters.end(), "no xmin defined");
        ss << it->second << ",";
        it = parameters.find("Ymin");
        ASSERTL0(it != parameters.end(), "no ymin defined");
        ss << it->second << ",";
        it = parameters.find("Xmax");
        ASSERTL0(it != parameters.end(), "no xmax defined");
        ss << it->second << ",";
        it = parameters.find("Ymax");
        ASSERTL0(it != parameters.end(), "no zmax defined");
        ss << it->second << ",";
        it = parameters.find("AOA");
        ASSERTL0(it != parameters.end(), "no aoa defined");
        ss << it->second;

        m_nacadomain = ss.str();
    }

    set<string>::iterator sit;
    sit        = boolparameters.find("SurfaceOptimiser");
    m_surfopti = sit != boolparameters.end();
    sit        = boolparameters.find("WriteOctree");
    m_woct     = sit != boolparameters.end();
    sit        = boolparameters.find("VariationalOptimiser");
    m_varopti  = sit != boolparameters.end();

    m_refine = refinement.size() > 0;
    if (m_refine)
    {
        stringstream ss;
        for (sit = refinement.begin(); sit != refinement.end(); sit++)
        {
            ss << *sit;
            ss << ":";
        }
        m_refinement = ss.str();
        m_refinement.erase(m_refinement.end() - 1);
    }
}

void Nektar::Utilities::InputMCF::Process ( )

virtual

Implements Nektar::NekMeshUtils::Module.

Definition at line 258 of file InputMCF.cpp.

References Nektar::LibUtilities::NekFactory< tKey, tBase, >::CreateInstance(), Nektar::FieldUtils::eProcessModule, Nektar::FieldUtils::GetModuleFactory(), m_2D, m_bllayers, m_blprog, m_blsurfs, m_blthick, m_cadfile, Nektar::NekMeshUtils::Module::m_config, m_eps, m_makeBL, m_manifold, m_maxDelta, Nektar::NekMeshUtils::Module::m_mesh, m_minDelta, m_naca, m_nacadomain, m_order, m_refine, m_refinement, m_splitBL, m_surfopti, m_varopti, m_woct, ParseFile(), Nektar::NekMeshUtils::Module::ProcessComposites(), Nektar::NekMeshUtils::Module::ProcessEdges(), Nektar::NekMeshUtils::Module::ProcessElements(), Nektar::NekMeshUtils::Module::ProcessFaces(), and Nektar::NekMeshUtils::Module::ProcessVertices().

{
    ParseFile(m_config["infile"].as<string>());

    m_mesh->m_expDim   = 3;
    m_mesh->m_spaceDim = 3;
    m_mesh->m_nummode  = boost::lexical_cast<int>(m_order) + 1;

    ModuleSharedPtr module;

    ////**** CAD ****////
    module = GetModuleFactory().CreateInstance(
        ModuleKey(eProcessModule, "loadcad"), m_mesh);
    module->RegisterConfig("filename", m_cadfile);

    if (m_2D)
    {
        module->RegisterConfig("2D", "");
    }
    if (m_naca)
    {
        module->RegisterConfig("NACA", m_nacadomain);
    }

    module->SetDefaults();
    module->Process();

    ////**** OCTREE ****////
    module = GetModuleFactory().CreateInstance(
        ModuleKey(eProcessModule, "loadoctree"), m_mesh);
    module->RegisterConfig("mindel", m_minDelta);
    module->RegisterConfig("maxdel", m_maxDelta);
    module->RegisterConfig("eps", m_eps);
    if (m_refine)
    {
        module->RegisterConfig("refinement", m_refinement);
    }
    if (m_woct)
    {
        module->RegisterConfig("writeoctree", "");
    }

    module->SetDefaults();
    module->Process();

    ////**** LINEAR MESHING ****////
    if (m_2D)
    {
        m_mesh->m_expDim   = 2;
        m_mesh->m_spaceDim = 2;
        module             = GetModuleFactory().CreateInstance(
            ModuleKey(eProcessModule, "2dgenerator"), m_mesh);
        if (m_makeBL)
        {
            module->RegisterConfig("blcurves", m_blsurfs);
            module->RegisterConfig("blthick", m_blthick);
        }

        try
        {
            module->SetDefaults();
            module->Process();
        }
        catch (runtime_error &e)
        {
            cout << "2D linear mesh generator failed with message:" << endl;
            cout << e.what() << endl;
            cout << "No mesh file has been created" << endl;
            abort();
        }
    }
    else
    {
        ////**** SurfaceMesh ****////
        module = GetModuleFactory().CreateInstance(
            ModuleKey(eProcessModule, "surfacemesh"), m_mesh);

        try
        {
            module->SetDefaults();
            module->Process();
        }
        catch (runtime_error &e)
        {
            cout << "Surface meshing has failed with message:" << endl;
            cout << e.what() << endl;
            cout << "Any surfaces which were succsessfully meshed will be "
                    "dumped as a manifold mesh"
                 << endl;
            m_mesh->m_expDim = 2;
            ProcessVertices();
            ProcessEdges();
            ProcessFaces();
            ProcessElements();
            ProcessComposites();
            return;
        }

        if(m_manifold)
        {
            //dont want to volume mesh
            m_mesh->m_expDim = 2;
        }
        else
        {
            ////**** VolumeMesh ****////
            module = GetModuleFactory().CreateInstance(
                ModuleKey(eProcessModule, "volumemesh"), m_mesh);
            if (m_makeBL)
            {
                module->RegisterConfig("blsurfs", m_blsurfs);
                module->RegisterConfig("blthick", m_blthick);
                module->RegisterConfig("bllayers", m_bllayers);
                module->RegisterConfig("blprog", m_blprog);
            }

            try
            {
                module->SetDefaults();
                module->Process();
            }
            catch (runtime_error &e)
            {
                cout << "Volume meshing has failed with message:" << endl;
                cout << e.what() << endl;
                cout << "The linear surface mesh be dumped as a manifold mesh"
                     << endl;
                m_mesh->m_expDim = 2;
                m_mesh->m_element[3].clear();
                ProcessVertices();
                ProcessEdges();
                ProcessFaces();
                ProcessElements();
                ProcessComposites();
                return;
            }
        }
    }

    ////**** HOSurface ****////
    module = GetModuleFactory().CreateInstance(
        ModuleKey(eProcessModule, "hosurface"), m_mesh);
    if (m_surfopti)
    {
        module->RegisterConfig("opti", "");
    }

    try
    {
        module->SetDefaults();
        module->Process();
    }
    catch (runtime_error &e)
    {
        cout << "High-order surface meshing has failed with message:" << endl;
        cout << e.what() << endl;
        cout << "The mesh will be written as normal but the incomplete surface "
                "will remain faceted"
             << endl;
        return;
    }

    ////*** VARIATIONAL OPTIMISATION ****////
    if (m_varopti)
    {
        unsigned int np = boost::thread::physical_concurrency();
        if (m_mesh->m_verbose)
        {
            cout << "Detecting 4 cores, will attempt to run in parrallel"
                 << endl;
        }
        module = GetModuleFactory().CreateInstance(
            ModuleKey(eProcessModule, "varopti"), m_mesh);
        module->RegisterConfig("hyperelastic", "");
        module->RegisterConfig("maxiter", "10");
        module->RegisterConfig("numthreads",
                                    boost::lexical_cast<string>(np));

        try
        {
            module->SetDefaults();
            module->Process();
        }
        catch (runtime_error &e)
        {
            cout << "Variational optimisation has failed with message:" << endl;
            cout << e.what() << endl;
            cout << "The mesh will be written as is, it may be invalid" << endl;
            return;
        }
    }

    ////**** SPLIT BL ****////
    if (m_splitBL)
    {
        module = GetModuleFactory().CreateInstance(
            ModuleKey(eProcessModule, "bl"), m_mesh);
        module->RegisterConfig("layers", m_bllayers);
        module->RegisterConfig("surf", m_blsurfs);
        module->RegisterConfig(
            "nq", boost::lexical_cast<string>(m_mesh->m_nummode));
        module->RegisterConfig("r", m_blprog);

        try
        {
            module->SetDefaults();
            module->Process();
        }
        catch (runtime_error &e)
        {
            cout << "Boundary layer splitting has failed with message:" << endl;
            cout << e.what() << endl;
            cout << "The mesh will be written as is, it may be invalid" << endl;
            return;
        }
    }


}

Member Data Documentation

ModuleKey Nektar::Utilities::InputMCF::className

static

Initial value:

= GetModuleFactory().RegisterCreatorFunction(
    ModuleKey(eInputModule, "mcf"), InputMCF::create,
    "Reads mesh configuration and will generate the mesh file.")

ModuleKey for class.

Definition at line 59 of file InputMCF.h.

bool Nektar::Utilities::InputMCF::m_2D

private

Definition at line 67 of file InputMCF.h.

Referenced by ParseFile(), and Process().

bool Nektar::Utilities::InputMCF::m_adjust

private

Definition at line 67 of file InputMCF.h.

bool Nektar::Utilities::InputMCF::m_adjustall

private

Definition at line 67 of file InputMCF.h.

std::string Nektar::Utilities::InputMCF::m_bllayers

private

Definition at line 64 of file InputMCF.h.

Referenced by ParseFile(), and Process().

std::string Nektar::Utilities::InputMCF::m_blprog

private

Definition at line 64 of file InputMCF.h.

Referenced by ParseFile(), and Process().

std::string Nektar::Utilities::InputMCF::m_blsurfs

private

Definition at line 64 of file InputMCF.h.

Referenced by ParseFile(), and Process().

std::string Nektar::Utilities::InputMCF::m_blthick

private

Definition at line 64 of file InputMCF.h.

Referenced by ParseFile(), and Process().

std::string Nektar::Utilities::InputMCF::m_cadfile

private

Definition at line 64 of file InputMCF.h.

Referenced by ParseFile(), and Process().

std::string Nektar::Utilities::InputMCF::m_eps

private

Definition at line 64 of file InputMCF.h.

Referenced by ParseFile(), and Process().

bool Nektar::Utilities::InputMCF::m_makeBL

private

Definition at line 67 of file InputMCF.h.

Referenced by ParseFile(), and Process().

bool Nektar::Utilities::InputMCF::m_manifold

private

Definition at line 67 of file InputMCF.h.

Referenced by ParseFile(), and Process().

std::string Nektar::Utilities::InputMCF::m_maxDelta

private

Definition at line 64 of file InputMCF.h.

Referenced by ParseFile(), and Process().

std::string Nektar::Utilities::InputMCF::m_minDelta

private

Definition at line 64 of file InputMCF.h.

Referenced by ParseFile(), and Process().

bool Nektar::Utilities::InputMCF::m_naca

private

Definition at line 67 of file InputMCF.h.

Referenced by ParseFile(), and Process().

std::string Nektar::Utilities::InputMCF::m_nacadomain

private

Definition at line 64 of file InputMCF.h.

Referenced by ParseFile(), and Process().

std::string Nektar::Utilities::InputMCF::m_order

private

Definition at line 64 of file InputMCF.h.

Referenced by ParseFile(), and Process().

bool Nektar::Utilities::InputMCF::m_refine

private

Definition at line 67 of file InputMCF.h.

Referenced by ParseFile(), and Process().

std::string Nektar::Utilities::InputMCF::m_refinement

private

Definition at line 64 of file InputMCF.h.

Referenced by ParseFile(), and Process().

bool Nektar::Utilities::InputMCF::m_splitBL

private

Definition at line 67 of file InputMCF.h.

Referenced by ParseFile(), and Process().

bool Nektar::Utilities::InputMCF::m_surfopti

private

Definition at line 67 of file InputMCF.h.

Referenced by ParseFile(), and Process().

bool Nektar::Utilities::InputMCF::m_varopti

private

Definition at line 67 of file InputMCF.h.

Referenced by ParseFile(), and Process().

bool Nektar::Utilities::InputMCF::m_woct

private

Definition at line 67 of file InputMCF.h.

Referenced by ParseFile(), and Process().