PETSc + HYPRE

[Márcio Ricardo Pivello]

Hi, I want to use hypre preconditioners coupled with PETSc, but so far I
have not succeeded. Here's what I've done:

Firstly I create the preconditioner:


           Mat A_Par(NSubSteps)
           Vec Unk_Par(NSubSteps)
           Vec B_Load_Par(NSubSteps)
           KSP KspSolv
--->      PC  precond

******************************

Later in the code I set the preconditioner type and create the Krylov
solver:

---->    call PCSetType(precond,'hypre',iError)
---->    call PCHYPRESetType(precond,'boomeramg',iError)
---->    call KSPCreate (PETSC_COMM_WORLD, KspSolv, iError)
---->    call KSPSetFromOptions (KspSolv, iError)
          call KSPSetOperators (KspSolv, A_Par(nstp), A_Par(nstp),
SAME_NONZERO_PATTERN, iError)
          call KSPSolve (KspSolv, B_Load_Par(nstp), Unk_Par(nstp), iError)


***************************

Then, when I run the program I put the following options in the command
line:

mpirun -np 2 /home/mpivello/bin/SolverGP.x -pc_type hypre -pc_hypre_type
boomeramg -pc_hypre_boomeramg_sweep_all true -pc_hypre_boomeramg_grid_sweeps
1 -pc_hypre_boomeramg_strong_threshold 0.9 -pc_hypre_boomeramg_max_iter 5
-pc_hypre_boomeramg_coarsen_type modifiedRuge-Stueben -f0 dummy.tmp 2>&1
-ksp_gmres_restart 200 -ksp_max_it 3000 -ksp_rtol 1.0e-10 -ksp_atol 1.0e-15
-ksp_monitor -log_summary < /dev/null > run.parallel.log &

But this proceeding is not working. What am I doing wrong?

Thanks in advance

Márcio Ricardo
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20080415/1fb62c95/attachment.htm>