16 #ifndef dealii_sparse_decomposition_h 17 #define dealii_sparse_decomposition_h 19 #include <deal.II/base/config.h> 21 #include <deal.II/lac/sparse_matrix.h> 25 DEAL_II_NAMESPACE_OPEN
109 template <
typename number>
213 template <
typename somenumber>
248 template <
class OutVector,
class InVector>
250 vmult_add(OutVector &dst,
const InVector &src)
const;
259 template <
class OutVector,
class InVector>
261 Tvmult_add(OutVector &dst,
const InVector &src)
const;
280 <<
"The strengthening parameter " << arg1
281 <<
" is not greater or equal than zero!");
288 template <
typename somenumber>
349 template <
typename number>
360 template <
typename number>
368 template <
typename number>
376 template <
typename number>
387 template <
typename number>
388 template <
class OutVector,
class InVector>
391 const InVector &src)
const 395 this->
vmult(tmp, src);
403 template <
typename number>
404 template <
class OutVector,
class InVector>
407 const InVector &src)
const 418 template <
typename number>
420 const double strengthen_diag,
421 const unsigned int extra_off_diag,
422 const bool use_prev_sparsity,
433 DEAL_II_NAMESPACE_CLOSE
435 #endif // dealii_sparse_decomposition_h void Tvmult(OutVector &dst, const InVector &src) const
void vmult(OutVector &dst, const InVector &src) const
AdditionalData(const double strengthen_diagonal=0, const unsigned int extra_off_diagonals=0, const bool use_previous_sparsity=false, const SparsityPattern *use_this_sparsity=nullptr)
void Tvmult_add(OutVector &dst, const InVector &src) const
static::ExceptionBase & ExcInvalidStrengthening(double arg1)
void copy_from(const SparseMatrix< somenumber > &matrix)
SparsityPattern * own_sparsity
bool use_previous_sparsity
double strengthen_diagonal
const SparsityPattern * use_this_sparsity
#define DeclException1(Exception1, type1, outsequence)
virtual std::size_t memory_consumption() const
virtual void clear() override
typename SparseMatrix< number >::size_type size_type
void vmult_add(OutVector &dst, const InVector &src) const
virtual ~SparseLUDecomposition() override=0
std::vector< const size_type * > prebuilt_lower_bound
virtual number get_strengthen_diagonal(const number rowsum, const size_type row) const
unsigned int extra_off_diagonals
types::global_dof_index size_type
virtual void strengthen_diagonal_impl()
double strengthen_diagonal
void initialize(const SparseMatrix< somenumber > &matrix, const AdditionalData parameters)
void prebuild_lower_bound()