time_loop.h

Go to the documentation of this file.

//
// SPDX-License-Identifier: Apache-2.0 WITH LLVM-exception
// Copyright (C) 2020 - 2024 by the ryujin authors
//
 
#pragma once
 
#include <compile_time_options.h>
 
#include "discretization.h"
#include "hyperbolic_module.h"
#include "initial_values.h"
#include "mesh_adaptor.h"
#include "mpi_ensemble.h"
#include "mpi_ensemble_container.h"
#include "offline_data.h"
#include "parabolic_module.h"
#include "postprocessor.h"
#include "quantities.h"
#include "solution_transfer.h"
#include "time_integrator.h"
#include "vtu_output.h"
 
#include <deal.II/base/parameter_acceptor.h>
#include <deal.II/base/timer.h>
 
#include <fstream>
 
namespace ryujin
{
  template <typename Description, int dim, typename Number = double>
  class TimeLoop final : public dealii::ParameterAcceptor
  {
  public:
 
    using HyperbolicSystem = typename Description::HyperbolicSystem;
 
    using View =
        typename Description::template HyperbolicSystemView<dim, Number>;
 
    using ParabolicSystem = typename Description::ParabolicSystem;
 
    using ParabolicSolver =
        typename Description::template ParabolicSolver<dim, Number>;
 
    using ScalarNumber = typename View::ScalarNumber;
 
    static constexpr auto problem_dimension = View::problem_dimension;
 
    static constexpr auto n_precomputed_values = View::n_precomputed_values;
 
    using StateVector = typename View::StateVector;
 
    using ScalarVector = Vectors::ScalarVector<Number>;
 
 
 
    TimeLoop(const MPI_Comm &mpi_comm);
 
    void run();
 
  protected:
 
    template <typename Callable>
    void read_checkpoint(StateVector &state_vector,
                         const std::string &base_name,
                         Number &t,
                         unsigned int &output_cycle,
                         const Callable &prepare_compute_kernels);
 
    void write_checkpoint(const StateVector &state_vector,
                          const std::string &base_name,
                          const Number &t,
                          const unsigned int &output_cycle);
 
    template <typename Callable>
    void adapt_mesh_and_transfer_state_vector(
        StateVector &state_vector, const Callable &prepare_compute_kernels);
 
    void compute_error(StateVector &state_vector, Number t);
 
    void output(StateVector &state_vector,
                const std::string &name,
                const Number t,
                const unsigned int cycle);
 
    void print_parameters(std::ostream &stream);
    void print_mpi_partition(std::ostream &stream);
 
    void print_info(const std::string &header);
 
    void print_head(const std::string &header,
                    const std::string &secondary,
                    std::ostream &stream);
 
    void print_information(unsigned int output_cycle,
                           Number last_checkpoint,
                           std::ostream &stream,
                           bool final_time = false);
    void print_memory_statistics(std::ostream &stream);
    void print_timers(std::ostream &stream);
    void print_throughput(unsigned int cycle,
                          Number t,
                          std::ostream &stream,
                          bool final_time = false);
 
    void print_cycle_statistics(unsigned int cycle,
                                Number t,
                                unsigned int output_cycle,
                                Number last_checkpoint,
                                bool write_to_logfile = false,
                                bool final_time = false);
 
  private:
 
    std::string base_name_;
    std::string base_name_ensemble_;
 
    std::string debug_filename_;
 
    Number t_final_;
    bool enforce_t_final_;
    Number timer_granularity_;
 
    bool enable_output_full_;
    bool enable_output_levelsets_;
    bool enable_compute_error_;
    bool enable_compute_quantities_;
    bool enable_mesh_adaptivity_;
 
    unsigned int timer_output_full_multiplier_;
    unsigned int timer_output_levelsets_multiplier_;
    unsigned int timer_compute_quantities_multiplier_;
 
    std::vector<std::string> error_quantities_;
    bool error_normalize_;
 
    bool resume_;
    bool resume_at_time_zero_;
 
    Number terminal_update_interval_;
    bool terminal_show_rank_throughput_;
 
    Number checkpoint_update_interval_;
 
 
 
    MPIEnsemble mpi_ensemble_;
 
    std::map<std::string, dealii::Timer> computing_timer_;
 
    MPIEnsembleContainer<HyperbolicSystem> hyperbolic_system_;
    MPIEnsembleContainer<ParabolicSystem> parabolic_system_;
    Discretization<dim> discretization_;
    OfflineData<dim, Number> offline_data_;
    MPIEnsembleContainer<InitialValues<Description, dim, Number>>
        initial_values_;
    HyperbolicModule<Description, dim, Number> hyperbolic_module_;
    ParabolicModule<Description, dim, Number> parabolic_module_;
    TimeIntegrator<Description, dim, Number> time_integrator_;
    MeshAdaptor<Description, dim, Number> mesh_adaptor_;
    SolutionTransfer<Description, dim, Number> solution_transfer_;
    Postprocessor<Description, dim, Number> postprocessor_;
    VTUOutput<Description, dim, Number> vtu_output_;
    Quantities<Description, dim, Number> quantities_;
 
    dealii::types::global_dof_index n_global_dofs_;
 
    std::ofstream logfile_; /* log file */
 
  };
 
} // namespace ryujin