17 #ifndef dealii_sparse_ilu_h 18 #define dealii_sparse_ilu_h 58 template <
typename number>
98 template <
typename somenumber>
109 template <
typename somenumber>
111 vmult(Vector<somenumber> &dst,
const Vector<somenumber> &src)
const;
120 template <
typename somenumber>
122 Tvmult(Vector<somenumber> &dst,
const Vector<somenumber> &src)
const;
142 <<
"The strengthening parameter " << arg1
143 <<
" is not greater or equal than zero!");
149 <<
"While computing the ILU decomposition, the algorithm " 150 "found a zero pivot on the diagonal of row " 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.");
164 #endif // dealii_sparse_ilu_h Contents is actually a matrix.
static ::ExceptionBase & ExcInvalidStrengthening(double arg1)
typename SparseLUDecomposition< number >::AdditionalData AdditionalData
#define DeclException1(Exception1, type1, outsequence)
typename SparseMatrix< number >::size_type size_type
void vmult(Vector< somenumber > &dst, const Vector< somenumber > &src) const
#define DEAL_II_NAMESPACE_CLOSE
std::size_t memory_consumption() const override
static ::ExceptionBase & ExcZeroPivot(size_type arg1)
void initialize(const SparseMatrix< somenumber > &matrix, const AdditionalData ¶meters=AdditionalData())
#define DEAL_II_NAMESPACE_OPEN
void Tvmult(Vector< somenumber > &dst, const Vector< somenumber > &src) const