On Fri, Mar 29, 2013 at 10:44 AM, John Mousel <elafint.john at gmail.com>wrote:
> I switched the arguments to > > CALL > MatXAIJSetPreallocation(Ap,1,0,dnnz,0,onnz,PETSC_NULL_OBJECT,PETSC_NULL_OBJECT,ierr) > You're passing too many arguments. http://www.mcs.anl.gov/petsc/petsc-dev/docs/manualpages/Mat/MatXAIJSetPreallocation.html When I first introduced MatXAIJSetPreallocation, it had more arguments, but we trimmed it down early on. > > but it doesn't seem to like that; > > [0]PETSC ERROR: MatSetValues() line 1075 in > /home/jmousel/NumericalLibraries/petsc-hg/src/mat/interface/matrix.c > [0]PETSC ERROR: --------------------- Error Message > ------------------------------------ > [0]PETSC ERROR: Object is in wrong state! > [0]PETSC ERROR: Must call MatXXXSetPreallocation() or MatSetUp() on > argument 1 "mat" before MatSetValues()! > > > > On Fri, Mar 29, 2013 at 10:42 AM, Jed Brown <jedbrown at mcs.anl.gov> wrote: > >> On Fri, Mar 29, 2013 at 10:31 AM, John Mousel <elafint.john at >> gmail.com>wrote: >> >>> So this >>> >>> Ap => proj%A >>> CALL MatCreate(solver%comm,Ap,ierr) >>> CALL MatSetSizes(Ap,nlp,nlp,ngp,ngp,ierr) >>> CALL MatSetType(Ap,MATMPIAIJ,ierr) >>> CALL MatMPIAIJSetPreallocation(Ap,0,dnnz,0,onnz,ierr) >>> >>> should go to this >>> >>> Ap => proj%A >>> CALL MatCreate(solver%comm,Ap,ierr) >>> CALL MatSetSizes(Ap,nlp,nlp,ngp,ngp,ierr) >>> CALL MatSetType(Ap,MATAIJ,ierr) >>> >> >> Yes >> >> >>> CALL MatXAIJSetPreallocation(Ap,0,dnnz,0,onnz,ierr) >>> >> >> Almost, see the actual arguments to MatXAIJSetPreallocation. >> >> >>> >>> ? >>> >>> >>> >>> >>> On Fri, Mar 29, 2013 at 10:25 AM, Jed Brown <jedbrown at mcs.anl.gov>wrote: >>> >>>> On Fri, Mar 29, 2013 at 10:22 AM, John Mousel <john.mousel at >>>> gmail.com>wrote: >>>> >>>>> I get the following error when I try to use GAMG on 1 core. I just >>>>> pulled dev and it still exists. >>>>> >>>> >>>> Are you explicitly setting your matrix type to MPIAIJ? Can you use >>>> MATAIJ instead (it will choose the correct type based on size)? >>>> >>>> >>>>> >>>>> John >>>>> >>>>> [0]PETSC ERROR: --------------------- Error Message >>>>> ------------------------------------ >>>>> [0]PETSC ERROR: Arguments are incompatible! >>>>> [0]PETSC ERROR: MatMatMult requires A, mpiaij, to be compatible with >>>>> B, seqaij! >>>>> [0]PETSC ERROR: >>>>> ------------------------------------------------------------------------ >>>>> [0]PETSC ERROR: Petsc Development GIT revision: unknown GIT Date: >>>>> unknown >>>>> [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: ../src/MUSSEL_exe on a intel-debug named >>>>> wv.iihr.uiowa.edu by jmousel Fri Mar 29 10:19:44 2013 >>>>> [0]PETSC ERROR: Libraries linked from >>>>> /home/jmousel/NumericalLibraries/petsc-hg/intel-debug/lib >>>>> [0]PETSC ERROR: Configure run at Fri Mar 29 10:00:56 2013 >>>>> [0]PETSC ERROR: Configure options --download-blacs=1 >>>>> --download-metis=1 --download-ml=1 --download-mpich=1 >>>>> --download-parmetis=/home/jmousel/NumericalLibraries/petsc-hg/parmetis-4.0.2-p3.tar.gz >>>>> --download-f-blas-lapack=1 --with-cc=icc --download-cmake=1 >>>>> --with-cxx=icpc >>>>> --with-fc=ifort >>>>> [0]PETSC ERROR: >>>>> ------------------------------------------------------------------------ >>>>> [0]PETSC ERROR: MatMatMult() line 8567 in >>>>> /home/jmousel/NumericalLibraries/petsc-hg/src/mat/interface/matrix.c >>>>> [0]PETSC ERROR: PCGAMGOptprol_AGG() line 1335 in >>>>> /home/jmousel/NumericalLibraries/petsc-hg/src/ksp/pc/impls/gamg/agg.c >>>>> [0]PETSC ERROR: PCSetUp_GAMG() line 680 in >>>>> /home/jmousel/NumericalLibraries/petsc-hg/src/ksp/pc/impls/gamg/gamg.c >>>>> [0]PETSC ERROR: PCSetUp() line 890 in >>>>> /home/jmousel/NumericalLibraries/petsc-hg/src/ksp/pc/interface/precon.c >>>>> [0]PETSC ERROR: KSPSetUp() line 278 in >>>>> /home/jmousel/NumericalLibraries/petsc-hg/src/ksp/ksp/interface/itfunc.c >>>>> [0]PETSC ERROR: KSPSolve() line 399 in >>>>> /home/jmousel/NumericalLibraries/petsc-hg/src/ksp/ksp/interface/itfunc.c >>>>> >>>>> >>>> >>> >> > -------------- next part -------------- An HTML attachment was scrubbed... URL: <http://lists.mcs.anl.gov/pipermail/petsc-dev/attachments/20130329/427445e9/attachment-0001.html>
