* argument.
*/
SparseMIC (const SparsityPattern &sparsity);
-
+
+ typedef SparseLUDecomposition<number>::AdditionalData AdditionalData;
/**
* Reinitialize the object but
*/
void reinit (const SparsityPattern &sparsity);
+ /**
+ * Same as @p{decompose}.
+ */
+ template <typename somenumber>
+ void initialize (const SparseMatrix<somenumber> &matrix,
+ const AdditionalData parameters);
+
/**
* Perform the incomplete LU
* factorization of the given
/**
* Exception
*/
- DeclException0 (ExcInternal);
+ DeclException0 (ExcStrengthenDiagonalTooSmall);
/**
* Exception
*/
+template <typename number>
+template <typename somenumber>
+inline
+void SparseMIC<number>::initialize (const SparseMatrix<somenumber> &matrix,
+ const AdditionalData data)
+{
+ decompose(matrix, data.strengthen_diagonal);
+}
+
+
+
#endif // __deal2__