data_out_rotation.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_rotation_h
#define dealii_data_out_rotation_h
 
 
#include <deal.II/base/config.h>
 
#include <deal.II/numerics/data_out_dof_data.h>
 
#include <string>
#include <vector>
 
DEAL_II_NAMESPACE_OPEN
 
 
namespace internal
{
  namespace DataOutRotationImplementation
  {
    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 unsigned int               n_patches_per_circle,
                   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);
 
      const unsigned int n_patches_per_circle;
 
      std::vector<Point<spacedim>> patch_evaluation_points;
    };
  } // namespace DataOutRotationImplementation
} // namespace internal
 
 

template <int dim, int spacedim = dim>
class DataOutRotation
  : public DataOut_DoFData<dim, dim + 1, spacedim, spacedim + 1>
{
  static_assert(dim == spacedim, "Not implemented for dim != spacedim.");
 
public:
  static constexpr int patch_dim      = dim + 1;
  static constexpr int patch_spacedim = spacedim + 1;

  using cell_iterator =
    typename DataOut_DoFData<dim, patch_dim, spacedim, patch_spacedim>::
      cell_iterator;

  virtual void
  build_patches(const unsigned int n_patches_per_circle,
                const unsigned int n_subdivisions = 0);

  virtual cell_iterator
  first_cell();

  virtual cell_iterator
  next_cell(const cell_iterator &cell);

  DeclException1(ExcRadialVariableHasNegativeValues,
                 double,
                 << "You are attempting to use this class on a triangulation "
                    "in which some vertices have a negative radial coordinate "
                    "value of "
                 << arg1
                 << ". If you rotate such a triangulation around an "
                    "axis, you will get (dim+1)-dimensional meshes "
                    "that are not likely what you hoped to see.");
 
private:
  void
  build_one_patch(
    const cell_iterator                                                  *cell,
    internal::DataOutRotationImplementation::ParallelData<dim, spacedim> &data,
    std::vector<DataOutBase::Patch<patch_dim, patch_spacedim>> &my_patches);
};
 
 
DEAL_II_NAMESPACE_CLOSE
 
#endif