PreconditionJacobi::initialize (const SparseMatrix &matrix,
const AdditionalData &additional_data)
{
+ preconditioner.release();
preconditioner = Teuchos::rcp (Ifpack().Create ("point relaxation", &*matrix.matrix, 0));
Assert (&*preconditioner != 0, ExcMessage ("Trilinos could not create this "
PreconditionSSOR::initialize (const SparseMatrix &matrix,
const AdditionalData &additional_data)
{
+ preconditioner.release();
preconditioner = Teuchos::rcp (Ifpack().Create ("point relaxation",
&*matrix.matrix,
PreconditionSOR::initialize (const SparseMatrix &matrix,
const AdditionalData &additional_data)
{
+ preconditioner.release();
preconditioner = Teuchos::rcp (Ifpack().Create ("point relaxation",
&*matrix.matrix,
PreconditionIC::initialize (const SparseMatrix &matrix,
const AdditionalData &additional_data)
{
+ preconditioner.release();
+
preconditioner = Teuchos::rcp (Ifpack().Create ("IC", &*matrix.matrix,
additional_data.overlap));
Assert (&*preconditioner != 0, ExcMessage ("Trilinos could not create this "
PreconditionILU::initialize (const SparseMatrix &matrix,
const AdditionalData &additional_data)
{
+ preconditioner.release();
+
preconditioner = Teuchos::rcp (Ifpack().Create ("ILU", &*matrix.matrix,
additional_data.overlap));
Assert (&*preconditioner != 0, ExcMessage ("Trilinos could not create this "