data_out_rotation.h

// ---------------------------------------------------------------------
//
// Copyright (C) 2000 - 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_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, typename DoFHandlerType = DoFHandler<dim>>
class DataOutRotation
  : public DataOut_DoFData<DoFHandlerType, DoFHandlerType::dimension + 1>
{
public:
  static const unsigned int dimension = DoFHandlerType::dimension;

  static const unsigned int space_dimension = DoFHandlerType::space_dimension;

  using cell_iterator =
    typename DataOut_DoFData<DoFHandlerType, dimension + 1>::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<dimension,
                                                          space_dimension>
      &data,
    std::vector<DataOutBase::Patch<dimension + 1, space_dimension + 1>>
      &my_patches);
};


DEAL_II_NAMESPACE_CLOSE

#endif