src/expde/grid/listbo/listbo1.h

Go to the documentation of this file.

//    expde: expression templates for partial differential equations.
//    Copyright (C) 2000  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.

// ------------------------------------------------------------
// listbo1.h
//
// ------------------------------------------------------------


#define initialize_bo true

// boundary points
class P_Bo2p {
 public:
  P_Bo2p(Index3D I, dir_3D d, P_Bo2p* Next) : 
    ind(I), direction(d), next(Next) {};

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

  // idem
  double *varM(const Grid* gitter) const;
  double *varP_slow(const Grid* gitter, int ebene, Ort1D i,Ort1D j,Ort1D k,
                    int level) const;


  // physical coordinates
  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_Bo2p* Next() { return next; };

  // give informationen
  dir_3D d() const { return direction; }
  Index3D Ind() const { return ind; }

  BoCeData* Give_BoData(dir_sons i) const;

  // label for boundary points
  bool Give_Label(int num, const Grid* gitter) const;
  void Put_Label(bool lab, int num, const Grid* gitter) const;
 private:
  Index3D ind;
  dir_3D  direction;
  P_Bo2p *next;

  double* pointer_varM;
  BoCeData* bocedatas[8]; 
};


inline void P_Bo2p::initialize(const Grid* grid, int l) {
  pointer_varM = grid->Give_variable(ind,direction);

  if(direction==Ndir) {
    bocedatas[WNDd]=grid->Give_Bo_cell(ind.next(WNDd,l+1));
    bocedatas[ENDd]=grid->Give_Bo_cell(ind.next(ENDd,l+1));
    bocedatas[WNTd]=grid->Give_Bo_cell(ind.next(WNTd,l+1));
    bocedatas[ENTd]=grid->Give_Bo_cell(ind.next(ENTd,l+1));

    bocedatas[WSDd]=NULL;
    bocedatas[ESDd]=NULL;
    bocedatas[WSTd]=NULL;
    bocedatas[ESTd]=NULL;
  }
  else if(direction==Sdir) {
    bocedatas[WSDd]=grid->Give_Bo_cell(ind.next(WSDd,l+1));
    bocedatas[ESDd]=grid->Give_Bo_cell(ind.next(ESDd,l+1));
    bocedatas[WSTd]=grid->Give_Bo_cell(ind.next(WSTd,l+1));
    bocedatas[ESTd]=grid->Give_Bo_cell(ind.next(ESTd,l+1));

    bocedatas[WNDd]=NULL; 
    bocedatas[ENDd]=NULL;
    bocedatas[WNTd]=NULL;
    bocedatas[ENTd]=NULL;
  }
  else if(direction==Wdir) {
    bocedatas[WSDd]=grid->Give_Bo_cell(ind.next(WSDd,l+1));
    bocedatas[WNDd]=grid->Give_Bo_cell(ind.next(WNDd,l+1));
    bocedatas[WSTd]=grid->Give_Bo_cell(ind.next(WSTd,l+1));
    bocedatas[WNTd]=grid->Give_Bo_cell(ind.next(WNTd,l+1));

    bocedatas[ESDd]=NULL; 
    bocedatas[ENDd]=NULL;
    bocedatas[ESTd]=NULL;
    bocedatas[ENTd]=NULL;
  }
  else if(direction==Edir) {
    bocedatas[ESDd]=grid->Give_Bo_cell(ind.next(ESDd,l+1));
    bocedatas[ENDd]=grid->Give_Bo_cell(ind.next(ENDd,l+1));
    bocedatas[ESTd]=grid->Give_Bo_cell(ind.next(ESTd,l+1));
    bocedatas[ENTd]=grid->Give_Bo_cell(ind.next(ENTd,l+1));

    bocedatas[WSDd]=NULL; 
    bocedatas[WNDd]=NULL;
    bocedatas[WSTd]=NULL;
    bocedatas[WNTd]=NULL;
  }
  else if(direction==Tdir) {
    bocedatas[WSTd]=grid->Give_Bo_cell(ind.next(WSTd,l+1));
    bocedatas[WNTd]=grid->Give_Bo_cell(ind.next(WNTd,l+1));
    bocedatas[ESTd]=grid->Give_Bo_cell(ind.next(ESTd,l+1));
    bocedatas[ENTd]=grid->Give_Bo_cell(ind.next(ENTd,l+1));

    bocedatas[WSDd]=NULL;
    bocedatas[WNDd]=NULL;
    bocedatas[ESDd]=NULL;
    bocedatas[ENDd]=NULL;
  }
  else {
    bocedatas[WSDd]=grid->Give_Bo_cell(ind.next(WSDd,l+1));
    bocedatas[WNDd]=grid->Give_Bo_cell(ind.next(WNDd,l+1));
    bocedatas[ESDd]=grid->Give_Bo_cell(ind.next(ESDd,l+1));
    bocedatas[ENDd]=grid->Give_Bo_cell(ind.next(ENDd,l+1));

    bocedatas[WSTd]=NULL;
    bocedatas[WNTd]=NULL;
    bocedatas[ESTd]=NULL;
    bocedatas[ENTd]=NULL;
  }
}


inline D3vector P_Bo2p::coordinate(const Grid* gitter) const {
  return gitter->transform_coord(ind,direction);
};
inline double P_Bo2p::coordinateX(const Grid* gitter) const {
  return gitter->transform_coord(ind,direction).x;
};
inline double P_Bo2p::coordinateY(const Grid* gitter) const {
  return gitter->transform_coord(ind,direction).y;
};
inline double P_Bo2p::coordinateZ(const Grid* gitter) const {
  return gitter->transform_coord(ind,direction).z;
};

inline bool P_Bo2p::Give_Label(int num, const Grid* gitter) const {
  return gitter->Give_label_bo(ind,direction,num);
};
inline void P_Bo2p::Put_Label(bool lab, int num, const Grid* gitter) const {
  gitter->Put_label_bo(lab,ind,direction,num);
};



inline double *P_Bo2p::varM(const Grid* gitter) const { 
  return pointer_varM;
};


inline double *P_Bo2p::varP_slow(const Grid* gitter,int var_ebene,
                             Ort1D i, Ort1D j, Ort1D k, int level) const {
  return gitter->Give_variable_slow(ind.next(i,j,k,level),var_ebene);
};

inline BoCeData *P_Bo2p::Give_BoData(dir_sons i) const {
  return bocedatas[i];
};

---