]> https://gitweb.dealii.org/ - dealii.git/commitdiff
Fix a bug in PreconditionIdentity.
authorBruno Turcksin <bruno.turcksin@gmail.com>
Wed, 14 Aug 2013 15:52:05 +0000 (15:52 +0000)
committerBruno Turcksin <bruno.turcksin@gmail.com>
Wed, 14 Aug 2013 15:52:05 +0000 (15:52 +0000)
git-svn-id: https://svn.dealii.org/trunk@30309 0785d39b-7218-0410-832d-ea1e28bc413d

deal.II/source/lac/trilinos_solver.cc

index 896d5d8499073eeda7c3480a84b5ab0d9b5c7cb5..a4580ef245b9f8a56fcfabea948a8a2f8f9963d9 100644 (file)
@@ -87,6 +87,7 @@ namespace TrilinosWrappers
                               &x.trilinos_vector(),
                               const_cast<Epetra_MultiVector *>(&b.trilinos_vector())));
 
+    std::cout<<"ewww"<<std::endl;
     execute_solve(preconditioner);
   }
 
@@ -209,13 +210,16 @@ namespace TrilinosWrappers
       }
 
     // Introduce the preconditioner, 
-    // if not the identity preconditioner ...
+    // if the identity preconditioner is used,
+    // the precondioner is set to none, ...
     if (preconditioner.preconditioner.use_count()!=0)
     {
       ierr = solver.SetPrecOperator (const_cast<Epetra_Operator *>
                                      (preconditioner.preconditioner.get()));
       AssertThrow (ierr == 0, ExcTrilinosError(ierr));
     }
+    else
+      solver.SetAztecOption(AZ_precond,AZ_none);
 
     // ... set some options, ...
     solver.SetAztecOption (AZ_output, additional_data.output_solver_details ?

In the beginning the Universe was created. This has made a lot of people very angry and has been widely regarded as a bad move.

Douglas Adams


Typeset in Trocchi and Trocchi Bold Sans Serif.