amroc/amr/Interfaces/SchemeBoundaryConditions.h

Go to the documentation of this file.

// -*- C++ -*-

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

#ifndef SCHEME_BOUNDARY_CONDITIONS_H
#define SCHEME_BOUNDARY_CONDITIONS_H

#include "BoundaryConditions.h"

template <class DataType>
class BCSide : public controlable {
public:
  BCSide(const int maxaux) : Type(-1), Scaling(0), NAux(0), MaxAux(maxaux) 
  { aux = new DataType[MaxAux]; }

  ~BCSide() { delete [] aux; }

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

template <class SchemeType, int dim>
class SchemeBoundaryConditions : public BoundaryConditions<typename SchemeType::VectorType,dim> {
  typedef BoundaryConditions<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 SchemeType::SideName SideName;

  SchemeBoundaryConditions(SchemeType &scheme) : base(), _scheme(scheme) {   
    for (int d=0; d<2*dim; d++)
      side[d] = new BCSide<DataType>(Scheme().NEquations());
  }

  ~SchemeBoundaryConditions() {
    for (int d=0; d<2*dim; d++)
      delete side[d];
  }

  virtual void register_at(ControlDevice& Ctrl, const std::string& prefix) {
    base::LocCtrl = Ctrl.getSubDevice(prefix+"BoundaryConditions");
    side[0]->register_at(base::LocCtrl,"Left");
    side[1]->register_at(base::LocCtrl,"Right");
    if (dim>1) {
      side[2]->register_at(base::LocCtrl,"Bottom"); 
      side[3]->register_at(base::LocCtrl,"Top"); 
    }
    if (dim>2) {
      side[4]->register_at(base::LocCtrl,"Front"); 
      side[5]->register_at(base::LocCtrl,"Back"); 
    }
  }

  virtual void SetBndry(vec_grid_data_type &fvec, const int& level, const BBox &bb, 
                        const int &dir, const double& time) {    
    
    BBox wholebbox = base::GH().wholebbox();
    if (dir==SchemeType::Left && bb.upper(0)+bb.stepsize(0)==wholebbox.lower(0))
      Scheme().BCStandard(fvec,bb,side[0]->Type,dir,side[0]->aux,side[0]->NAux,side[0]->Scaling); 
    else if (dir==SchemeType::Right && bb.lower(0)==wholebbox.upper(0)+wholebbox.stepsize(0))
      Scheme().BCStandard(fvec,bb,side[1]->Type,dir,side[1]->aux,side[1]->NAux,side[1]->Scaling); 
    else if (dir==SchemeType::Bottom && bb.upper(1)+bb.stepsize(1)==wholebbox.lower(1))
      Scheme().BCStandard(fvec,bb,side[2]->Type,dir,side[2]->aux,side[2]->NAux,side[2]->Scaling);
    else if (dir==SchemeType::Top && bb.lower(1)==wholebbox.upper(1)+wholebbox.stepsize(1))
      Scheme().BCStandard(fvec,bb,side[3]->Type,dir,side[3]->aux,side[3]->NAux,side[3]->Scaling);
    else if (dir==SchemeType::Front && bb.upper(2)+bb.stepsize(2)==wholebbox.lower(2))
      Scheme().BCStandard(fvec,bb,side[4]->Type,dir,side[4]->aux,side[4]->NAux,side[4]->Scaling);
    else if (dir==SchemeType::Back && bb.lower(2)==wholebbox.upper(2)+wholebbox.stepsize(2))
      Scheme().BCStandard(fvec,bb,side[5]->Type,dir,side[5]->aux,side[5]->NAux,side[5]->Scaling);
  }

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

  inline const BCSide<DataType>& Side(const int d) const { return *side[d]; } 
  
protected:
  SchemeType &_scheme;
  BCSide<DataType> *side[2*dim];
};
  

#endif

---