amroc/amr/Interfaces/SchemeGFMBoundary.h

Go to the documentation of this file.

// -*- C++ -*-

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

#ifndef SCHEME_GFMBOUNDARY_H
#define SCHEME_GFMBOUNDARY_H

#include "GFMBoundary.h"

template <class SchemeType, int dim>
class SchemeGFMBoundary : public GFMBoundary<typename SchemeType::VectorType,dim> {
  typedef GFMBoundary<typename SchemeType::VectorType,dim> base; 

public:
  typedef typename SchemeType::VectorType VectorType;
  typedef typename VectorType::InternalDataType DataType;
  typedef typename base::vec_grid_fct_type vec_grid_fct_type;  
  typedef typename base::vec_grid_data_type vec_grid_data_type;
  typedef typename base::grid_fct_type grid_fct_type;  
  typedef typename base::grid_data_type grid_data_type;
  typedef typename base::point_type point_type;

public:
  SchemeGFMBoundary(SchemeType &scheme) : base(), _scheme(scheme), Type(0), Scaling(0) {} 

  virtual ~SchemeGFMBoundary() {}
  
  virtual void register_at(ControlDevice& Ctrl, const std::string& prefix) {
    base::register_at(Ctrl,prefix);
    RegisterAt(base::LocCtrl,"Type",Type);
    RegisterAt(base::LocCtrl,"Scaling",Scaling);
  }

  virtual void SetupData(GridHierarchy* gh, const int& ghosts) {
    if (Type==0) 
      base::_bndtype = EXTRAPOLATION;
    base::SetupData(gh, ghosts);
  }

  virtual void SetGrid(vec_grid_data_type& gdu, grid_data_type& gdphi, const BBox& bb,
                       const int& Level, double t, const int& nc, 
                       const int* idx, const point_type* xc, 
                       DataType* distance, point_type* normal) {

    Coords ex = gdu.extents();
    VectorType* u = (VectorType*) 0;
    DataType* aux = (DataType*) 0;
    DCoords dx = base::GH().worldStep(gdu.stepsize());
    if (base::NAux()>0)
      aux = new DataType[nc*base::NAux()];

    if (base::NormalFactor()>0) {
      if (UseTransform()) {
        vec_grid_data_type gdhelp(gdu.bbox());  
        Transform(gdu,gdhelp);
        base::Extrapolation(gdhelp,gdphi,u,Level,nc,idx,xc,
                            distance,normal,base::NormalFactor());
        if (base::NAux()>0) {
          START_WATCH
            SetBndryAux(gdhelp,gdphi,u,aux,Level,t,nc,idx,xc,distance,normal);
          END_WATCH(GFM_AUXILIARY_VALUES)
        }
      }
      else {
        base::Extrapolation(gdu,gdphi,u,Level,nc,idx,xc,
                            distance,normal,base::NormalFactor());
        if (base::NAux()>0) {
          START_WATCH
            SetBndryAux(gdu,gdphi,u,aux,Level,t,nc,idx,xc,distance,normal);
          END_WATCH(GFM_AUXILIARY_VALUES)
        }
      }
    }
    else if (base::NAux()>0) {
      START_WATCH
        SetBndryAux(gdu,gdphi,u,aux,Level,t,nc,idx,xc,distance,normal);
      END_WATCH(GFM_AUXILIARY_VALUES)
    }

    START_WATCH
      SetBndry(gdu,nc,idx,u,xc,distance,normal,aux,base::NAux(),t);
    END_WATCH(GFM_SETBNDRY)
    
    if (u) delete [] u;
    if (aux) delete [] aux;
  }

  virtual bool UseTransform()
  { return Scheme().GFMUseTransform(); }

  virtual void Transform(vec_grid_data_type& gdu, vec_grid_data_type& gdhelp)
  { Scheme().GFMTransform(gdu,gdhelp); }
     
  virtual void SetBndry(vec_grid_data_type& gdu, const int& nc, const int* idx,
                        const VectorType* u, const point_type* xc, const DataType* distance, 
                        const point_type* normal, DataType* aux, const int naux, const double t)  
  { Scheme().GFMBCStandard(gdu,Type,nc,idx,u,xc,distance,normal,aux,naux,Scaling); }

  virtual void SetBndryAux(vec_grid_data_type& gdu, grid_data_type& gdphi, 
                           const VectorType* u, DataType* aux, const int& Level, 
                           double t, const int& nc, const int* idx, 
                           const point_type* xc, const DataType* distance, 
                           const point_type* normal) {}

  inline SchemeType& Scheme() { return _scheme; }
  inline const SchemeType& Scheme() const { return _scheme; }
  
protected:
  SchemeType &_scheme;
  int Type, Scaling;
};

#endif

---