Blockstructured Adaptive Mesh Refinement in object-oriented C++

---

       

---

ClpExactSolution.h

Go to the documentation of this file.

#ifndef AMROC_CLP_EXACT_SOLUTION_H
#define AMROC_CLP_EXACT_SOLUTION_H

#include "ExactSolution.h"

#ifndef ClpExactSolutionName
#define ClpExactSolution(dim)      name2(ClpExactSolution,dim)
#define ClpExactSolutionName
#endif

#if DIM == 1
# define f_exact FORTRAN_NAME(exact_, EXACT, exact)
extern "C" {
  void f_exact(const INTEGER& maxmx, 
               const INTEGER& meqn, const INTEGER& mbc,
               const INTEGER& mx, 
               const DOUBLE x[], 
               const DOUBLE& dx, 
               const DOUBLE& t, Vector_type* q);
}
#elif DIM == 2
# define f_exact FORTRAN_NAME(exact_, EXACT, exact)
extern "C" {
  void f_exact(const INTEGER& maxmx, const INTEGER& maxmy, 
               const INTEGER& meqn, const INTEGER& mbc,
               const INTEGER& mx, const INTEGER& my, 
               const DOUBLE x[], const DOUBLE y[], 
               const DOUBLE& dx, const DOUBLE& dy, 
               const DOUBLE& t, Vector_type* q);
}
#elif DIM == 3  
# define f_exact FORTRAN_NAME(exact_, EXACT, exact)
extern "C" {
  void f_exact(const INTEGER& maxmx, const INTEGER& maxmy, const INTEGER& maxmz, 
               const INTEGER& meqn, const INTEGER& mbc,
               const INTEGER& mx, const INTEGER& my, const INTEGER& mz, 
               const DOUBLE x[], const DOUBLE y[], const DOUBLE z[], 
               const DOUBLE& dx, const DOUBLE& dy, const DOUBLE& dz, 
               const DOUBLE& t, Vector_type* q);
}
#endif

template <class VectorType>
class ClpExactSolution(DIM) : public ExactSolution(DIM)<VectorType> {
  typedef typename VectorType::InternalDataType DataType;
  typedef GridData(DIM)<VectorType> vec_grid_data_type;

public:
  ClpExactSolution(DIM)(const int gh) : ExactSolution(DIM)<VectorType>(gh) {}

  virtual void SetGrid(vec_grid_data_type& gd, double t, 
                       DCoords& dx) {
    Coords ex = gd.extents();
    DCoords lbcorner = GH().worldCoords(gd.lower(), gd.stepsize());
    int maxm[DIM], mx[DIM], d;
    DataType* x[DIM];
    for (d=0; d<DIM; d++) {
      maxm[d] = ex(d);
      mx[d] = ex(d)-2*NGhosts();
      x[d] = new DataType[maxm[d]];
      for (int i=0; i<maxm[d]; i++) x[d][i] = (i+0.5)*dx(d)+lbcorner(d);
    }

    f_exact(AA(DIM,mx), NEquations(), NGhosts(), AA(DIM,mx), 
            AA(DIM,x), AA(DIM,dx()), t, gd.data()); 

    for (d=0; d<DIM; d++) 
      delete [] x[d];
  }
};


#endif

---

Quickstart     Users Guide     Programmers Reference     Installation      Examples     Download

---

AMROC Main      Home      Contact

last update: 06/01/04