Dear all,
Summer greetings,
I am back with an other query.
Before I successfully implement the projection preconditioner which was simply
the coarse grid correction operator P = I - A*(P*A_H*R);
I implemented this simply keeping both pre and post smoothing dummy in PCMG
setup.
Now I need to revise this and re-implement this where I replace A and coarse
operator A_H with
preconditioned one i.e. M^-1 A and M^-1 A_H respectively. Thus new projection
reads as
P_new = I - (M^-1 A) {P*(M_H^-1 A_H)*R}
Any suggestion to implement this in Petsc would be appreciated.
Thanking in anticipation.
with regards, Abdul
-------------- next part --------------
An HTML attachment was scrubbed...
URL:
<http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20120608/59e11e14/attachment.html>
