Reference documentation for deal.II version Git 4fbb5374f0 2021-01-22 15:02:13 -0500
\(\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\}}\)
sparse_ilu.h
Go to the documentation of this file.
1 // ---------------------------------------------------------------------
2 //
3 // Copyright (C) 1999 - 2018 by the deal.II authors
4 //
5 // This file is part of the deal.II library.
6 //
7 // The deal.II library is free software; you can use it, redistribute
8 // it, and/or modify it under the terms of the GNU Lesser General
9 // Public License as published by the Free Software Foundation; either
10 // version 2.1 of the License, or (at your option) any later version.
11 // The full text of the license can be found in the file LICENSE.md at
12 // the top level directory of deal.II.
13 //
14 // ---------------------------------------------------------------------
15 
16 
17 #ifndef dealii_sparse_ilu_h
18 #define dealii_sparse_ilu_h
19 
20 
21 #include <deal.II/base/config.h>
22 
23 #include <deal.II/lac/exceptions.h>
26 
28 
58 template <typename number>
59 class SparseILU : public SparseLUDecomposition<number>
60 {
61 public:
66 
73  SparseILU() = default;
74 
80 
98  template <typename somenumber>
99  void
101  const AdditionalData & parameters = AdditionalData());
102 
109  template <typename somenumber>
110  void
111  vmult(Vector<somenumber> &dst, const Vector<somenumber> &src) const;
112 
113 
120  template <typename somenumber>
121  void
122  Tvmult(Vector<somenumber> &dst, const Vector<somenumber> &src) const;
123 
124 
129  std::size_t
130  memory_consumption() const override;
131 
141  double,
142  << "The strengthening parameter " << arg1
143  << " is not greater or equal than zero!");
148  size_type,
149  << "While computing the ILU decomposition, the algorithm "
150  "found a zero pivot on the diagonal of row "
151  << arg1
152  << ". This must stop the ILU algorithm because it means "
153  "that the matrix for which you try to compute a "
154  "decomposition is singular.");
156 };
157 
159 //---------------------------------------------------------------------------
160 
161 
163 
164 #endif // dealii_sparse_ilu_h
Contents is actually a matrix.
static ::ExceptionBase & ExcInvalidStrengthening(double arg1)
SparseILU()=default
typename SparseLUDecomposition< number >::AdditionalData AdditionalData
Definition: sparse_ilu.h:79
#define DeclException1(Exception1, type1, outsequence)
Definition: exceptions.h:515
typename SparseMatrix< number >::size_type size_type
void vmult(Vector< somenumber > &dst, const Vector< somenumber > &src) const
#define DEAL_II_NAMESPACE_CLOSE
Definition: config.h:380
std::size_t memory_consumption() const override
static ::ExceptionBase & ExcZeroPivot(size_type arg1)
void initialize(const SparseMatrix< somenumber > &matrix, const AdditionalData &parameters=AdditionalData())
#define DEAL_II_NAMESPACE_OPEN
Definition: config.h:379
void Tvmult(Vector< somenumber > &dst, const Vector< somenumber > &src) const