Roy and others,

I've successfully been able to use HYPRE/BoomerAMG with libMesh.  Here is what I did.

I compiled libMesh with petsc-2.3.2-p7 and hypre-1.1.19b.  My first attempt at using
HYPRE/BoomerAMG ("-pc_type hypre -pc_hypre_type boomeramg") ran into the same hurdle
as found by Roy.   Also, using HYPRE/Euclid ("-pc_hypre_type euclid") worked fine.

The problem was that there is a slightly different initialization process for BoomerAMG within
HYPRE and the current set of calls in petsc_linear_solver.C works for HYPRE/Euclid but not for HYPRE/BoomerAMG.  In particular, it was the init(...) and solve(...) calls that needed minor reworking. 
What I've done so far is created a new init(...) that takes the matrix (a PetscMatrix) as an argument
so that I can call KSPSetOperators in init(...). 

It is presently called in solve(...) which is too late for HYPRE/BoomerAMG.  HYPRE/BoomerAMG
needs KSPSetOperators(...) called before KSPSetFromOptions(...).  This is all done within init(...) in petsc_linear_solver.C. I then commented out KSPSetOperators(...)  within solve(...) since it was
already called.

It is probably best to completely replace the init() call with an init(PetscMatrix <T>* matrix) call. 
Any thoughts on the cleanest way to fix this up.

Cheers,
Travis

----- Original Message ----
From: Travis Austin <austint73@yahoo.com>
To: Roy Stogner <roystgnr@ices.utexas.edu>; Travis Austin <austint73@gmail.com>
Cc: libMesh maillist <Libmesh-users@lists.sourceforge.net>
Sent: Friday, December 15, 2006 8:49:03 AM
Subject: Re: [Libmesh-users] Multilevel and multigrid on libMesh

Roy,

That would figure.  It looks like I still need to get hypre installed with my version of petsc.  As this is
something that I want to do for my research (get AMG working) I'll probably be working on it in the
next several weeks.  If anyone else gets it going it would be nice to hear if there are known hurdles
to get over.

Anyway I assumed that BoomerAMG came installed with hypre by default.  I'll check later.

Travis



----- Original Message ----
From: Roy Stogner <roystgnr@ices.utexas.edu>
To: Travis Austin <austint73@gmail.com>
Cc: libMesh maillist <Libmesh-users@lists.sourceforge.net>
Sent: Friday, December 15, 2006 8:13:17 AM
Subject: Re: [Libmesh-users] Multilevel and multigrid on libMesh

On Thu, 14 Dec 2006, Travis Austin wrote:

> Sounds like a great idea and is the direction that we (at the
> University of Auckland) want to take libMesh.  But I was figuring on
> using Algebraic Multigrid (AMG) for all of our multilevel computing
> needs.  BoomerAMG from LLNL is available within Petsc (may need to
> be downloaded separately) and easy to use.  Is there any particular
> reason that you are *not* using AMG other than better efficiency?  I
> just wanted to make sure you were aware of using AMG as an optimal
> solver.

I actually didn't realize that was available; unfortunately even now
it doesn't seem to be working on my PETSc installation.
"-pc_type hypre -pc_hypre_type euclid" works fine, and boomeramg is in
the "-help" list of arguments for "-pc_hypre_type", but "-pc_type
hypre -pc_hypre_type boomeramg" gives me a segfault.  Perhaps I don't
have the BoomerAMG code installed, and HYPRE or PETSc just isn't being
smart about giving me an error message to that effect.
---
Roy

-------------------------------------------------------------------------
Take Surveys. Earn Cash. Influence the Future of IT
Join SourceForge.net's Techsay panel and you'll get the chance to share your
opinions on IT & business topics through brief surveys - and earn cash
http://www.techsay.com/default.php?page=join.php&p=sourceforge&CID=DEVDEV
_______________________________________________
Libmesh-users mailing list
Libmesh-users@lists.sourceforge.net
https://lists.sourceforge.net/lists/listinfo/libmesh-users






____________________________________________________________________________________
Do you Yahoo!?
Everyone is raving about the all-new Yahoo! Mail beta.
http://new.mail.yahoo.com

-------------------------------------------------------------------------
Take Surveys. Earn Cash. Influence the Future of IT
Join SourceForge.net's Techsay panel and you'll get the chance to share your
opinions on IT & business topics through brief surveys - and earn cash
http://www.techsay.com/default.php?page=join.php&p=sourceforge&CID=DEVDEV
_______________________________________________
Libmesh-users mailing list
Libmesh-users@lists.sourceforge.net
https://lists.sourceforge.net/lists/listinfo/libmesh-users


Don't get soaked. Take a quick peak at the forecast
with theYahoo! Search weather shortcut.