Tecplot output class, specifically for binary field output. More...

#include <OutputTecplot.h>

Inheritance diagram for Nektar::FieldUtils::OutputTecplotBinary:

Collaboration diagram for Nektar::FieldUtils::OutputTecplotBinary:

Public Member Functions
	OutputTecplotBinary (FieldSharedPtr f)

virtual	~OutputTecplotBinary ()

Public Member Functions inherited from Nektar::FieldUtils::OutputTecplot
	OutputTecplot (FieldSharedPtr f)

virtual	~OutputTecplot ()

virtual void	Process (po::variables_map &vm)
	Set up member variables to dump Tecplot format output. More...

Public Member Functions inherited from Nektar::FieldUtils::OutputModule
	OutputModule (FieldSharedPtr p_f)

FIELD_UTILS_EXPORT void	OpenStream ()
	Open a file for output. More...

Public Member Functions inherited from Nektar::FieldUtils::Module
FIELD_UTILS_EXPORT	Module (FieldSharedPtr p_f)

FIELD_UTILS_EXPORT void	RegisterConfig (string key, string value)
	Register a configuration option with a module. More...

FIELD_UTILS_EXPORT void	PrintConfig ()
	Print out all configuration options for a module. More...

FIELD_UTILS_EXPORT void	SetDefaults ()
	Sets default configuration options for those which have not been set. More...

FIELD_UTILS_EXPORT bool	GetRequireEquiSpaced (void)

FIELD_UTILS_EXPORT void	SetRequireEquiSpaced (bool pVal)

FIELD_UTILS_EXPORT void	EvaluateTriFieldAtEquiSpacedPts (LocalRegions::ExpansionSharedPtr &exp, const Array< OneD, const NekDouble > &infield, Array< OneD, NekDouble > &outfield)

Static Public Member Functions
static boost::shared_ptr< Module >	create (FieldSharedPtr f)
	Creates an instance of this class. More...

Static Public Member Functions inherited from Nektar::FieldUtils::OutputTecplot
static boost::shared_ptr< Module >	create (FieldSharedPtr f)
	Creates an instance of this class. More...

Static Public Attributes
static ModuleKey	m_className

Static Public Attributes inherited from Nektar::FieldUtils::OutputTecplot
static ModuleKey	m_className

Protected Member Functions
void	WriteDoubleOrFloat (std::ofstream &outfile, Array< OneD, NekDouble > &data)
	Write either double-precision or single-precision output of field data. More...

virtual void	WriteTecplotHeader (std::ofstream &outfile, std::vector< std::string > &var)
	Write Tecplot files header in binary format. More...

virtual void	WriteTecplotZone (std::ofstream &outfile)

virtual void	WriteTecplotConnectivity (std::ofstream &outfile)
	Write Tecplot connectivity information (ASCII) More...

Protected Member Functions inherited from Nektar::FieldUtils::OutputTecplot
int	GetNumTecplotBlocks ()
	Calculate number of Tecplot blocks. More...

void	CalculateConnectivity ()
	Calculate connectivity information for each expansion dimension. More...

virtual std::string	GetModuleName ()
	Returns this module's name. More...

Protected Member Functions inherited from Nektar::FieldUtils::Module
	Module ()

Additional Inherited Members
Protected Attributes inherited from Nektar::FieldUtils::OutputTecplot
bool	m_binary
	True if writing binary field output. More...

bool	m_oneOutputFile
	True if writing a single output file. More...

TecplotZoneType	m_zoneType
	Tecplot zone type of output. More...

vector< int >	m_numPoints
	Number of points per block in Tecplot file. More...

int	m_numBlocks
	Number of blocks in Tecplot file. More...

int	m_coordim
	Coordinate dimension of output. More...

int	m_totConn
	Total number of connectivity entries. More...

vector< Array< OneD, int > >	m_conn
	Connectivty for each block: one per element. More...

Array< OneD, int >	m_rankFieldSizes
	Each rank's field sizes. More...

Array< OneD, int >	m_rankConnSizes
	Each rank's connectivity sizes. More...

Array< OneD, Array< OneD, NekDouble > >	m_fields
	Field data to output. More...

Protected Attributes inherited from Nektar::FieldUtils::OutputModule
ofstream	m_fldFile
	Output stream. More...

Protected Attributes inherited from Nektar::FieldUtils::Module
FieldSharedPtr	m_f
	Field object. More...

map< string, ConfigOption >	m_config
	List of configuration values. More...

bool	m_requireEquiSpaced

Detailed Description

Tecplot output class, specifically for binary field output.

Definition at line 118 of file OutputTecplot.h.

Constructor & Destructor Documentation

Nektar::FieldUtils::OutputTecplotBinary::OutputTecplotBinary ( FieldSharedPtr f )

inline

Definition at line 128 of file OutputTecplot.h.

References Nektar::FieldUtils::OutputTecplot::m_binary, and Nektar::FieldUtils::Module::m_config.

                                           : OutputTecplot(f)
     {
         m_binary = true;
         m_config["double"] =
             ConfigOption(true, "0", "Write double-precision data: more "
                                     "accurate but more disk space required");
     }

