face_info.h

Go to the documentation of this file.

// ------------------------------------------------------------------------
//
// SPDX-License-Identifier: LGPL-2.1-or-later
// Copyright (C) 2018 - 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_matrix_free_face_info_h
#define dealii_matrix_free_face_info_h
 
 
#include <deal.II/base/config.h>
 
#include <deal.II/base/exceptions.h>
#include <deal.II/base/memory_consumption.h>
#include <deal.II/base/table.h>
 
 
DEAL_II_NAMESPACE_OPEN
 
 
 
namespace internal
{
  namespace MatrixFreeFunctions
  {
    template <int vectorization_width>
    struct FaceToCellTopology
    {
      std::array<unsigned int, vectorization_width> cells_interior;
 
      std::array<unsigned int, vectorization_width> cells_exterior;
 
      types::boundary_id exterior_face_no;
 
      unsigned char interior_face_no;
 
      unsigned char subface_index;
 
      unsigned char face_orientation;
 
      unsigned char face_type;
 
      std::size_t
      memory_consumption() const
      {
        return sizeof(*this);
      }
    };
 
 
 
    template <int vectorization_width>
    struct FaceInfo
    {
      void
      clear()
      {
        faces.clear();
        cell_and_face_to_plain_faces.reinit(TableIndices<3>(0, 0, 0));
        cell_and_face_boundary_id.reinit(TableIndices<3>(0, 0, 0));
      }
 
      std::size_t
      memory_consumption() const
      {
        return sizeof(faces) +
               cell_and_face_to_plain_faces.memory_consumption() +
               cell_and_face_boundary_id.memory_consumption();
      }
 
      std::vector<FaceToCellTopology<vectorization_width>> faces;
 
      ::Table<3, unsigned int> cell_and_face_to_plain_faces;
 
      ::Table<3, types::boundary_id> cell_and_face_boundary_id;
    };
  } // end of namespace MatrixFreeFunctions
} // end of namespace internal
 
DEAL_II_NAMESPACE_CLOSE
 
#endif