tria_object.h

// ---------------------------------------------------------------------
//
// Copyright (C) 1998 - 2018 by the deal.II authors
//
// This file is part of the deal.II library.
//
// The deal.II library is free software; you can use it, redistribute
// it, and/or modify it under the terms of the GNU Lesser General
// Public License as published by the Free Software Foundation; either
// version 2.1 of the License, or (at your option) any later version.
// The full text of the license can be found in the file LICENSE.md at
// the top level directory of deal.II.
//
// ---------------------------------------------------------------------

#ifndef dealii_tria_object_h
#define dealii_tria_object_h


#include <deal.II/base/config.h>

#include <deal.II/base/exceptions.h>
#include <deal.II/base/geometry_info.h>

DEAL_II_NAMESPACE_OPEN

namespace internal
{
  namespace TriangulationImplementation
  {
    template <int structdim>
    class TriaObject
    {
    public:
      static const unsigned int dimension = structdim;

      TriaObject();

      TriaObject(const int i0, const int i1);

      TriaObject(const int i0, const int i1, const int i2, const int i3);

      TriaObject(const int i0,
                 const int i1,
                 const int i2,
                 const int i3,
                 const int i4,
                 const int i5);


      int
      face(const unsigned int i) const;

      void
      set_face(const unsigned int i, const int index);

      static std::size_t
      memory_consumption();

      template <class Archive>
      void
      serialize(Archive &ar, const unsigned int version);

    protected:
      int faces[GeometryInfo<structdim>::faces_per_cell];
    };

    //----------------------------------------------------------------------//

    template <int structdim>
    inline TriaObject<structdim>::TriaObject()
    {
      for (unsigned int i = 0; i < GeometryInfo<structdim>::faces_per_cell; ++i)
        faces[i] = -1;
    }


    template <int structdim>
    inline TriaObject<structdim>::TriaObject(const int i0, const int i1)
    {
      Assert(structdim == 1, ExcImpossibleInDim(structdim));
      faces[0] = i0;
      faces[1] = i1;
    }


    template <int structdim>
    inline TriaObject<structdim>::TriaObject(const int i0,
                                             const int i1,
                                             const int i2,
                                             const int i3)
    {
      Assert(structdim == 2, ExcImpossibleInDim(structdim));
      faces[0] = i0;
      faces[1] = i1;
      faces[2] = i2;
      faces[3] = i3;
    }


    template <int structdim>
    inline TriaObject<structdim>::TriaObject(const int i0,
                                             const int i1,
                                             const int i2,
                                             const int i3,
                                             const int i4,
                                             const int i5)
    {
      Assert(structdim == 3, ExcImpossibleInDim(structdim));
      faces[0] = i0;
      faces[1] = i1;
      faces[2] = i2;
      faces[3] = i3;
      faces[4] = i4;
      faces[5] = i5;
    }


    template <int structdim>
    inline int
    TriaObject<structdim>::face(const unsigned int i) const
    {
      Assert(i < GeometryInfo<structdim>::faces_per_cell,
             ExcIndexRange(i, 0, GeometryInfo<structdim>::faces_per_cell));
      return faces[i];
    }



    template <int structdim>
    inline void
    TriaObject<structdim>::set_face(const unsigned int i, const int index)
    {
      Assert(i < GeometryInfo<structdim>::faces_per_cell,
             ExcIndexRange(i, 0, GeometryInfo<structdim>::faces_per_cell));
      faces[i] = index;
    }



    template <int structdim>
    inline std::size_t
    TriaObject<structdim>::memory_consumption()
    {
      return sizeof(TriaObject<structdim>);
    }


    template <int structdim>
    template <class Archive>
    void
    TriaObject<structdim>::serialize(Archive &ar, const unsigned int)
    {
      ar &faces;
    }


  } // namespace TriangulationImplementation
} // namespace internal


DEAL_II_NAMESPACE_CLOSE

#endif