src/expde/grid/listint/listint4.h

Go to the documentation of this file.

//    expde: expression templates for partial differential equations.
//    Copyright (C) 2001  Christoph Pflaum
//    This program is free software; you can redistribute it and/or modify
//    it under the terms of the GNU General Public License as published by
//    the Free Software Foundation; either version 2 of the License, or
//    (at your option) any later version.
//
//    This program is distributed in the hope that it will be useful,
//    but WITHOUT ANY WARRANTY; without even the implied warranty of
//    MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
//    GNU General Public License for more details.
//
//                 SEE  Notice1.doc made by 
//                 LAWRENCE LIVERMORE NATIONAL LABORATORY
//

// ------------------------------------------------------------
// listint3.h
//
// ------------------------------------------------------------

#define initialize_interior true


class P_interior {
 public:
  P_interior(Index3D I, P_interior* Next) : ind(I), next(Next) {};

  // initialize
  void initialize(const Grid* gitter, int var_ebene);

  // same level
  // --------------
  // idem:
  double *varM(const Grid* gitter,int level) const;
  // next:
  double *varN(const Grid* gitter, int level) const;
  double *varS(const Grid* gitter, int level) const;
  double *varW(const Grid* gitter, int level) const;
  double *varE(const Grid* gitter, int level) const;
  double *varT(const Grid* gitter, int level) const;
  double *varD(const Grid* gitter, int level) const;
  // double next:
  double *varND(const Grid* gitter, int level) const;
  double *varWN(const Grid* gitter, int level) const;
  double *varWT(const Grid* gitter, int level) const;
  double *varED(const Grid* gitter, int level) const;
  double *varST(const Grid* gitter, int level) const;
  double *varES(const Grid* gitter, int level) const;
  // third next:
  double *varEST(const Grid* gitter, int level) const;
  double *varWND(const Grid* gitter, int level) const;

  // double next special
  double *varWD(const Grid* gitter, int level) const;
  double *varET(const Grid* gitter, int level) const;

  // cell variable:
  double *cell_varWSD(const Grid* gitter,int lev) const;
  double *cell_varWST(const Grid* gitter,int lev) const;
  double *cell_varWND(const Grid* gitter,int lev) const;
  double *cell_varWNT(const Grid* gitter,int lev) const;
  double *cell_varESD(const Grid* gitter,int lev) const;
  double *cell_varEST(const Grid* gitter,int lev) const;
  double *cell_varEND(const Grid* gitter,int lev) const;
  double *cell_varENT(const Grid* gitter,int lev) const;

  // finer level
  // --------------
  // idem:
  double *var_fine_M(const Grid* gitter, int level) const;
  // next:
  double *var_fine_N(const Grid* gitter, int level) const;
  double *var_fine_S(const Grid* gitter, int level) const;
  double *var_fine_W(const Grid* gitter, int level) const;
  double *var_fine_E(const Grid* gitter, int level) const;
  double *var_fine_T(const Grid* gitter, int level) const;
  double *var_fine_D(const Grid* gitter, int level) const;
  // double next:
  double *var_fine_ND(const Grid* gitter, int level) const;
  double *var_fine_WN(const Grid* gitter, int level) const;
  double *var_fine_WT(const Grid* gitter, int level) const;
  double *var_fine_ED(const Grid* gitter, int level) const;
  double *var_fine_ST(const Grid* gitter, int level) const;
  double *var_fine_ES(const Grid* gitter, int level) const;
  // third next:
  double *var_fine_EST(const Grid* gitter, int level) const;
  double *var_fine_WND(const Grid* gitter, int level) const;

  // coarse level
  // --------------
  // idem:
  double *var_coarse_M(const Grid* gitter, int level) const;
  // next:
  double *var_coarse_N(const Grid* gitter, int level) const;
  double *var_coarse_S(const Grid* gitter, int level) const;
  double *var_coarse_W(const Grid* gitter, int level) const;
  double *var_coarse_E(const Grid* gitter, int level) const;
  double *var_coarse_T(const Grid* gitter, int level) const;
  double *var_coarse_D(const Grid* gitter, int level) const;
  // double next:
  double *var_coarse_ND(const Grid* gitter, int level) const;
  double *var_coarse_WN(const Grid* gitter, int level) const;
  double *var_coarse_WT(const Grid* gitter, int level) const;
  double *var_coarse_ED(const Grid* gitter, int level) const;
  double *var_coarse_ST(const Grid* gitter, int level) const;
  double *var_coarse_ES(const Grid* gitter, int level) const;
  // third next:
  double *var_coarse_EST(const Grid* gitter, int level) const;
  double *var_coarse_WND(const Grid* gitter, int level) const;

