grid_impl.hpp

Go to the documentation of this file.

/*
 * Copyright (c) 2009-2015:  G-CSC, Goethe University Frankfurt
 * Author: Sebastian Reiter
 * 
 * This file is part of UG4.
 * 
 * UG4 is free software: you can redistribute it and/or modify it under the
 * terms of the GNU Lesser General Public License version 3 (as published by the
 * Free Software Foundation) with the following additional attribution
 * requirements (according to LGPL/GPL v3 §7):
 * 
 * (1) The following notice must be displayed in the Appropriate Legal Notices
 * of covered and combined works: "Based on UG4 (www.ug4.org/license)".
 * 
 * (2) The following notice must be displayed at a prominent place in the
 * terminal output of covered works: "Based on UG4 (www.ug4.org/license)".
 * 
 * (3) The following bibliography is recommended for citation and must be
 * preserved in all covered files:
 * "Reiter, S., Vogel, A., Heppner, I., Rupp, M., and Wittum, G. A massively
 *   parallel geometric multigrid solver on hierarchically distributed grids.
 *   Computing and visualization in science 16, 4 (2013), 151-164"
 * "Vogel, A., Reiter, S., Rupp, M., Nägel, A., and Wittum, G. UG4 -- a novel
 *   flexible software system for simulating pde based models on high performance
 *   computers. Computing and visualization in science 16, 4 (2013), 165-179"
 * 
 * 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 Lesser General Public License for more details.
 */
 
#ifndef __H__LIB_GRID__GRID_IMPLEMENTATION__
#define __H__LIB_GRID__GRID_IMPLEMENTATION__
 
//#include <cassert>
#include "common/common.h"
#include "common/static_assert.h"
#include "grid_util.h"
#include "grid.h"
 
