scratch_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/discretization/discretization.h>
 
#include <deal.II/base/quadrature_lib.h>
#include <deal.II/hp/fe_values.h>
#include <deal.II/lac/full_matrix.h>
 
namespace ryujin
{
 
  template <int dim>
  class AssemblyScratchData
  {
  public:
    AssemblyScratchData(const AssemblyScratchData<dim> &assembly_scratch_data)
        : AssemblyScratchData(assembly_scratch_data.discretization_)
    {
    }
 
 
    AssemblyScratchData(const Discretization<dim> &discretization)
        : discretization_(discretization)
        , hp_fe_values_(discretization_.mapping(),
                        discretization_.finite_element(),
                        discretization_.quadrature(),
                        dealii::update_values | dealii::update_gradients |
                            dealii::update_quadrature_points |
                            dealii::update_JxW_values)
        , hp_fe_face_values_(
              discretization_.mapping(),
              discretization_.finite_element(),
              discretization_.face_quadrature(),
              dealii::update_values | dealii::update_quadrature_points |
                  dealii::update_JxW_values | dealii::update_normal_vectors)
        , hp_fe_face_values_nodal_(discretization_.mapping(),
                                   discretization_.finite_element(),
                                   discretization_.face_nodal_quadrature(),
                                   dealii::update_values |
                                       dealii::update_quadrature_points)
        , hp_fe_neighbor_face_values_(discretization_.mapping(),
                                      discretization_.finite_element(),
                                      discretization_.face_quadrature(),
                                      dealii::update_values)
        , hp_fe_neighbor_face_values_nodal_(
              discretization_.mapping(),
              discretization_.finite_element(),
              discretization_.face_nodal_quadrature(),
              dealii::update_values)
    {
    }
 
    const Discretization<dim> &discretization_;
 
    dealii::hp::FEValues<dim> hp_fe_values_;
    dealii::hp::FEFaceValues<dim> hp_fe_face_values_;
    dealii::hp::FEFaceValues<dim> hp_fe_face_values_nodal_;
    dealii::hp::FEFaceValues<dim> hp_fe_neighbor_face_values_;
    dealii::hp::FEFaceValues<dim> hp_fe_neighbor_face_values_nodal_;
  };
 
  template <int dim, typename Number = double>
  class AssemblyCopyData
  {
  public:
    bool is_locally_owned_;
    std::vector<dealii::types::global_dof_index> local_dof_indices_;
    dealii::FullMatrix<Number> cell_mass_matrix_;
    dealii::FullMatrix<Number> cell_mass_matrix_inverse_;
    dealii::FullMatrix<Number> cell_betaij_matrix_;
    std::array<dealii::FullMatrix<Number>, dim> cell_cij_matrix_;
    Number cell_measure_;
 
    static constexpr unsigned int n_faces = 2 * dim;
    std::array<std::vector<dealii::types::global_dof_index>, n_faces>
        neighbor_local_dof_indices_;
    std::array<std::array<dealii::FullMatrix<Number>, dim>, n_faces>
        interface_cij_matrix_;
    std::array<dealii::FullMatrix<Number>, n_faces> interface_incidence_matrix_;
  };
 
} // namespace ryujin