GridGenerator.C

Go to the documentation of this file.

#include "Testing.H"
#include "Simulation.H"
#include "OperatorKernels.H"
#include "Stencil.H"
#include "PCPPCommUtil.H"
#include "PCPPReport.H"
#include "PCPPIntervalUtils.H"
#include "PCPPHDF5.H"
#include "Report.H"
#include "PC2IO.H"

namespace gridgen {

  template<int numDim>
  int RectilinearGrid1(const std::vector<size_t> &ijk,const std::vector<double> &scale,std::vector<double> &xyz)
  {
    //    double scale[] = {101,51,21};
    for(int iDim = 0;iDim < numDim;iDim++){
      double iVar = ((double)(ijk[iDim]+1))/scale[iDim];
      xyz[iDim] = (iDim+1)*std::pow(iVar,2.0);
    }
    return(0);
  };

  template<int numDim>
  int RectilinearMetric1(const std::vector<size_t> &ijk,
                         const std::vector<double> &scale,
                         std::vector<double> &xyz,double &jacm1)
  {
    //    double scale[] = {101,51,21};
    double dxDXi  = 2.0*(ijk[0]+1)/(scale[0]*scale[0]);
    double dyDEta = 4.0*(ijk[1]+1)/(scale[1]*scale[1]);
    xyz[0] = dyDEta;
    xyz[1] = dxDXi;
    jacm1 = xyz[0]*xyz[1];
    if (numDim == 3){
      double dzDZeta = 6.0*(ijk[2]+1)/(scale[2]*scale[2]);
      xyz[0] = dyDEta*dzDZeta;
      xyz[1] = dxDXi*dzDZeta;
      xyz[2] = dxDXi*dyDEta;
      jacm1 = xyz[0]*xyz[1]*xyz[2];
    }
    return(0);
  };

}