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