#include <source/offline_data.h>

Inheritance diagram for ryujin::OfflineData< dim, Number >:

Collaboration diagram for ryujin::OfflineData< dim, Number >:

Public Types
using	ScalarVector = Vectors::ScalarVector< Number >

using	ScalarVectorFloat = Vectors::ScalarVector< float >

using	BoundaryDescription = std::tuple< unsigned int, dealii::Tensor< 1, dim, Number >, Number, Number, dealii::types::boundary_id, dealii::Point< dim > >

using	CouplingDescription = std::tuple< unsigned int, unsigned int, unsigned int >

Public Member Functions
	OfflineData (const MPI_Comm &mpi_communicator, const Discretization< dim > &discretization, const std::string &subsection="/OfflineData")

void	prepare (const unsigned int problem_dimension, const unsigned int n_precomputed_values)

auto &	dof_handler () const

auto &	affine_constraints () const

const auto &	scalar_partitioner () const

const auto &	hyperbolic_vector_partitioner () const

const auto &	precomputed_vector_partitioner () const

auto &	n_export_indices () const

auto &	n_locally_internal () const

auto &	n_locally_owned () const

auto &	n_locally_relevant () const

auto &	boundary_map () const

auto &	coupling_boundary_pairs () const

auto &	level_boundary_map () const

auto &	sparsity_pattern () const

auto &	sparsity_pattern_simd () const

auto &	mass_matrix () const

auto &	mass_matrix_inverse () const

auto &	lumped_mass_matrix () const

auto &	lumped_mass_matrix_inverse () const

auto &	level_lumped_mass_matrix () const

auto &	cij_matrix () const

auto &	incidence_matrix () const

auto &	measure_of_omega () const

auto &	discretization () const

template<typename ITERATOR1 , typename ITERATOR2 >
auto	construct_boundary_map (const ITERATOR1 &begin, const ITERATOR2 &end, const Utilities::MPI::Partitioner &partitioner) const -> BoundaryMap

template<typename ITERATOR1 , typename ITERATOR2 >
auto	collect_coupling_boundary_pairs (const ITERATOR1 &begin, const ITERATOR2 &end, const Utilities::MPI::Partitioner &partitioner) const -> CouplingBoundaryPairs

Detailed Description

template<int dim, typename Number = double>
class ryujin::OfflineData< dim, Number >

A class to store all data that can be precomputed offline.

This class takes a reference to a Discretization object (that itself holds a Triangulation, FiniteElement, Mapping, and Quadrature object).

Most notably this class sets up a DoFHandler, the SparsityPattern, various IndexSet objects to hold locally owned and locally relevant indices, and precomputes all matrices (mass matrix, lumped mass matrix, $c_{ij}$ matrices, and $n_{ij}$ matrices).

After prepare() is called, all getter functions return valid references.

Note: The offline data precomputed in this class is problem independent, it only depends on the chosen geometry and ansatz stored in the Discretization class.

Definition at line 47 of file offline_data.h.

Member Typedef Documentation

◆ ScalarVector

template<int dim, typename Number = double>

using ryujin::OfflineData< dim, Number >::ScalarVector = Vectors::ScalarVector<Number>

Definition at line 53 of file offline_data.h.

◆ ScalarVectorFloat

template<int dim, typename Number = double>

using ryujin::OfflineData< dim, Number >::ScalarVectorFloat = Vectors::ScalarVector<float>

Scalar vector storing single-precision floats

Definition at line 58 of file offline_data.h.

◆ BoundaryDescription

template<int dim, typename Number = double>

using ryujin::OfflineData< dim, Number >::BoundaryDescription = std::tuple<unsigned int , dealii::Tensor<1, dim, Number> , Number , Number , dealii::types::boundary_id , dealii::Point<dim> >

A tuple describing (local) dof index, boundary normal, normal mass, boundary mass, boundary id, and position of the boundary degree of freedom.

Definition at line 65 of file offline_data.h.

◆ CouplingDescription

template<int dim, typename Number = double>

using ryujin::OfflineData< dim, Number >::CouplingDescription = std::tuple<unsigned int , unsigned int , unsigned int >

A tuple describing coupling boundary degrees of freedom on directly enforced boundaries for which we have to symmetrize the d_ij matrix.

