16 #ifndef dealii__sparse_mic_h 17 #define dealii__sparse_mic_h 19 #include <deal.II/lac/sparse_matrix.h> 20 #include <deal.II/lac/sparse_decomposition.h> 22 DEAL_II_NAMESPACE_OPEN
46 template <
typename number>
97 template <
typename somenumber>
107 template <
typename somenumber>
120 template <
typename somenumber>
144 <<
"The strengthening parameter " << arg1
145 <<
" is not greater or equal than zero!");
150 <<
"The diagonal element (" <<arg1<<
","<<arg1<<
") is " 151 << arg2 <<
", but must be positive");
174 number
get_rowsum (
const size_type row)
const;
179 DEAL_II_NAMESPACE_CLOSE
#define DeclException2(Exception2, type1, type2, outsequence)
std::vector< number > inner_sums
static::ExceptionBase & ExcDecompositionNotStable(int arg1, double arg2)
void initialize(const SparseMatrix< somenumber > &matrix, const AdditionalData ¶meters=AdditionalData())
std::vector< number > diag
#define DeclException1(Exception1, type1, outsequence)
unsigned int global_dof_index
std::size_t memory_consumption() const
#define DeclException0(Exception0)
std::vector< number > inv_diag
types::global_dof_index size_type
void Tvmult(Vector< somenumber > &dst, const Vector< somenumber > &src) const
number get_rowsum(const size_type row) const
void vmult(Vector< somenumber > &dst, const Vector< somenumber > &src) const
static::ExceptionBase & ExcInvalidStrengthening(double arg1)
static::ExceptionBase & ExcStrengthenDiagonalTooSmall()
SparseLUDecomposition< number >::AdditionalData AdditionalData