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 ()
 
- Public Member Functions inherited from Nektar::FieldUtils::OutputFileBase
 OutputFileBase (FieldSharedPtr f)
 
virtual ~OutputFileBase ()
 
- 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
 
void Process (po::variables_map &vm)
 
std::string GetModuleName ()
 
std::string GetModuleDescription ()
 
const ConfigOptionGetConfigOption (const std::string &key) const
 
ModulePriority GetModulePriority ()
 
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 v_WriteTecplotHeader (std::ofstream &outfile, std::vector< std::string > &var) override
 Write Tecplot files header in binary format. More...
 
virtual void v_WriteTecplotZone (std::ofstream &outfile) override
 
virtual void v_WriteTecplotConnectivity (std::ofstream &outfile) override
 Write Tecplot connectivity information (ASCII) More...
 
- Protected Member Functions inherited from Nektar::FieldUtils::OutputTecplot
virtual void v_Process (po::variables_map &vm) override
 Write fld to output file. More...
 
virtual std::string v_GetModuleName () override
 
virtual void v_OutputFromPts (po::variables_map &vm) override
 Write from pts to output file. More...
 
virtual void v_OutputFromExp (po::variables_map &vm) override
 Write from m_exp to output file. More...
 
virtual void v_OutputFromData (po::variables_map &vm) override
 Write from data to output file. More...
 
virtual fs::path v_GetPath (std::string &filename, po::variables_map &vm) override
 
virtual fs::path v_GetFullOutName (std::string &filename, po::variables_map &vm) override
 
void WriteTecplotHeader (std::ofstream &outfile, std::vector< std::string > &var)
 
void WriteTecplotZone (std::ofstream &outfile)
 
void WriteTecplotConnectivity (std::ofstream &outfile)
 
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::OutputFileBase
virtual std::string v_GetModuleDescription () override
 
virtual ModulePriority v_GetModulePriority () override
 
fs::path GetPath (std::string &filename, po::variables_map &vm)
 
fs::path GetFullOutName (std::string &filename, po::variables_map &vm)
 
- 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 149 of file OutputTecplot.h.

Constructor & Destructor Documentation

◆ OutputTecplotBinary()

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

Definition at line 159 of file OutputTecplot.h.

159  : OutputTecplot(f)
160  {
161  m_binary = true;
162  }
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 164 of file OutputTecplot.h.

165  {
166  }

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 153 of file OutputTecplot.h.

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

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

◆ v_WriteTecplotConnectivity()

void Nektar::FieldUtils::OutputTecplotBinary::v_WriteTecplotConnectivity ( std::ofstream &  outfile)
overrideprotectedvirtual

Write Tecplot connectivity information (ASCII)

Parameters
outfileOutput file

Reimplemented from Nektar::FieldUtils::OutputTecplot.

Definition at line 905 of file OutputTecplot.cpp.

906 {
907  if (m_oneOutputFile && m_f->m_comm->GetSpaceComm()->GetRank() > 0)
908  {
909  // Need to amalgamate connectivity information
910  Array<OneD, int> conn(m_totConn);
911  for (int i = 0, cnt = 0; i < m_conn.size(); ++i)
912  {
913  Vmath::Vcopy(m_conn[i].size(), &m_conn[i][0], 1, &conn[cnt], 1);
914  cnt += m_conn[i].size();
915  }
916  m_f->m_comm->GetSpaceComm()->Send(0, conn);
917  }
918  else
919  {
920  for (int i = 0; i < m_conn.size(); ++i)
921  {
922  WriteStream(outfile, m_conn[i]);
923  }
924 
925  if (m_oneOutputFile && m_f->m_comm->GetSpaceComm()->GetRank() == 0)
926  {
927  int offset = m_rankFieldSizes[0];
928 
929  for (int n = 1; n < m_f->m_comm->GetSpaceComm()->GetSize(); ++n)
930  {
931  Array<OneD, int> conn(m_rankConnSizes[n]);
932  m_f->m_comm->GetSpaceComm()->Recv(n, conn);
933 
934  for (int j = 0; j < conn.size(); ++j)
935  {
936  conn[j] += offset;
937  }
938 
939  WriteStream(outfile, conn);
940  offset += m_rankFieldSizes[n];
941  }
942  }
943  }
944 }
FieldSharedPtr m_f
Field object.
Definition: Module.h:234
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 WriteStream(std::ostream &outfile, T data)
Helper function to write binary data to stream.
void Vcopy(int n, const T *x, const int incx, T *y, const int incy)
Definition: Vmath.cpp:1255

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

◆ v_WriteTecplotHeader()

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

Write Tecplot files header in binary format.

Parameters
outfileOutput file name
varVariables names

Reimplemented from Nektar::FieldUtils::OutputTecplot.

Definition at line 470 of file OutputTecplot.cpp.

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

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

◆ v_WriteTecplotZone()

void Nektar::FieldUtils::OutputTecplotBinary::v_WriteTecplotZone ( std::ofstream &  outfile)
overrideprotectedvirtual

Write Tecplot zone output in binary

Parameters
outfileOutput file name.
expansionExpansion that is considered

Reimplemented from Nektar::FieldUtils::OutputTecplot.

Definition at line 700 of file OutputTecplot.cpp.

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

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

◆ 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 673 of file OutputTecplot.cpp.

675 {
676  // Data format: either double or single depending on user options
677  bool useDoubles = m_config["double"].as<bool>();
678 
679  if (useDoubles)
680  {
681  // For doubles, we can just write data.
682  WriteStream(outfile, data);
683  }
684  else
685  {
686  // For single precision, needs typecast first.
687  int nPts = data.size();
688  std::vector<float> tmp(data.size());
689  std::copy(&data[0], &data[0] + nPts, &tmp[0]);
690  WriteStream(outfile, tmp);
691  }
692 }
def copy(self)
Definition: pycml.py:2663

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

Referenced by v_WriteTecplotZone().

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:198
std::pair< ModuleType, std::string > ModuleKey
Definition: Module.h:317
ModuleFactory & GetModuleFactory()
Definition: Module.cpp:49

Definition at line 158 of file OutputTecplot.h.