Calling hypre

[Ben Tay]

Sorry to trouble you ppl again.

The 1st step worked but subsequently I got the error msg:

 --------------------- Error Message ------------------------------------
[0]PETSC ERROR: Operation done in wrong order!
[0]PETSC ERROR: Cannot set the HYPRE preconditioner type once it has been
set!

I have used :

call PCHYPRESetType(pc,'boomeramg',ierr)

call PCSetType(pc,'hypre',ierr)


I changed the 2 order but it still can't work.

Thanks again






On 5/22/07, Ben Tay <zonexo at gmail.com> wrote:
>
> Ops.... how careless. Thanks again.
>
> On 5/22/07, Barry Smith <bsmith at mcs.anl.gov> wrote:
> >
> >
> > pcsettype
> >
> >
> > On Tue, 22 May 2007, Ben Tay wrote:
> >
> > > Sorry I got this error during linking:
> > >
> > > petsc_sub.o(.text+0x20c): In function
> > `petsc_sub_mp_petsc_solver_pois_':
> > > /nfs/home/enduser/g0306332/ns2d_c/petsc_sub.F:125: undefined reference
> > to
> > > `pcsetype_'
> > >
> > > Did I miss out something?
> > >
> > > Thanks
> > >
> > >
> > > On 5/22/07, Satish Balay <balay at mcs.anl.gov> wrote:
> > > >
> > > > use:
> > > >
> > > > call PCSetType(pc,'hypre',ierr)
> > > > call PCHYPRESetType(pc,'boomeramg',ierr)
> > > >
> > > > Satish
> > > >
> > > > On Tue, 22 May 2007, Ben Tay wrote:
> > > >
> > > > > Hi,
> > > > >
> > > > > presently I am using ./a.out -pc_type hypre -pc_hypre_type
> > boomeramg to
> > > > use
> > > > > hypre.
> > > > >
> > > > > However, I have 2 matrix to solve and now I only want one of them
> > to use
> > > > > hypre, while the other use jacobi + bcgs.
> > > > >
> > > > > for the solver using hypre,
> > > > >
> > > > > initially it's just:
> > > > >
> > > > > call PCSetType(pc,PCLU,ierr) - direct solver
> > > > >
> > > > > I now changed to
> > > > >
> > > > > call PCSetType(pc,hypre,ierr)
> > > > >
> > > > > call PCHYPRESetType(pc,boomeramg,ierr)
> > > > >
> > > > > However, I got the error msg:
> > > > >
> > > > > fortcom: Error: petsc_sub.F, line 125: This name does not have a
> > type,
> > > > and
> > > > > must have an explicit type.   [HYPRE]
> > > > >        call PCSetType(pc,hypre,ierr)
> > > > >
> > > > > Sorry I don't really understand the eg. hyppilut.c since my
> > knowledge of
> > > > the
> > > > > c language is very basic.
> > > > >
> > > > > So how should I set the options to use hypre?
> > > > >
> > > > > Thank you.
> > > > >
> > > >
> > > >
> > >
> >
> >
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20070522/6946a3a2/attachment.htm>