[petsc-users] Argument out of range error in MatPermute

[Eda Oktay]

I attached whole program I wrote where the problem is in line 285. One of
the matrices I used was airfoil1_binary, included in the folder. Also, I
included makefile. Is that what you want?

Matthew Knepley <knepley at gmail.com>, 25 Mar 2019 Pzt, 13:41 tarihinde şunu
yazdı:

> That should not happen. Can you send in a small example that we can debug.
>
>   Thanks,
>
>     Matt
>
> On Mon, Mar 25, 2019 at 12:38 AM Eda Oktay via petsc-users <
> petsc-users at mcs.anl.gov> wrote:
>
>> Hello,
>>
>> I am trying to permute a vector A using following lines:
>>
>> ierr =
>> ISCreateGeneral(PETSC_COMM_SELF,siz,idx,PETSC_COPY_VALUES,&is);CHKERRQ(ierr);
>>   ierr = ISSetPermutation(is);CHKERRQ(ierr);
>>   ierr = ISDuplicate(is,&newIS);CHKERRQ(ierr);
>>   ierr =
>> MatSetOption(A,MAT_NEW_NONZERO_ALLOCATION_ERR,PETSC_FALSE);CHKERRQ(ierr);
>>   ierr = MatPermute(A,is,newIS,&PL);CHKERRQ(ierr);
>>
>> However, in MatPermute line, I get the following error even if I used
>> MatSetOption before this line:
>>
>> [0]PETSC ERROR: --------------------- Error Message
>> --------------------------------------------------------------
>> [0]PETSC ERROR: Argument out of range
>> [0]PETSC ERROR: New nonzero at (0,485) caused a malloc
>> Use MatSetOption(A, MAT_NEW_NONZERO_ALLOCATION_ERR, PETSC_FALSE) to turn
>> off this check
>> [0]PETSC ERROR: See http://www.mcs.anl.gov/petsc/documentation/faq.html
>> for trouble shooting.
>> [0]PETSC ERROR: Petsc Release Version 3.10.3, Dec, 18, 2018
>> [0]PETSC ERROR: ./DENEME_TEMIZ_ENYENI_FINAL on a arch-linux2-c-debug
>> named 1232.wls.metu.edu.tr by edaoktay Mon Mar 25 12:15:14 2019
>> [0]PETSC ERROR: Configure options --with-cc=gcc --with-cxx=g++
>> --with-fc=gfortran --with-cxx-dialect=C++11 --download-openblas
>> --download-metis --download-parmetis --download-superlu_dist
>> --download-slepc --download-mpich
>> [0]PETSC ERROR: #1 MatSetValues_MPIAIJ() line 579 in
>> /home/edaoktay/petsc-3.10.3/src/mat/impls/aij/mpi/mpiaij.c
>> [0]PETSC ERROR: #2 MatAssemblyEnd_MPIAIJ() line 807 in
>> /home/edaoktay/petsc-3.10.3/src/mat/impls/aij/mpi/mpiaij.c
>> [0]PETSC ERROR: #3 MatAssemblyEnd() line 5340 in
>> /home/edaoktay/petsc-3.10.3/src/mat/interface/matrix.c
>> [0]PETSC ERROR: #4 MatPermute_MPIAIJ() line 1723 in
>> /home/edaoktay/petsc-3.10.3/src/mat/impls/aij/mpi/mpiaij.c
>> [0]PETSC ERROR: #5 MatPermute() line 4997 in
>> /home/edaoktay/petsc-3.10.3/src/mat/interface/matrix.c
>> [0]PETSC ERROR: #6 main() line 285 in
>> /home/edaoktay/petsc-3.10.3/arch-linux2-c-debug/share/slepc/examples/src/eda/DENEME_TEMIZ_ENYENI_FINAL.c
>> [0]PETSC ERROR: PETSc Option Table entries:
>> [0]PETSC ERROR: -f
>> /home/edaoktay/petsc-3.10.3/share/petsc/datafiles/matrices/binary_files/airfoil1_binary
>> [0]PETSC ERROR: -mat_partitioning_type parmetis
>> [0]PETSC ERROR: -weighted
>> [0]PETSC ERROR: ----------------End of Error Message -------send entire
>> error message to petsc-maint at mcs.anl.gov----------
>>
>> I'll be glad if you can help me.
>>
>> Thanks!
>>
>> Eda
>>
>
>
> --
> 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
>
> https://www.cse.buffalo.edu/~knepley/
> <http://www.cse.buffalo.edu/~knepley/>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20190325/dda7faa6/attachment-0001.html>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: example.zip
Type: application/zip
Size: 42504 bytes
Desc: not available
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20190325/dda7faa6/attachment-0001.zip>