On Tue, Nov 13, 2012 at 8:16 AM, Rob Kunz <span dir="ltr"><<a href="mailto:rfk102@arl.psu.edu" target="_blank">rfk102@arl.psu.edu</a>></span> wrote:<br><div class="gmail_extra"><div class="gmail_quote"><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">






<div lang="EN-US" link="blue" vlink="purple">
<div>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1f497d">I have successfully implemented the GAMG and I am getting convergence now as expected from a multigrid solver!<u></u><u></u></span></p>

<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1f497d">Two issues:<u></u><u></u></span></p>
<p><u></u><span style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1f497d"><span>1)<span style="font:7.0pt "Times New Roman"">     
</span></span></span><u></u><span style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1f497d">Every iteration if the CFD solver the pressure equation takes more solution time. Iteration 1 the code executes at 5e-5 wall secs/cell*iteration.
 Iteration 1000 the code executes at 3e-4 wall secs/cell*iteration. Steady increase. Why does this happen?</span></p></div></div></blockquote><div>The most common cause is a memory leak. Maybe you are creating an object and not freeing it? Run a small model with -malloc_dump.</div>
<div> </div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div lang="EN-US" link="blue" vlink="purple"><div><p><span style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1f497d"><u></u><u></u></span></p>

<p><u></u><span style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1f497d"><span>2)<span style="font:7.0pt "Times New Roman"">     
</span></span></span><u></u><span style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1f497d">I would like to see some documentation on how to control the GAMG, the online docs are too limited for my modest experience. Also there is only one
 example and that has very basic implementation. I am particularly interested in trying to get the solver to act only in the strong coefficient direction, with multiple cell agglomeration – i.e., try and make the solver act like a block correction multigrid.</span></p>
</div></div></blockquote><div>I'm afraid there aren't comprehensive docs at this time. If you are familiar with multigrid (see, e.g., Trottenberg's Multigrid book), the multigrid options should be easy to understand. You'll have to read about smoothed aggregation to interpret some of the -pc_gamg_* options.<br>
<br>We are refactoring some of the components in GAMG (should not affect the high-level interface) and will write better documentation once that is complete.</div></div></div>