Nektar::LibUtilities::PtsIO Class Reference

#include <PtsIO.h>

Public Member Functions
	PtsIO ()
void	Import (const string &inFile, LibUtilities::PtsFieldSharedPtr &ptsField)
	Import a pts field from file. More...
void	Write (const string &outFile, const LibUtilities::PtsFieldSharedPtr &ptsField)
	Save a pts field to a file. More...

Detailed Description

Definition at line 64 of file PtsIO.h.

Constructor & Destructor Documentation

Nektar::LibUtilities::PtsIO::PtsIO ( )

inline

Definition at line 68 of file PtsIO.h.

        {
        };

Member Function Documentation

void Nektar::LibUtilities::PtsIO::Import	(	const string &	inFile,
		LibUtilities::PtsFieldSharedPtr &	ptsField
	)

Import a pts field from file.

Parameters

inFile	filename of the file to read
ptsField	the resulting pts field.

Definition at line 64 of file PtsIO.cpp.

References Nektar::MemoryManager< DataType >::AllocateSharedPtr(), ASSERTL0, Nektar::ParseUtils::GenerateOrderedStringVector(), and npts.

Referenced by Nektar::LibUtilities::Import().

{
    TiXmlDocument docInput;
    ASSERTL0(docInput.LoadFile(inFile), "Unable to open file '" + inFile + "'.");

    TiXmlElement *nektar = docInput.FirstChildElement("NEKTAR");
    TiXmlElement *points = nektar->FirstChildElement("POINTS");
    int dim;
    int err = points->QueryIntAttribute("DIM", &dim);

    ASSERTL0(err == TIXML_SUCCESS, "Unable to read attribute DIM.");

    std::string fields = points->Attribute("FIELDS");

    vector<string> fieldNames;
    if (!fields.empty())
    {
        bool valid = ParseUtils::GenerateOrderedStringVector(
                         fields.c_str(), fieldNames);
        ASSERTL0(valid,
                 "Unable to process list of field variable in  FIELDS attribute:  " + fields);
    }

    int nfields = fieldNames.size();
    int totvars = dim + nfields;

    TiXmlNode *pointsBody = points->FirstChild();

    std::istringstream pointsDataStrm(pointsBody->ToText()->Value());

    vector<NekDouble> ptsSerial;
    Array<OneD,  Array<OneD,  NekDouble> > pts(totvars);

    try
    {
        NekDouble      ptsStream;
        while (!pointsDataStrm.fail())
        {
            pointsDataStrm >> ptsStream;

            ptsSerial.push_back(ptsStream);
        }
    }
    catch (...)
    {
        ASSERTL0(false, "Unable to read Points data.");
    }

    int npts = ptsSerial.size() / totvars;

    for (int i = 0; i < totvars; ++i)
    {
        pts[i] = Array<OneD, NekDouble>(npts);
    }

    for (int i = 0; i < npts; ++i)
    {
        for (int j = 0; j < totvars; ++j)
        {
            pts[j][i] = ptsSerial[i * totvars + j];
        }
    }

    ptsField = MemoryManager<PtsField>::AllocateSharedPtr(dim, fieldNames, pts);
}

void Nektar::LibUtilities::PtsIO::Write	(	const string &	outFile,
		const LibUtilities::PtsFieldSharedPtr &	ptsField
	)

Save a pts field to a file.

Parameters

inFile	filename of the file
ptsField	the pts field

Definition at line 137 of file PtsIO.cpp.

References ASSERTL0.

Referenced by Nektar::LibUtilities::Write().

{

    ASSERTL0(false, "Not implemented yet");
}