Nektar::Utilities::InputCAD Class Reference

#include <InputCAD.h>

Inheritance diagram for Nektar::Utilities::InputCAD:

Collaboration diagram for Nektar::Utilities::InputCAD:

Public Member Functions
	InputCAD (MeshSharedPtr m)
	Set up InputCAD object. More...
virtual	~InputCAD ()
virtual void	Process ()
Public Member Functions inherited from Nektar::Utilities::InputModule
	InputModule (FieldSharedPtr p_m)
void	AddFile (string fileType, string fileName)
	InputModule (MeshSharedPtr p_m)
void	OpenStream ()
	Open a file for input. More...
Public Member Functions inherited from Nektar::Utilities::Module
	Module (FieldSharedPtr p_f)
virtual void	Process (po::variables_map &vm)=0
virtual std::string	GetModuleName ()=0
void	RegisterConfig (string key, string value)
	Register a configuration option with a module. More...
void	PrintConfig ()
	Print out all configuration options for a module. More...
void	SetDefaults ()
	Sets default configuration options for those which have not been set. More...
bool	GetRequireEquiSpaced (void)
void	SetRequireEquiSpaced (bool pVal)
void	EvaluateTriFieldAtEquiSpacedPts (LocalRegions::ExpansionSharedPtr &exp, const Array< OneD, const NekDouble > &infield, Array< OneD, NekDouble > &outfield)
	Module (MeshSharedPtr p_m)
void	RegisterConfig (std::string key, std::string value)
void	PrintConfig ()
void	SetDefaults ()
MeshSharedPtr	GetMesh ()
virtual void	ProcessVertices ()
	Extract element vertices. More...
virtual void	ProcessEdges (bool ReprocessEdges=true)
	Extract element edges. More...
virtual void	ProcessFaces (bool ReprocessFaces=true)
	Extract element faces. More...
virtual void	ProcessElements ()
	Generate element IDs. More...
virtual void	ProcessComposites ()
	Generate composites. More...
virtual void	ClearElementLinks ()

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

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

Private Attributes
NekDouble	m_minDelta
NekDouble	m_maxDelta
NekDouble	m_eps
NekDouble	m_blthick
int	m_order
string	m_CADName
string	m_udsName
bool	m_makeBL
bool	m_writeoctree

Additional Inherited Members
Protected Member Functions inherited from Nektar::Utilities::InputModule
void	PrintSummary ()
	Print summary of elements. More...
void	PrintSummary ()
	Print summary of elements. More...
Protected Member Functions inherited from Nektar::Utilities::Module
	Module ()
void	ReorderPrisms (PerMap &perFaces)
	Reorder node IDs so that prisms and tetrahedra are aligned correctly. More...
void	PrismLines (int prism, PerMap &perFaces, std::set< int > &prismsDone, std::vector< ElementSharedPtr > &line)
Protected Attributes inherited from Nektar::Utilities::InputModule
set< string >	m_allowedFiles
std::ifstream	m_mshFile
	Input stream. More...
Protected Attributes inherited from Nektar::Utilities::Module
FieldSharedPtr	m_f
	Field object. More...
map< string, ConfigOption >	m_config
	List of configuration values. More...
bool	m_requireEquiSpaced
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 InputCAD.h.

Constructor & Destructor Documentation

Nektar::Utilities::InputCAD::InputCAD

(

MeshSharedPtr

m

)

Set up InputCAD object.

Definition at line 67 of file InputCAD.cpp.

                                  : InputModule(m)
{
}

Nektar::Utilities::InputCAD::~InputCAD ( )

virtual

Definition at line 71 of file InputCAD.cpp.

{
}

Member Function Documentation

static ModuleSharedPtr Nektar::Utilities::InputCAD::create ( MeshSharedPtr  m )

inlinestatic

Creates an instance of this class.

Definition at line 54 of file InputCAD.h.

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

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

void Nektar::Utilities::InputCAD::Process ( )

virtual

Implements Nektar::Utilities::Module.

Definition at line 75 of file InputCAD.cpp.

