cageInit.h

Go to the documentation of this file.

#ifndef CAGEINIT_H
#define CAGEINIT_H

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


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

class cageInit {
  
 public:
  cageInit();

  void solve(int numberOfCage, 
             int cellNumberPerPatch, 
             int dimension, 
             // for u
             real ****u, 
             int **pointNumberMeshU, 
             int  ****uCagePoint, 
             // for v
             real ****v, 
             int **pointNumberMeshV, 
             int  ****vCagePoint, 
             // for w
             real ****w, 
             int **pointNumberMeshW, 
             int  ****wCagePoint
             );

  ~cageInit();
  
 private: 

  void setNVariableCage(int &nVariableCage,
                        int ** pointNumberMeshVariable,
                        int **** variableCagePoint);
  void setIndex(int **** indexVariableCage,
                int * indexIVariableCage,
                int * indexJVariableCage,
                int * indexKVariableCage,
                //outer class 
                int ** pointNumberMeshVariable,
                int **** variableCagePoint);

  int * allocate1dArray(int size);
  void delete1dArray(int * x);

  int **** allocate4dIntArrayCage(int nMeshes, int **pointNumberMesh);
  void delete4dArrayCage(int ****x, int nMeshes, int **pointNumberMesh);

  int i, j, k, m; // dummy index
 
  // number of pressure points and vortivity points
  int nCage;
  int nMeshes;
  int nCellPerPatch;
  int dim;

  int **** indexUCage;
  int * indexIUCage;
  int * indexJUCage;
  int * indexKUCage;
  int nUCage;

  int **** indexVCage;
  int * indexIVCage;
  int * indexJVCage;
  int * indexKVCage;
  int nVCage;

  int **** indexWCage;
  int * indexIWCage;
  int * indexJWCage;
  int * indexKWCage;
  int nWCage;


};

#endif

---