Output2D.h

// =======================================================================
// @(#)Output2D.h        1.3 11/15/99
// 
// Class:       TOutput2D
// Purpose:     store given data and realize output
//
// Author:      Gunar Matthies (21.08.1998)
//
// History:     start of implementation 21.08.1998 (Gunar Matthies)
//        :     parallel vtk output 26.09.09 (Sashikumaar Ganesan)
// =======================================================================

#ifndef __OUTPUT2D__
#define __OUTPUT2D__

#include <FEVectFunct2D.h>
#include <Domain.h>

#ifdef _MPI
#  include "mpi.h"
#endif

class TOutput2D
{
  protected:
    TCollection *Coll;

    int N_FESpaces;

    int MaxN_FESpaces;

    TFESpace2D **FESpaceArray;

    int N_ScalarVar;

    int MaxN_ScalarVar;

    TFEFunction2D **FEFunctionArray;

    int N_VectorVar;

    int MaxN_VectorVar;

    TFEVectFunct2D **FEVectFunctArray;

    int N_Parameters;

    int MaxN_Parameters;

    double *ParameterValues;

    const char **ParameterDescription;

    TDomain *Domain;

    int AddFESpace(TFESpace2D *fespace);

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

      double **FEFuncValues; // first two for coords

      enum CELLTYPE {TRIA=2, QUAD} Type;

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

    TOutputData *Data;

  protected:
    void ComputeOutputData();
    void ComputeFEValues();

  public:
    TOutput2D(int maxn_fespaces, int maxn_scalar, int maxn_vect, 
              int maxn_parameters, TDomain *domain);

    ~TOutput2D();

    int AddFEFunction(TFEFunction2D *fefunction);

    int AddFEVectFunct(TFEVectFunct2D *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 WriteGnuplot(const char *name);

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

    int WriteMatlab(const char *name);

    int WriteMatlabOld(const char *name);

    int WriteGNU_iso(const char *name, int scalar);

    int WriteGMV(const char *name);    

    int Write_ParVTK(
#ifdef _MPI
                          MPI_Comm comm,
#endif
                          int img, char *subID);
    int WriteAsciiPlt(const char *filename);
    int WriteBinaryPlt(const char *filename);
};

#endif