virtual Nektar::FieldUtils::OutputTecplotBinary::~OutputTecplotBinary ( )

inlinevirtual

Definition at line 136 of file OutputTecplot.h.

137 {

138 }

Member Function Documentation

static boost::shared_ptr<Module> Nektar::FieldUtils::OutputTecplotBinary::create ( FieldSharedPtr f )

inlinestatic

Creates an instance of this class.

Definition at line 122 of file OutputTecplot.h.

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

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

void Nektar::FieldUtils::OutputTecplotBinary::WriteDoubleOrFloat	(	std::ofstream &	outfile,
		Array< OneD, NekDouble > &	data
	)

protected

Write either double-precision or single-precision output of field data.

Parameters

outfile	Output file name.
expansion	Expansion that is considered

Definition at line 616 of file OutputTecplot.cpp.

References CellMLToNektar.pycml::copy(), Nektar::FieldUtils::Module::m_config, and Nektar::FieldUtils::WriteStream().

Referenced by WriteTecplotZone().

 {
     // Data format: either double or single depending on user options
     bool useDoubles = m_config["double"].m_beenSet;
 
     if (useDoubles)
     {
         // For doubles, we can just write data.
         WriteStream(outfile, data);
     }
     else
     {
         // For single precision, needs typecast first.
         int nPts = data.num_elements();
         vector<float> tmp(data.num_elements());
         std::copy(&data[0], &data[0] + nPts, &tmp[0]);
         WriteStream(outfile, tmp);
     }
 }

void Nektar::FieldUtils::OutputTecplotBinary::WriteTecplotConnectivity ( std::ofstream & outfile )

protectedvirtual

Write Tecplot connectivity information (ASCII)

Parameters

outfile Output file

Reimplemented from Nektar::FieldUtils::OutputTecplot.

Definition at line 835 of file OutputTecplot.cpp.

References Nektar::FieldUtils::OutputTecplot::m_conn, Nektar::FieldUtils::Module::m_f, Nektar::FieldUtils::OutputTecplot::m_oneOutputFile, Nektar::FieldUtils::OutputTecplot::m_rankConnSizes, Nektar::FieldUtils::OutputTecplot::m_rankFieldSizes, Nektar::FieldUtils::OutputTecplot::m_totConn, Vmath::Vcopy(), and Nektar::FieldUtils::WriteStream().

 {
     if (m_oneOutputFile && m_f->m_comm->GetRank() > 0)
     {
         // Need to amalgamate connectivity information
         Array<OneD, int> conn(m_totConn);
         for (int i = 0, cnt = 0; i < m_conn.size(); ++i)
         {
             Vmath::Vcopy(m_conn[i].num_elements(), &m_conn[i][0], 1,
                          &conn[cnt], 1);
             cnt += m_conn[i].num_elements();
         }
         m_f->m_comm->Send(0, conn);
     }
     else
     {
         for (int i = 0; i < m_conn.size(); ++i)
         {
             WriteStream(outfile, m_conn[i]);
         }
 
         if (m_oneOutputFile && m_f->m_comm->GetRank() == 0)
         {
             int offset = m_rankFieldSizes[0];
 
             for (int n = 1; n < m_f->m_comm->GetSize(); ++n)
             {
                 Array<OneD, int> conn(m_rankConnSizes[n]);
                 m_f->m_comm->Recv(n, conn);
 
                 for (int j = 0; j < conn.num_elements(); ++j)
                 {
                     conn[j] += offset;
                 }
 
                 WriteStream(outfile, conn);
                 offset += m_rankFieldSizes[n];
             }
         }
     }
 }

void Nektar::FieldUtils::OutputTecplotBinary::WriteTecplotHeader	(	std::ofstream &	outfile,
		std::vector< std::string > &	var
	)

protectedvirtual

Write Tecplot files header in binary format.

Parameters

outfile	Output file name
var	Variables names

Reimplemented from Nektar::FieldUtils::OutputTecplot.

Definition at line 476 of file OutputTecplot.cpp.

References Nektar::FieldUtils::Module::m_f, Nektar::FieldUtils::OutputTecplot::m_oneOutputFile, CG_Iterations::title, and Nektar::FieldUtils::WriteStream().

 {
     if (m_oneOutputFile && m_f->m_comm->GetRank() > 0)
     {
         return;
     }
 
     // Version number
     outfile << "#!TDV112";
 
     // Int value of 1 for endian check
     WriteStream(outfile, 1);
 
     // We'll probably write a full solution field
     WriteStream(outfile, 0);
 
     // Title
     std::string title = "";
     WriteStream(outfile, title);
 
     // Number of variables
     WriteStream(outfile, (int)var.size());
 
     for (int i = 0; i < var.size(); ++i)
     {
         WriteStream(outfile, var[i]);
     }
 }

void Nektar::FieldUtils::OutputTecplotBinary::WriteTecplotZone ( std::ofstream & outfile )

protectedvirtual

Write Tecplot zone output in binary

Parameters

outfile	Output file name.
expansion	Expansion that is considered

Reimplemented from Nektar::FieldUtils::OutputTecplot.

