vecdoubleghost.h

Go to the documentation of this file.

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

class GhostBC
{
  OrthoMesh &m_mesh;
  std::vector<int> m_vIndices;
  unsigned _off;
  unsigned _numBC;
  GhostBC(GhostBC &ghost);
  GhostBC&  operator = (GhostBC &ghost);
public:
  GhostBC(OrthoMesh &mesh)
    :m_mesh(mesh),m_vIndices(numGhostCells(),0);
  {
    _off = m_mesh.numCells();m_vIndices[ghostIndex-_off];
  }
  bool isGhostIndex(Index ghostIndex);
  bool isMirrorBC(Index ghostIndex);
  Index getPrescribedBCIndex(Index ghostIndex);
  void setGhostCellAsMirror(unsigned ghostIndex,unsigned cellIndex);
  void setGhostBC(Function3D &f);
  void setGhostBCVector(Function3D &f,VecDouble &v);
  

};



};

  
  void setGhostCellAsValue(unsigned ghostIndex)
  {
    assert(ghostIndex >= m_mesh.numCells());
    unsigned index = ghostIndex-m_mesh.numCells();
    m_vIndices[index]=

  }
    
}
class VecDoubleGhost 
{
 protected:
  std::vector<int> m_vIndices;
  std::vector<double> m_vFBC;
  VecDouble *m_vRealData;
 public:
  VecDoubleGhost(OrthoMesh &mesh)
    :m_mesh(mesh)
  {
    m_vRealData=new VecDouble(mesh.numCells());
  }
  
 VecDoubleGhost(VecDouble &v,OrthoMesh &mesh)
    :m_mesh(mesh)
  {
    assert(v.size() == mesh.numRawCells());
    m_vRealData=&v;
  }

  void setGhostCellAsMirror(unsigned ghostIndex,unsigned cellIndex)
  {
    assert(ghostIndex >= m_mesh.numCells());
    assert(cellIndex <  m_mesh.numCells());
    unsigned index = ghostIndex-m_mesh.numCells();
    assert(index < m_vIndices.size());
    m_vIndices[index]=cellIndex;
  }
  
  

  void setGhostCells(Function3D &f)
  {
    OrthoMesh::Raw_Face_It face = m_mesh.begin_raw_face();
      
  }
                



  ~VecDoubleGhost();

};

#endif

---