Blockstructured Adaptive Mesh Refinement in object-oriented C++

---

       

---

ClpInitialCondition.h

Go to the documentation of this file.

#ifndef AMROC_CLP_INITIALCONDITION_H
#define AMROC_CLP_INITIALCONDITION_H

#include "InitialCondition.h"

#ifndef ClpInitialConditionName
#define ClpInitialCondition(dim)      name2(ClpInitialCondition,dim)
#define ClpInitialConditionName
#endif

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

public:
  ClpInitialCondition(DIM)(const int gh) :
    InitialCondition(DIM)<VectorType>(gh) {}

  virtual void SetGrid(vec_grid_data_type& gd, 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_initial(AA(DIM,mx), NEquations(), NGhosts(), AA(DIM,mx), 
              AA(DIM,x), AA(DIM,dx()), 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