[petsc-users] with hypre

[Matthew Knepley]

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

> I am not sure which PETSc makefile to use as there are many in PETSc
> directory. Also do you mean PETSc makefile just to build this example or to
> build hypre from scratch? In the makefile I sent you I added the path of
> libfblas and libflapack which PETSc had downloaded and installed as an
> external package but it didn't work either.
>

There is an entire chapter in the manual on makefiles, including a section
on making user makefiles. Have
you read this? You can see our full link line using "make getlinklibs", but
that should be used by advanced
users of make.

    Matt


> On Fri, Aug 3, 2012 at 8:24 AM, Matthew Knepley <knepley at gmail.com> wrote:
>
>> 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
>>
>
>
>
> --
> 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/7614bd60/attachment.html>