F77BoundaryConditions.h

Go to the documentation of this file.

// -*- C++ -*-

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

#ifndef AMROC_F77_BOUNDARYCONDITIONS_H
#define AMROC_F77_BOUNDARYCONDITIONS_H

#include "BoundaryConditions.h"

template <class VectorType, int dim>
class F77BoundaryConditions : public BoundaryConditions<VectorType,dim> {
  typedef BoundaryConditions<VectorType,dim> base;
public:
  typedef typename base::vec_grid_data_type vec_grid_data_type;

  typedef generic_fortran_func generic_func_type;

  typedef void (*boundary_1_func_type) ( FI(1,VectorType), BI, const DOUBLE xc[], 
                                         const DOUBLE dx[], const INTEGER& dir,
                                         const DOUBLE *bnd, const INTEGER& mb,  
                                         const DOUBLE& time, const INTEGER& meqn);
  typedef void (*boundary_2_func_type) ( FI(2,VectorType), BI, const DOUBLE xc[], 
                                         const DOUBLE dx[], const INTEGER& dir,
                                         const DOUBLE *bnd, const INTEGER& mb,  
                                         const DOUBLE& time, const INTEGER& meqn);
  typedef void (*boundary_3_func_type) ( FI(3,VectorType), BI, const DOUBLE xc[], 
                                         const DOUBLE dx[], const INTEGER& dir,
                                         const DOUBLE *bnd, const INTEGER& mb,  
                                         const DOUBLE& time, const INTEGER& meqn);

  F77BoundaryConditions() : base(), f_bndry(0) {}
  F77BoundaryConditions(generic_func_type bndry) : base(), f_bndry(bndry) {}
  
  virtual ~F77BoundaryConditions() {}
  
  virtual void SetBndry(vec_grid_data_type &gd, const int& level, const BBox &bb, 
                        const int &dir, const double& time) {    

    assert(f_bndry != (generic_func_type) 0);
    DCoords lbc = base::GH().worldCoords(gd.lower(), gd.stepsize());
    DCoords dx = base::GH().worldStep(gd.stepsize());

    if (dim == 1) 
      ((boundary_1_func_type) f_bndry)(FORTRAN_ARGS(1,gd),BOUNDING_BOX(bb),lbc(),dx(),dir,
                                       base::GH().wholebndry(),base::GH().nbndry(),time,
                                       base::NEquations());
    else if (dim == 2) 
      ((boundary_2_func_type) f_bndry)(FORTRAN_ARGS(2,gd),BOUNDING_BOX(bb),lbc(),dx(),dir,
                                       base::GH().wholebndry(),base::GH().nbndry(),time,
                                       base::NEquations());
    else if (dim == 3) 
      ((boundary_3_func_type) f_bndry)(FORTRAN_ARGS(3,gd),BOUNDING_BOX(bb),lbc(),dx(),dir,
                                       base::GH().wholebndry(),base::GH().nbndry(),time,
                                       base::NEquations());
  }

  inline void SetFunc(generic_func_type bndry) { f_bndry = bndry; }
  generic_func_type GetFunc() const { return f_bndry; }

protected:
  generic_func_type f_bndry;
};


#endif

---