src/expde/grid/listce/listce1.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
//

// ------------------------------------------------------------
// listce1.h
//
// ------------------------------------------------------------


// interior cells

class P_interior_cell {
 public:
  P_interior_cell(Index3D I, P_interior_cell* Next, double* var) : 
    pointer_varM(var), ind(I), next(Next) {};

  // idem
  double *varM(const Grid* gitter) const;

  // corner point variable:
  double *varWSD(const Grid* gitter,int lev) const;
  double *varWST(const Grid* gitter,int lev) const;
  double *varWND(const Grid* gitter,int lev) const;
  double *varWNT(const Grid* gitter,int lev) const;
  double *varESD(const Grid* gitter,int lev) const;
  double *varEST(const Grid* gitter,int lev) const;
  double *varEND(const Grid* gitter,int lev) const;
  double *varENT(const Grid* gitter,int lev) 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_interior_cell* Next() { return next; };

  ~P_interior_cell() {
    if(next!=NULL)         delete next;
    if(pointer_varM!=NULL) delete pointer_varM;
  }

 private:
  double* pointer_varM;

  Index3D ind;
  P_interior_cell *next;  
};


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

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

inline double *P_interior_cell::varENT(const Grid* gitter,int lev) const {
  return gitter->Give_variable(ind.neighbour(ENTd),lev);
};
inline double *P_interior_cell::varEND(const Grid* gitter,int lev) const {
  return gitter->Give_variable(ind.neighbour(ENDd),lev);
};
inline double *P_interior_cell::varEST(const Grid* gitter,int lev) const {
  return gitter->Give_variable(ind.neighbour(ESTd),lev);
};
inline double *P_interior_cell::varESD(const Grid* gitter,int lev) const {
  return gitter->Give_variable(ind.neighbour(ESDd),lev);
};
inline double *P_interior_cell::varWNT(const Grid* gitter,int lev) const {
  return gitter->Give_variable(ind.neighbour(WNTd),lev);
};
inline double *P_interior_cell::varWND(const Grid* gitter,int lev) const {
  return gitter->Give_variable(ind.neighbour(WNDd),lev);
};
inline double *P_interior_cell::varWST(const Grid* gitter,int lev) const {
  return gitter->Give_variable(ind.neighbour(WSTd),lev);
};
inline double *P_interior_cell::varWSD(const Grid* gitter,int lev) const {
  return gitter->Give_variable(ind.neighbour(WSDd),lev);
};

// boundary tets

class P_boundary_tet {
 public:
  P_boundary_tet(D3vector coord, P_boundary_tet* Next, double* var,
                 Tetraeder_storage* ts, BoCell *bc) : 
    pointer_varM(var), tet(ts), in_boundary_cell(bc),
    middle_points_coord(coord), next(Next) {};

  // idem
  double *varM(const Grid* gitter) const; // middle of tet

  double *var0(const Grid* gitter) const; // corner of tet
  double *var1(const Grid* gitter) const; // corner of tet
  double *var2(const Grid* gitter) const; // corner of tet
  double *var3(const Grid* gitter) const; // corner of tet

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

  ~P_boundary_tet() {
    if(next!=NULL)         delete next;
    if(pointer_varM!=NULL) delete pointer_varM;
  }
 private:
  double* pointer_varM;      // middle of tet

  Tetraeder_storage* tet;    // pointer to tet
  BoCell *in_boundary_cell;  // boundary cell in which tet is contained

  D3vector middle_points_coord;
  //  Index3D ind;
  P_boundary_tet *next;  
};


// coordinates
inline D3vector P_boundary_tet::coordinate(const Grid* gitter) const {
  return middle_points_coord;
};
inline double P_boundary_tet::coordinateX(const Grid* gitter) const {
  return middle_points_coord.x;
};
inline double P_boundary_tet::coordinateY(const Grid* gitter) const {
  return middle_points_coord.y;
};
inline double P_boundary_tet::coordinateZ(const Grid* gitter) const {
  return middle_points_coord.z;
};

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

inline double *P_boundary_tet::var0(const Grid* gitter) const {
  return in_boundary_cell->vars[tet->N0()];
};
inline double *P_boundary_tet::var1(const Grid* gitter) const {
  return in_boundary_cell->vars[tet->N1()];
};
inline double *P_boundary_tet::var2(const Grid* gitter) const {
  return in_boundary_cell->vars[tet->N2()];
};
inline double *P_boundary_tet::var3(const Grid* gitter) const {
  return in_boundary_cell->vars[tet->N3()];
};

---