Matt, It's not technically an error. It's an inefficient preallocation. I'm assuming Mark doesn't want this to be happening.
[0] PCSetUp(): Setting up new PC [0] PetscCommDuplicate(): Using internal PETSc communicator 1140850689 -2080374777 [2] PetscCommDuplicate(): Using internal PETSc communicator 1140850689 -2080374780 [1] PetscCommDuplicate(): Using internal PETSc communicator 1140850689 -2080374780 [3] PetscCommDuplicate(): Using internal PETSc communicator 1140850689 -2080374780 [0] PetscCommDuplicate(): Using internal PETSc communicator 1140850689 -2080374777 [2] PetscCommDuplicate(): Using internal PETSc communicator 1140850689 -2080374780 [1] PetscCommDuplicate(): Using internal PETSc communicator 1140850689 -2080374780 [3] PetscCommDuplicate(): Using internal PETSc communicator 1140850689 -2080374780 [1] MatAssemblyBegin_MPIAIJ(): Stash has 1830 entries, uses 0 mallocs. [3] MatAssemblyBegin_MPIAIJ(): Stash has 914 entries, uses 0 mallocs. [2] MatAssemblyBegin_MPIAIJ(): Stash has 1991 entries, uses 0 mallocs. [0] MatStashScatterBegin_Private(): No of messages: 3 [0] MatStashScatterBegin_Private(): Mesg_to: 1: size: 16200 [0] MatStashScatterBegin_Private(): Mesg_to: 2: size: 4192 [0] MatStashScatterBegin_Private(): Mesg_to: 3: size: 1800 [0] MatAssemblyBegin_MPIAIJ(): Stash has 2771 entries, uses 0 mallocs. [2] MatAssemblyEnd_SeqAIJ(): Matrix size: 10840 X 10840; storage space: 0 unneeded,119860 used [2] MatAssemblyEnd_SeqAIJ(): Number of mallocs during MatSetValues() is 0 [2] MatAssemblyEnd_SeqAIJ(): Maximum nonzeros in any row is 26 [2] Mat_CheckInode(): Found 10840 nodes out of 10840 rows. Not using Inode routines [3] MatAssemblyEnd_SeqAIJ(): Matrix size: 10862 X 10862; storage space: 0 unneeded,117330 used [3] MatAssemblyEnd_SeqAIJ(): Number of mallocs during MatSetValues() is 0 [3] MatAssemblyEnd_SeqAIJ(): Maximum nonzeros in any row is 26 [3] Mat_CheckInode(): Found 10862 nodes out of 10862 rows. Not using Inode routines [1] MatAssemblyEnd_SeqAIJ(): Matrix size: 12138 X 12138; storage space: 0 unneeded,116975 used [1] MatAssemblyEnd_SeqAIJ(): Number of mallocs during MatSetValues() is 0 [1] MatAssemblyEnd_SeqAIJ(): Maximum nonzeros in any row is 26 [1] Mat_CheckInode(): Found 12138 nodes out of 12138 rows. Not using Inode routines [0] MatAssemblyEnd_SeqAIJ(): Matrix size: 10765 X 10765; storage space: 0 unneeded,110555 used [0] MatAssemblyEnd_SeqAIJ(): Number of mallocs during MatSetValues() is 0 [0] MatAssemblyEnd_SeqAIJ(): Maximum nonzeros in any row is 25 [0] Mat_CheckInode(): Found 10765 nodes out of 10765 rows. Not using Inode routines [3] PetscCommDuplicate(): Using internal PETSc communicator 1140850689 -2080374780 [2] PetscCommDuplicate(): Using internal PETSc communicator 1140850689 -2080374780 [1] PetscCommDuplicate(): Using internal PETSc communicator 1140850689 -2080374780 [0] PetscCommDuplicate(): Using internal PETSc communicator 1140850689 -2080374777 [0] PetscCommDuplicate(): Using internal PETSc communicator 1140850689 -2080374777 [2] PetscCommDuplicate(): Using internal PETSc communicator 1140850689 -2080374780 [3] PetscCommDuplicate(): Using internal PETSc communicator 1140850689 -2080374780 [1] PetscCommDuplicate(): Using internal PETSc communicator 1140850689 -2080374780 [0] VecScatterCreateCommon_PtoS(): Using blocksize 1 scatter [0] VecScatterCreate(): General case: MPI to Seq [2] MatAssemblyEnd_SeqAIJ(): Matrix size: 10840 X 566; storage space: 118775 unneeded,1100 used [2] MatAssemblyEnd_SeqAIJ(): Number of mallocs during MatSetValues() is 1 [2] MatAssemblyEnd_SeqAIJ(): Maximum nonzeros in any row is 10 [3] MatAssemblyEnd_SeqAIJ(): Matrix size: 10862 X 259; storage space: 116725 unneeded,650 used [3] MatAssemblyEnd_SeqAIJ(): Number of mallocs during MatSetValues() is 3 [3] MatAssemblyEnd_SeqAIJ(): Maximum nonzeros in any row is 11 [0] MatAssemblyEnd_SeqAIJ(): Matrix size: 10765 X 1295; storage space: 107825 unneeded,3165 used [0] MatAssemblyEnd_SeqAIJ(): Number of mallocs during MatSetValues() is 29 [0] MatAssemblyEnd_SeqAIJ(): Maximum nonzeros in any row is 13 [1] MatAssemblyEnd_SeqAIJ(): Matrix size: 12138 X 1131; storage space: 114609 unneeded,2591 used [1] MatAssemblyEnd_SeqAIJ(): Number of mallocs during MatSetValues() is 15 [1] MatAssemblyEnd_SeqAIJ(): Maximum nonzeros in any row is 13 John On Tue, Sep 25, 2012 at 5:41 AM, Matthew Knepley <knepley at gmail.com> wrote: > On Mon, Sep 24, 2012 at 9:48 PM, John Mousel <john.mousel at gmail.com>wrote: > >> I'm having a problem with GAMG causing a MatSetValues() is [some non-zero >> number] after >> > > I do not see any errors in your output. > > Matt > > >> [0] PCSetUp(): Setting up new PC >> >> I don't get this with either ILU or HYPRE/BoomerAMG. I've attached the >> output of -info. I see there is a recent change to some pre-allocation for >> GAMG here: >> >> http://petsc.cs.iit.edu/petsc/petsc-dev/rev/fe806009181b >> >> I'm not sure if that has anything to do with it, but my blocksize is 1. >> I'm running with the following options. >> >> -ksp_type bcgsl -pc_type gamg -pc_gamg_threshold 0.01 -mg_levels_ksp_type >> richardson -mg_levels_pc_type sor -mg_coarse_ksp_type richardson >> -mg_coarse_pc_type sor -mg_coarse_pc_sor_its 4 -pc_gamg_type agg >> -pc_gamg_agg_nsmooths 1 -pc_gamg_sym_graph true -info >> -ksp_monitor_true_residual >> >> Thanks, >> John >> >> >> >> >> >> >> > > > -- > 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 > -------------- next part -------------- An HTML attachment was scrubbed... URL: <http://lists.mcs.anl.gov/pipermail/petsc-dev/attachments/20120925/558fdcc2/attachment.html>
