Nektar++
Public Member Functions | Static Public Member Functions | Static Public Attributes | Protected Member Functions | List of all members
Nektar::FieldUtils::OutputTecplotBinary Class Reference

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

#include <OutputTecplot.h>

Inheritance diagram for Nektar::FieldUtils::OutputTecplotBinary:
[legend]

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)
 Write fld to output file. More...
 
virtual std::string GetModuleName ()
 
- Public Member Functions inherited from Nektar::FieldUtils::OutputFileBase
 OutputFileBase (FieldSharedPtr f)
 
virtual ~OutputFileBase ()
 
virtual std::string GetModuleDescription ()
 
virtual ModulePriority GetModulePriority ()
 
- 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)
 
virtual ~Module ()=default
 
const ConfigOptionGetConfigOption (const std::string &key) const
 
FIELD_UTILS_EXPORT void RegisterConfig (std::string key, std::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 void AddFile (std::string fileType, std::string fileName)
 
FIELD_UTILS_EXPORT void EvaluateTriFieldAtEquiSpacedPts (LocalRegions::ExpansionSharedPtr &exp, const Array< OneD, const NekDouble > &infield, Array< OneD, NekDouble > &outfield)
 

Static Public Member Functions

static std::shared_ptr< Modulecreate (FieldSharedPtr f)
 Creates an instance of this class. More...
 
- Static Public Member Functions inherited from Nektar::FieldUtils::OutputTecplot
static std::shared_ptr< Modulecreate (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
virtual void OutputFromPts (po::variables_map &vm)
 Write from pts to output file. More...
 
virtual void OutputFromExp (po::variables_map &vm)
 Write from m_exp to output file. More...
 
virtual void OutputFromData (po::variables_map &vm)
 Write from data to output file. More...
 
virtual fs::path GetPath (std::string &filename, po::variables_map &vm)
 
virtual fs::path GetFullOutName (std::string &filename, po::variables_map &vm)
 
void WriteTecplotFile (po::variables_map &vm)
 
int GetNumTecplotBlocks ()
 Calculate number of Tecplot blocks. More...
 
void CalculateConnectivity ()
 Calculate connectivity information for each expansion dimension. More...
 
- Protected Member Functions inherited from Nektar::FieldUtils::Module
 Module ()
 

Additional Inherited Members

- Public Attributes inherited from Nektar::FieldUtils::Module
FieldSharedPtr m_f
 Field object. More...
 
- 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...
 
bool m_writeHeader
 True if writing header. More...
 
TecplotZoneType m_zoneType
 Tecplot zone type of output. More...
 
std::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...
 
std::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::OutputFileBase
bool m_requireEquiSpaced
 
- Protected Attributes inherited from Nektar::FieldUtils::OutputModule
std::ofstream m_fldFile
 Output stream. More...
 
- Protected Attributes inherited from Nektar::FieldUtils::Module
std::map< std::string, ConfigOptionm_config
 List of configuration values. More...
 
std::set< std::string > m_allowedFiles
 List of allowed file formats. More...
 

Detailed Description

Tecplot output class, specifically for binary field output.

Definition at line 137 of file OutputTecplot.h.

Constructor & Destructor Documentation

◆ OutputTecplotBinary()

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

Definition at line 147 of file OutputTecplot.h.

147  : OutputTecplot(f)
148  {
149  m_binary = true;
150  }
bool m_binary
True if writing binary field output.
Definition: OutputTecplot.h:98

References Nektar::FieldUtils::OutputTecplot::m_binary.

◆ ~OutputTecplotBinary()

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

Definition at line 152 of file OutputTecplot.h.

153  {
154  }

Member Function Documentation

◆ create()

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

Creates an instance of this class.

Definition at line 141 of file OutputTecplot.h.

142  {
144  }
static std::shared_ptr< DataType > AllocateSharedPtr(const Args &...args)
Allocate a shared pointer from the memory pool.

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

◆ WriteDoubleOrFloat()

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
outfileOutput file name.
expansionExpansion that is considered

Definition at line 675 of file OutputTecplot.cpp.

677 {
678  // Data format: either double or single depending on user options
679  bool useDoubles = m_config["double"].as<bool>();
680 
681  if (useDoubles)
682  {
683  // For doubles, we can just write data.
684  WriteStream(outfile, data);
685  }
686  else
687  {
688  // For single precision, needs typecast first.
689  int nPts = data.size();
690  std::vector<float> tmp(data.size());
691  std::copy(&data[0], &data[0] + nPts, &tmp[0]);
692  WriteStream(outfile, tmp);
693  }
694 }
std::map< std::string, ConfigOption > m_config
List of configuration values.
Definition: Module.h:233
def copy(self)
Definition: pycml.py:2663
void WriteStream(std::ostream &outfile, T data)
Helper function to write binary data to stream.

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

Referenced by WriteTecplotZone().

◆ WriteTecplotConnectivity()

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

Write Tecplot connectivity information (ASCII)

Parameters
outfileOutput file

Reimplemented from Nektar::FieldUtils::OutputTecplot.

Definition at line 903 of file OutputTecplot.cpp.

904 {
905  if (m_oneOutputFile && m_f->m_comm->GetRank() > 0)
906  {
907  // Need to amalgamate connectivity information
908  Array<OneD, int> conn(m_totConn);
909  for (int i = 0, cnt = 0; i < m_conn.size(); ++i)
910  {
911  Vmath::Vcopy(m_conn[i].size(), &m_conn[i][0], 1,
912  &conn[cnt], 1);
913  cnt += m_conn[i].size();
914  }
915  m_f->m_comm->Send(0, conn);
916  }
917  else
918  {
919  for (int i = 0; i < m_conn.size(); ++i)
920  {
921  WriteStream(outfile, m_conn[i]);
922  }
923 
924  if (m_oneOutputFile && m_f->m_comm->GetRank() == 0)
925  {
926  int offset = m_rankFieldSizes[0];
927 
928  for (int n = 1; n < m_f->m_comm->GetSize(); ++n)
929  {
930  Array<OneD, int> conn(m_rankConnSizes[n]);
931  m_f->m_comm->Recv(n, conn);
932 
933  for (int j = 0; j < conn.size(); ++j)
934  {
935  conn[j] += offset;
936  }
937 
938  WriteStream(outfile, conn);
939  offset += m_rankFieldSizes[n];
940  }
941  }
942  }
943 }
FieldSharedPtr m_f
Field object.
Definition: Module.h:230
int m_totConn
Total number of connectivity entries.
std::vector< Array< OneD, int > > m_conn
Connectivty for each block: one per element.
Array< OneD, int > m_rankConnSizes
Each rank's connectivity sizes.
bool m_oneOutputFile
True if writing a single output file.
Array< OneD, int > m_rankFieldSizes
Each rank's field sizes.
void Vcopy(int n, const T *x, const int incx, T *y, const int incy)
Definition: Vmath.cpp:1199

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

◆ WriteTecplotHeader()

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

Write Tecplot files header in binary format.

Parameters
outfileOutput file name
varVariables names

Reimplemented from Nektar::FieldUtils::OutputTecplot.

Definition at line 477 of file OutputTecplot.cpp.

479 {
480  if (m_oneOutputFile && m_f->m_comm->GetRank() > 0)
481  {
482  return;
483  }
484 
485  // Version number
486  outfile << "#!TDV112";
487 
488  // Int value of 1 for endian check
489  WriteStream(outfile, 1);
490 
491  // We'll probably write a full solution field
492  WriteStream(outfile, 0);
493 
494  // Title
495  std::string title = "";
496  WriteStream(outfile, title);
497 
498  // Number of variables
499  WriteStream(outfile, (int)var.size());
500 
501  for (int i = 0; i < var.size(); ++i)
502  {
503  WriteStream(outfile, var[i]);
504  }
505 }

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

◆ WriteTecplotZone()

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

Write Tecplot zone output in binary

Parameters
outfileOutput file name.
expansionExpansion that is considered

Reimplemented from Nektar::FieldUtils::OutputTecplot.

Definition at line 702 of file OutputTecplot.cpp.

703 {
704  Array<OneD, NekDouble> fieldMin(m_fields.size());
705  Array<OneD, NekDouble> fieldMax(m_fields.size());
706 
707  // Data format: either double or single depending on user options
708  bool useDoubles = m_config["double"].as<bool>();
709 
710  if ((m_oneOutputFile && m_f->m_comm->GetRank() == 0) || !m_oneOutputFile)
711  {
712  // Don't bother naming zone
713  WriteStream(outfile, 299.0f); // Zone marker
714 
715  // Write same name as preplot
716  int rank = m_f->m_comm->GetRank();
717  string zonename = "ZONE " + boost::lexical_cast<string>(rank);
718  WriteStream(outfile, zonename);
719 
720  WriteStream(outfile, -1); // No parent zone
721  WriteStream(outfile, -1); // No strand ID
722  WriteStream(outfile, 0.0); // Solution time
723  WriteStream(outfile, -1); // Unused, set to -1
724 
725  // Zone type: 1 = lineseg, 3 = quad, 5 = brick
726  WriteStream(outfile, (int)m_zoneType);
727 
728  WriteStream(outfile, 0); // Data at nodes
729  WriteStream(outfile, 0); // No 1-1 face neighbours
730  WriteStream(outfile, 0); // No user-defined connections
731 
732  if (m_zoneType == eOrdered)
733  {
734  for (int i = 0; i < m_numPoints.size(); ++i)
735  {
736  WriteStream(outfile, m_numPoints[i]);
737  }
738 
739  for (int i = m_numPoints.size(); i < 3; ++i)
740  {
741  WriteStream(outfile, 0);
742  }
743  }
744  else
745  {
746  // Number of points in zone
747  int nPoints = m_oneOutputFile ?
748  Vmath::Vsum(m_f->m_comm->GetSize(), m_rankFieldSizes, 1) :
749  m_fields[0].size();
750 
751  WriteStream(outfile, nPoints); // Total number of points
752  WriteStream(outfile, m_numBlocks); // Number of blocks
753  WriteStream(outfile, 0); // Unused
754  WriteStream(outfile, 0); // Unused
755  WriteStream(outfile, 0); // Unused
756  }
757 
758  WriteStream(outfile, 0); // No auxiliary data names
759 
760  // Finalise header
761  WriteStream(outfile, 357.0f);
762 
763  // Now start to write data section so that we can dump geometry
764  // information
765 
766  // Data marker
767  WriteStream(outfile, 299.0f);
768 
769  for (int j = 0; j < m_fields.size(); ++j)
770  {
771  WriteStream(outfile, useDoubles ? 2 : 1);
772  }
773 
774  // No passive variables or variable sharing, no zone connectivity
775  // sharing (we only dump one zone)
776  WriteStream(outfile, 0);
777  WriteStream(outfile, 0);
778  WriteStream(outfile, -1);
779  }
780 
781  for (int i = 0; i < m_fields.size(); ++i)
782  {
783  fieldMin[i] = Vmath::Vmin(m_fields[i].size(), m_fields[i], 1);
784  fieldMax[i] = Vmath::Vmax(m_fields[i].size(), m_fields[i], 1);
785  }
786 
787  m_f->m_comm->AllReduce(fieldMin, LibUtilities::ReduceMin);
788  m_f->m_comm->AllReduce(fieldMax, LibUtilities::ReduceMax);
789 
790  // Write out min/max of field data
791  if ((m_oneOutputFile && m_f->m_comm->GetRank() == 0) || !m_oneOutputFile)
792  {
793  for (int i = 0; i < m_fields.size(); ++i)
794  {
795  WriteStream(outfile, fieldMin[i]);
796  WriteStream(outfile, fieldMax[i]);
797  }
798  }
799 
800  if (m_oneOutputFile && m_f->m_comm->GetRank() == 0)
801  {
802  for (int i = 0; i < m_fields.size(); ++i)
803  {
804  WriteDoubleOrFloat(outfile, m_fields[i]);
805 
806  for (int n = 1; n < m_f->m_comm->GetSize(); ++n)
807  {
808  Array<OneD, NekDouble> tmp(m_rankFieldSizes[n]);
809  m_f->m_comm->Recv(n, tmp);
810  WriteDoubleOrFloat(outfile, tmp);
811  }
812  }
813  }
814  else if (m_oneOutputFile && m_f->m_comm->GetRank() > 0)
815  {
816  for (int i = 0; i < m_fields.size(); ++i)
817  {
818  m_f->m_comm->Send(0, m_fields[i]);
819  }
820  }
821  else
822  {
823  for (int i = 0; i < m_fields.size(); ++i)
824  {
825  WriteDoubleOrFloat(outfile, m_fields[i]);
826  }
827  }
828 }
void WriteDoubleOrFloat(std::ofstream &outfile, Array< OneD, NekDouble > &data)
Write either double-precision or single-precision output of field data.
TecplotZoneType m_zoneType
Tecplot zone type of output.
std::vector< int > m_numPoints
Number of points per block in Tecplot file.
int m_numBlocks
Number of blocks in Tecplot file.
Array< OneD, Array< OneD, NekDouble > > m_fields
Field data to output.
T Vmin(int n, const T *x, const int incx)
Return the minimum element in x - called vmin to avoid conflict with min.
Definition: Vmath.cpp:992
T Vsum(int n, const T *x, const int incx)
Subtract return sum(x)
Definition: Vmath.cpp:846
T Vmax(int n, const T *x, const int incx)
Return the maximum element in x – called vmax to avoid conflict with max.
Definition: Vmath.cpp:892

References Nektar::FieldUtils::eOrdered, Nektar::FieldUtils::Module::m_config, Nektar::FieldUtils::Module::m_f, Nektar::FieldUtils::OutputTecplot::m_fields, Nektar::FieldUtils::OutputTecplot::m_numBlocks, Nektar::FieldUtils::OutputTecplot::m_numPoints, Nektar::FieldUtils::OutputTecplot::m_oneOutputFile, Nektar::FieldUtils::OutputTecplot::m_rankFieldSizes, Nektar::FieldUtils::OutputTecplot::m_zoneType, Nektar::LibUtilities::ReduceMax, Nektar::LibUtilities::ReduceMin, Vmath::Vmax(), Vmath::Vmin(), Vmath::Vsum(), WriteDoubleOrFloat(), and Nektar::FieldUtils::WriteStream().

Member Data Documentation

◆ m_className

ModuleKey Nektar::FieldUtils::OutputTecplotBinary::m_className
static
Initial value:
=
"Writes a Tecplot file in binary plt format.")
static std::shared_ptr< Module > create(FieldSharedPtr f)
Creates an instance of this class.
tKey RegisterCreatorFunction(tKey idKey, CreatorFunction classCreator, std::string pDesc="")
Register a class with the factory.
Definition: NekFactory.hpp:200
std::pair< ModuleType, std::string > ModuleKey
Definition: Module.h:290
ModuleFactory & GetModuleFactory()
Definition: Module.cpp:49

Definition at line 146 of file OutputTecplot.h.