  // physical coordinates of the point
  D3vector coordinate(const Grid* gitter) const;
  double coordinateX(const Grid* gitter) const;
  double coordinateY(const Grid* gitter) const;
  double coordinateZ(const Grid* gitter) const;

  // for iteration:
  P_interior* Next() { return next; };

  // index of the point
  Index3D Ind() const { return ind; };
 private:
  // same level:
  double* pointer_varM;

  double* pointer_varN;
  double* pointer_varS;
  double* pointer_varE;
  double* pointer_varW;
  double* pointer_varT;
  double* pointer_varD;

  double* pointer_varND;
  double* pointer_varWN;
  double* pointer_varWT;
  double* pointer_varED;
  double* pointer_varST;
  double* pointer_varES;

  double* pointer_varEST;
  double* pointer_varWND;

  // finer level:
  double* pointer_fine_varM;

  double* pointer_fine_varN;
  double* pointer_fine_varS;
  double* pointer_fine_varE;
  double* pointer_fine_varW;
  double* pointer_fine_varT;
  double* pointer_fine_varD;

  double* pointer_fine_varND;
  double* pointer_fine_varWN;
  double* pointer_fine_varWT;
  double* pointer_fine_varED;
  double* pointer_fine_varST;
  double* pointer_fine_varES;

  double* pointer_fine_varEST;
  double* pointer_fine_varWND;

  Index3D ind;
  P_interior *next;
};



// coordinates
inline D3vector P_interior::coordinate(const Grid* gitter) const {
  return gitter->transform_coord(ind);
};
inline double P_interior::coordinateX(const Grid* gitter) const {
  return gitter->transform_coordX(ind);
};
inline double P_interior::coordinateY(const Grid* gitter) const {
  return gitter->transform_coordY(ind);
};
inline double P_interior::coordinateZ(const Grid* gitter) const {
  return gitter->transform_coordZ(ind);
};


inline void P_interior::initialize(const Grid* gitter, int level) {
  // same level:
  pointer_varM = gitter->Give_variable(ind,level);

  pointer_varN = gitter->Give_variable(ind.next_NS(Rd,level),level);
  pointer_varS = gitter->Give_variable(ind.next_NS(Ld,level),level);
  pointer_varE = gitter->Give_variable(ind.next_EW(Rd,level),level);
  pointer_varW = gitter->Give_variable(ind.next_EW(Ld,level),level);
  pointer_varT = gitter->Give_variable(ind.next_TD(Rd,level),level);
  pointer_varD = gitter->Give_variable(ind.next_TD(Ld,level),level);
  pointer_varND = 
    gitter->Give_variable(ind.next_NS(Rd,level).next_TD(Ld,level),level);
  pointer_varWN = 
    gitter->Give_variable(ind.next_NS(Rd,level).next_EW(Ld,level),level);
  pointer_varWT = 
    gitter->Give_variable(ind.next_TD(Rd,level).next_EW(Ld,level),level);
  pointer_varED =
    gitter->Give_variable(ind.next_TD(Ld,level).next_EW(Rd,level),level);
  pointer_varST =
    gitter->Give_variable(ind.next_TD(Rd,level).next_NS(Ld,level),level);
  pointer_varES =
    gitter->Give_variable(ind.next_EW(Rd,level).next_NS(Ld,level),level);
  pointer_varEST =
    gitter->Give_variable(
          ind.next_EW(Rd,level).next_NS(Ld,level).next_TD(Rd,level),level);
  pointer_varWND =
    gitter->Give_variable(
          ind.next_EW(Ld,level).next_NS(Rd,level).next_TD(Ld,level),level);

  // finer level
  if(level<gitter->Max_level()) {
    pointer_fine_varM = gitter->Give_variable(ind,level+1);
    
    pointer_fine_varN = gitter->Give_variable(ind.next_NS(Rd,level+1),level+1);
    pointer_fine_varS = gitter->Give_variable(ind.next_NS(Ld,level+1),level+1);
    pointer_fine_varE = gitter->Give_variable(ind.next_EW(Rd,level+1),level+1);
    pointer_fine_varW = gitter->Give_variable(ind.next_EW(Ld,level+1),level+1);
    pointer_fine_varT = gitter->Give_variable(ind.next_TD(Rd,level+1),level+1);
    pointer_fine_varD = gitter->Give_variable(ind.next_TD(Ld,level+1),level+1);
    pointer_fine_varND = 
      gitter->Give_variable(ind.next_NS(Rd,level+1).next_TD(Ld,level+1),level+1);
    pointer_fine_varWN = 
      gitter->Give_variable(ind.next_NS(Rd,level+1).next_EW(Ld,level+1),level+1);
    pointer_fine_varWT = 
      gitter->Give_variable(ind.next_TD(Rd,level+1).next_EW(Ld,level+1),level+1);
    pointer_fine_varED =
      gitter->Give_variable(ind.next_TD(Ld,level+1).next_EW(Rd,level+1),level+1);
    pointer_fine_varST =
      gitter->Give_variable(ind.next_TD(Rd,level+1).next_NS(Ld,level+1),level+1);
    pointer_fine_varES =
      gitter->Give_variable(ind.next_EW(Rd,level+1).next_NS(Ld,level+1),level+1);
    pointer_fine_varEST =
      gitter->Give_variable(
        ind.next_EW(Rd,level+1).next_NS(Ld,level+1).next_TD(Rd,level+1),level+1);
    pointer_fine_varWND =
      gitter->Give_variable(
        ind.next_EW(Ld,level+1).next_NS(Rd,level+1).next_TD(Ld,level+1),level+1);
  }
};

