[petsc-dev] GAMG MatMatMult error on 1 core

[Jed Brown]

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.html>