amroc/amr/F77Interfaces/F77UpdateLevelTransfer.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@cacr.caltech.edu

#ifndef AMROC_F77_UPDATE_LEVELTRANSFER_H
#define AMROC_F77_UPDATE_LEVELTRANSFER_H

#include "F77Interfaces/F77LevelTransfer.h"

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

  typedef void (*update_1_func_type) ( FI(1,VectorType), BI, const INTEGER& meqn );
  typedef void (*update_2_func_type) ( FI(2,VectorType), BI, const INTEGER& meqn );
  typedef void (*update_3_func_type) ( FI(3,VectorType), BI, const INTEGER& meqn );

  F77UpdateLevelTransfer() : base(), f_update(0) {}
  F77UpdateLevelTransfer(generic_func_type prolong, generic_func_type restrct,
                         generic_func_type update) : base(prolong,restrct), f_update(update) {}
  
  virtual ~F77UpdateLevelTransfer() {}
  
  virtual void SetupData(GridHierarchy* gh, const int& ghosts) {  
    base::SetAdaptiveBoundaryType(DAGHAdaptBoundaryBoth);
    base::SetupData(gh,ghosts);
  }

  void Update(vec_grid_data_type &target, const BBox &bb) {    
    if (!f_update) return;
    if (dim == 1) 
      ((update_1_func_type) f_update)(FORTRAN_ARGS(1,target),BOUNDING_BOX(bb),base::NEquations()); 
    else if (dim == 2) 
      ((update_2_func_type) f_update)(FORTRAN_ARGS(2,target),BOUNDING_BOX(bb),base::NEquations()); 
    else if (dim == 3) 
      ((update_3_func_type) f_update)(FORTRAN_ARGS(3,target),BOUNDING_BOX(bb),base::NEquations()); 
  }

  virtual void Prolong(vec_grid_data_type &source, const int& source_level, 
                       vec_grid_data_type &target, const int& target_level, const BBox &bb) {
    base::Prolong(source,source_level,target,target_level,bb);
    Update(target,bb);
  }    

  virtual void Restrict(vec_grid_data_type &source, const int& source_level, 
                        vec_grid_data_type &target, const int& target_level, const BBox &bb) { 
    base::Restrict(source,source_level,target,target_level,bb);
    Update(target,bb);
  }

  virtual void SetAdaptBndry(vec_grid_data_type &target, vec_grid_data_type &target_help_previous, 
                             vec_grid_data_type &target_help_next, vec_grid_data_type &source_next, const double &frac,
                             vec_grid_data_type &source_previous, const double &oneminusfrac, 
                             const int &target_level, const BBox &bb, const int side) {
    base::SetAdaptBndry(target,target_help_previous,target_help_next,source_next,frac,
                        source_previous,oneminusfrac,target_level,bb,side);
    Update(target,bb);    
  }


  inline void SetUpdateFunc(generic_func_type update) { f_update = update; }
  generic_func_type GetUpdateFunc() const { return f_update; }

protected:
  generic_func_type f_update;
};


#endif

---