Isn't this just putting M and M_H as the preconditioning operator for the smoother?<br><br><div class="gmail_quote">On Fri, Jun 8, 2012 at 11:58 AM, Abdul Hanan Sheikh <span dir="ltr"><<a href="mailto:hanangul12@yahoo.co.uk" target="_blank">hanangul12@yahoo.co.uk</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div><div style="font-size:10pt;font-family:verdana,helvetica,sans-serif"><div><span>Dear all, <br></span></div><div><span>Summer greetings, <br>
</span></div><div><span>I am back with an other query.</span></div><div><span>Before I successfully implement the projection preconditioner which was simply <br></span></div><div><span>the coarse grid correction operator P = I - A*(P*A_H*R); <br>
</span></div><div><span>I implemented this simply keeping both pre and post smoothing dummy in PCMG setup. <br></span></div><div><span>Now I need to revise this and re-implement this where I replace A and coarse operator A_H with <br>
</span></div><div><span>preconditioned one i.e. M^-1 A and M^-1 A_H respectively. Thus new projection reads as</span></div><div><span> <br></span></div><div><br><span></span></div><div><span>P_new = I - (M^-1 A) {P*(M_H^-1 A_H)*R}</span></div>
<div><br><span></span></div><div><span>Any suggestion to
 implement this in Petsc would be appreciated. <br></span></div><div><br><span></span></div><div><span>Thanking in anticipation. <br></span></div><div><span>with regards, Abdul <br></span></div><div><span style="color:rgb(0,0,255)"><span style="font-family:verdana,helvetica,sans-serif"><span style="font-family:lucida console,sans-serif"></span></span></span><span style="color:rgb(0,0,255)"><span style="font-family:verdana,helvetica,sans-serif"></span></span></div>
</div></div></blockquote></div><br>