ghostPoint.h

Go to the documentation of this file.

#ifndef GHOSTPOINT_H
#define GHOSTPOINT_H

#include <iostream>
#include <string>
#include <cmath>
using namespace std;

#ifdef SINGLE
#define real float
#else
#define real double
#endif

class ghostPoint {

 public:
  ghostPoint();

  void initialize(int numberOfPatches, 
                  real t,  
                  real dt, 
                  int cellNumberPerPatch, 
                  int ** StartIndex,    
                  int ** pointNumberMeshU, 
                  int ** pointNumberMeshV, 
                  int ** pointNumberMeshW, 
                  int dimension, 
                  real ***xU, 
                  real ***xV, 
                  real ***xW, 
                  real *h, 
                  real ****u, 
                  real ****v, 
                  real ****w  
                  );

  ~ghostPoint();
  
 private: 

  void uCoarseGhostPoint(real ****u);

  void vCoarseGhostPoint(real ****v);

  void wCoarseGhostPoint(real ****w);
  
  void uFineGhostPoint(real iy, // 0 for odd; 1 for even in y direction
                       real iz, // 0 for odd; 1 for even in z direction
                       real ****u);

  void vFineGhostPoint(real ix, // 0 for odd; 1 for even in x direction
                       real iz, // 0 for odd; 1 for even in z direction
                       real ****v);

  void wFineGhostPoint(real ix, // 0 for odd; 1 for even in x direction
                       real iy, // 0 for odd; 1 for even in y direction
                       real ****w);
    
  int nPatch; // each patch carries a square fine mesh
  int nMeshes;   // including the coarse mesh
  // the coarse mesh is the zeroth mesh and fine meshes are from 1st to the last
  int nCellPerPatch; // cell number in the coarse mesh that a patch uses

  int dim;


  // iRHS for index in matrix A = totalUnknowns x totalUnknowns
  // jRHS
  int i, j, k;

  int iFine, jFine, kFine;
  int iCoarse, jCoarse, kCoarse;

  // m for mesh iteration: meshes for patchs and coarse grid
  int m;


};

#endif

---