Using -pc_jacobi_rowmax or -pc_jacobi_rowsum

[Matthew Knepley]

On Wed, Feb 11, 2009 at 9:52 AM, Nguyen, Hung V ERDC-ITL-MS
<Hung.V.Nguyen at usace.army.mil> wrote:
> All,
>
> I tried to use the jacobi preconditioner with maximum entry in each row as
> the diagonal preconditioner and got the error message below. The matrix was
> built with mpiaij type. Does it require different type of matrix?

1) We need the entire error message, which will include a stack trace.

2) What version of PETSc are you running?

3) I cannot reproduce this problem in PETSc 3.0.0

  cd $PETSC_DIR/src/ksp/ksp/examples/tutorials
  make ex2
  $MPIEXEC -n 2 ./ex2 -pc_type jacobi -pc_jacobi_rowmax

  Thanks,

     Matt

> Thank you.
>
> -Hung
>
>  hvnguyen:jade20% aprun -n 2 ./test_read_petsc -ksp_type cg -pc_type jacobi
> -ksp_rtol 1.0e-12 -ksp_max_it 50000
> Time in PETSc solver: 4.725237 seconds
> The number of iteration       = 670
> The solution residual error = 3.221739e-09
> 2 norm 4.813996e-04
>  infinity norm 5.369620e-06
>  1 norm 1.255338e-01
> Application 555708 resources: utime 486, stime 12
> hvnguyen:jade20%
>
> hvnguyen:jade20%  aprun -n 2  ./test_read_petsc  -ksp_type cg -pc_type jacobi
> -ksp_rtol 1.0e-12 -ksp_max_it 50000 -pc_jacobi_rowmax
> [1]PETSC ERROR: [0]PETSC ERROR: --------------------- Error Message
> ------------------------------------
> --------------------- Error Message ------------------------------------
> [1]PETSC ERROR: [0]PETSC ERROR: No support for this operation for this object
> type!
> [1]PETSC ERROR: Mat type mpiaij!
> [1]PETSC ERROR:
> ------------------------------------------------------------------------
> [1]PETSC ERROR: No support for this operation for this object type!
> [0]PETSC ERROR: Mat type mpiaij!
> [0]PETSC ERROR: Petsc Release Version 2.3.3, Patch 15, Tue Sep 23 10:02:49
> CDT 2008 HG revision:
> 31306062cd1a6f6a2496fccb4878f485c9b91760-------------------------------------
> -----------------------------------
> [1]PETSC ERROR: See docs/changes/index.html for recent updates.
> [1]PETSC ERROR: See docs/faq.html for hints about trouble shooting.
> [1]PETSC ERROR: See docs/index.html for manual pages.
>



-- 
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