Definition at line 643 of file OutputTecplot.cpp.

 {
     Array<OneD, NekDouble> fieldMin(m_fields.num_elements());
     Array<OneD, NekDouble> fieldMax(m_fields.num_elements());
 
     // Data format: either double or single depending on user options
     bool useDoubles = m_config["double"].m_beenSet;
 
     if ((m_oneOutputFile && m_f->m_comm->GetRank() == 0) || !m_oneOutputFile)
     {
         // Don't bother naming zone
         WriteStream(outfile, 299.0f); // Zone marker
 
         // Write same name as preplot
         int rank   = m_f->m_comm->GetRank();
         string zonename = "ZONE " + boost::lexical_cast<string>(rank);
         WriteStream(outfile, zonename);
 
         WriteStream(outfile, -1); // No parent zone
         WriteStream(outfile, -1); // No strand ID
         WriteStream(outfile, 0.0); // Solution time
         WriteStream(outfile, -1); // Unused, set to -1
 
         // Zone type: 1 = lineseg, 3 = quad, 5 = brick
         WriteStream(outfile, (int)m_zoneType);
 
         WriteStream(outfile, 0); // Data at nodes
         WriteStream(outfile, 0); // No 1-1 face neighbours
         WriteStream(outfile, 0); // No user-defined connections
 
         if (m_zoneType == eOrdered)
         {
             for (int i = 0; i < m_numPoints.size(); ++i)
             {
                 WriteStream(outfile, m_numPoints[i]);
             }
 
             for (int i = m_numPoints.size(); i < 3; ++i)
             {
                 WriteStream(outfile, 0);
             }
         }
         else
         {
             // Number of points in zone
             int nPoints = m_oneOutputFile ?
                 Vmath::Vsum(m_f->m_comm->GetSize(), m_rankFieldSizes, 1) :
                 m_fields[0].num_elements();
 
             WriteStream(outfile, nPoints); // Total number of points
             WriteStream(outfile, m_numBlocks); // Number of blocks
             WriteStream(outfile, 0); // Unused
             WriteStream(outfile, 0); // Unused
             WriteStream(outfile, 0); // Unused
         }
 
         WriteStream(outfile, 0); // No auxiliary data names
 
         // Finalise header
         WriteStream(outfile, 357.0f);
 
         // Now start to write data section so that we can dump geometry
         // information
 
         // Data marker
         WriteStream(outfile, 299.0f);
 
         for (int j = 0; j < m_fields.num_elements(); ++j)
         {
             WriteStream(outfile, useDoubles ? 2 : 1);
         }
 
         // No passive variables or variable sharing, no zone connectivity
         // sharing (we only dump one zone)
         WriteStream(outfile, 0);
         WriteStream(outfile, 0);
         WriteStream(outfile, -1);
     }
 
     for (int i = 0; i < m_fields.num_elements(); ++i)
     {
         fieldMin[i] = Vmath::Vmin(m_fields[i].num_elements(), m_fields[i], 1);
         fieldMax[i] = Vmath::Vmax(m_fields[i].num_elements(), m_fields[i], 1);
     }
 
     m_f->m_comm->AllReduce(fieldMin, LibUtilities::ReduceMin);
     m_f->m_comm->AllReduce(fieldMax, LibUtilities::ReduceMax);
 
     // Write out min/max of field data
     if ((m_oneOutputFile && m_f->m_comm->GetRank() == 0) || !m_oneOutputFile)
     {
         for (int i = 0; i < m_fields.num_elements(); ++i)
         {
             WriteStream(outfile, fieldMin[i]);
             WriteStream(outfile, fieldMax[i]);
         }
     }
 
     if (m_oneOutputFile && m_f->m_comm->GetRank() == 0)
     {
         for (int i = 0; i < m_fields.num_elements(); ++i)
         {
             WriteDoubleOrFloat(outfile, m_fields[i]);
 
             for (int n = 1; n < m_f->m_comm->GetSize(); ++n)
             {
                 Array<OneD, NekDouble> tmp(m_rankFieldSizes[n]);
                 m_f->m_comm->Recv(n, tmp);
                 WriteDoubleOrFloat(outfile, tmp);
             }
         }
     }
     else if (m_oneOutputFile && m_f->m_comm->GetRank() > 0)
     {
         for (int i = 0; i < m_fields.num_elements(); ++i)
         {
             m_f->m_comm->Send(0, m_fields[i]);
         }
     }
     else
     {
         for (int i = 0; i < m_fields.num_elements(); ++i)
         {
             WriteDoubleOrFloat(outfile, m_fields[i]);
         }
     }
 }

Member Data Documentation

ModuleKey Nektar::FieldUtils::OutputTecplotBinary::m_className

static

Initial value:

=
    GetModuleFactory().RegisterCreatorFunction(
        ModuleKey(eOutputModule, "plt"),
        OutputTecplotBinary::create,
        "Writes a Tecplot file in binary plt format.")

Definition at line 127 of file OutputTecplot.h.

Public Member Functions

Static Public Member Functions

Static Public Attributes

Protected Member Functions

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation