orthocellaccessorwithholes.cpp

Go to the documentation of this file.

#include "orthomesh.h"


 OrthoCellAccessorWithHoles::OrthoCellAccessorWithHoles(const OrthoMesh &mesh,unsigned rawIndex) 
   :OrthoCellAccessor(mesh,rawIndex)
{
  _delIt = mesh._delCellLst.findLUEntry(rawIndex);
  _delItFcX=mesh._delFaceLst.begin();
  _delItFcY=_delItFcX;
  _delItFcZ=_delItFcX;
  _delItVert=mesh._delVertLst.begin();
  _newIndex = mesh._delCellLst.adjustI(_index,_delIt);
}

void  OrthoCellAccessorWithHoles::operator++(int) 
{
  OrthoCellAccessor &rawCell = *this;
  assert(rawCell.index() < _delIt->index());
  OrthoCellAccessor::operator++(0);
      

  if (rawCell.index() <  _delIt->index())
  {
    _newIndex = _delIt->adjust(rawCell.index());
    return;
  }
  else //we know that delIt->index() == rawCell.index()
  {
    do {
      OrthoCellAccessor::operator++(0);
      _delIt++;
        } while(rawCell.index() == _delIt->index());
      _newIndex = _delIt->adjust(rawCell.index());
  }
}






unsigned OrthoCellAccessorWithHoles::vertex_index(VertexDirection3D dir) 
{
  return getMesh()._delVertLst.adjust(OrthoCellAccessor::vertex_index(dir),_delItVert);
}


unsigned OrthoCellAccessorWithHoles::neighbor_index(unsigned int dir) const 
{
  unsigned i=OrthoCellAccessor::neighbor_index(dir);
  if (i != OrthoMesh::invalidIndex())
    return _pMesh->_delCellLst.adjustI(OrthoCellAccessor::neighbor_index(dir),_delIt);
  else
    return OrthoMesh::invalidIndex();
}



unsigned OrthoCellAccessorWithHoles::face_index(FaceDirection3D dir) 
{
  return _pMesh->getFaceIndexFromCell(dir,_i,_j,_k,_delItFcX,_delItFcY,_delItFcZ);
}

unsigned OrthoCellAccessorWithHoles::face_index_up() 
{
  return _pMesh->_delFaceLst.adjust(_i + _j*_pMesh->_fzYStride + (_k +1)*_pMesh->_fzZStride + _pMesh->_fzOff,_delItFcZ);
}

unsigned OrthoCellAccessorWithHoles::face_index_down() 
{
  return  _pMesh->_delFaceLst.adjust(_i + _j*_pMesh->_fzYStride + _k*_pMesh->_fzZStride + _pMesh->_fzOff,_delItFcZ);
}

unsigned OrthoCellAccessorWithHoles::face_index_left() 
{
  return _pMesh->_delFaceLst.adjust(_i+_j*_pMesh->_fxYStride + _k*_pMesh->_fxZStride,_delItFcX);
}

unsigned OrthoCellAccessorWithHoles::face_index_right() 
{
  return _pMesh->_delFaceLst.adjust(_i + 1  + _j*_pMesh->_fxYStride + _k*_pMesh->_fxZStride,_delItFcX);
}

unsigned OrthoCellAccessorWithHoles::face_index_front() 
{
  return _pMesh->_delFaceLst.adjust(_i + _j*_pMesh->_fyYStride + _k*_pMesh->_fyZStride + _pMesh->_fyOff,_delItFcY);

}

unsigned OrthoCellAccessorWithHoles::face_index_back() 
{
  return _pMesh->_delFaceLst.adjust(_i + (_j+1)*_pMesh->_fyYStride + _k*_pMesh->_fyZStride + _pMesh->_fyOff,_delItFcY);
}





bool OrthoCellAccessorWithHoles::at_boundary() 
{
  
  return (neighbor_index(LEFT_CELL  ) == OrthoMesh::invalidIndex() ||
          neighbor_index(RIGHT_CELL ) == OrthoMesh::invalidIndex() ||
          neighbor_index(UP_CELL    ) == OrthoMesh::invalidIndex() ||
          neighbor_index(BOTTOM_CELL) == OrthoMesh::invalidIndex() ||
          neighbor_index(FRONT_CELL ) == OrthoMesh::invalidIndex() ||
          neighbor_index(BACK_CELL  ) == OrthoMesh::invalidIndex() );
}






OrthoMesh::Face_It OrthoCellAccessorWithHoles::face(FaceDirection3D dir) 
{
  OrthoFaceAccessorWithHoles faceAcc(*this,dir);
  return OrthoMesh::Face_It(faceAcc);
}


OrthoMesh::Cell_It OrthoCellAccessorWithHoles::neighbor(CellDirection3D dir) 
{
  OrthoCellAccessorWithHoles cell = *this;
  if (cell.advance(dir))
    return OrthoMesh::Cell_It(cell);
  else
    return _pMesh->end_cell();
}


bool OrthoCellAccessorWithHoles::advance(CellDirection3D dir) 
{
  OrthoCellAccessorWithHoles backTrack = *this;
  if (OrthoCellAccessor::advance(dir))
  {
    //We advance the cell as there was no holes.
    //Now calculate the index in the mesh with holes
    _newIndex=_pMesh->_delCellLst.adjustI(raw_index(),_delIt);

    //If _newIndex is invalid it means the neighbor cell is not a valid
    //cell because it is inside the well. So backtrack and return false
    if (_newIndex == OrthoMesh::INVALID_INDEX)
    {
      *this=backTrack;
      return false;
    }
    else
      return true;
  }
  else //Neighbor is out of the domain
  {
    return false;
  }
}







unsigned OrthoCellAccessorWithHoles::advance(CellDirection3D dir,unsigned stps) 
{
  unsigned i;
  for (i=0;i<stps;i++)
  {
    if (!this->advance(dir))
    {
      return i;
    }
  }
  return i;
}


unsigned OrthoCellAccessorWithHoles::index_up() 
{
  assert(this->isValid());

  if (_k == _pMesh->_LastK)
    return OrthoMesh::INVALID_INDEX;
  else
    return _pMesh->_delCellLst.adjustI(_index+_pMesh->_ZStride,_delIt);

}

unsigned OrthoCellAccessorWithHoles::index_down() 
{
  assert(this->isValid());
  if (_k == 0)
    return OrthoMesh::INVALID_INDEX;
  else
     return _pMesh->_delCellLst.adjustI(_index-_pMesh->_ZStride,_delIt);

}

unsigned OrthoCellAccessorWithHoles::index_left() 
{
  assert(this->isValid());
  if (_i == 0)
    return OrthoMesh::INVALID_INDEX;
  else
    return _pMesh->_delCellLst.adjustI(_index-1,_delIt);
}

unsigned OrthoCellAccessorWithHoles::index_right() 
{
  assert(this->isValid());
  if (_i == _pMesh->_LastI)
    return OrthoMesh::INVALID_INDEX;
  else
    return  _pMesh->_delCellLst.adjustI(_index+1,_delIt);

}



unsigned OrthoCellAccessorWithHoles::index_front() 
{
  assert(this->isValid());
  if (_j == 0)
    return OrthoMesh::INVALID_INDEX;
  else
    return _pMesh->_delCellLst.adjustI(_index-_pMesh->_YStride,_delIt);
}




unsigned OrthoCellAccessorWithHoles::index_back() 
{
  assert(this->isValid());
  if (_j == _pMesh->_LastJ)
    return OrthoMesh::INVALID_INDEX;
  else
    return _pMesh->_delCellLst.adjustI(_index+_pMesh->_YStride,_delIt);

}

---