Mark, I have a non-symmetric matrix. I am running with the following options.
-pc_type gamg -pc_gamg_sym_graph -ksp_monitor_true_residual and with the inclusion of -pc_gamg_sym_graph, I get a new malloc error: 0]PETSC ERROR: --------------------- Error Message ------------------------------------ [0]PETSC ERROR: Argument out of range! [0]PETSC ERROR: New nonzero at (5150,9319) caused a malloc! [0]PETSC ERROR: ------------------------------------------------------------------------ [0]PETSC ERROR: Petsc Development HG revision: 587b25035091aaa309c87c90ac64c13408ecf34e HG Date: Wed Mar 14 09:22:54 2012 -0500 [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: ../JohnRepo/VFOLD_exe on a linux-deb named wv.iihr.uiowa.eduby jmousel Wed Mar 14 11:51:35 2012 [0]PETSC ERROR: Libraries linked from /home/jmousel/NumericalLibraries/petsc-hg/petsc-dev/linux-debug/lib [0]PETSC ERROR: Configure run at Wed Mar 14 09:46:39 2012 [0]PETSC ERROR: Configure options --download-blacs=1 --download-hypre=1 --download-metis=1 --download-ml=1 --download-mpich=1 --download-parmetis=1 --download-scalapack=1 --with-blas-lapack-dir=/opt/intel11/mkl/lib/em64t --with-cc=gcc --with-cmake=/usr/local/bin/cmake --with-cxx=g++ --with-fc=ifort PETSC_ARCH=linux-debug [0]PETSC ERROR: ------------------------------------------------------------------------ [0]PETSC ERROR: MatSetValues_MPIAIJ() line 506 in /home/jmousel/NumericalLibraries/petsc-hg/petsc-dev/src/mat/impls/aij/mpi/mpiaij.c [0]PETSC ERROR: MatSetValues() line 1141 in /home/jmousel/NumericalLibraries/petsc-hg/petsc-dev/src/mat/interface/matrix.c [0]PETSC ERROR: scaleFilterGraph() line 155 in /home/jmousel/NumericalLibraries/petsc-hg/petsc-dev/src/ksp/pc/impls/gamg/tools.c [0]PETSC ERROR: PCGAMGgraph_AGG() line 865 in /home/jmousel/NumericalLibraries/petsc-hg/petsc-dev/src/ksp/pc/impls/gamg/agg.c [0]PETSC ERROR: PCSetUp_GAMG() line 516 in /home/jmousel/NumericalLibraries/petsc-hg/petsc-dev/src/ksp/pc/impls/gamg/gamg.c [0]PETSC ERROR: PCSetUp() line 832 in /home/jmousel/NumericalLibraries/petsc-hg/petsc-dev/src/ksp/pc/interface/precon.c [0]PETSC ERROR: KSPSetUp() line 261 in /home/jmousel/NumericalLibraries/petsc-hg/petsc-dev/src/ksp/ksp/interface/itfunc.c [0]PETSC ERROR: KSPSolve() line 385 in /home/jmousel/NumericalLibraries/petsc-hg/petsc-dev/src/ksp/ksp/interface/itfunc.c John On Wed, Mar 14, 2012 at 11:27 AM, Mark F. Adams <mark.adams at columbia.edu>wrote: > > On Mar 14, 2012, at 11:56 AM, John Mousel wrote: > > Mark, > > The matrix is asymmetric. Does this require the setting of an option? > > > Yes: -pc_gamg_sym_graph > > Mark > > I pulled petsc-dev this morning, so I should have (at least close to) the > latest code. > > John > > On Wed, Mar 14, 2012 at 10:54 AM, Mark F. Adams <mark.adams at > columbia.edu>wrote: > >> >> On Mar 14, 2012, at 11:08 AM, John Mousel wrote: >> >> I'm getting the following error when using GAMG. >> >> petsc-dev/src/ksp/pc/impls/gamg/agg.c:508: smoothAggs: Assertion >> `sgid==-1' failed. >> >> >> Is it possible that your matrix is structurally asymmetric? >> >> This code is evolving fast and so you will need to move to the dev >> version if you are not already using it. (I think I fixed a bug that hit >> this assert). >> >> >> When I try to alter the type of aggregation at the command line using >> -pc_gamg_type pa, I'm getting >> >> [0]PETSC ERROR: [1]PETSC ERROR: --------------------- Error Message >> ------------------------------------ >> [1]PETSC ERROR: Unknown type. Check for miss-spelling or missing external >> package needed for type: >> see http://www.mcs.anl.gov/petsc/documentation/installation.html#external >> ! >> [1]PETSC ERROR: Unknown GAMG type pa given! >> >> Has there been a change in the aggregation options? I just pulled >> petsc-dev this morning. >> >> >> Yes, this option is gone now. You can use -pc_gamg_type agg for now. >> >> Mark >> >> John >> >> >> > > -------------- next part -------------- An HTML attachment was scrubbed... URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20120314/26635f09/attachment.htm>
