Nektar::Utilities::InputDat Class Reference

Input module for Xml files. More...

#include <InputDat.h>

Inheritance diagram for Nektar::Utilities::InputDat:

Collaboration diagram for Nektar::Utilities::InputDat:

Public Member Functions
	InputDat (FieldSharedPtr f)
	Set up InputDat object. More...
virtual	~InputDat ()
virtual void	Process (po::variables_map &vm)
virtual std::string	GetModuleName ()
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)
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)
virtual void	Process ()=0
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 (FieldSharedPtr f)
	Creates an instance of this class. More...

Static Public Attributes
static ModuleKey	m_className []
	ModuleKey for class. More...

Private Member Functions
void	ReadTecplotFEBlockZone (std::ifstream &datFile, string &line, Array< OneD, Array< OneD, NekDouble > > &pts, vector< Array< OneD, int > > &ptsConn)

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

Input module for Xml files.

Definition at line 47 of file InputDat.h.

Constructor & Destructor Documentation

Nektar::Utilities::InputDat::InputDat

(

FieldSharedPtr

f

)

Set up InputDat object.

Definition at line 63 of file InputDat.cpp.

References Nektar::Utilities::InputModule::m_allowedFiles.

                                   : InputModule(f)
{
    m_allowedFiles.insert("dat");
}

Nektar::Utilities::InputDat::~InputDat ( )

virtual

Definition at line 72 of file InputDat.cpp.

{
}

Member Function Documentation

static ModuleSharedPtr Nektar::Utilities::InputDat::create ( FieldSharedPtr  f )

inlinestatic

Creates an instance of this class.

Definition at line 55 of file InputDat.h.

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

        {
            return MemoryManager<InputDat>::AllocateSharedPtr(f);
        }

virtual std::string Nektar::Utilities::InputDat::GetModuleName ( )

inlinevirtual

Implements Nektar::Utilities::Module.

Definition at line 62 of file InputDat.h.

        {
            return "InputDat";
        }

void Nektar::Utilities::InputDat::Process ( po::variables_map &  vm )

virtual

Implements Nektar::Utilities::Module.

Definition at line 80 of file InputDat.cpp.

References Nektar::MemoryManager< DataType >::AllocateSharedPtr(), ASSERTL0, Nektar::LibUtilities::ePtsTriBlock, Nektar::ParseUtils::GenerateOrderedStringVector(), Nektar::Utilities::Module::m_f, and ReadTecplotFEBlockZone().

{

    if(m_f->m_verbose)
    {
        if(m_f->m_comm->TreatAsRankZero())
        {
            cout << "Processing input dat file" << endl;
        }
    }

    string      line, word, tag;
    std::ifstream datFile;
    stringstream s;

    // Open the file stream.
    string fname = m_f->m_inputfiles["dat"][0];


    datFile.open(fname.c_str());
    if (!datFile.good())
    {
        cerr << "Error opening file: " << fname << endl;
        abort();
    }

    // read variables
    // currently assum there are x y and z coordinates
    int dim = 3;
    vector<string> fieldNames;
    while (!datFile.eof())
    {
        getline(datFile, line);

        if(line.find("VARIABLES") != string::npos)
        {
            std::size_t  pos = line.find('=');
            pos++;

            // note this expects a comma separated list but
            // does not work for white space separated lists!
            bool valid = ParseUtils::GenerateOrderedStringVector(
                                line.substr(pos).c_str(), fieldNames);
            ASSERTL0(valid,"Unable to process list of field variable in "
                     " VARIABLES list:  "+ line.substr(pos));

            // remove coordinates from fieldNames
            fieldNames.erase(fieldNames.begin(), fieldNames.begin() + dim);

            break;
        }
    }

    // set up basic parameters
    int nfields = fieldNames.size();
    int totvars = dim + nfields;
    Array<OneD, Array<OneD, NekDouble> > pts(totvars);
    vector<Array<OneD, int> > ptsConn;

    // read zones
    while (!datFile.eof())
    {
        getline(datFile, line);

        if((line.find("ZONE") != string::npos)||
           (line.find("Zone") != string::npos)||
           (line.find("zone") != string::npos))
        {
            ReadTecplotFEBlockZone(datFile, line, pts, ptsConn);
        }
    }

    datFile.close();

    m_f->m_fieldPts =
        MemoryManager<LibUtilities::PtsField>::AllocateSharedPtr(
            dim, fieldNames, pts);
    m_f->m_fieldPts->SetPtsType(LibUtilities::ePtsTriBlock);
    m_f->m_fieldPts->SetConnectivity(ptsConn);
}

void Nektar::Utilities::InputDat::ReadTecplotFEBlockZone ( std::ifstream &  datFile, string &  line, Array< OneD, Array< OneD, NekDouble > > &  pts, vector< Array< OneD, int > > &  ptsConn  )

private

Definition at line 165 of file InputDat.cpp.

References ASSERTL0.

Referenced by Process().

{
    ASSERTL0(line.find("FEBlock") != string::npos,
             "Routine only set up for FEBLock format");
    ASSERTL0(line.find("ET") != string::npos,
             "Routine only set up TRIANLES");

    // read the number of nodes

    stringstream s;
    string tag;
    int start,end;

    s.clear();
    s.str(line);
    tag = s.str();

    // read the number of vertices
    start = tag.find("N=");
    end   = tag.find_first_of(',',start);
    int nvert = atoi(tag.substr(start+2,end).c_str());

    // read the number of elements
    start = tag.find("E=");
    end   = tag.find_first_of(',',start);
    int nelmt = atoi(tag.substr(start+2,end).c_str());

    // set-up or extend m_pts array;
    int norigpts  = pts[0].num_elements();
    int totfields = pts.num_elements();
    Array<OneD, Array<OneD, NekDouble> > origpts(totfields);
    for(int i = 0; i < totfields; ++i)
    {
        origpts[i] = pts[i];
        pts[i] = Array<OneD, NekDouble>(norigpts + nvert);
    }

    NekDouble value;
    for(int n = 0; n < totfields; ++n)
    {

        for(int i = 0; i < norigpts; ++i)
        {
            pts[n][i] = origpts[n][i];
        }
        for(int i = 0; i < nvert; ++i)
        {
            datFile >> value;
            pts[n][norigpts+i] = value;
        }
    }

    // read connectivity and add to list
    int intvalue;
    Array<OneD, int> conn(3*nelmt);
    for(int i = 0; i < 3*nelmt; ++i)
    {
        datFile >> intvalue;
        intvalue -=1; // decrement intvalue by 1 for c array convention
        conn[i] = norigpts + intvalue;
    }
    ptsConn.push_back(conn);

    getline(datFile, line);
}

Member Data Documentation

ModuleKey Nektar::Utilities::InputDat::m_className

static

Initial value:

= {
    GetModuleFactory().RegisterCreatorFunction(
            ModuleKey(eInputModule, "dat"),
            InputDat::create,
            "Reads Tecplot dat file for FE block triangular format."),
}

ModuleKey for class.

Definition at line 60 of file InputDat.h.