WENOGFMBoundary.h

Go to the documentation of this file.

// -*- C++ -*-

// Copyright (C) 2003-2007 California Institute of Technology
// Ralf Deiterding, ralf@amroc.net

#ifndef AMROC_WENO_GFMBOUNDARY_H
#define AMROC_WENO_GFMBOUNDARY_H

#if DIM == 1
#elif DIM == 2

#include "GFMBoundary2.h"

#define INTERNALBND_PARAMETER_LIST ( \
        FI(2,Vector_type), const INTEGER &, const INTEGER &, \
        const INTEGER *, const DOUBLE *, const DOUBLE *, \
        const DOUBLE *, const DOUBLE *, const DOUBLE *, \
        const DOUBLE &)
extern "C" {
  typedef void (*intbnd_func_type) INTERNALBND_PARAMETER_LIST;
}
#define INTTRANSBND_PARAMETER_LIST ( \
        const INTEGER& mx, const INTEGER& my, \
        const INTEGER& meqn, Vector_type* q, Vector_type* qt)
extern "C" {
  typedef void (*inttrans_func_type) INTTRANSBND_PARAMETER_LIST;
}
#elif DIM == 3  

#include "GFMBoundary3.h"

#define INTERNALBND_PARAMETER_LIST ( \
        FI(3,Vector_type), const INTEGER &, const INTEGER &, \
        const INTEGER *, const DOUBLE *, const DOUBLE *, \
        const DOUBLE *, const DOUBLE *, const DOUBLE *, \
        const DOUBLE &)
extern "C" {
  typedef void (*intbnd_func_type) INTERNALBND_PARAMETER_LIST;
}
#define INTTRANSBND_PARAMETER_LIST ( \
        const INTEGER& mx, const INTEGER& my, const INTEGER& mz, \
        const INTEGER& meqn, Vector_type* q, Vector_type* qt)
extern "C" {
  typedef void (*inttrans_func_type) INTTRANSBND_PARAMETER_LIST;
}

#endif

#ifndef WENOGFMBoundaryName
#define WENOGFMBoundary(dim)      name2(WENOGFMBoundary,dim)
#define WENOGFMBoundaryName
#endif

template <class VectorType>
class WENOGFMBoundary(DIM) : public GFMBoundary(DIM)<VectorType> {
  typedef typename VectorType::InternalDataType DataType;
  typedef GridData(DIM)<DataType>   grid_data_type;   
  typedef GridData(DIM)<VectorType> vec_grid_data_type;  
  typedef GFMBoundary(DIM)<VectorType> base; 
  typedef Vector(DIM)<DataType>  point_type;

public:
  WENOGFMBoundary(DIM)(const int gh, intbnd_func_type bf, inttrans_func_type tf) : 
    base(gh), bnd_func(bf), trans_func(tf) {} 

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

    Coords ex = gdu.extents();
    VectorType* u = (VectorType*) 0;

    if (base::NormalFactor()>0)
      if (trans_func) {
        vec_grid_data_type gdhelp(gdu.bbox());  
        (*trans_func)(AA(DIM,ex()),base::NEquations(),gdu.data(),gdhelp.data());
        base::Extrapolation(gdhelp,gdphi,u,Level,dx,Ncells,indices,xc,
                            distance,normal,base::NormalFactor());
      }
      else 
        base::Extrapolation(gdu,gdphi,u,Level,dx,Ncells,indices,xc,
                            distance,normal,base::NormalFactor());

    if (bnd_func)
      (*bnd_func)(FORTRAN_ARGS(DIM,gdu),base::NEquations(),Ncells,
                  indices,u[0].data(),xc[0].data(),distance,
                  normal[0].data(),dx(),t);
    
    if (u) delete [] u;
  }

protected:
  intbnd_func_type bnd_func;
  inttrans_func_type trans_func;
};


#endif