vtu_output.h

Go to the documentation of this file.

//
// SPDX-License-Identifier: Apache-2.0 WITH LLVM-exception
// Copyright (C) 2020 - 2026 by the ryujin authors
//
 
#pragma once
 
#include <ryujin/base/compile_time_options.h>
#include <ryujin/base/mpi_ensemble.h>
#include <ryujin/base/observer_pointer.h>
#include <ryujin/discretization/offline_data.h>
#include <ryujin/interface/hyperbolic_system.h>
#include <ryujin/interface/postprocessor.h>
#include <ryujin/interface/simulation.h>
#include <ryujin/interface/task.h>
 
#include <deal.II/base/parameter_acceptor.h>
#include <deal.II/grid/intergrid_map.h>
#include <deal.II/multigrid/mg_transfer_matrix_free.h>
 
namespace ryujin
{
 
  template <int dim, typename Number = double>
  class VTUOutput final : public dealii::ParameterAcceptor,
                          public Interface::Task<Number>
  {
  public:
 
    VTUOutput(const Interface::Simulation<dim, Number> &simulation,
              bool levelsets = false,
              const std::string &subsection = "/VTUOutput");
 
 
 
    void prepare(TaskContext<Number> &context) override;
 
    void execute(TaskContext<Number> &context) override;
 
 
  private:
 
    bool use_mpi_io_;
    bool output_analytic_solution_;
 
    Number beta_;
    bool recompute_bounds_;
 
    std::vector<std::string> manifolds_;
 
    std::vector<std::string> vtu_output_quantities_;
 
 
 
    const bool levelsets_;
    Trigger<Number> trigger_;
 
    const Interface::Simulation<dim, Number> &simulation_;
    const MPIEnsemble &mpi_ensemble_;
 
    dealii::ObserverPointer<const OfflineData<dim, Number>> offline_data_;
 
    const Interface::Postprocessor<dim, Number> &postprocessor_;
 
    typename Interface::Postprocessor<dim, Number>::Selection selection_;
 
    std::vector<std::pair<Number, Number>> bounds_;
 
 
 
    void schedule_output(const StateVector<Number> &state_vector,
                         std::string name,
                         Number t,
                         unsigned int cycle);
 
  };
 
} // namespace ryujin