deal.II version GIT relicensing-2017-g7ae82fad8b 2024-10-22 19:20:00+00:00
\(\newcommand{\dealvcentcolon}{\mathrel{\mathop{:}}}\) \(\newcommand{\dealcoloneq}{\dealvcentcolon\mathrel{\mkern-1.2mu}=}\) \(\newcommand{\jump}[1]{\left[\!\left[ #1 \right]\!\right]}\) \(\newcommand{\average}[1]{\left\{\!\left\{ #1 \right\}\!\right\}}\)
Loading...
Searching...
No Matches
precondition_selector.h
Go to the documentation of this file.
1// ------------------------------------------------------------------------
2//
3// SPDX-License-Identifier: LGPL-2.1-or-later
4// Copyright (C) 1999 - 2024 by the deal.II authors
5//
6// This file is part of the deal.II library.
7//
8// Part of the source code is dual licensed under Apache-2.0 WITH
9// LLVM-exception OR LGPL-2.1-or-later. Detailed license information
10// governing the source code and code contributions can be found in
11// LICENSE.md and CONTRIBUTING.md at the top level directory of deal.II.
12//
13// ------------------------------------------------------------------------
14
15#ifndef dealii_precondition_selector_h
16#define dealii_precondition_selector_h
17
18
19#include <deal.II/base/config.h>
20
23
24#include <string>
25
27
28// Forward declarations
29#ifndef DOXYGEN
30template <class number>
31class Vector;
32template <class number>
33class SparseMatrix;
34#endif
35
36
98template <typename MatrixType = SparseMatrix<double>,
99 typename VectorType = ::Vector<double>>
101{
102public:
106 using size_type = typename MatrixType::size_type;
107
112 PreconditionSelector(const std::string &preconditioning,
113 const typename VectorType::value_type &omega = 1.);
114
118 virtual ~PreconditionSelector() override;
119
124 void
125 use_matrix(const MatrixType &M);
126
132 m() const;
133
139 n() const;
140
145 virtual void
146 vmult(VectorType &dst, const VectorType &src) const;
147
152 virtual void
153 Tvmult(VectorType &dst, const VectorType &src) const;
154
165 static std::string
167
178
180protected:
184 std::string preconditioning;
185
186private:
191 ObserverPointer<const MatrixType,
194
198 const typename VectorType::value_type omega;
199};
200
202/* --------------------- Inline and template functions ------------------- */
203
204
205template <typename MatrixType, typename VectorType>
207 const std::string &preconditioning,
208 const typename VectorType::value_type &omega)
209 : preconditioning(preconditioning)
210 , omega(omega)
211{}
212
213
214template <typename MatrixType, typename VectorType>
216{
217 // release the matrix A
218 A = nullptr;
219}
220
221
222template <typename MatrixType, typename VectorType>
223void
225{
226 A = &M;
227}
228
229
230template <typename MatrixType, typename VectorType>
233{
234 Assert(A != nullptr, ExcNoMatrixGivenToUse());
235 return A->m();
236}
237
238
239template <typename MatrixType, typename VectorType>
242{
243 Assert(A != nullptr, ExcNoMatrixGivenToUse());
244 return A->n();
245}
246
247
248
249template <typename MatrixType, typename VectorType>
250void
252 const VectorType &src) const
253{
254 if (preconditioning == "none")
255 {
256 dst = src;
257 }
258 else
259 {
260 Assert(A != nullptr, ExcNoMatrixGivenToUse());
261
262 if (preconditioning == "jacobi")
263 {
264 A->precondition_Jacobi(dst, src, omega);
265 }
266 else if (preconditioning == "sor")
267 {
268 A->precondition_SOR(dst, src, omega);
269 }
270 else if (preconditioning == "ssor")
271 {
272 A->precondition_SSOR(dst, src, omega);
273 }
274 else
276 }
277}
278
279
280template <typename MatrixType, typename VectorType>
281void
283 VectorType &dst,
284 const VectorType &src) const
285{
286 if (preconditioning == "none")
287 {
288 dst = src;
289 }
290 else
291 {
292 Assert(A != nullptr, ExcNoMatrixGivenToUse());
293
294 if (preconditioning == "jacobi")
295 {
296 A->precondition_Jacobi(dst, src, omega); // Symmetric operation
297 }
298 else if (preconditioning == "sor")
299 {
300 A->precondition_TSOR(dst, src, omega);
301 }
302 else if (preconditioning == "ssor")
303 {
304 A->precondition_SSOR(dst, src, omega); // Symmetric operation
305 }
306 else
308 }
309}
310
311
312template <typename MatrixType, typename VectorType>
313std::string
315{
316 return "none|jacobi|sor|ssor";
317}
318
319
321
322#endif
virtual void vmult(VectorType &dst, const VectorType &src) const
virtual void Tvmult(VectorType &dst, const VectorType &src) const
typename MatrixType::size_type size_type
const VectorType::value_type omega
static std::string get_precondition_names()
void use_matrix(const MatrixType &M)
virtual ~PreconditionSelector() override
PreconditionSelector(const std::string &preconditioning, const typename VectorType::value_type &omega=1.)
ObserverPointer< const MatrixType, PreconditionSelector< MatrixType, VectorType > > A
#define DEAL_II_NAMESPACE_OPEN
Definition config.h:498
#define DEAL_II_NAMESPACE_CLOSE
Definition config.h:499
static ::ExceptionBase & ExcNoMatrixGivenToUse()
#define DeclException0(Exception0)
Definition exceptions.h:466
#define Assert(cond, exc)
#define DEAL_II_NOT_IMPLEMENTED()