Definition at line 77 of file offline_data.h.

Constructor & Destructor Documentation

◆ OfflineData()

template<int dim, typename Number >

ryujin::OfflineData< dim, Number >::OfflineData	(	const MPI_Comm &	mpi_communicator,
		const Discretization< dim > &	discretization,
		const std::string &	subsection = `"/OfflineData< dim, Number >"`
	)

Constructor

Definition at line 38 of file offline_data.template.h.

Member Function Documentation

◆ prepare()

template<int dim, typename Number = double>

void ryujin::OfflineData< dim, Number >::prepare	(	const unsigned int	problem_dimension,
		const unsigned int	n_precomputed_values
	)

inline

Prepare offline data. A call to prepare() internally calls setup() and assemble().

The problem_dimension and n_precomputed_values parameters is used to set up appropriately sized vector partitioners for the state and precomputed MultiComponentVector.

Definition at line 96 of file offline_data.h.

◆ dof_handler()

template<int dim, typename Number = double>

auto & ryujin::OfflineData< dim, Number >::dof_handler ( ) const

inline

The DofHandler for our (scalar) CG ansatz space in (deal.II typical) global numbering.

Definition at line 108 of file offline_data.h.

Referenced by ryujin::Checkpointing::load_state_vector(), and ryujin::Checkpointing::write_checkpoint().

◆ affine_constraints()

template<int dim, typename Number = double>

auto & ryujin::OfflineData< dim, Number >::affine_constraints ( ) const

inline

An AffineConstraints object storing constraints in (Deal.II typical) global numbering.

Definition at line 114 of file offline_data.h.

◆ scalar_partitioner()

template<int dim, typename Number = double>

const auto & ryujin::OfflineData< dim, Number >::scalar_partitioner ( ) const

inline

An MPI partitioner for the (scalar) Vector storing a scalar-valued quantity.

Definition at line 120 of file offline_data.h.

Referenced by ryujin::Checkpointing::load_state_vector(), and ryujin::Checkpointing::write_checkpoint().

◆ hyperbolic_vector_partitioner()

template<int dim, typename Number = double>

const auto & ryujin::OfflineData< dim, Number >::hyperbolic_vector_partitioner ( ) const

inline

An MPI partitioner for the MultiComponentVector storing a vector-valued quantity of size HyperbolicSystem::problem_dimension.

Definition at line 126 of file offline_data.h.

Referenced by ryujin::Vectors::reinit_state_vector().

◆ precomputed_vector_partitioner()

template<int dim, typename Number = double>

const auto & ryujin::OfflineData< dim, Number >::precomputed_vector_partitioner ( ) const

inline

An MPI partitioner for the MultiComponentVector storing a vector-valued quantity of size HyperbolicSystem::problem_dimension.

Definition at line 132 of file offline_data.h.

Referenced by ryujin::Vectors::reinit_state_vector().

◆ n_export_indices()

template<int dim, typename Number = double>

auto & ryujin::OfflineData< dim, Number >::n_export_indices ( ) const

inline

The subinterval $[0,\texttt{n_export_indices()})$ contains all (SIMD-vectorized) indices of the interval $[0,\texttt{n_locally_internal()})$ that are exported to neighboring MPI ranks.

Note: The interval $[\texttt{n_locally_internal()}, \texttt{n_locally_relevant()})$ (consisting of non-SIMD-vectorized indices) contains additional degrees of freedom that might have to be exported to neighboring MPI ranks.

Definition at line 145 of file offline_data.h.

◆ n_locally_internal()

template<int dim, typename Number = double>

auto & ryujin::OfflineData< dim, Number >::n_locally_internal ( ) const

inline

Number of locally owned internal degrees of freedom: In (MPI rank) local numbering all indices in the half open interval [0, n_locally_internal_) are owned by this processor, have standard connectivity, and are not situated at a boundary.

Definition at line 153 of file offline_data.h.

◆ n_locally_owned()

template<int dim, typename Number = double>

auto & ryujin::OfflineData< dim, Number >::n_locally_owned ( ) const

inline

Number of locally owned degrees of freedom: In (MPI rank) local numbering all indices in the half open interval [0, n_locally_owned_) are owned by this processor.