// same level
// --------------
// idem:
inline double *P_interior::varM(const Grid* gitter, int level) const {
  return pointer_varM;
};
// next:
inline double *P_interior::varN(const Grid* gitter, int level)
     const {
  return pointer_varN;
};
inline double *P_interior::varS(const Grid* gitter, int level)
     const {
  return pointer_varS;
};
inline double *P_interior::varE(const Grid* gitter, int level)
     const {
  return pointer_varE;
};
inline double *P_interior::varW(const Grid* gitter, int level)
     const {
  return pointer_varW;
};
inline double *P_interior::varT(const Grid* gitter, int level)
     const {
  return pointer_varT;
};
inline double *P_interior::varD(const Grid* gitter, int level)
     const {
  return pointer_varD;
};
// double next:
inline double *P_interior::varND(const Grid* gitter, int level)
     const {
  return pointer_varND;
};
inline double *P_interior::varWN(const Grid* gitter, int level)
     const {
  return pointer_varWN;
};
inline double *P_interior::varWT(const Grid* gitter, int level)
     const {
  return pointer_varWT;
};
inline double *P_interior::varED(const Grid* gitter, int level)
     const {
  return pointer_varED;
};
inline double *P_interior::varST(const Grid* gitter, int level)
     const {
  return pointer_varST;
};
inline double *P_interior::varES(const Grid* gitter, int level)
     const {
  return pointer_varES;
};

// double next special:
inline double *P_interior::varWD(const Grid* gitter, int level)
     const {
  return gitter->Give_variable(ind.next_WD(level),level);
};
inline double *P_interior::varET(const Grid* gitter, int level)
     const {
  return gitter->Give_variable(ind.next_ET(level),level);
};

// third next:
inline double *P_interior::varEST(const Grid* gitter, int level)
     const {
  return pointer_varEST;
};
inline double *P_interior::varWND(const Grid* gitter, int level)
     const {
  return pointer_varWND;
};

// cell variable
inline double *P_interior::cell_varENT(const Grid* gitter,int lev) const {
  return gitter->Give_cell_variable(ind.next_ENT(lev+1));
};
inline double *P_interior::cell_varEND(const Grid* gitter,int lev) const {
  return gitter->Give_cell_variable(ind.next_END(lev+1));
};
inline double *P_interior::cell_varEST(const Grid* gitter,int lev) const {
  return gitter->Give_cell_variable(ind.next_EST(lev+1));
};
inline double *P_interior::cell_varESD(const Grid* gitter,int lev) const {
  return gitter->Give_cell_variable(ind.next_ESD(lev+1));
};
inline double *P_interior::cell_varWNT(const Grid* gitter,int lev) const {
  return gitter->Give_cell_variable(ind.next_WNT(lev+1));
};
inline double *P_interior::cell_varWND(const Grid* gitter,int lev) const {
  return gitter->Give_cell_variable(ind.next_WND(lev+1));
};
inline double *P_interior::cell_varWST(const Grid* gitter,int lev) const {
  return gitter->Give_cell_variable(ind.next_WST(lev+1));
};
inline double *P_interior::cell_varWSD(const Grid* gitter,int lev) const {
  return gitter->Give_cell_variable(ind.next_WSD(lev+1));
};


