ExaDG
Loading...
Searching...
No Matches
include
exadg
operators
inverse_mass_parameters.h
1
/* ______________________________________________________________________
2
*
3
* ExaDG - High-Order Discontinuous Galerkin for the Exa-Scale
4
*
5
* Copyright (C) 2021 by the ExaDG authors
6
*
7
* This program is free software: you can redistribute it and/or modify
8
* it under the terms of the GNU General Public License as published by
9
* the Free Software Foundation, either version 3 of the License, or
10
* (at your option) any later version.
11
*
12
* This program is distributed in the hope that it will be useful,
13
* but WITHOUT ANY WARRANTY; without even the implied warranty of
14
* MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
15
* GNU General Public License for more details.
16
*
17
* You should have received a copy of the GNU General Public License
18
* along with this program. If not, see <https://www.gnu.org/licenses/>.
19
* ______________________________________________________________________
20
*/
21
22
#ifndef INCLUDE_EXADG_OPERATORS_INVERSE_MASS_PARAMETERS_H_
23
#define INCLUDE_EXADG_OPERATORS_INVERSE_MASS_PARAMETERS_H_
24
25
// ExaDG
26
#include <exadg/solvers_and_preconditioners/solvers/solver_data.h>
27
28
namespace
ExaDG
29
{
30
enum class
InverseMassType
31
{
32
MatrixfreeOperator,
// currently only available via deal.II for Hypercube elements with n_nodes_1d
33
// = n_q_points_1d
34
ElementwiseKrylovSolver,
35
BlockMatrices
36
};
37
38
enum class
PreconditionerMass
39
{
40
None,
41
PointJacobi
42
};
43
48
struct
InverseMassParameters
49
{
50
InverseMassParameters
()
51
: implementation_type(InverseMassType::MatrixfreeOperator),
52
preconditioner(PreconditionerMass::PointJacobi),
53
solver_data(
SolverData
(1000, 1e-12, 1e-12))
54
{
55
}
56
57
// The implementation type used to invert the mass operator.
58
InverseMassType implementation_type;
59
60
// This parameter is only relevant if the mass operator is inverted by an iterative solver with
61
// matrix-free implementation, InverseMassType::ElementwiseKrylovSolver.
62
PreconditionerMass preconditioner;
63
64
// solver data for iterative solver in case of implementation type
65
// InverseMassType::ElementwiseKrylovSolver.
66
SolverData
solver_data;
67
};
68
74
struct
InverseMassParametersHdiv
75
{
76
InverseMassParametersHdiv
()
77
: preconditioner(PreconditionerMass::PointJacobi), solver_data(
SolverData
(1000, 1e-12, 1e-12))
78
{
79
}
80
81
// The preconditioner used to iteratively solve the global mass problem
82
PreconditionerMass preconditioner;
83
84
// solver data for iterative solver
85
SolverData
solver_data;
86
};
87
}
// namespace ExaDG
88
89
90
#endif
/* INCLUDE_EXADG_OPERATORS_INVERSE_MASS_PARAMETERS_H_ */
ExaDG
Definition
driver.cpp:33
ExaDG::InverseMassParametersHdiv
Definition
inverse_mass_parameters.h:75
ExaDG::InverseMassParameters
Definition
inverse_mass_parameters.h:49
ExaDG::SolverData
Definition
solver_data.h:34
Generated by
1.10.0