Definition at line 160 of file offline_data.h.

◆ n_locally_relevant()

template<int dim, typename Number = double>

auto & ryujin::OfflineData< dim, Number >::n_locally_relevant ( ) const

inline

Number of locally relevant degrees of freedom: This number is the toal number of degrees of freedom we store locally on this MPI rank. I.e., we can access the half open interval [0, n_locally_relevant_) on this machine.

Definition at line 168 of file offline_data.h.

◆ boundary_map()

template<int dim, typename Number = double>

auto & ryujin::OfflineData< dim, Number >::boundary_map ( ) const

inline

The boundary map. Local numbering.

For every degree of freedom that has nonzero support at the boundary we record the global degree of freedom index along with a weighted boundary normal, the associated boundary id, and position.

This map is later used in OfflineData to handle boundary degrees of freedom after every time step (for example to implement reflective boundary conditions).

Definition at line 181 of file offline_data.h.

Referenced by ryujin::NavierStokes::VelocityMatrix< dim, Number, Number2 >::vmult(), and ryujin::NavierStokes::EnergyMatrix< dim, Number, Number2 >::vmult().

◆ coupling_boundary_pairs()

template<int dim, typename Number = double>

auto & ryujin::OfflineData< dim, Number >::coupling_boundary_pairs ( ) const

inline

A vector of tuples describing coupling degrees of freedom i and j where both degrees of freedom are collocated at the boundary (and hence the d_ij matrix has to be symmetrized). The function returns a reference to a vector of tuples consisting of (i, col_idx, j).

Definition at line 189 of file offline_data.h.

◆ level_boundary_map()

template<int dim, typename Number = double>

auto & ryujin::OfflineData< dim, Number >::level_boundary_map ( ) const

inline

The boundary map on all levels of the grid in case multilevel support was enabled.

Definition at line 195 of file offline_data.h.

Referenced by ryujin::NavierStokes::EnergyMatrix< dim, Number, Number2 >::compute_diagonal(), ryujin::NavierStokes::VelocityMatrix< dim, Number, Number2 >::compute_diagonal(), ryujin::NavierStokes::VelocityMatrix< dim, Number, Number2 >::vmult(), and ryujin::NavierStokes::EnergyMatrix< dim, Number, Number2 >::vmult().

◆ sparsity_pattern()

template<int dim, typename Number = double>

auto & ryujin::OfflineData< dim, Number >::sparsity_pattern ( ) const

inline

A sparsity pattern for (standard deal.II) matrices storing indices in (Deal.II typical) global numbering.

Definition at line 201 of file offline_data.h.

◆ sparsity_pattern_simd()

template<int dim, typename Number = double>

auto & ryujin::OfflineData< dim, Number >::sparsity_pattern_simd ( ) const

inline

A sparsity pattern for matrices in vectorized format. Local numbering.

Definition at line 207 of file offline_data.h.

◆ mass_matrix()

template<int dim, typename Number = double>

auto & ryujin::OfflineData< dim, Number >::mass_matrix ( ) const

inline

The mass matrix. (SIMD storage, local numbering)

Definition at line 212 of file offline_data.h.

◆ mass_matrix_inverse()

template<int dim, typename Number = double>

auto & ryujin::OfflineData< dim, Number >::mass_matrix_inverse ( ) const

inline

The inverse mass matrix. (SIMD storage, local numbering)

This matrix is only available for a discontinuous finite Element ansatz.

Definition at line 220 of file offline_data.h.

◆ lumped_mass_matrix()

template<int dim, typename Number = double>

auto & ryujin::OfflineData< dim, Number >::lumped_mass_matrix ( ) const

inline

The lumped mass matrix. (stored as vector, local numbering)

Definition at line 225 of file offline_data.h.

Referenced by ryujin::NavierStokes::VelocityMatrix< dim, Number, Number2 >::vmult(), and ryujin::NavierStokes::EnergyMatrix< dim, Number, Number2 >::vmult().

◆ lumped_mass_matrix_inverse()

template<int dim, typename Number = double>

auto & ryujin::OfflineData< dim, Number >::lumped_mass_matrix_inverse ( ) const

inline

The inverse of the lumped mass matrix. (stored as vector, local numbering)

