[petsc-dev] Deflated Krylov solvers for PETSc

[Jie Chen]

My personal experience is to use approximate eigenvectors to deflate conjugate gradient. In that regard, it is kinda like DGMRES, except that the iterative method is different. A lot of situations I have seen use eigenvectors, but just as what Vuik does, you never know when someone will propose something unconventional. Deflation being a general framework, it does not hurt (from a user's point of view) to add something to the library and let the user play.

Jie



----- Original Message -----
From: "Matthew Knepley" <knepley at gmail.com>
To: "For users of the development version of PETSc" <petsc-dev at mcs.anl.gov>
Sent: Wednesday, February 27, 2013 4:05:47 AM
Subject: Re: [petsc-dev] Deflated Krylov solvers for PETSc



I think the problem with that approach to deflation is preconditioning, unless you use injection as Vuik does. You can easily get 
this kind of projection by using MatShell. 


I tried Ronald Morgan's stuff when I was in grad school and it never made a bit of difference. Is there any situation where deflation 
is useful and cannot be phrased as our current MG? 


Matt