- The deprecated constructor of SparseILU.
- SparseILU::apply_decomposition.
- The deprecated constructor of SparseMIC.
- - SparseMIC::decompose.
+ - SparseMIC::decompose and SparseILU::decompose.
- The compress() functions without argument in the various vector
classes. You should use the versions with a VectorOperation
argument instead.
void initialize (const SparseMatrix<somenumber> &matrix,
const AdditionalData ¶meters = AdditionalData());
- /**
- * @deprecated This method is deprecated, and left for backward compability.
- * It will be removed in later versions.
- */
- template <typename somenumber>
- void decompose (const SparseMatrix<somenumber> &matrix,
- const double strengthen_diagonal=0.) DEAL_II_DEPRECATED;
-
/**
* Apply the incomplete decomposition, i.e. do one forward-backward step
* $dst=(LU)^{-1}src$.
{
SparseLUDecomposition<number>::initialize(matrix, data);
- decompose(matrix, data.strengthen_diagonal);
-}
-
-
-
-template <typename number>
-template <typename somenumber>
-void SparseILU<number>::decompose (const SparseMatrix<somenumber> &matrix,
- const double strengthen_diagonal)
-{
Assert (matrix.m()==matrix.n(), ExcNotQuadratic ());
Assert (this->m()==this->n(), ExcNotQuadratic ());
Assert (matrix.m()==this->m(), ExcDimensionMismatch(matrix.m(), this->m()));
- Assert (strengthen_diagonal>=0,
- ExcInvalidStrengthening (strengthen_diagonal));
+ Assert (data.strengthen_diagonal>=0,
+ ExcInvalidStrengthening (data.strengthen_diagonal));
- SparseLUDecomposition<number>::decompose (matrix, strengthen_diagonal);
+ SparseLUDecomposition<number>::decompose (matrix, data.strengthen_diagonal);
- if (strengthen_diagonal>0)
+ if (data.strengthen_diagonal>0)
this->strengthen_diagonal_impl();
// in the following, we implement algorithm 10.4 in the book by Saad by
-
template <typename number>
template <typename somenumber>
void SparseILU<number>::vmult (Vector<somenumber> &dst,
// ---------------------------------------------------------------------
//
-// Copyright (C) 1999 - 2014 by the deal.II authors
+// Copyright (C) 1999 - 2015 by the deal.II authors
//
// This file is part of the deal.II library.
//
template class SparseILU<double>;
template void SparseILU<double>::initialize<double> (const SparseMatrix<double> &,
const AdditionalData &data);
-template void SparseILU<double>::decompose<double> (const SparseMatrix<double> &,
- const double);
template void SparseILU<double>::vmult <double> (Vector<double> &,
const Vector<double> &) const;
template void SparseILU<double>::Tvmult <double> (Vector<double> &,
const Vector<double> &) const;
template void SparseILU<double>::initialize<float> (const SparseMatrix<float> &,
const AdditionalData &data);
-template void SparseILU<double>::decompose<float> (const SparseMatrix<float> &,
- const double);
template void SparseILU<double>::vmult<float> (Vector<float> &,
const Vector<float> &) const;
template void SparseILU<double>::Tvmult<float> (Vector<float> &,
template class SparseILU<float>;
template void SparseILU<float>::initialize<double> (const SparseMatrix<double> &,
const AdditionalData &data);
-template void SparseILU<float>::decompose<double> (const SparseMatrix<double> &,
- const double);
template void SparseILU<float>::vmult<double> (Vector<double> &,
const Vector<double> &) const;
template void SparseILU<float>::Tvmult<double> (Vector<double> &,
const Vector<double> &) const;
template void SparseILU<float>::initialize<float> (const SparseMatrix<float> &,
const AdditionalData &data);
-template void SparseILU<float>::decompose<float> (const SparseMatrix<float> &,
- const double);
template void SparseILU<float>::vmult<float> (Vector<float> &,
const Vector<float> &) const;
template void SparseILU<float>::Tvmult<float> (Vector<float> &,