data_out_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) 2000 - 2023 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_data_out_faces_h
#define dealii_data_out_faces_h
 
 
#include <deal.II/base/config.h>
 
#include <deal.II/numerics/data_out.h>
 
#include <string>
#include <vector>
 
DEAL_II_NAMESPACE_OPEN
 
 
namespace internal
{
  namespace DataOutFacesImplementation
  {
    template <int dim, int spacedim>
    struct ParallelData
      : public internal::DataOutImplementation::ParallelDataBase<dim, spacedim>
    {
      ParallelData(const unsigned int               n_datasets,
                   const unsigned int               n_subdivisions,
                   const std::vector<unsigned int> &n_postprocessor_outputs,
                   const Mapping<dim, spacedim>    &mapping,
                   const std::vector<
                     std::shared_ptr<::hp::FECollection<dim, spacedim>>>
                                    &finite_elements,
                   const UpdateFlags update_flags);
 
      std::vector<Point<spacedim>> patch_evaluation_points;
    };
  } // namespace DataOutFacesImplementation
} // namespace internal
 
 
template <int dim, int spacedim = dim>
class DataOutFaces : public DataOut_DoFData<dim, dim - 1, spacedim, spacedim>
{
  static_assert(dim == spacedim, "Not implemented for dim != spacedim.");
 
public:
  static constexpr int patch_dim      = dim - 1;
  static constexpr int patch_spacedim = spacedim;
 
  using cell_iterator =
    typename DataOut_DoFData<dim, patch_dim, spacedim, patch_spacedim>::
      cell_iterator;
 
  DataOutFaces(const bool surface_only = true);
 
  virtual void
  build_patches(const unsigned int n_subdivisions = 0);
 
  virtual void
  build_patches(const Mapping<dim, spacedim> &mapping,
                const unsigned int            n_subdivisions = 0);
 
  using FaceDescriptor = typename std::pair<cell_iterator, unsigned int>;
 
 
  virtual FaceDescriptor
  first_face();
 
  virtual FaceDescriptor
  next_face(const FaceDescriptor &face);
 
private:
  const bool surface_only;
 
  void
  build_one_patch(
    const FaceDescriptor *cell_and_face,
    internal::DataOutFacesImplementation::ParallelData<dim, spacedim> &data,
    DataOutBase::Patch<patch_dim, patch_spacedim>                     &patch);
};
 
 
DEAL_II_NAMESPACE_CLOSE
 
#endif