Definition at line 231 of file offline_data.h.

◆ level_lumped_mass_matrix()

template<int dim, typename Number = double>

auto & ryujin::OfflineData< dim, Number >::level_lumped_mass_matrix ( ) const

inline

The lumped mass matrix on all levels of the grid in case multilevel support was enabled.

Definition at line 237 of file offline_data.h.

Referenced by ryujin::NavierStokes::EnergyMatrix< dim, Number, Number2 >::compute_diagonal(), ryujin::NavierStokes::VelocityMatrix< dim, Number, Number2 >::compute_diagonal(), ryujin::NavierStokes::VelocityMatrix< dim, Number, Number2 >::vmult(), and ryujin::NavierStokes::EnergyMatrix< dim, Number, Number2 >::vmult().

◆ cij_matrix()

template<int dim, typename Number = double>

auto & ryujin::OfflineData< dim, Number >::cij_matrix ( ) const

inline

The $(c_{ij})$ matrix. (SIMD storage, local numbering)

Definition at line 242 of file offline_data.h.

◆ incidence_matrix()

template<int dim, typename Number = double>

auto & ryujin::OfflineData< dim, Number >::incidence_matrix ( ) const

inline

The incidence matrix $(beta_{ij})$: 1 for coupling face degrees of freedom that share the same support point coordinate, 0 otherwise.

(SIMD storage, local numbering)

This matrix is only available for a discontinuous finite Element ansatz.

Definition at line 253 of file offline_data.h.

◆ measure_of_omega()

template<int dim, typename Number = double>

auto & ryujin::OfflineData< dim, Number >::measure_of_omega ( ) const

inline

Size of computational domain.

Definition at line 258 of file offline_data.h.

◆ discretization()

template<int dim, typename Number = double>

auto & ryujin::OfflineData< dim, Number >::discretization ( ) const

inline

Returns a reference of the underlying Discretization object.

Definition at line 263 of file offline_data.h.

Referenced by ryujin::Checkpointing::write_checkpoint().

◆ construct_boundary_map()

template<int dim, typename Number = double>

template<typename ITERATOR1 , typename ITERATOR2 >

auto ryujin::OfflineData< dim, Number >::construct_boundary_map	(	const ITERATOR1 &	begin,
		const ITERATOR2 &	end,
		const Utilities::MPI::Partitioner &	partitioner
	)		const -> BoundaryMap

Definition at line 1173 of file offline_data.template.h.

References ryujin::no_slip, ryujin::periodic, and ryujin::slip.

◆ collect_coupling_boundary_pairs()

template<int dim, typename Number = double>

template<typename ITERATOR1 , typename ITERATOR2 >

auto ryujin::OfflineData< dim, Number >::collect_coupling_boundary_pairs	(	const ITERATOR1 &	begin,
		const ITERATOR2 &	end,
		const Utilities::MPI::Partitioner &	partitioner
	)		const -> CouplingBoundaryPairs

Definition at line 1371 of file offline_data.template.h.

References ryujin::periodic.

The documentation for this class was generated from the following files:

source/offline_data.h
source/offline_data.template.h

Public Types

Public Member Functions

Detailed Description

Member Typedef Documentation

◆ ScalarVector

◆ ScalarVectorFloat

◆ BoundaryDescription

◆ CouplingDescription

Constructor & Destructor Documentation

◆ OfflineData()

Member Function Documentation

◆ prepare()

◆ dof_handler()

◆ affine_constraints()

◆ scalar_partitioner()

◆ hyperbolic_vector_partitioner()

◆ precomputed_vector_partitioner()

◆ n_export_indices()

◆ n_locally_internal()

◆ n_locally_owned()

◆ n_locally_relevant()

◆ boundary_map()

◆ coupling_boundary_pairs()

◆ level_boundary_map()

◆ sparsity_pattern()

◆ sparsity_pattern_simd()

◆ mass_matrix()

◆ mass_matrix_inverse()

◆ lumped_mass_matrix()

◆ lumped_mass_matrix_inverse()

◆ level_lumped_mass_matrix()

◆ cij_matrix()

◆ incidence_matrix()

◆ measure_of_omega()

◆ discretization()

◆ construct_boundary_map()

◆ collect_coupling_boundary_pairs()