namespace ug
{
//  parallelism
bool Grid::
is_parallel() const
{
  return m_distGridMgr != NULL;
}
 
DistributedGridManager* Grid::
distributed_grid_manager()
{
  return m_distGridMgr;
}
 
const DistributedGridManager* Grid::
distributed_grid_manager() const
{
  return m_distGridMgr;
}
 
 
//  create functions
template<class TGeomObj>
typename geometry_traits<TGeomObj>::iterator
Grid::create(GridObject* pParent)
{
  STATIC_ASSERT(geometry_traits<TGeomObj>::CONTAINER_SECTION != -1
    &&  geometry_traits<TGeomObj>::BASE_OBJECT_ID != -1,
    invalid_geometry_type);
 
  TGeomObj* geomObj = new TGeomObj;
//  int baseObjectType = geometry_traits<GeomObjType>::base_object_type();
//  geomObj->m_elemHandle = m_elementStorage[baseObjectType].m_sectionContainer.insert_element(geomObj, geometry_traits<GeomObjType>::container_section());
//  m_elementStorage[baseObjectType].m_attachmentPipe.register_element(geomObj);
 
  register_element(geomObj, pParent);
 
  return iterator_cast<typename geometry_traits<TGeomObj>::iterator>(get_iterator(geomObj));
}
 
template <class TGeomObj>
typename geometry_traits<TGeomObj>::iterator
Grid::create(const typename geometry_traits<TGeomObj>::Descriptor& descriptor,
      GridObject* pParent)
{
  STATIC_ASSERT(geometry_traits<TGeomObj>::CONTAINER_SECTION != -1
      &&  geometry_traits<TGeomObj>::BASE_OBJECT_ID != -1,
      invalid_geometry_type);
 
  TGeomObj* geomObj = new TGeomObj(descriptor);
 
//  int baseObjectType = geometry_traits<TGeomObj>::base_object_type();
//  geomObj->m_elemHandle = m_elementStorage[baseObjectType].m_sectionContainer.insert_element(geomObj, geometry_traits<GeomObjType>::container_section());
//  m_elementStorage[baseObjectType].m_attachmentPipe.register_element(geomObj);
 
  register_element(geomObj, pParent);
 
  return iterator_cast<typename geometry_traits<TGeomObj>::iterator>(get_iterator(geomObj));
}
 
template<class TGeomObj>
typename geometry_traits<TGeomObj>::iterator
Grid::create_and_replace(typename geometry_traits<TGeomObj>::grid_base_object* pReplaceMe)
{
  STATIC_ASSERT(geometry_traits<TGeomObj>::CONTAINER_SECTION != -1
    &&  geometry_traits<TGeomObj>::BASE_OBJECT_ID != -1,
    invalid_geometry_type);
 
  TGeomObj* geomObj = new TGeomObj;
 
  if(geomObj->reference_object_id() == pReplaceMe->reference_object_id())
  {
    register_and_replace_element(geomObj, pReplaceMe);
    return iterator_cast<typename geometry_traits<TGeomObj>::iterator>(get_iterator(geomObj));
  }
  else
  {
    LOG("ERROR in Grid::create_and_replace(...): reference objects do not match!");
    assert(!"ERROR in Grid::create_and_replace(...): reference objects do not match!");
    delete geomObj;
    return end<TGeomObj>();
  }
}
 
template <class TGeomObj>
void Grid::reserve(size_t num)
{
  STATIC_ASSERT(geometry_traits<TGeomObj>::BASE_OBJECT_ID != -1,
        invalid_geometry_type);
 
  element_storage<TGeomObj>().m_attachmentPipe.reserve(num);
}
 
//  erase
template <class GeomObjIter>
void Grid::erase(const GeomObjIter& iterBegin, const GeomObjIter& iterEnd)
{
  GeomObjIter iter = iterBegin;
  while(iter != iterEnd)
  {
    GeomObjIter tmpIter = iter;
    ++iter;
    erase(*tmpIter);
  }
}
 
template <class TGeomObj>
void Grid::clear()
{
  while(begin<TGeomObj>() != end<TGeomObj>())
    erase(*begin<TGeomObj>());
}
 
//  Iterators
template <class TGeomObj>
typename geometry_traits<TGeomObj>::iterator
Grid::begin()
{
  STATIC_ASSERT(geometry_traits<TGeomObj>::BASE_OBJECT_ID != -1,
    invalid_GeomObj);
 
  return iterator_cast<typename geometry_traits<TGeomObj>::iterator>
    (element_storage<TGeomObj>().m_sectionContainer.section_begin(geometry_traits<TGeomObj>::CONTAINER_SECTION));
}
 
template <class TGeomObj>
typename geometry_traits<TGeomObj>::iterator
Grid::end()
{
  STATIC_ASSERT(geometry_traits<TGeomObj>::BASE_OBJECT_ID != -1,
    invalid_GeomObj);
 
  return iterator_cast<typename geometry_traits<TGeomObj>::iterator>
    (element_storage<TGeomObj>().m_sectionContainer.section_end(geometry_traits<TGeomObj>::CONTAINER_SECTION));
}
 
template <class TGeomObj>
typename geometry_traits<TGeomObj>::const_iterator
Grid::begin() const
{
  STATIC_ASSERT(geometry_traits<TGeomObj>::BASE_OBJECT_ID != -1,
    invalid_GeomObj);
 
  return iterator_cast<typename geometry_traits<TGeomObj>::const_iterator>
    (element_storage<TGeomObj>().m_sectionContainer.section_begin(geometry_traits<TGeomObj>::CONTAINER_SECTION));
}
 
template <class TGeomObj>
typename geometry_traits<TGeomObj>::const_iterator
Grid::end() const
{
  STATIC_ASSERT(geometry_traits<TGeomObj>::BASE_OBJECT_ID != -1,
    invalid_GeomObj);
 
  return iterator_cast<typename geometry_traits<TGeomObj>::const_iterator>
    (element_storage<TGeomObj>().m_sectionContainer.section_end(geometry_traits<TGeomObj>::CONTAINER_SECTION));
}
 
template <class TGeomObj>
TGeomObj*
Grid::front()
{
  STATIC_ASSERT(geometry_traits<TGeomObj>::BASE_OBJECT_ID != -1,
    invalid_GeomObj);
 
  return static_cast<TGeomObj*>(*element_storage<TGeomObj>().m_sectionContainer.
                    front(geometry_traits<TGeomObj>::CONTAINER_SECTION));
}
 
template <class TGeomObj>
TGeomObj*
Grid::back()
{
  STATIC_ASSERT(geometry_traits<TGeomObj>::BASE_OBJECT_ID != -1,
    invalid_GeomObj);
 
  return static_cast<TGeomObj*>(*element_storage<TGeomObj>().m_sectionContainer.
                    back(geometry_traits<TGeomObj>::CONTAINER_SECTION));
}
//  element numbers
template <class TGeomObj>
size_t Grid::num() const
{
  STATIC_ASSERT(geometry_traits<TGeomObj>::BASE_OBJECT_ID != -1,
    invalid_GeomObj);
 
  int secIndex = geometry_traits<TGeomObj>::CONTAINER_SECTION;
 
  if(secIndex == -1)
    return element_storage<TGeomObj>().m_sectionContainer.num_elements();
 
  return element_storage<TGeomObj>().m_sectionContainer.num_elements(secIndex);
}
 
inline void Grid::
objects_will_be_merged(Vertex* target, Vertex* elem1,
            Vertex* elem2)
{
  for(Grid::ObserverContainer::iterator iter = m_vertexObservers.begin();
    iter != m_vertexObservers.end(); iter++)
  {
    (*iter)->vertices_to_be_merged(this, target, elem1, elem2);
  }
}
 
inline void Grid::
objects_will_be_merged(Edge* target, Edge* elem1,
            Edge* elem2)
{
  for(Grid::ObserverContainer::iterator iter = m_edgeObservers.begin();
    iter != m_edgeObservers.end(); iter++)
  {
    (*iter)->edges_to_be_merged(this, target, elem1, elem2);
  }
}
 
inline void Grid::
objects_will_be_merged(Face* target, Face* elem1,
            Face* elem2)
{
  for(Grid::ObserverContainer::iterator iter = m_faceObservers.begin();
    iter != m_faceObservers.end(); iter++)
  {
    (*iter)->faces_to_be_merged(this, target, elem1, elem2);
  }
}
 
inline void Grid::
objects_will_be_merged(Volume* target, Volume* elem1,
            Volume* elem2)
{
  for(Grid::ObserverContainer::iterator iter = m_volumeObservers.begin();
    iter != m_volumeObservers.end(); iter++)
  {
    (*iter)->volumes_to_be_merged(this, target, elem1, elem2);
  }
}
 
template <class TGeomObj>
size_t Grid::attachment_container_size() const
{
  return element_storage<TGeomObj>().m_attachmentPipe.num_data_entries();
}
 
//  attachment handling
template <class TGeomObjClass>
void Grid::attach_to(IAttachment& attachment, bool passOnValues)
{
  STATIC_ASSERT(geometry_traits<TGeomObjClass>::BASE_OBJECT_ID != -1,
      invalid_GeomObjClass);
 
//  setup the options for this attachment.
  int options = 0;
  if(passOnValues)
    options = 1;
 
  element_storage<TGeomObjClass>().m_attachmentPipe.attach(attachment, options);
}
 
inline void Grid::attach_to_all(IAttachment& attachment, bool passOnValues)
{
  attach_to<Vertex>(attachment, passOnValues);
  attach_to<Edge>(attachment, passOnValues);
  attach_to<Face>(attachment, passOnValues);
  attach_to<Volume>(attachment, passOnValues);
}
 
inline void Grid::attach_to_all(IAttachment& attachment)
{
  attach_to<Vertex>(attachment);
  attach_to<Edge>(attachment);
  attach_to<Face>(attachment);
  attach_to<Volume>(attachment);
}
 
template <class TGeomObjClass, class TAttachment>
void Grid::attach_to_dv(TAttachment& attachment, const typename TAttachment::ValueType& defaultValue)
{
  attach_to_dv<TGeomObjClass, TAttachment>(attachment, defaultValue, attachment.default_pass_on_behaviour());
}
 
template <class TAttachment>
inline void Grid::
attach_to_all_dv(TAttachment& attachment,
         const typename TAttachment::ValueType& defaultValue)
{
  attach_to_dv<Vertex>(attachment, defaultValue);
  attach_to_dv<Edge>(attachment, defaultValue);
  attach_to_dv<Face>(attachment, defaultValue);
  attach_to_dv<Volume>(attachment, defaultValue);
}
 
template <class TGeomObjClass, class TAttachment>
void Grid::attach_to_dv(TAttachment& attachment, const typename TAttachment::ValueType& defaultValue, bool passOnValues)
{
  STATIC_ASSERT(geometry_traits<TGeomObjClass>::BASE_OBJECT_ID != -1,
      invalid_GeomObjClass);
 
//  setup the options for this attachment.
  int options = 0;
  if(passOnValues)
    options = 1;
 
  element_storage<TGeomObjClass>().m_attachmentPipe.attach(attachment, defaultValue, options);
}
 
template <class TAttachment>
inline void Grid::
attach_to_all_dv(TAttachment& attachment,
         const typename TAttachment::ValueType& defaultValue,
         bool passOnValues)
{
  attach_to_dv<Vertex>(attachment, defaultValue, passOnValues);
  attach_to_dv<Edge>(attachment, defaultValue, passOnValues);
  attach_to_dv<Face>(attachment, defaultValue, passOnValues);
  attach_to_dv<Volume>(attachment, defaultValue, passOnValues);
}
 
template <class TGeomObjClass>
void Grid::detach_from(IAttachment& attachment)
{
  STATIC_ASSERT(geometry_traits<TGeomObjClass>::BASE_OBJECT_ID != -1,
        invalid_GeomObjClass);
 
  element_storage<TGeomObjClass>().m_attachmentPipe.detach(attachment);
}
 
inline void Grid::detach_from_all(IAttachment& attachment)
{
  detach_from<Vertex>(attachment);
  detach_from<Edge>(attachment);
  detach_from<Face>(attachment);
  detach_from<Volume>(attachment);
}
 
/*
template <class TGeomObjClass>
util::IAttachmentDataContainer* Grid::get_data_container(util::IAttachment& attachment)
{
  assert(geometry_traits<TGeomObjClass>::base_object_type() != -1
      && "ERROR in Grid::get_data_container(...). Invalid base_object_type of GeomObjClass!");
 
  return m_elementStorage[geometry_traits<TGeomObjClass>::base_object_type()].
    m_attachmentPipe.get_data_container(attachment);
}
*/
 
template <class TGeomObj, class TAttachment>
typename TAttachment::ContainerType*
Grid::get_attachment_data_container(TAttachment& attachment)
{
  STATIC_ASSERT(geometry_traits<TGeomObj>::BASE_OBJECT_ID != -1,
      invalid_GeomObj);
 
  return element_storage<TGeomObj>().m_attachmentPipe.get_data_container(attachment);
}
 
template <class TGeomObj>
typename Grid::traits<TGeomObj>::AttachmentPipe&
Grid::get_attachment_pipe()
{
  STATIC_ASSERT(geometry_traits<TGeomObj>::BASE_OBJECT_ID != -1,
      invalid_GeomObj);
 
  return element_storage<TGeomObj>().m_attachmentPipe;
}
 
template <class TGeomObj>
uint
Grid::get_attachment_data_index(TGeomObj* pObj) const
{
  typedef typename geometry_traits<TGeomObj>::grid_base_object BaseObj;
  return attachment_traits<BaseObj*, ElementStorage<BaseObj> >::
      get_data_index(&element_storage<TGeomObj>(), pObj);
}
 
inline void
Grid::autoenable_option(uint option, const char* caller, const char* optionName)
{
  if(!option_is_enabled(option))
  {
    LOG("WARNING in " << caller << ": auto-enabling " << optionName << "." << std::endl);
    enable_options(option);
  }
}
 
//  marks
template <class TIterator>
void Grid::mark(TIterator begin, TIterator end)
{
  for(TIterator iter = begin; iter != end; ++iter)
    mark(*iter);
}
 
template <class TIterator>
void Grid::unmark(TIterator begin, TIterator end)
{
  for(TIterator iter = begin; iter != end; ++iter)
    unmark(*iter);
}
 
 
template <class TContainer>
void Grid::get_associated(TContainer& container, GridObject* o)
{
  switch(o->base_object_id()){
    case VERTEX: return get_associated(container, static_cast<Vertex*>(o));
    case EDGE: return get_associated(container, static_cast<Edge*>(o));
    case FACE: return get_associated(container, static_cast<Face*>(o));
    case VOLUME: return get_associated(container, static_cast<Volume*>(o));
  }
}
 
template <class TElem>
void Grid::associated_elements(traits<Vertex>::secure_container& elemsOut, TElem* e)
{
  get_associated(elemsOut, e);
}
 
template <class TElem>
void Grid::associated_elements(traits<Edge>::secure_container& elemsOut, TElem* e)
{
  get_associated(elemsOut, e);
}
 
template <class TElem>
void Grid::associated_elements(traits<Face>::secure_container& elemsOut, TElem* e)
{
  get_associated(elemsOut, e);
}
 
template <class TElem>
void Grid::associated_elements(traits<Volume>::secure_container& elemsOut, TElem* e)
{
  get_associated(elemsOut, e);
}
 
template <class TElem>
void Grid::get_associated(typename traits<typename TElem::grid_base_object>
              ::secure_container& elems, TElem* e)
{
//  we have to retrieve a valid pointer on the element pointer. The only way
//  to receive it, is to return the pointer to the entry in which e is stored in
//  the element storage.
  elems.set_external_array(
    element_storage<typename TElem::grid_base_object>().m_sectionContainer.
      get_container().get_pointer_to_element(e),
    1);
}
 
template <class TElem>
void Grid::associated_elements_sorted(traits<Vertex>::secure_container& elemsOut, TElem* e)
{
  get_associated_sorted(elemsOut, e);
}
 
template <class TElem>
void Grid::associated_elements_sorted(traits<Edge>::secure_container& elemsOut, TElem* e)
{
  get_associated_sorted(elemsOut, e);
}
 
template <class TElem>
void Grid::associated_elements_sorted(traits<Face>::secure_container& elemsOut, TElem* e)
{
  get_associated_sorted(elemsOut, e);
}
 
template <class TElem>
void Grid::associated_elements_sorted(traits<Volume>::secure_container& elemsOut, TElem* e)
{
  get_associated_sorted(elemsOut, e);
}
 
 
template <class TElem>
void Grid::get_associated_sorted(typename traits<typename TElem::grid_base_object>
                 ::secure_container& elems, TElem* e)
{
//  we have to retrieve a valid pointer on the element pointer. The only way
//  to receive it, is to return the pointer to the entry in which e is stored in
//  the element storage.
  elems.set_external_array(
    element_storage<typename TElem::grid_base_object>().m_sectionContainer.
      get_container().get_pointer_to_element(e),
    1);
}
 
 
//  neighbourhood access
template <class TGeomObj>
Edge* Grid::find_edge_in_associated_edges(TGeomObj* obj, const EdgeVertices& ev)
{
  GRID_PROFILE_FUNC();
 
  AssociatedEdgeIterator iterEnd = associated_edges_end(obj);
  for(AssociatedEdgeIterator iter = associated_edges_begin(obj);
    iter != iterEnd; ++iter)
  {
    Edge* e = *iter;
    if(CompareVertices(e, &ev))
      return e;
  }
  return NULL;
}
 
                    
template <class TGeomObj>
Face* Grid::find_face_in_associated_faces(TGeomObj* obj, const FaceVertices& fv)
{
  GRID_PROFILE_FUNC();
 
  unsigned long key = hash_key(&fv);
  AssociatedFaceIterator iterEnd = associated_faces_end(obj);
  for(AssociatedFaceIterator iter = associated_faces_begin(obj);
    iter != iterEnd; ++iter)
  {
    Face* f = *iter;
    if(key == hash_key(f))
    {
      if(CompareVertices(f, &fv))
        return f;
    }
  }
  
  return NULL;
}
                    
template <class TGeomObj>
Volume* Grid::find_volume_in_associated_volumes(TGeomObj* obj,
                        const VolumeVertices& vv)
{
  GRID_PROFILE_FUNC();
 
  unsigned long key = hash_key(&vv);
  AssociatedVolumeIterator iterEnd = associated_volumes_end(obj);
  for(AssociatedVolumeIterator iter = associated_volumes_begin(obj);
    iter != iterEnd; ++iter)
  {
    Volume* v = *iter;
    if(key == hash_key(v))
    {
      if(CompareVertices(v, &vv))
        return v;
    }
  }
  
  return NULL;
}
 
 
//  implementation of Grids AttachmentAccessors
 
//  AttachmentAccessor
template <class TElem, class TAttachment>
Grid::AttachmentAccessor<TElem, TAttachment>::
AttachmentAccessor() :
ug::AttachmentAccessor<TElem*, TAttachment, ElementStorage<TElem> >()
{
}
 
template <class TElem, class TAttachment>
Grid::AttachmentAccessor<TElem, TAttachment>::
AttachmentAccessor(const AttachmentAccessor& aa) :
ug::AttachmentAccessor<TElem*, TAttachment, ElementStorage<TElem> >(aa)
{
}
 
template <class TElem, class TAttachment>
Grid::AttachmentAccessor<TElem, TAttachment>::
AttachmentAccessor(Grid& grid, TAttachment& a) :
ug::AttachmentAccessor<typename TElem::grid_base_object*, TAttachment,
             typename traits<TElem>::ElementStorage>
  (grid.get_attachment_pipe<TElem>(), a)
{
}
 
template <class TElem, class TAttachment>
Grid::AttachmentAccessor<TElem, TAttachment>::
AttachmentAccessor(Grid& grid, TAttachment& a, bool autoAttach) :
ug::AttachmentAccessor<typename TElem::grid_base_object*, TAttachment,
             typename traits<TElem>::ElementStorage>()
{
  if(autoAttach){
    if(!grid.has_attachment<TElem>(a))
      grid.attach_to<TElem>(a);
  }
 
  access(grid, a);
}
 
//  VertexAttachmentAccessor
template <class TAttachment>
Grid::VertexAttachmentAccessor<TAttachment>::
VertexAttachmentAccessor() :
  Grid::AttachmentAccessor<Vertex, TAttachment>()
{
}
 
template <class TAttachment>
Grid::VertexAttachmentAccessor<TAttachment>::
VertexAttachmentAccessor(const VertexAttachmentAccessor& aa) :
  Grid::AttachmentAccessor<Vertex, TAttachment>(aa)
{
}
 
template <class TAttachment>
Grid::VertexAttachmentAccessor<TAttachment>::
VertexAttachmentAccessor(Grid& grid, TAttachment& a) :
  Grid::AttachmentAccessor<Vertex, TAttachment>(grid, a)
{
}
 
template <class TAttachment>
Grid::VertexAttachmentAccessor<TAttachment>::
VertexAttachmentAccessor(Grid& grid, TAttachment& a, bool autoAttach) :
  Grid::AttachmentAccessor<Vertex, TAttachment>(grid, a, autoAttach)
{
}
 
 
//  EdgeAttachmentAccessor
template <class TAttachment>
Grid::EdgeAttachmentAccessor<TAttachment>::
EdgeAttachmentAccessor() :
  Grid::AttachmentAccessor<Edge, TAttachment>()
{
}
 
template <class TAttachment>
Grid::EdgeAttachmentAccessor<TAttachment>::
EdgeAttachmentAccessor(const EdgeAttachmentAccessor& aa) :
  Grid::AttachmentAccessor<Edge, TAttachment>(aa)
{
}
 
template <class TAttachment>
Grid::EdgeAttachmentAccessor<TAttachment>::
EdgeAttachmentAccessor(Grid& grid, TAttachment& a) :
  Grid::AttachmentAccessor<Edge, TAttachment>(grid, a)
{
}
 
template <class TAttachment>
Grid::EdgeAttachmentAccessor<TAttachment>::
EdgeAttachmentAccessor(Grid& grid, TAttachment& a, bool autoAttach) :
  Grid::AttachmentAccessor<Edge, TAttachment>(grid, a, autoAttach)
{
}
 
 
//  FaceAttachmentAccessor
template <class TAttachment>
Grid::FaceAttachmentAccessor<TAttachment>::
FaceAttachmentAccessor() :
  Grid::AttachmentAccessor<Face, TAttachment>()
{
}
 
template <class TAttachment>
Grid::FaceAttachmentAccessor<TAttachment>::
FaceAttachmentAccessor(const FaceAttachmentAccessor& aa) :
  Grid::AttachmentAccessor<Face, TAttachment>(aa)
{
}
 
template <class TAttachment>
Grid::FaceAttachmentAccessor<TAttachment>::
FaceAttachmentAccessor(Grid& grid, TAttachment& a) :
  Grid::AttachmentAccessor<Face, TAttachment>(grid, a)
{
}
 
template <class TAttachment>
Grid::FaceAttachmentAccessor<TAttachment>::
FaceAttachmentAccessor(Grid& grid, TAttachment& a, bool autoAttach) :
  Grid::AttachmentAccessor<Face, TAttachment>(grid, a, autoAttach)
{
}
 
 
//  FaceAttachmentAccessor
template <class TAttachment>
Grid::VolumeAttachmentAccessor<TAttachment>::
VolumeAttachmentAccessor() :
  Grid::AttachmentAccessor<Volume, TAttachment>()
{
}
 
template <class TAttachment>
Grid::VolumeAttachmentAccessor<TAttachment>::
VolumeAttachmentAccessor(const VolumeAttachmentAccessor& aa) :
  Grid::AttachmentAccessor<Volume, TAttachment>(aa)
{
}
 
template <class TAttachment>
Grid::VolumeAttachmentAccessor<TAttachment>::
VolumeAttachmentAccessor(Grid& grid, TAttachment& a) :
  Grid::AttachmentAccessor<Volume, TAttachment>(grid, a)
{
}
 
template <class TAttachment>
Grid::VolumeAttachmentAccessor<TAttachment>::
VolumeAttachmentAccessor(Grid& grid, TAttachment& a, bool autoAttach) :
  Grid::AttachmentAccessor<Volume, TAttachment>(grid, a, autoAttach)
{
}
 
 
//  marks
inline void Grid::mark(GridObject* obj)
{
  const int typeID = obj->base_object_id();
  switch(typeID){
    case VERTEX: mark(static_cast<Vertex*>(obj)); break;
    case EDGE: mark(static_cast<Edge*>(obj)); break;
    case FACE: mark(static_cast<Face*>(obj)); break;
    case VOLUME: mark(static_cast<Volume*>(obj)); break;
  }
}
 
inline void Grid::mark(Vertex* obj)
{
  assert(m_bMarking && "ERROR: Grid::mark may only be called between calls to Grid::begin_marking and Grid::end_marking.");
  m_aaMarkVRT[obj] = m_currentMark;
}
 
inline void Grid::mark(Edge* obj)
{
  assert(m_bMarking && "ERROR: Grid::mark may only be called between calls to Grid::begin_marking and Grid::end_marking.");
  m_aaMarkEDGE[obj] = m_currentMark;
}
 
inline void Grid::mark(Face* obj)
{
  assert(m_bMarking && "ERROR: Grid::mark may only be called between calls to Grid::begin_marking and Grid::end_marking.");
  m_aaMarkFACE[obj] = m_currentMark;
}
 
inline void Grid::mark(Volume* obj)
{
  assert(m_bMarking && "ERROR: Grid::mark may only be called between calls to Grid::begin_marking and Grid::end_marking.");
  m_aaMarkVOL[obj] = m_currentMark;
}
 
inline void Grid::unmark(GridObject* obj)
{
  const int typeID = obj->base_object_id();
  switch(typeID){
    case VERTEX: unmark(static_cast<Vertex*>(obj)); break;
    case EDGE: unmark(static_cast<Edge*>(obj)); break;
    case FACE: unmark(static_cast<Face*>(obj)); break;
    case VOLUME: unmark(static_cast<Volume*>(obj)); break;
  }
}
 
inline void Grid::unmark(Vertex* obj)
{
  assert(m_bMarking && "ERROR: Grid::unmark may only be called between calls to Grid::begin_marking and Grid::end_marking.");
  m_aaMarkVRT[obj] = 0;
}
 
inline void Grid::unmark(Edge* obj)
{
  assert(m_bMarking && "ERROR: Grid::unmark may only be called between calls to Grid::begin_marking and Grid::end_marking.");
  m_aaMarkEDGE[obj] = 0;
}
 
inline void Grid::unmark(Face* obj)
{
  assert(m_bMarking && "ERROR: Grid::unmark may only be called between calls to Grid::begin_marking and Grid::end_marking.");
  m_aaMarkFACE[obj] = 0;
}
 
inline void Grid::unmark(Volume* obj)
{
  assert(m_bMarking && "ERROR: Grid::unmark may only be called between calls to Grid::begin_marking and Grid::end_marking.");
  m_aaMarkVOL[obj] = 0;
}
 
inline bool Grid::is_marked(GridObject* obj) const
{
  const int typeID = obj->base_object_id();
  switch(typeID){
    case VERTEX: return is_marked(static_cast<Vertex*>(obj));
    case EDGE: return is_marked(static_cast<Edge*>(obj));
    case FACE: return is_marked(static_cast<Face*>(obj));
    case VOLUME: return is_marked(static_cast<Volume*>(obj));
    default: return false;
  }
}
 
inline bool Grid::is_marked(Vertex* obj) const
{
  if(m_currentMark == 0)
    return false;
  return m_aaMarkVRT[obj] == m_currentMark;
}
 
inline bool Grid::is_marked(Edge* obj) const
{
  if(m_currentMark == 0)
    return false;
  return m_aaMarkEDGE[obj] == m_currentMark;
}
 
inline bool Grid::is_marked(Face* obj) const
{
  if(m_currentMark == 0)
    return false;
  return m_aaMarkFACE[obj] == m_currentMark;
}
 
inline bool Grid::is_marked(Volume* obj) const
{
  if(m_currentMark == 0)
    return false;
  return m_aaMarkVOL[obj] == m_currentMark;
}
 
}// end of namespace libGrid
#endif