face_info.h

// ---------------------------------------------------------------------
//
// Copyright (C) 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_matrix_free_face_info_h
#define dealii_matrix_free_face_info_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
    {
      unsigned int cells_interior[vectorization_width];

      unsigned int cells_exterior[vectorization_width];

      unsigned char interior_face_no;

      unsigned char exterior_face_no;

      unsigned char subface_index;

      unsigned char face_orientation;

      std::size_t
      memory_consumption() const
      {
        return sizeof(*this);
      }
    };



    template <int vectorization_width>
    struct FaceInfo
    {
      void
      clear()
      {
        faces = std::vector<FaceToCellTopology<vectorization_width>>();
        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