References Nektar::MemoryManager< DataType >::AllocateSharedPtr(), ASSERTL0, Nektar::Utilities::Module::ClearElementLinks(), Nektar::LibUtilities::SessionReader::CreateInstance(), Nektar::LibUtilities::NekFactory< tKey, tBase, >::CreateInstance(), Nektar::Utilities::eOutputModule, Nektar::StdRegions::find(), Nektar::ParseUtils::GenerateSeqVector(), Nektar::Utilities::GetModuleFactory(), m_blthick, m_CADName, Nektar::Utilities::Module::m_config, m_eps, m_makeBL, m_maxDelta, Nektar::Utilities::Module::m_mesh, m_minDelta, m_order, m_udsName, m_writeoctree, Nektar::Utilities::Module::ProcessComposites(), Nektar::Utilities::Module::ProcessEdges(), Nektar::Utilities::Module::ProcessElements(), Nektar::Utilities::Module::ProcessFaces(), and Nektar::Utilities::Module::ProcessVertices().

{
    vector<string> filename;
    filename.push_back(m_config["infile"].as<string>());
    string fn = filename[0].substr(0, filename[0].find("."));

    LibUtilities::SessionReaderSharedPtr pSession =
        LibUtilities::SessionReader::CreateInstance(0, NULL, filename);

    // these parameters must be defined for any mesh generation to work
    pSession->LoadParameter("MinDelta", m_minDelta);
    pSession->LoadParameter("MaxDelta", m_maxDelta);
    pSession->LoadParameter("EPS", m_eps);
    pSession->LoadParameter("Order", m_order);
    m_CADName = pSession->GetSolverInfo("CADFile");

    if (pSession->DefinesSolverInfo("MeshType"))
    {
        if (pSession->GetSolverInfo("MeshType") == "BL")
        {
            m_makeBL = true;
            pSession->LoadParameter("BLThick", m_blthick);
        }
        else
        {
            m_makeBL = false;
        }
    }
    else
    {
        m_makeBL = false;
    }

    if (pSession->DefinesSolverInfo("WriteOctree"))
    {
        m_writeoctree = pSession->GetSolverInfo("WriteOctree") == "TRUE";
    }

    vector<unsigned int> symsurfs;
    vector<unsigned int> blsurfs;
    if (m_makeBL)
    {
        string sym, bl;
        sym = pSession->GetSolverInfo("SymPlane");
        bl = pSession->GetSolverInfo("BLSurfs");
        ParseUtils::GenerateSeqVector(sym.c_str(), symsurfs);
        ParseUtils::GenerateSeqVector(bl.c_str(), blsurfs);
        sort(symsurfs.begin(), symsurfs.end());
        sort(blsurfs.begin(), blsurfs.end());
        ASSERTL0(blsurfs.size() > 0,
                 "No surfaces selected to make boundary layer on");
    }

    if (pSession->DefinesSolverInfo("UserDefinedSpacing"))
    {
        m_udsName = pSession->GetSolverInfo("UserDefinedSpacing");
        ASSERTL0(boost::filesystem::exists(m_udsName.c_str()),
                 "UserDefinedSpacing file does not exist");
    }
    else
    {
        m_udsName = "N";
    }

    CADSystemSharedPtr m_cad =
        MemoryManager<CADSystem>::AllocateSharedPtr(m_CADName);

    if (m_mesh->m_verbose)
    {
        cout << "Building mesh for: " << m_CADName << endl;
    }

    ASSERTL0(m_cad->LoadCAD(), "Failed to load CAD");

    if (m_mesh->m_verbose)
    {
        cout << "With parameters:" << endl;
        cout << "\tmin delta: " << m_minDelta << endl
             << "\tmax delta: " << m_maxDelta << endl
             << "\tesp: " << m_eps << endl
             << "\torder: " << m_order << endl;
        m_cad->Report();
    }

    if (m_makeBL && m_mesh->m_verbose)
    {

        cout << "\tWill make boundary layers on surfs: ";
        for (int i = 0; i < blsurfs.size(); i++)
        {
            cout << blsurfs[i] << " ";
        }
        cout << endl << "\tWith the symmetry planes: ";
        for (int i = 0; i < symsurfs.size(); i++)
        {
            cout << symsurfs[i] << " ";
        }
        cout << endl << "\tWith thickness " << m_blthick << endl;
    }

    // create octree
    OctreeSharedPtr m_octree = MemoryManager<Octree>::AllocateSharedPtr(
        m_cad, m_mesh->m_verbose, m_minDelta, m_maxDelta, m_eps, m_udsName);

    m_octree->Build();

    if (m_writeoctree)
    {
        MeshSharedPtr oct = boost::shared_ptr<Mesh>(new Mesh());
        oct->m_expDim     = 3;
        oct->m_spaceDim   = 3;
        oct->m_nummode    = 2;

        m_octree->GetOctreeMesh(oct);

        ModuleSharedPtr mod = GetModuleFactory().CreateInstance(
            ModuleKey(eOutputModule, "xml"), oct);
        mod->RegisterConfig("outfile", fn + "_oct.xml");
        mod->ProcessVertices();
        mod->ProcessEdges();
        mod->ProcessFaces();
        mod->ProcessElements();
        mod->ProcessComposites();
        mod->Process();
    }

    m_mesh->m_expDim   = 3;
    m_mesh->m_spaceDim = 3;
    m_mesh->m_nummode = m_order + 1;
    if (m_makeBL)
    {
        m_mesh->m_numcomp = 2; // prisms and tets
    }
    else
    {
        m_mesh->m_numcomp = 1; // just tets
    }
    // m_mesh->m_nummode = 2;

    // create surface mesh
    m_mesh->m_expDim--; // just to make it easier to surface mesh for now
    SurfaceMeshSharedPtr m_surfacemesh =
        MemoryManager<SurfaceMesh>::AllocateSharedPtr(
            m_mesh, m_cad, m_octree, symsurfs, m_blthick);

    m_surfacemesh->Mesh();

    ProcessVertices();
    ProcessEdges();
    ProcessFaces();
    ProcessElements();
    ProcessComposites();

    m_surfacemesh->Report();

    m_mesh->m_expDim = 3;
    m_mesh->m_fields.push_back("u");
    m_mesh->m_fields.push_back("v");
    m_mesh->m_fields.push_back("w");
    m_mesh->m_fields.push_back("p");

    map<int, FaceSharedPtr> surftopriface;
    // map of surface element id to opposite prism
    // face for psudo surface in tetmesh

    TetMeshSharedPtr m_tet;
    if (m_makeBL)
    {
        BLMeshSharedPtr m_blmesh = MemoryManager<BLMesh>::AllocateSharedPtr(
            m_mesh, blsurfs, symsurfs, m_blthick);

        m_blmesh->Mesh();

        // create tet mesh
        m_tet = MemoryManager<TetMesh>::AllocateSharedPtr(
            m_mesh, m_octree, m_blmesh);
    }
    else
    {
        m_tet = MemoryManager<TetMesh>::AllocateSharedPtr(m_mesh, m_octree);
    }

    m_tet->Mesh();

    ClearElementLinks();
    ProcessVertices();
    ProcessEdges();
    ProcessFaces();
    ProcessElements();
    ProcessComposites();

    m_surfacemesh->HOSurf();

    if (m_mesh->m_verbose)
    {
        cout << endl;
        cout << m_mesh->m_element[3].size() << endl;
    }
}

