On Fri, Aug 3, 2012 at 7:28 AM, Amneet Bhalla <span dir="ltr"><<a href="mailto:mail2amneet@gmail.com" target="_blank">mail2amneet@gmail.com</a>></span> wrote:<br><div class="gmail_quote"><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
Hi All, <div><br></div><div>I am trying to solve system of PDEs on AMR grid. I using PETSc 's Krylov solver for that. This requires me</div><div>to use a FAC preconditioner with a non-Krylov solver for the level 0. I actually tried PETSc's Krylov Solver</div>



<div>for level 0 (by building  MPIAIJ Mat and Vec), but it does not really improve the convergence rate. So I believe</div><div>I will have to use a multigrid solver for level 0. Can I use my already built PETSc Mat and Vec and call multigrid solvers from PETSc?  </div>
</blockquote><div><br></div><div>Use AMG to start here since it is easy. You can either use ML or GAMG (which is builtin) to start, or Hypre if</div><div>you straighten out the link (likely broken from not using the makefiles).</div>
<div> </div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div>If thats not possible, I am planning to use hypre's "SysPFMG" to write a multigrid solver for level 0. </div>
<div>I have installed PETSc with hypre by using --download-hypre=1 option.  I can build and run many </div>


<div>examples in PETSC_ARCH/externalpackages/hypre-2.7.0b/src/examples/ but not ex9; which is the</div><div>closest I could find for my purpose.</div><div><a href="http://acts.nersc.gov/hypre/example09/ex9.c.html" target="_blank">http://acts.nersc.gov/hypre/example09/ex9.c.html</a></div>


<div><br></div><div> I am getting some linkage error like</div></blockquote><div><br></div><div>Always send the entire output, including the command that was run. I have no idea what file you are trying to</div><div>compile here. Does the example work?</div>
<div><br></div><div>   Matt</div><div> </div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div><div>$HOME/SOFTWARES/LIBRARIES/PETSc/linux/lib/libHYPRE.a(par_gsmg.o): In function `hypre_BoomerAMGFitVectors':</div>
<div>$HOME/SOFTWARES/LIBRARIES/PETSc/externalpackages/hypre-2.7.0b/src/parcsr_ls/par_gsmg.c:735: undefined reference to `dgels_'</div>

<div>$HOME/SOFTWARES/LIBRARIES/PETSc/linux/lib/libHYPRE.a(par_relax.o): In function `hypre_BoomerAMGRelax':</div><div>$HOME/SOFTWARES/LIBRARIES/PETSc/externalpackages/hypre-2.7.0b/src/parcsr_ls/par_relax.c:3135: undefined reference to `dgetrf_'</div>


<div>$HOME/SOFTWARES/LIBRARIES/PETSc/externalpackages/hypre-2.7.0b/src/parcsr_ls/par_relax.c:3146: undefined reference to `dgetrs_'</div><div>$HOME/SOFTWARES/LIBRARIES/PETSc/linux/lib/libHYPRE.a(schwarz.o): In function `hypre_AMGNodalSchwarzSmoother':</div>


<div>$HOME/SOFTWARES/LIBRARIES/PETSc/externalpackages/hypre-2.7.0b/src/parcsr_ls/schwarz.c:312: undefined reference to `dpotrf_'</div></div><div><br></div><div><div>Any help will be appreciated.</div><div><br></div><div>


Thanks,</div><span class="HOEnZb"><font color="#888888">-- <br><div>Amneet <br><br></div><div><br></div><div><br></div><br>
</font></span></div>
</blockquote></div><br><br clear="all"><div><br></div>-- <br>What most experimenters take for granted before they begin their experiments is infinitely more interesting than any results to which their experiments lead.<br>
-- Norbert Wiener<br>