src/expde/grid/listcf/listcf2.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
//

// ------------------------------------------------------------
// listcf2.h
//
// ------------------------------------------------------------

#define initialize_cellpoi true


// points of Randzellfreiheitsgrad
class P_cellpoi {
 public:
  P_cellpoi(Index3D I, P_cellpoi* Next) : 
    ind(I), next(Next) {};

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

  // idem
  double *varM(const Grid* gitter) 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;

  // index of cell
  Index3D Ind() const { return ind; };

  // for iteration:
  P_cellpoi* Next() { return next; };
  
 private:
  Index3D ind;
  P_cellpoi *next;

  double* pointer_varM;
};


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


inline D3vector P_cellpoi::coordinate(const Grid* gitter) const {
  return gitter->transform_coord_cellpoi(ind);
};
inline double P_cellpoi::coordinateX(const Grid* gitter) const {
  return gitter->transform_coordX_cellpoi(ind);
};
inline double P_cellpoi::coordinateY(const Grid* gitter) const {
  return gitter->transform_coordY_cellpoi(ind);
};
inline double P_cellpoi::coordinateZ(const Grid* gitter) const {
  return gitter->transform_coordZ_cellpoi(ind);
};

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

---