Member Data Documentation

ModuleKey Nektar::Utilities::InputCAD::className

static

Initial value:

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

ModuleKey for class.

Definition at line 59 of file InputCAD.h.

NekDouble Nektar::Utilities::InputCAD::m_blthick

private

Definition at line 62 of file InputCAD.h.

Referenced by Process().

string Nektar::Utilities::InputCAD::m_CADName

private

Definition at line 64 of file InputCAD.h.

Referenced by Process().

NekDouble Nektar::Utilities::InputCAD::m_eps

private

Definition at line 62 of file InputCAD.h.

Referenced by Process().

bool Nektar::Utilities::InputCAD::m_makeBL

private

Definition at line 65 of file InputCAD.h.

Referenced by Process().

NekDouble Nektar::Utilities::InputCAD::m_maxDelta

private

Definition at line 62 of file InputCAD.h.

Referenced by Process().

NekDouble Nektar::Utilities::InputCAD::m_minDelta

private

Definition at line 62 of file InputCAD.h.

Referenced by Process().

int Nektar::Utilities::InputCAD::m_order

private

Definition at line 63 of file InputCAD.h.

Referenced by Process().

string Nektar::Utilities::InputCAD::m_udsName

private

Definition at line 64 of file InputCAD.h.

Referenced by Process().

bool Nektar::Utilities::InputCAD::m_writeoctree

private

Definition at line 65 of file InputCAD.h.

Referenced by Process().