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>
