Output3D.h

// =======================================================================
// %W% %G%
// 
// Class:       TOutput3D
// Purpose:     store given data and realize output
//
// Author:      Gunar Matthies (24.07.2000)
//
// History:     start of implementation 24.07.2000 (Gunar Matthies)
//
// =======================================================================
#ifdef _MPI
#  include "mpi.h"
#endif


#ifndef __OUTPUT3D__
#define __OUTPUT3D__

#include <FEVectFunct3D.h>
#include <Domain.h>

class TOutput3D
{
  protected:
    TCollection *Coll;

    int N_FESpaces;

    int MaxN_FESpaces;

    TFESpace3D **FESpaceArray;

    int N_ScalarVar;

    int MaxN_ScalarVar;

    TFEFunction3D **FEFunctionArray;

    int N_VectorVar;

    int MaxN_VectorVar;

    TFEVectFunct3D **FEVectFunctArray;

    int N_Parameters;

    int MaxN_Parameters;

    double *ParameterValues;

    const char **ParameterDescription;

    TDomain *Domain;

    int AddFESpace(TFESpace3D *fespace);

    struct TOutputData
    {
      int N_Nodes;
      int N_Data;
      TVertex **Nodes;
      int *ConList;

      double **FEFuncValues;

      enum CELLTYPE {TETRAHEDRON=4, BRICK} Type;

      TOutputData() : Nodes(0), ConList(0), N_Nodes(0),
              N_Data(0), FEFuncValues(0) {};
      ~TOutputData();
    };

    TOutputData *Data;

    char *Name;
  
  protected:
    void ComputeOutputData();
    void ComputeFEValues();

  public:
    TOutput3D(int maxn_fespaces, int maxn_scalar, int maxn_vect, 
              int maxn_parameters, TDomain *domain, TCollection *coll=NULL,
          const char *name=NULL);

    ~TOutput3D();

    int AddFEFunction(TFEFunction3D *fefunction);

    int AddFEVectFunct(TFEVectFunct3D *fevectfunct);

    int AddParameter(double value, const char *descr);

    int Write(std::string basename, int i=1, double t=0.);
    
    int WriteGrape(const char *name);

    int WriteTecplot(const char *name);

    int WriteGMV(const char *name);

    int WriteAmira(const char *name);

    int WriteVtk(const char *name);
    
    void WriteVtkDiscontinuous(const char *fileName, 
                               int N_LocVertices, TVertex **Vertices);

    int Write_ParVTK(
#ifdef _MPI
                                MPI_Comm comm,
#endif
                               int img, char *subID);

    int WriteBinaryPlt(const char *filename);

};

#endif