amroc/amr/Interfaces/SchemeFileInput.h

Go to the documentation of this file.

// -*- C++ -*-

// Copyright (C) 2013 Ralf Deiterding
// German Aerospace Center (DLR)

#ifndef SCHEME_FILEINPUT_H
#define SCHEME_FILEINPUT_H

#include "FileInput.h"
#include "Interfaces/SchemeFileOutput.h"

template <class SchemeType, int dim>
class SchemeInBase {
  typedef typename SchemeType::VectorType VectorType;
  typedef typename VectorType::InternalDataType DataType;
  typedef GridFunction<VectorType,dim> vec_grid_fct_type;  
  typedef GridFunction<DataType,dim> grid_fct_type;  

public:
  SchemeInBase(SchemeType &scheme) : _scheme(scheme) {}

  virtual ~SchemeInBase() {}

  virtual void Transform(grid_fct_type& work, vec_grid_fct_type& u, const int Time, const int& Level, 
                         const int cnt, const double& t) {
    forall(u,Time,Level,c)
      Scheme().Input(u(Time,Level,c),work(Time,Level,c),cnt);
    end_forall 
  }

  inline SchemeType& Scheme() { return _scheme; }
  inline const SchemeType& Scheme() const { return _scheme; }

protected:
  SchemeType &_scheme;
};


template <class SchemeType, int dim>
class SchemeFileInput : public SchemeOutBase<SchemeType,dim>, public SchemeInBase<SchemeType,dim>,
                        public FileInput<typename SchemeType::VectorType,dim> {
  typedef typename SchemeType::VectorType VectorType;
  typedef typename VectorType::InternalDataType DataType;
  typedef FileInput<VectorType,dim> base;
  typedef SchemeOutBase<SchemeType,dim> out_base;
  typedef SchemeInBase<SchemeType,dim> in_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;

  SchemeFileInput(SchemeType &scheme) : out_base(scheme), in_base(scheme), base(), _scheme(scheme) {}

  virtual ~SchemeFileInput() {}
  
  virtual void ReadIn(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;
      int lev;
      double t;
      for (lev=0; lev<=FineLevel(base::GH()); lev++) {
        int Time = CurrentTime(base::GH(),lev); 
        t = GetPhysicalTime(u,Time,lev);
        ((out_base*) this)->Transform(u,IOfunc,Time,lev,cnt,t);
      }
      ReadIn(IOfunc,base::CompName[cnt-1].c_str());
      for (lev=0; lev<=FineLevel(base::GH()); lev++) {
        int Time = CurrentTime(base::GH(),lev); 
        t = GetPhysicalTime(u,Time,lev);
        ((in_base*) this)->Transform(IOfunc,u,Time,lev,cnt,t);
      }
#ifndef DAGH_NO_MPI
      MPI_Barrier(comm_service::comm()); 
#endif
    } 
  }   

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

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

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

  inline SchemeType& Scheme() { return _scheme; }
  inline const SchemeType& Scheme() const { return _scheme; }

protected:
  SchemeType &_scheme;
};

#endif      

---