amroc/amr/F77Interfaces/F77LevelTransfer.h

Go to the documentation of this file.

// -*- C++ -*-

// Copyright (C) 2002 Ralf Deiterding
// Brandenburgische Universitaet Cottbus

#ifndef AMROC_F77_LEVELTRANSFER_H
#define AMROC_F77_LEVELTRANSFER_H

#include "LevelTransfer.h"

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

  typedef generic_fortran_func generic_func_type;

  typedef void (*transfer_1_func_type) ( FI(1,VectorType), FI(1,VectorType), BI, 
                                         const INTEGER& meqn, const INTEGER& mbc);
  typedef void (*transfer_2_func_type) ( FI(2,VectorType), FI(2,VectorType), BI, 
                                         const INTEGER& meqn, const INTEGER& mbc);
  typedef void (*transfer_3_func_type) ( FI(3,VectorType), FI(3,VectorType), BI, 
                                         const INTEGER& meqn, const INTEGER& mbc);

  F77LevelTransfer() : base(), f_prolong(0), f_restrict(0) {}
  F77LevelTransfer(generic_func_type prolong, generic_func_type restrct) : 
    base(), f_prolong(prolong), f_restrict(restrct) {}
  
  virtual ~F77LevelTransfer() {}
  
  virtual void Prolong(vec_grid_data_type &source, const int& source_level, 
                       vec_grid_data_type &target, const int& target_level, const BBox &bb) { 
    if (!f_prolong) return;
    if (dim == 1) 
      ((transfer_1_func_type) f_prolong)(FORTRAN_ARGS(1,source),FORTRAN_ARGS(1,target),BOUNDING_BOX(bb),
                                         base::NEquations(),base::NGhosts()); 
    else if (dim == 2) 
      ((transfer_2_func_type) f_prolong)(FORTRAN_ARGS(2,source),FORTRAN_ARGS(2,target),BOUNDING_BOX(bb),
                                         base::NEquations(),base::NGhosts()); 
    else if (dim == 3) 
      ((transfer_3_func_type) f_prolong)(FORTRAN_ARGS(3,source),FORTRAN_ARGS(3,target),BOUNDING_BOX(bb),
                                         base::NEquations(),base::NGhosts()); 
  }    

  virtual void Restrict(vec_grid_data_type &source, const int& source_level, 
                        vec_grid_data_type &target, const int& target_level, const BBox &bb) { 
    if (!f_restrict) return;
    if (dim == 1) 
      ((transfer_1_func_type) f_restrict)(FORTRAN_ARGS(1,source),FORTRAN_ARGS(1,target),BOUNDING_BOX(bb),
                                          base::NEquations(),base::NGhosts()); 
    else if (dim == 2) 
      ((transfer_2_func_type) f_restrict)(FORTRAN_ARGS(2,source),FORTRAN_ARGS(2,target),BOUNDING_BOX(bb),
                                          base::NEquations(),base::NGhosts()); 
    else if (dim == 3) 
      ((transfer_3_func_type) f_restrict)(FORTRAN_ARGS(3,source),FORTRAN_ARGS(3,target),BOUNDING_BOX(bb),
                                          base::NEquations(),base::NGhosts()); 
  }

  inline void SetProlongFunc(generic_func_type prolong) { f_prolong = prolong; }
  generic_func_type GetProlongFunc() const { return f_prolong; }
  inline void SetRestrictFunc(generic_func_type restrct) { f_restrict = restrct; }
  generic_func_type GetRestrictFunc() const { return f_restrict; }

protected:
  generic_func_type f_prolong, f_restrict;
};


#endif

---