[petsc-users] PCGAMG

Wed Mar 27 18:30:34 CDT 2013

Sorry, the codes attached were not the most recent versions.
Here is the code you can reproduce the error with.

Cheers,
Mohamad

On Wed, Mar 27, 2013 at 4:08 PM, Mohamad M. Nasr-Azadani
<mmnasr at gmail.com>wrote:

> Thanks Jed,
> How did you create these matrices?
>     ierr = DMCreateMatrix(DA2d,  MATMPIAIJ, &LHS_mat); CHKERRQ(ierr);
>
> (The DA2d is a DMDA in two dimensions, with 1 point stencil width, 1 dof,
> STAR_STENCIL).
> It is a 5-point based discretization of the laplace equation on uniform
> grid).
>
>  Always send the FULL error message, including the stack trace.
> [0]PCSetData_AGG bs=1 MM=11000
> [0]PETSC ERROR: --------------------- Error Message
> ------------------------------------
> [0]PETSC ERROR: No support for this operation for this object type!
> [0]PETSC ERROR: Not implemented for A=mpiaij and P=seqaij!
> [0]PETSC ERROR:
> ------------------------------------------------------------------------
> [0]PETSC ERROR: Petsc Release Version 3.3.0, Patch 6, Mon Feb 11 12:26:34
> CST 2013
> [0]PETSC ERROR: See docs/changes/index.html for recent updates.
> [0]PETSC ERROR: See docs/faq.html for hints about trouble shooting.
> [0]PETSC ERROR: See docs/index.html for manual pages.
> [0]PETSC ERROR:
> ------------------------------------------------------------------------
> [0]PETSC ERROR: ./ex1 on a arch-linu named yoda.mee.ucsb.edu by mmnasr
> Wed Mar 27 13:12:28 2013
> [0]PETSC ERROR: Libraries linked from
> /home/mmnasr/Softwares/petsc-3.3-p6/arch-linux2-c-debug/lib
> [0]PETSC ERROR: Configure run at Tue Mar 26 14:35:30 2013
> [0]PETSC ERROR: Configure options --with-cc=gcc --with-fc=gfortran
> --download-f-blas-lapack --download-mpich --download-hypre --download-mumps
> --download-blacs --download-scalapack --download-ptscotch --with-debuging=0
> [0]PETSC ERROR:
> ------------------------------------------------------------------------
> [0]PETSC ERROR: MatPtAPSymbolic_MPIAIJ() line 81 in
> src/mat/impls/aij/mpi/mpiptap.c
> [0]PETSC ERROR: MatPtAPSymbolic() line 8362 in src/mat/interface/matrix.c
> [0]PETSC ERROR: MatPtAP_Basic() line 12 in src/mat/utils/ptap.c
> [0]PETSC ERROR: MatPtAP() line 8256 in src/mat/interface/matrix.c
> [0]PETSC ERROR: createLevel() line 152 in src/ksp/pc/impls/gamg/gamg.c
> [0]PETSC ERROR: PCSetUp_GAMG() line 722 in src/ksp/pc/impls/gamg/gamg.c
> [0]PETSC ERROR: PCSetUp() line 832 in src/ksp/pc/interface/precon.c
> [0]PETSC ERROR: KSPSetUp() line 278 in src/ksp/ksp/interface/itfunc.c
> [0]PETSC ERROR: KSPSolve() line 402 in src/ksp/ksp/interface/itfunc.c
> [0]PETSC ERROR: main() line 158 in ex1.c
> application called MPI_Abort(MPI_COMM_WORLD, 56) - process 0
> [cli_0]: aborting job:
> application called MPI_Abort(MPI_COMM_WORLD, 56) - process 0
>
>
>
> Can you set up an example that can reproduce?
> See attached.
>
> Thanks,
> Mohamad
>
>
>
>
> On Wed, Mar 27, 2013 at 3:54 PM, Jed Brown <jedbrown at mcs.anl.gov> wrote:
>
>> How did you create these matrices? Always send the FULL error message,
>> including the stack trace. Can you set up an example that can reproduce?
>>
>>
>> On Wed, Mar 27, 2013 at 3:18 PM, Mohamad M. Nasr-Azadani <
>> mmnasr at gmail.com> wrote:
>>
>>> Hi,
>>>
>>> I am trying to use PCGAMG.
>>> The matrix I have is based on a 2-D DMDA, Laplace equation.
>>> This is the part of the code I have to setup the solver:
>>>
>>>  KSP solver;
>>> PC pc;
>>> ierr = KSPCreate(PETSC_COMM_WORLD, &solver);
>>> ierr = KSPGetPC(solver, &pc);CHKERRQ(ierr);
>>> ierr = PCSetType(pc, PCGAMG);CHKERRQ(ierr);
>>> ierr = KSPSetFromOptions(solver); CHKERRQ(ierr);
>>> ierr = KSPSetOperators(solver, LHS_mat, LHS_mat, SAME_NONZERO_PATTERN);
>>> CHKERRQ(ierr);
>>>
>>> When I run the code in parallel, it is fine.
>>> But when I run it in serial, this is the error I get:
>>>
>>> [0]PCSetData_AGG bs=1 MM=11000
>>> [0]PETSC ERROR: --------------------- Error Message
>>> ------------------------------------
>>> [0]PETSC ERROR: No support for this operation for this object type!
>>> [0]PETSC ERROR: Not implemented for A=mpiaij and P=seqaij!
>>> [0]PETSC ERROR:
>>> ------------------------------------------------------------------------
>>> [0]PETSC ERROR: Petsc Release Version 3.3.0, Patch 6, Mon Feb 11
>>> 12:26:34 CST 2013
>>> [0]PETSC ERROR: See docs/changes/index.html for recent updates.
>>> [0]PETSC ERROR: See docs/faq.html for hints about trouble shooting.
>>> [0]PETSC ERROR: See docs/index.html for manual pages.
>>> [0]PETSC ERROR:
>>> ------------------------------------------------------------------------
>>>
>>>
>>> Any thoughts?
>>>
>>> Mohamad
>>>
>>>
>>>
>>>
>>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20130327/ff934b8b/attachment-0001.html>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: ex1.c
Type: text/x-csrc
Size: 18759 bytes
Desc: not available
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20130327/ff934b8b/attachment-0001.c>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: ex1.h
Type: text/x-chdr
Size: 787 bytes
Desc: not available
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20130327/ff934b8b/attachment-0001.h>