mappingorthomesh.h

Go to the documentation of this file.

#ifndef  _MY_MappingOrthoMesh_
#define _MY_MappingOrthoMesh_
#include "orthomesh.h"


class MappingOrthoMesh 
{
 protected:
  Point3D _X0;
  double _j;
  Matrix _invGradT,_GradT;
  double _dx,_dy,_dz;
  const OrthoMesh *_pMesh;
 public:
  

  void reinit(const OrthoMesh &mesh) 
  {
    _j=mesh.cellVolume();
    _invGradT.reinit(3,3);
    _GradT.reinit(3,3);
    _invGradT=0.0;
    _GradT=0.0;

    _dx=mesh.get_dx();
    _dy=mesh.get_dy();
    _dz=mesh.get_dz();
    
    _invGradT(0,0)=1.0/_dx;
    _invGradT(1,1)=1.0/_dy;
    _invGradT(2,2)=1.0/_dz;

    _GradT(0,0)=_dx;
    _GradT(1,1)=_dy;
    _GradT(2,2)=_dz;
    _pMesh = &mesh;
  }



  void reinit(OrthoMesh::Cell_It &cell) 
  {
    //Assert that the cell has as mesh,
    //the mesh passed in the reinit(OrthoMesh&) method
    assert(&(cell->getMesh()) == _pMesh); 
    cell->vertex(VERTEX_000,_X0);
  }

virtual double J(Point3D &p)
  {
    return _j;
  }


virtual void  Inv_T_Grad(const VecDouble &xi,Matrix &InvTGrad)
  {
    assert(xi.size() ==3);
    InvTGrad = _invGradT;
  }


  virtual void T_grad(const VecDouble &xi,Matrix &gradT)
  {
    assert(xi.size() == 3);
    gradT=_GradT;
  }


  virtual void T(const VecDouble &p,VecDouble &resp)
  {
    assert(p.size() == 3);
    assert(resp.size() == 3);
    resp(0)=_X0[0] + _dx*p(0);
    resp(1)=_X0[1] + _dy*p(1);
    resp(2)=_X0[2] + _dz*p(2);

  }

};

#endif

---