On Fri, Mar 29, 2013 at 10:31 AM, John Mousel <[email protected]>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 <[email protected]> wrote:
On Fri, Mar 29, 2013 at 10:22 AM, John Mousel <[email protected]>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