offline_data.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/convenience_macros.h>
#include <ryujin/base/gpu.h>
#include <ryujin/base/mpi_ensemble.h>
#include <ryujin/discretization/discretization.h>
#include <ryujin/interface/offline_data.h>
#include <ryujin/linear_algebra/sparse_matrix.h>
#include <ryujin/linear_algebra/sparsity_pattern.h>
#include <ryujin/linear_algebra/state_vector.h>
 
#include <deal.II/base/parameter_acceptor.h>
#include <deal.II/dofs/dof_handler.h>
 
namespace ryujin
{
  template <int dim, typename Number = double>
  class OfflineData final : public Interface::OfflineData<dim, Number>
  {
  public:
 
    using typename Interface::OfflineData<dim, Number>::BoundaryDescription;
 
    using typename Interface::OfflineData<dim, Number>::CouplingDescription;
 
    using typename Interface::OfflineData<dim, Number>::BoundaryMap;
 
    using typename Interface::OfflineData<dim, Number>::CouplingBoundaryPairs;
 
 
 
    OfflineData(const MPIEnsemble &mpi_ensemble,
                const Interface::Discretization<dim> &discretization,
                const std::string &subsection = "/OfflineData");
 
    void prepare(const unsigned int problem_dimension,
                 const unsigned int n_precomputed_values) override;
 
 
  private:
 
    bool treat_fe_nothing_as_boundary_;
    double incidence_relaxation_even_;
    double incidence_relaxation_odd_;
 
 
 
    const MPIEnsemble &mpi_ensemble_;
 
 
 
    dealii::DoFHandler<dim> &writable_dof_handler()
    {
      if (this->discretization_->have_discontinuous_ansatz()) {
        Assert(this->dof_handler_dg_, dealii::ExcInternalError());
        return *this->dof_handler_dg_;
      } else {
        Assert(this->dof_handler_cg_, dealii::ExcInternalError());
        return *this->dof_handler_cg_;
      }
    }
 
    void create_dof_handlers();
 
    void renumber_for_simd();
 
    void create_constraints_and_sparsity_pattern();
 
    void ensure_simd_stride_consistency();
 
    void create_partitioner_and_simd_sparsity(
        const unsigned int problem_dimension,
        const unsigned int n_precomputed_values);
 
    void create_matrices();
 
    void create_multigrid_data();
 
    template <typename ITERATOR1, typename ITERATOR2>
    BoundaryMap construct_boundary_map(
        const ITERATOR1 &begin,
        const ITERATOR2 &end,
        const dealii::Utilities::MPI::Partitioner &partitioner) const;
 
    template <typename ITERATOR1, typename ITERATOR2>
    CouplingBoundaryPairs collect_coupling_boundary_pairs(
        const ITERATOR1 &begin,
        const ITERATOR2 &end,
        const dealii::Utilities::MPI::Partitioner &partitioner) const;
 
  };
 
} // namespace ryujin