amroc/amr/Interfaces/SchemeGFMFileOutput.h

Go to the documentation of this file.

// -*- C++ -*-

// Copyright (C) 2012 Oak Ridge National Laboratory
// Ralf Deiterding, deiterdingr@ornl.gov

#ifndef SCHEME_GFM_FILEOUTPUT_H
#define SCHEME_GFM_FILEOUTPUT_H

#include "GFMFileOutput.h"

template <class SchemeType, class FixupType, class FlagType, int dim>
class SchemeGFMFileOutput : public SchemeOutBase<SchemeType,dim>,
  public GFMFileOutput<typename SchemeType::VectorType,FixupType,FlagType,dim> {
  typedef typename SchemeType::VectorType VectorType;
  typedef typename VectorType::InternalDataType DataType;
  typedef GFMFileOutput<VectorType,FixupType,FlagType,dim> base;
  typedef SchemeOutBase<SchemeType,dim> out_base;
public:
  typedef typename base::vec_grid_fct_type vec_grid_fct_type;  
  typedef typename base::grid_fct_type grid_fct_type;
  typedef typename base::grid_data_type grid_data_type;
  typedef typename base::gfm_solver_type gfm_solver_type;

  SchemeGFMFileOutput(gfm_solver_type& solver, SchemeType &lbm) : out_base(lbm), base(solver) {}

  virtual ~SchemeGFMFileOutput() {}
  
  virtual void WriteOut(vec_grid_fct_type& u, grid_fct_type& IOfunc) {
    for (int cnt=1; cnt<=base::Ncnt(); cnt++) { 
      if (base::CompName[cnt-1].c_str()[0] == '-') 
        continue;
      for (int lev=0; lev<=FineLevel(base::GH()); lev++) {
        int Time = CurrentTime(base::GH(),lev); 
        double t = GetPhysicalTime(u,Time,lev);
        this->Transform(u,IOfunc,Time,lev,cnt,t);
        base::WriteOut(IOfunc,base::CompName[cnt-1].c_str(),Time,lev,t);
      }
#ifndef DAGH_NO_MPI
      MPI_Barrier(comm_service::comm()); 
#endif
    }
  }

  virtual void WriteOut(grid_fct_type& IOfunc, const char* name) 
  { base::WriteOut(IOfunc,name); }  

  virtual void WriteOut(grid_fct_type& IOfunc, const char* name, 
                        const int& Time, const int& Level, const double& t) 
  { base::WriteOut(IOfunc,name,Time,Level,t); }

  virtual void WriteOut(grid_data_type& IOdata, const char* name, 
                        const int& Time, const int& Level, const double& t) {
    base::WriteOut(IOdata,name,Time,Level,t);
  }
};

#endif      

---