tria_faces.h

Go to the documentation of this file.

// ------------------------------------------------------------------------
//
// SPDX-License-Identifier: LGPL-2.1-or-later
// Copyright (C) 2006 - 2023 by the deal.II authors
//
// This file is part of the deal.II library.
//
// Part of the source code is dual licensed under Apache-2.0 WITH
// LLVM-exception OR LGPL-2.1-or-later. Detailed license information
// governing the source code and code contributions can be found in
// LICENSE.md and CONTRIBUTING.md at the top level directory of deal.II.
//
// ------------------------------------------------------------------------
 
#ifndef dealii_tria_faces_h
#define dealii_tria_faces_h
 
#include <deal.II/base/config.h>
 
#include <deal.II/grid/reference_cell.h>
#include <deal.II/grid/tria_objects.h>
 
 
DEAL_II_NAMESPACE_OPEN
 
namespace internal
{
  namespace TriangulationImplementation
  {
    class TriaFaces
    {
    public:
      TriaFaces(const unsigned int dim);
 
      TriaFaces() = default;
 
      unsigned int dim;
 
      TriaObjects quads;
 
      std::vector<bool> quads_line_orientations;
 
      std::vector<bool> quad_is_quadrilateral;
 
      ReferenceCell
      get_quad_type(const std::size_t index) const;
 
      void
      set_quad_type(const std::size_t index, const ReferenceCell face_type);
 
      TriaObjects lines;
 
      std::size_t
      memory_consumption() const;
 
      template <class Archive>
      void
      serialize(Archive &ar, const unsigned int version);
    };
 
 
 
    inline ReferenceCell
    TriaFaces::get_quad_type(const std::size_t index) const
    {
      AssertIndexRange(index, quad_is_quadrilateral.size());
      return quad_is_quadrilateral[index] ? ReferenceCells::Quadrilateral :
                                            ReferenceCells::Triangle;
    }
 
 
 
    inline void
    TriaFaces::set_quad_type(const std::size_t   index,
                             const ReferenceCell face_type)
    {
      AssertIndexRange(index, quad_is_quadrilateral.size());
      Assert(face_type == ReferenceCells::Quadrilateral ||
               face_type == ReferenceCells::Triangle,
             ExcInternalError());
      if (face_type == ReferenceCells::Quadrilateral)
        quad_is_quadrilateral[index] = true;
      else
        quad_is_quadrilateral[index] = false;
    }
 
 
 
    template <class Archive>
    void
    TriaFaces::serialize(Archive &ar, const unsigned int)
    {
      ar                                        &dim;
      ar &quads &lines &quads_line_orientations &quad_is_quadrilateral;
    }
  } // namespace TriangulationImplementation
} // namespace internal
 
DEAL_II_NAMESPACE_CLOSE
 
#endif