src/expde/grid/listbo/listbo2.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
//

// ------------------------------------------------------------
// listbo2.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; }

  // 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;
};


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


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);
};

---