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) CALL MatXAIJSetPreallocation(Ap,0,dnnz,0,onnz,ierr) ? 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/9194fa24/attachment-0001.html>
