hdf5_diagnostics.hpp

Go to the documentation of this file.

#ifndef HDF5_DIAG
#define HDF5_DIAG 1

#include <H5Cpp.h>
#include <string>
#include "hdsim2d.hpp"
#include "../../misc/int2str.hpp"
#include "diagnostics.hpp"
#include "../../tessellation/Delaunay.hpp"

class Snapshot
{
public:

  Snapshot(void);

  Snapshot(const Snapshot& source);

  Snapshot& operator=(const Snapshot& other);

  vector<Vector2D> mesh_points;

  vector<ComputationalCell> cells;

  double time;

  int cycle;

#ifdef RICH_MPI
  vector<Vector2D> proc_points;
#endif

  TracerStickerNames tracerstickernames;
};

Snapshot read_hdf5_snapshot
(const string& fname,bool mpioverride=false);

class DiagnosticAppendix
{
public:

  virtual vector<double> operator()(const hdsim& sim) const = 0;

  virtual string getName(void) const = 0;

  virtual ~DiagnosticAppendix(void);
};

void write_snapshot_to_hdf5(hdsim const& sim,string const& fname,
                            const vector<DiagnosticAppendix*>& appendices=vector<DiagnosticAppendix*>());
void read_hdf5_snapshot(ResetDump &dump,string const& fname,EquationOfState
        const* eos);
void ConvertHDF5toBinary(string const& input, string const& output);

void WriteDelaunay(Delaunay const& tri, string const& filename);

Snapshot ReDistributeData(string const& filename, Tessellation const& proctess, size_t snapshot_number);

Snapshot ReDistributeData2(string const& filename, Tessellation const& proctess, size_t snapshot_number,
        bool mpioverride = false);


void WriteTess(Tessellation const& tess, string const& filename);
#endif // HDF5_DIAG