tria_faces.h

Go to the documentation of this file.

// -----------------------------------------------------------------------------
//
// SPDX-License-Identifier: Apache-2.0 WITH LLVM-exception OR LGPL-2.1-or-later
// Copyright (C) 2006 - 2026 by the deal.II authors
//
// This file is part of the deal.II library.
//
// Detailed license information governing the source code and contributions
// can be found in LICENSE.md and CONTRIBUTING.md at the top level directory.
//
// -----------------------------------------------------------------------------
 
#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
  {
    template <int dim>
    class TriaFaces
    {
    public:
      TriaFaces(const unsigned int max_children_per_quad,
                const unsigned int max_lines_per_quad);
 
      TriaFaces() = default;
 
      void
      allocate(const std::size_t n_lines, const std::size_t n_quads);
 
      void
      allocate_end(const std::size_t new_quads_in_pairs,
                   const std::size_t new_quads_single);
 
      unsigned int lines_per_quad;
 
      TriaObjects quads;
 
      types::geometric_orientation
      get_line_orientation(const int index, const unsigned int line_no) const;
 
      void
      set_line_orientation(const int                    index,
                           const unsigned int           line_no,
                           types::geometric_orientation line_orientation);
 
      ReferenceCell<dim - 1>
      get_quad_type(const int index) const;
 
      void
      set_quad_type(const int index, const ReferenceCell<dim - 1> face_type);
 
      TriaObjects lines;
 
      std::size_t
      memory_consumption() const;
 
      template <class Archive>
      void
      serialize(Archive &ar, const unsigned int version);
 
    private:
      std::vector<bool> quad_is_quadrilateral;
 
      std::vector<bool> quads_line_orientations;
    };
 
 
 
    template <int dim>
    inline ReferenceCell<dim - 1>
    TriaFaces<dim>::get_quad_type(const int index) const
    {
      Assert(index >= 0, ExcInternalError());
      AssertIndexRange(index, quad_is_quadrilateral.size());
      if constexpr (dim == 3)
        return quad_is_quadrilateral[index] ? ReferenceCells::Quadrilateral :
                                              ReferenceCells::Triangle;
      else
        {
          DEAL_II_ASSERT_UNREACHABLE();
          return ReferenceCells::Invalid<dim - 1>;
        }
    }
 
 
 
    template <int dim>
    inline void
    TriaFaces<dim>::set_quad_type(const int                    index,
                                  const ReferenceCell<dim - 1> face_type)
    {
      Assert(index >= 0, ExcInternalError());
      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 <int dim>
    inline types::geometric_orientation
    TriaFaces<dim>::get_line_orientation(const int          index,
                                         const unsigned int line_no) const
    {
      Assert(index >= 0, ExcInternalError());
      AssertIndexRange(line_no, lines_per_quad);
      const std::size_t i = index * lines_per_quad + line_no;
      AssertIndexRange(i, quads_line_orientations.size());
      if constexpr (dim == 3)
        return quads_line_orientations[i] ?
                 numbers::default_geometric_orientation :
                 numbers::reverse_line_orientation;
      else
        {
          DEAL_II_ASSERT_UNREACHABLE();
          return numbers::invalid_geometric_orientation;
        }
    }
 
 
 
    template <int dim>
    inline void
    TriaFaces<dim>::set_line_orientation(
      const int                          index,
      const unsigned int                 line_no,
      const types::geometric_orientation line_orientation)
    {
      Assert(index >= 0, ExcInternalError());
      AssertIndexRange(line_no, lines_per_quad);
      const std::size_t i = index * lines_per_quad + line_no;
      AssertIndexRange(i, quads_line_orientations.size());
      Assert(line_orientation == numbers::default_geometric_orientation ||
               line_orientation == numbers::reverse_line_orientation,
             ExcInternalError());
      if constexpr (dim == 3)
        quads_line_orientations[i] =
          line_orientation == numbers::default_geometric_orientation;
      else
        DEAL_II_ASSERT_UNREACHABLE();
    }
 
 
 
    template <int dim>
    template <class Archive>
    void
    TriaFaces<dim>::serialize(Archive &ar, const unsigned int)
    {
      ar &lines_per_quad &quads &lines &quads_line_orientations
        &quad_is_quadrilateral;
    }
  } // namespace TriangulationImplementation
} // namespace internal
 
DEAL_II_NAMESPACE_CLOSE
 
#endif