15#ifndef dealii_sparse_decomposition_h
16#define dealii_sparse_decomposition_h
105template <
typename number>
209 template <
typename somenumber>
243 template <
class OutVector,
class InVector>
254 template <
class OutVector,
class InVector>
275 <<
"The strengthening parameter " <<
arg1
276 <<
" is not greater or equal than zero!");
283 template <
typename somenumber>
344template <
typename number>
348 const size_type )
const
350 return strengthen_diagonal;
355template <
typename number>
363template <
typename number>
371template <
typename number>
382template <
typename number>
383template <
class OutVector,
class InVector>
390 this->vmult(tmp, src);
398template <
typename number>
399template <
class OutVector,
class InVector>
406 this->Tvmult(tmp, src);
413template <
typename number>
unsigned int extra_off_diagonals
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)
double strengthen_diagonal
bool use_previous_sparsity
const SparsityPattern * use_this_sparsity
virtual number get_strengthen_diagonal(const number rowsum, const size_type row) const
virtual void strengthen_diagonal_impl()
void Tvmult_add(OutVector &dst, const InVector &src) const
typename SparseMatrix< number >::size_type size_type
double strengthen_diagonal
virtual ~SparseLUDecomposition() override=0
void vmult_add(OutVector &dst, const InVector &src) const
void copy_from(const SparseMatrix< somenumber > &matrix)
virtual std::size_t memory_consumption() const
std::vector< const size_type * > prebuilt_lower_bound
virtual void clear() override
void prebuild_lower_bound()
void initialize(const SparseMatrix< somenumber > &matrix, const AdditionalData parameters)
SparsityPattern * own_sparsity
#define DEAL_II_NAMESPACE_OPEN
#define DEAL_II_NAMESPACE_CLOSE
static ::ExceptionBase & ExcInvalidStrengthening(double arg1)
#define DeclException1(Exception1, type1, outsequence)