// finer level
// --------------
// idem:
inline double *P_interior::var_fine_M(const Grid* gitter, int level) const {
  return pointer_fine_varM;
};
// next:
inline double *P_interior::var_fine_N(const Grid* gitter, int level)
     const {
  return pointer_fine_varN;
};
inline double *P_interior::var_fine_S(const Grid* gitter, int level)
     const {
  return pointer_fine_varS;
};
inline double *P_interior::var_fine_E(const Grid* gitter, int level)
     const {
  return pointer_fine_varE;
};
inline double *P_interior::var_fine_W(const Grid* gitter, int level)
     const {
  return pointer_fine_varW;
};
inline double *P_interior::var_fine_T(const Grid* gitter, int level)
     const {
  return pointer_fine_varT;
};
inline double *P_interior::var_fine_D(const Grid* gitter, int level)
     const {
  return pointer_fine_varD;
};
// double next:
inline double *P_interior::var_fine_ND(const Grid* gitter, int level)
     const {
  return pointer_fine_varND;
};
inline double *P_interior::var_fine_WN(const Grid* gitter, int level)
     const {
  return pointer_fine_varWN;
};
inline double *P_interior::var_fine_WT(const Grid* gitter, int level)
     const {
  return pointer_fine_varWT;
};
inline double *P_interior::var_fine_ED(const Grid* gitter, int level)
     const {
  return pointer_fine_varED;
};
inline double *P_interior::var_fine_ST(const Grid* gitter, int level)
     const {
  return pointer_fine_varST;
};
inline double *P_interior::var_fine_ES(const Grid* gitter, int level)
     const {
  return pointer_fine_varES;
};
// third next:
inline double *P_interior::var_fine_EST(const Grid* gitter, int level)
     const {
  return pointer_fine_varEST;
};
inline double *P_interior::var_fine_WND(const Grid* gitter, int level)
     const {
  return pointer_fine_varWND;
};

// Coarse level
// --------------
// idem:
inline double *P_interior::var_coarse_M(const Grid* gitter, int level) const {
  return gitter->Give_variable(ind,level-1);
};
// next:
inline double *P_interior::var_coarse_N(const Grid* gitter, int level)
     const {
  return gitter->Give_variable(ind.next_NS(Rd,level),level-1);
};
inline double *P_interior::var_coarse_S(const Grid* gitter, int level)
     const {
  return gitter->Give_variable(ind.next_NS(Ld,level),level-1);
};
inline double *P_interior::var_coarse_E(const Grid* gitter, int level)
     const {
  return gitter->Give_variable(ind.next_EW(Rd,level),level-1);
};
inline double *P_interior::var_coarse_W(const Grid* gitter, int level)
     const {
  return gitter->Give_variable(ind.next_EW(Ld,level),level-1);
};
inline double *P_interior::var_coarse_T(const Grid* gitter, int level)
     const {
  return gitter->Give_variable(ind.next_TD(Rd,level),level-1);
};
inline double *P_interior::var_coarse_D(const Grid* gitter, int level)
     const {
  return gitter->Give_variable(ind.next_TD(Ld,level),level-1);
};
// double next:
inline double *P_interior::var_coarse_ND(const Grid* gitter, int level)
     const {
  return gitter->Give_variable(ind.next_ND(level),level-1);
};
inline double *P_interior::var_coarse_WN(const Grid* gitter, int level)
     const {
  return gitter->Give_variable(ind.next_WN(level),level-1);
};
inline double *P_interior::var_coarse_WT(const Grid* gitter, int level)
     const {
  return gitter->Give_variable(ind.next_WT(level),level-1);
};
inline double *P_interior::var_coarse_ED(const Grid* gitter, int level)
     const {
  return gitter->Give_variable(ind.next_ED(level),level-1);
};
inline double *P_interior::var_coarse_ST(const Grid* gitter, int level)
     const {
  return gitter->Give_variable(ind.next_ST(level),level-1);
};
inline double *P_interior::var_coarse_ES(const Grid* gitter, int level)
     const {
  return gitter->Give_variable(ind.next_ES(level),level-1);
};
// third next:
inline double *P_interior::var_coarse_EST(const Grid* gitter, int level)
     const {
  return gitter->Give_variable(ind.next_EST(level),level-1);
};
inline double *P_interior::var_coarse_WND(const Grid* gitter, int level)
     const {
  return gitter->Give_variable(ind.next_WND(level),level-1);
};

---