solution_transfer.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) 2024 - 2024 by the ryujin authors
//
 
#pragma once
 
#include <compile_time_options.h>
 
#include "mpi_ensemble.h"
#include "offline_data.h"
#include "state_vector.h"
 
#include <deal.II/base/parameter_acceptor.h>
 
namespace ryujin
{
  template <typename Description, int dim, typename Number = double>
  class SolutionTransfer
  {
  public:
 
    using HyperbolicSystem = typename Description::HyperbolicSystem;
    using ParabolicSystem = typename Description::ParabolicSystem;
 
    using View =
        typename Description::template HyperbolicSystemView<dim, Number>;
 
    static constexpr auto problem_dimension = View::problem_dimension;
 
    using state_type = typename View::state_type;
 
    using StateVector = typename View::StateVector;
    using HyperbolicVector = typename View::HyperbolicVector;
 
    using ScalarVector = Vectors::ScalarVector<Number>;
 
 
 
    SolutionTransfer(const MPIEnsemble &mpi_ensemble,
                     typename Discretization<dim>::Triangulation &triangulation,
                     const OfflineData<dim, Number> &offline_data,
                     const HyperbolicSystem &hyperbolic_system,
                     const ParabolicSystem &parabolic_system);
 
    ~SolutionTransfer() = default;
 
 
 
    void prepare_projection(const StateVector &old_state_vector);
 
    unsigned int get_handle() const
    {
      Assert(handle_ != dealii::numbers::invalid_unsigned_int,
             dealii::ExcMessage("Invalid handle: Cannot retrieve a valid "
                                "handle because get_handle() can only be "
                                "called after a call to prepare_projection()"));
      return handle_;
    }
 
    void set_handle(unsigned int handle)
    {
      Assert(handle_ == dealii::numbers::invalid_unsigned_int,
             dealii::ExcMessage(
                 "Invalid state: Cannot set handle because we already have a "
                 "valid handle due to a prior call to prepare_projection()."));
      handle_ = handle;
    }
 
    void project(StateVector &new_state_vector);
 
  private:
 
 
    const MPIEnsemble &mpi_ensemble_;
 
    dealii::SmartPointer<typename Discretization<dim>::Triangulation>
        triangulation_;
    dealii::SmartPointer<const OfflineData<dim, Number>> offline_data_;
    dealii::SmartPointer<const HyperbolicSystem> hyperbolic_system_;
    dealii::SmartPointer<const ParabolicSystem> parabolic_system_;
 
    unsigned int handle_;
 
    state_type get_tensor(const HyperbolicVector &U,
                          const dealii::types::global_dof_index global_i);
 
    void add_tensor(HyperbolicVector &U,
                    const state_type &U_i,
                    const dealii::types::global_dof_index global_i);
  };
} // namespace ryujin