#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

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.

64 {

65 }

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_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";
     if (it->second == "2DBndLayer")
     {
         m_makeBL = true;
         m_2D = true;
     }
     if(!m_makeBL && !m_2D)
     {
         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::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_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, and ParseFile().

 {
     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;
 
     vector<ModuleSharedPtr> mods;
 
     ////**** CAD ****////
     mods.push_back(GetModuleFactory().CreateInstance(
         ModuleKey(eProcessModule, "loadcad"), m_mesh));
     mods.back()->RegisterConfig("filename", m_cadfile);
 
     if(m_2D)
     {
         mods.back()->RegisterConfig("2D","");
     }
     if(m_naca)
     {
         mods.back()->RegisterConfig("NACA",m_nacadomain);
     }
 
     ////**** Octree ****////
     mods.push_back(GetModuleFactory().CreateInstance(
         ModuleKey(eProcessModule, "loadoctree"), m_mesh));
     mods.back()->RegisterConfig("mindel", m_minDelta);
     mods.back()->RegisterConfig("maxdel", m_maxDelta);
     mods.back()->RegisterConfig("eps", m_eps);
     if (m_refine)
     {
         mods.back()->RegisterConfig("refinement", m_refinement);
     }
     if (m_woct)
     {
         mods.back()->RegisterConfig("writeoctree", "");
     }
 
     if(m_2D)
     {
         m_mesh->m_expDim = 2;
         m_mesh->m_spaceDim = 2;
         mods.push_back(GetModuleFactory().CreateInstance(
             ModuleKey(eProcessModule, "2dgenerator"), m_mesh));
         if (m_makeBL)
         {
             mods.back()->RegisterConfig("blcurves", m_blsurfs);
             mods.back()->RegisterConfig("blthick", m_blthick);
         }
     }
     else
     {
         ////**** SurfaceMesh ****////
         mods.push_back(GetModuleFactory().CreateInstance(
             ModuleKey(eProcessModule, "surfacemesh"), m_mesh));
 
         ////**** VolumeMesh ****////
         mods.push_back(GetModuleFactory().CreateInstance(
             ModuleKey(eProcessModule, "volumemesh"), m_mesh));
         if(m_makeBL)
         {
             mods.back()->RegisterConfig("blsurfs",m_blsurfs);
             mods.back()->RegisterConfig("blthick",m_blthick);
             mods.back()->RegisterConfig("bllayers",m_bllayers);
             mods.back()->RegisterConfig("blprog",m_blprog);
         }
     }
 
     ////**** HOSurface ****////
     mods.push_back(GetModuleFactory().CreateInstance(
         ModuleKey(eProcessModule, "hosurface"), m_mesh));
     if (m_surfopti)
     {
         mods.back()->RegisterConfig("opti", "");
     }
 
     ////*** 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;
         }
         mods.push_back(GetModuleFactory().CreateInstance(
             ModuleKey(eProcessModule, "varopti"), m_mesh));
         mods.back()->RegisterConfig("hyperelastic","");
         mods.back()->RegisterConfig("maxiter","10");
         mods.back()->RegisterConfig("numthreads",boost::lexical_cast<string>(np));
     }
 
     ////**** SPLIT BL ****////
     if(m_splitBL)
     {
         mods.push_back(GetModuleFactory().CreateInstance(
             ModuleKey(eProcessModule, "bl"), m_mesh));
         mods.back()->RegisterConfig("layers",m_bllayers);
         mods.back()->RegisterConfig("surf",m_blsurfs);
         mods.back()->RegisterConfig("nq",boost::lexical_cast<string>(m_mesh->m_nummode));
         mods.back()->RegisterConfig("r",m_blprog);
     }
 
     for(int i = 0; i < mods.size(); i++)
     {
         mods[i]->SetDefaults();
         mods[i]->Process();
     }
 }

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().

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().

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().

Public Member Functions

Static Public Member Functions

Static Public Attributes

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation