amroc/amr/F77Interfaces/F77GFMInterpolation.h

Go to the documentation of this file.

// -*- C++ -*-

// Copyright (C) 2003-2007 California Institute of Technology
// Ralf Deiterding, ralf@cacr.caltech.edu

#ifndef AMROC_F77_GFM_INTERPOLATION_H
#define AMROC_F77_GFM_INTERPOLATION_H

#include "F77Interfaces/F77Interpolation.h"

template <class VectorType, int dim>
class F77GFMInterpolation : public F77Interpolation<VectorType,dim> {
  typedef F77Interpolation<VectorType,dim> base;
  typedef typename VectorType::InternalDataType DataType;
public:
  typedef typename base::vec_grid_data_type vec_grid_data_type;
  typedef typename base::grid_data_type grid_data_type;  
  typedef typename base::point_type point_type;
  typedef typename base::bool_grid_data_type bool_grid_data_type;

  typedef generic_fortran_func generic_func_type;

  typedef void (*intex_1_func_type) ( const DOUBLE q[], const INTEGER &mx, 
                                      const INTEGER &meqn, const DOUBLE phi[], const INTEGER &nc, 
                                      DOUBLE qint[], const DOUBLE xc[],
                                      const DOUBLE corn[], const DOUBLE dx[], 
                                      const DOUBLE& ev);

  typedef void (*intex_2_func_type) ( const DOUBLE q[], const INTEGER &mx, 
                                      const INTEGER &my, 
                                      const INTEGER &meqn, const DOUBLE phi[], const INTEGER &nc, 
                                      DOUBLE qint[], const DOUBLE xc[],
                                      const DOUBLE corn[], const DOUBLE dx[], 
                                      const DOUBLE& ev);

  typedef void (*intex_3_func_type) ( const DOUBLE q[], const INTEGER &mx, 
                                      const INTEGER &my, const INTEGER &mz, 
                                      const INTEGER &meqn, const DOUBLE phi[], const INTEGER &nc, 
                                      DOUBLE qint[], const DOUBLE xc[],
                                      const DOUBLE corn[], const DOUBLE dx[], 
                                      const DOUBLE& ev);

  F77GFMInterpolation() : base(), f_intexpol(0) {}
  F77GFMInterpolation(generic_func_type intexpol) : base(), f_intexpol(intexpol) {}
  F77GFMInterpolation(generic_func_type intpol, generic_func_type intexpol) : base(intpol), 
    f_intexpol(intexpol) {}
  
  virtual ~F77GFMInterpolation() {}
  
  virtual void Interpolate(GridHierarchy& GH, vec_grid_data_type& gdu, const int& nc, 
                           const point_type* xc, VectorType* uv, const DataType& ErrorValue) {
    base::Interpolate(GH,gdu,nc,xc,uv,ErrorValue);
  }    

  virtual void Interpolate(GridHierarchy& GH, vec_grid_data_type& gdu, 
                           grid_data_type& gdphi, const int& nc, const point_type* xc, 
                           VectorType* uv, const DataType& ErrorValue) {
    if (!f_intexpol) return;
    DCoords dx = GH.worldStep(gdu.stepsize());
    DCoords corn = GH.worldCoords(gdu.lower(),gdu.stepsize());
    int mx[3];
    for (register int d=0; d<dim; d++) 
      mx[d] = gdu.extents(d);

    if (dim == 1)
      ((intex_1_func_type) f_intexpol)((DataType*)(gdu.data()),AA(1,mx),base::NEquations(),
                                       gdphi.data(),nc,uv[0].data(),xc[0].data(),corn(),dx(),
                                       ErrorValue); 
    else if (dim == 2) 
      ((intex_2_func_type) f_intexpol)((DataType*)(gdu.data()),AA(2,mx),base::NEquations(),
                                       gdphi.data(),nc,uv[0].data(),xc[0].data(),corn(),dx(),
                                       ErrorValue); 
    else if (dim == 3) 
      ((intex_3_func_type) f_intexpol)((DataType*)(gdu.data()),AA(3,mx),base::NEquations(),
                                       gdphi.data(),nc,uv[0].data(),xc[0].data(),corn(),dx(),
                                       ErrorValue); 
  }    

  virtual void Interpolate(GridHierarchy& GH, vec_grid_data_type& gdu, 
                           bool_grid_data_type& gdflg, const int& nc, const point_type* xc, 
                           VectorType* uv, const DataType& ErrorValue) {
    if (!f_intexpol) return;
    grid_data_type phihelp(gdflg.bbox());
    base::SetPhiFromBool(phihelp,gdflg);
    Interpolate(GH,gdu,phihelp,nc,xc,uv,ErrorValue);
  }

  virtual void Interpolate(GridHierarchy& GH, grid_data_type& gd, const int& nc, 
                           const point_type* xc, DataType* u, const DataType& ErrorValue) {
    base::Interpolate(GH,gd,nc,xc,u,ErrorValue);
  }

  virtual void Interpolate(GridHierarchy& GH, grid_data_type& gd, 
                           grid_data_type& gdphi, const int& nc, const point_type* xc, 
                           DataType* u, const DataType& ErrorValue) {
    if (!f_intexpol) return;
    DCoords dx = GH.worldStep(gd.stepsize());
    DCoords corn = GH.worldCoords(gd.lower(),gd.stepsize());
    int mx[3];
    for (register int d=0; d<dim; d++) 
      mx[d] = gd.extents(d);
    int meqn = 1;

    if (dim == 1) 
      ((intex_1_func_type) f_intexpol)(gd.data(),AA(1,mx),meqn,gdphi.data(),nc,
                                       u,xc[0].data(),corn(),dx(),ErrorValue); 
    else if (dim == 2) 
      ((intex_2_func_type) f_intexpol)(gd.data(),AA(2,mx),meqn,gdphi.data(),nc,
                                       u,xc[0].data(),corn(),dx(),ErrorValue); 
    else if (dim == 3) 
      ((intex_3_func_type) f_intexpol)(gd.data(),AA(3,mx),meqn,gdphi.data(),nc,
                                       u,xc[0].data(),corn(),dx(),ErrorValue); 
  }

  virtual void Interpolate(GridHierarchy& GH, grid_data_type& gd, 
                           bool_grid_data_type& gdflg, const int& nc, const point_type* xc, 
                           DataType* u, const DataType& ErrorValue) {
    if (!f_intexpol) return;
    grid_data_type phihelp(gdflg.bbox());
    base::SetPhiFromBool(phihelp,gdflg);
    Interpolate(GH,gd,phihelp,nc,xc,u,ErrorValue);
  }

  inline void SetIntExPolFunc(generic_func_type intexpol) { f_intexpol = intexpol; }
  generic_func_type GetIntExPolFunc() const { return f_intexpol; }

protected:
  generic_func_type f_intexpol;
};


#endif

---