amroc/amr/Interfaces/SchemeInitialCondition.h

Go to the documentation of this file.

// -*- C++ -*-

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

#ifndef SCHEME_INITIAL_CONDITION_H
#define SCHEME_INITIAL_CONDITION_H

#include "InitialCondition.h"

template <class SchemeType, int dim>
class SchemeInitialCondition : public InitialCondition<typename SchemeType::VectorType,dim> {
  typedef InitialCondition<typename SchemeType::VectorType,dim> base;
public:    
  typedef typename SchemeType::VectorType VectorType;
  typedef typename VectorType::InternalDataType DataType;
  typedef typename base::vec_grid_data_type vec_grid_data_type;
  typedef typename base::grid_data_type grid_data_type;

  SchemeInitialCondition(SchemeType &scheme) : base(), _scheme(scheme), Type(0), Scaling(0), NAux(0) 
  { aux = new DataType[Scheme().NEquations()]; }
    
  ~SchemeInitialCondition() { delete [] aux; }

  virtual void register_at(ControlDevice& Ctrl, const std::string& prefix) {
    base::LocCtrl = Ctrl.getSubDevice(prefix+"InitialCondition");
    RegisterAt(base::LocCtrl,"Type",Type);
    RegisterAt(base::LocCtrl,"Scaling",Scaling);
    RegisterAt(base::LocCtrl,"NAux",NAux);
    char VariableName[16];
    for (int i=0; i<Scheme().NEquations(); i++) {
      aux[i]=DataType(0.);
      std::sprintf(VariableName,"Aux(%d)",i+1);
      RegisterAt(base::LocCtrl,VariableName,aux[i]);
    }
  }

  virtual void SetGrid(vec_grid_data_type& fvec, grid_data_type& workvec, const int& level) 
  { Scheme().ICStandard(fvec,Type,aux,NAux,Scaling); }

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

protected:
  SchemeType &_scheme;
  int Type, Scaling, NAux;
  DataType *aux;
};
  

#endif

---