[petsc-users] with hypre

[Matthew Knepley]

On Fri, Aug 3, 2012 at 8:12 AM, Amneet Bhalla <mail2amneet at gmail.com> wrote:

> Thanks Matt for a quick reply.
>
> Attached is the Makefile I used for building ex9. I used the same Makefile
> provided by  hypre. I just changed the
> path of hypre's directory. Also attached is the command and full output I
> got. Using the same Makefile I was able to
> build and run ex1, ex2, ex3 and ex4 of hypre.
>

You are missing the BLAS. Use the PETSc makefiles.


> Can you point me to PETSc's examples of AMG, ML or GAMG?
>

cd src/snes/examples/tutorials
make ex5
./ex5 -pc_type ml -snes_monitor

   Matt


> On Fri, Aug 3, 2012 at 7:47 AM, Matthew Knepley <knepley at gmail.com> wrote:
>
>> On Fri, Aug 3, 2012 at 7:28 AM, Amneet Bhalla <mail2amneet at gmail.com>wrote:
>>
>>> Hi All,
>>>
>>> I am trying to solve system of PDEs on AMR grid. I using PETSc 's Krylov
>>> solver for that. This requires me
>>> to use a FAC preconditioner with a non-Krylov solver for the level 0. I
>>> actually tried PETSc's Krylov Solver
>>> for level 0 (by building  MPIAIJ Mat and Vec), but it does not really
>>> improve the convergence rate. So I believe
>>> 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?
>>>
>>
>> Use AMG to start here since it is easy. You can either use ML or GAMG
>> (which is builtin) to start, or Hypre if
>> you straighten out the link (likely broken from not using the makefiles).
>>
>>
>>> If thats not possible, I am planning to use hypre's "SysPFMG" to write a
>>> multigrid solver for level 0.
>>> I have installed PETSc with hypre by using --download-hypre=1 option.  I
>>> can build and run many
>>> examples in PETSC_ARCH/externalpackages/hypre-2.7.0b/src/examples/ but
>>> not ex9; which is the
>>> closest I could find for my purpose.
>>> http://acts.nersc.gov/hypre/example09/ex9.c.html
>>>
>>>  I am getting some linkage error like
>>>
>>
>> Always send the entire output, including the command that was run. I have
>> no idea what file you are trying to
>> compile here. Does the example work?
>>
>>    Matt
>>
>>
>>> $HOME/SOFTWARES/LIBRARIES/PETSc/linux/lib/libHYPRE.a(par_gsmg.o): In
>>> function `hypre_BoomerAMGFitVectors':
>>> $HOME/SOFTWARES/LIBRARIES/PETSc/externalpackages/hypre-2.7.0b/src/parcsr_ls/par_gsmg.c:735:
>>> undefined reference to `dgels_'
>>> $HOME/SOFTWARES/LIBRARIES/PETSc/linux/lib/libHYPRE.a(par_relax.o): In
>>> function `hypre_BoomerAMGRelax':
>>> $HOME/SOFTWARES/LIBRARIES/PETSc/externalpackages/hypre-2.7.0b/src/parcsr_ls/par_relax.c:3135:
>>> undefined reference to `dgetrf_'
>>> $HOME/SOFTWARES/LIBRARIES/PETSc/externalpackages/hypre-2.7.0b/src/parcsr_ls/par_relax.c:3146:
>>> undefined reference to `dgetrs_'
>>> $HOME/SOFTWARES/LIBRARIES/PETSc/linux/lib/libHYPRE.a(schwarz.o): In
>>> function `hypre_AMGNodalSchwarzSmoother':
>>> $HOME/SOFTWARES/LIBRARIES/PETSc/externalpackages/hypre-2.7.0b/src/parcsr_ls/schwarz.c:312:
>>> undefined reference to `dpotrf_'
>>>
>>> Any help will be appreciated.
>>>
>>> Thanks,
>>> --
>>> Amneet
>>>
>>>
>>>
>>>
>>>
>>
>>
>> --
>> What most experimenters take for granted before they begin their
>> experiments is infinitely more interesting than any results to which their
>> experiments lead.
>> -- Norbert Wiener
>>
>
>
>
> --
> Amneet
>
>
>
>
>


-- 
What most experimenters take for granted before they begin their
experiments is infinitely more interesting than any results to which their
experiments lead.
-- Norbert Wiener
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20120803/2e8f1b6f/attachment-0001.html>