ryujin 2.1.1 revision 350e54cc11f3d21282bddcf3e6517944dbc508bf
offline_data.h
Go to the documentation of this file.
1//
2// SPDX-License-Identifier: Apache-2.0 WITH LLVM-exception
3// Copyright (C) 2020 - 2024 by the ryujin authors
4//
5
6#pragma once
7
8#include <compile_time_options.h>
9
10#include "convenience_macros.h"
11#include "discretization.h"
12#include "mpi_ensemble.h"
13#include "sparse_matrix_simd.h"
14#include "state_vector.h"
15
16#include <deal.II/base/parameter_acceptor.h>
17#include <deal.II/base/partitioner.h>
18#include <deal.II/base/smartpointer.h>
19#include <deal.II/dofs/dof_handler.h>
20#include <deal.II/lac/affine_constraints.h>
21#include <deal.II/lac/la_parallel_vector.h>
22#include <deal.II/lac/sparse_matrix.h>
23
24#include <deal.II/numerics/data_out.h>
25
26namespace ryujin
27{
51 template <int dim, typename Number = double>
52 class OfflineData : public dealii::ParameterAcceptor
53 {
54 public:
59
64
71 std::tuple<unsigned int /*i*/,
72 dealii::Tensor<1, dim, Number> /*normal*/,
73 Number /*normal mass*/,
74 Number /*boundary mass*/,
75 dealii::types::boundary_id /*id*/,
76 dealii::Point<dim>> /*position*/;
77
82 using CouplingDescription = std::tuple<unsigned int /*i*/, //
83 unsigned int /*col_idx*/,
84 unsigned int /*j*/>;
85
89 OfflineData(const MPIEnsemble &mpi_ensemble,
91 const std::string &subsection = "/OfflineData");
92
101 void prepare(const unsigned int problem_dimension,
102 const unsigned int n_precomputed_values,
103 const unsigned int n_parabolic_state_vectors)
104 {
105 setup(problem_dimension, n_precomputed_values);
106
107 create_matrices();
108
109 if (!dof_handler_->has_hp_capabilities())
110 create_multigrid_data();
111
112 n_parabolic_state_vectors_ = n_parabolic_state_vectors;
113 }
114
120
121
126
132
138
144
149
162
170
177
185
198
206
212
218
224
229
237
242
248
254
259
270
275
280
281 private:
286
291 void create_constraints_and_sparsity_pattern();
292
300 void setup(const unsigned int problem_dimension,
301 const unsigned int n_precomputed_values);
302
306 void create_matrices();
307
311 void create_multigrid_data();
312
314
318
319 const MPIEnsemble &mpi_ensemble_;
320
321 std::unique_ptr<dealii::DoFHandler<dim>> dof_handler_;
322
323 dealii::AffineConstraints<Number> affine_constraints_;
324
325 std::shared_ptr<const dealii::Utilities::MPI::Partitioner>
326 scalar_partitioner_;
327
328 std::shared_ptr<const dealii::Utilities::MPI::Partitioner>
329 hyperbolic_vector_partitioner_;
330
331 std::shared_ptr<const dealii::Utilities::MPI::Partitioner>
332 precomputed_vector_partitioner_;
333
334 unsigned int n_parabolic_state_vectors_;
335
336 unsigned int n_export_indices_;
337 unsigned int n_locally_internal_;
338 unsigned int n_locally_owned_;
339 unsigned int n_locally_relevant_;
340
341 using BoundaryMap = std::vector<BoundaryDescription>;
342 BoundaryMap boundary_map_;
343 std::vector<BoundaryMap> level_boundary_map_;
344
345 using CouplingBoundaryPairs = std::vector<CouplingDescription>;
346 CouplingBoundaryPairs coupling_boundary_pairs_;
347
348 dealii::DynamicSparsityPattern sparsity_pattern_;
349
350 SparsityPatternSIMD<dealii::VectorizedArray<Number>::size()>
351 sparsity_pattern_simd_;
352
353 SparseMatrixSIMD<Number> mass_matrix_;
354 SparseMatrixSIMD<Number> mass_matrix_inverse_;
355
356 ScalarVector lumped_mass_matrix_;
357 ScalarVector lumped_mass_matrix_inverse_;
358
359 std::vector<ScalarVectorFloat> level_lumped_mass_matrix_;
360
361 SparseMatrixSIMD<Number, dim> cij_matrix_;
362 SparseMatrixSIMD<Number> incidence_matrix_;
363
364 Number measure_of_omega_;
365
366 dealii::SmartPointer<const Discretization<dim>> discretization_;
367
371 template <typename ITERATOR1, typename ITERATOR2>
372 BoundaryMap construct_boundary_map(
373 const ITERATOR1 &begin,
374 const ITERATOR2 &end,
375 const dealii::Utilities::MPI::Partitioner &partitioner) const;
376
381 template <typename ITERATOR1, typename ITERATOR2>
382 CouplingBoundaryPairs collect_coupling_boundary_pairs(
383 const ITERATOR1 &begin,
384 const ITERATOR2 &end,
385 const dealii::Utilities::MPI::Partitioner &partitioner) const;
386
388
392
393 double incidence_relaxation_even_;
394 double incidence_relaxation_odd_;
395
397 };
398
399} /* namespace ryujin */
const auto & n_export_indices() const
Definition: offline_data.h:161
Vectors::ScalarVector< float > ScalarVectorFloat
Definition: offline_data.h:63
const auto & boundary_map() const
Definition: offline_data.h:197
const auto & n_parabolic_state_vectors() const
Definition: offline_data.h:148
const auto & level_lumped_mass_matrix() const
Definition: offline_data.h:253
const auto & lumped_mass_matrix_inverse() const
Definition: offline_data.h:247
const auto & level_boundary_map() const
Definition: offline_data.h:211
const auto & scalar_partitioner() const
Definition: offline_data.h:131
void prepare(const unsigned int problem_dimension, const unsigned int n_precomputed_values, const unsigned int n_parabolic_state_vectors)
Definition: offline_data.h:101
const auto & sparsity_pattern() const
Definition: offline_data.h:217
const auto & n_locally_owned() const
Definition: offline_data.h:176
const auto & lumped_mass_matrix() const
Definition: offline_data.h:241
const auto & affine_constraints() const
Definition: offline_data.h:125
const auto & hyperbolic_vector_partitioner() const
Definition: offline_data.h:137
const auto & precomputed_vector_partitioner() const
Definition: offline_data.h:143
std::tuple< unsigned int, unsigned int, unsigned int > CouplingDescription
Definition: offline_data.h:84
Vectors::ScalarVector< Number > ScalarVector
Definition: offline_data.h:58
std::tuple< unsigned int, dealii::Tensor< 1, dim, Number >, Number, Number, dealii::types::boundary_id, dealii::Point< dim > > BoundaryDescription
Definition: offline_data.h:76
const auto & mass_matrix_inverse() const
Definition: offline_data.h:236
const auto & mass_matrix() const
Definition: offline_data.h:228
const auto & discretization() const
Definition: offline_data.h:279
OfflineData(const MPIEnsemble &mpi_ensemble, const Discretization< dim > &discretization, const std::string &subsection="/OfflineData")
const auto & measure_of_omega() const
Definition: offline_data.h:274
const auto & n_locally_relevant() const
Definition: offline_data.h:184
const auto & coupling_boundary_pairs() const
Definition: offline_data.h:205
const auto & sparsity_pattern_simd() const
Definition: offline_data.h:223
const auto & dof_handler() const
Definition: offline_data.h:119
const auto & incidence_matrix() const
Definition: offline_data.h:269
const auto & n_locally_internal() const
Definition: offline_data.h:169
const auto & cij_matrix() const
Definition: offline_data.h:258
#define ACCESSOR_READ_ONLY(member)
#define ACCESSOR_READ_ONLY_NO_DEREFERENCE(member)
dealii::LinearAlgebra::distributed::Vector< Number > ScalarVector
Definition: state_vector.h:31