It seems MatFindOffBlockDiagonalEntries_MPIAIJ is broken for a special case. Can you send a MWE to reproduce problem so that we can fix it properly and add it to the testsuite?
Il giorno lun 22 giu 2020 alle ore 12:20 Eda Oktay <eda.ok...@metu.edu.tr> ha scritto: > I am so sorry for the inconvenience, the last error I sent you was a > different one, I forgot to change the codes. I am sending you the > actual error where I used MatFindOffBlockDiagonalEntries again in the > newest version. > > [1]PETSC ERROR: --------------------- Error Message > -------------------------------------------------------------- > [1]PETSC ERROR: Null argument, when expecting valid pointer > [2]PETSC ERROR: --------------------- Error Message > -------------------------------------------------------------- > [2]PETSC ERROR: Null argument, when expecting valid pointer > [2]PETSC ERROR: Null Object: Parameter # 1 > [2]PETSC ERROR: See > https://www.mcs.anl.gov/petsc/documentation/faq.html for trouble > shooting. > [2]PETSC ERROR: Petsc Release Version 3.13.2, Jun 02, 2020 > [2]PETSC ERROR: [3]PETSC ERROR: --------------------- Error Message > -------------------------------------------------------------- > [3]PETSC ERROR: Null argument, when expecting valid pointer > [3]PETSC ERROR: Null Object: Parameter # 1 > [3]PETSC ERROR: See > https://www.mcs.anl.gov/petsc/documentation/faq.html for trouble > shooting. > [3]PETSC ERROR: Petsc Release Version 3.13.2, Jun 02, 2020 > [3]PETSC ERROR: > ./approx_cut_deneme_clustering_son_final_edgecut_without_parmetis on a > arch-linux2-c-debug named b342.wls.metu.edu.tr by edaoktay Mon Jun 22 > 12:11:44 2020 > [1]PETSC ERROR: Null Object: Parameter # 1 > [1]PETSC ERROR: See > https://www.mcs.anl.gov/petsc/documentation/faq.html for trouble > shooting. > [1]PETSC ERROR: Petsc Release Version 3.13.2, Jun 02, 2020 > [1]PETSC ERROR: > ./approx_cut_deneme_clustering_son_final_edgecut_without_parmetis on a > arch-linux2-c-debug named b342.wls.metu.edu.tr by edaoktay Mon Jun 22 > 12:11:44 2020 > [1]PETSC ERROR: Configure options --with-cc=gcc --with-cxx=g++ > --download-mpich --download-openblas --download-slepc --download-metis > --download-parmetis --download-chaco --with-X=1 > [1]PETSC ERROR: #1 ISGetSize() line 1101 in > /home/edaoktay/petsc-3.13.2/src/vec/is/is/interface/index.c > ./approx_cut_deneme_clustering_son_final_edgecut_without_parmetis on a > arch-linux2-c-debug named b342.wls.metu.edu.tr by edaoktay Mon Jun 22 > 12:11:44 2020 > [2]PETSC ERROR: Configure options --with-cc=gcc --with-cxx=g++ > --download-mpich --download-openblas --download-slepc --download-metis > --download-parmetis --download-chaco --with-X=1 > [2]PETSC ERROR: #1 ISGetSize() line 1101 in > /home/edaoktay/petsc-3.13.2/src/vec/is/is/interface/index.c > [2]PETSC ERROR: [3]PETSC ERROR: Configure options --with-cc=gcc > --with-cxx=g++ --download-mpich --download-openblas --download-slepc > --download-metis --download-parmetis --download-chaco --with-X=1 > [3]PETSC ERROR: #1 ISGetSize() line 1101 in > /home/edaoktay/petsc-3.13.2/src/vec/is/is/interface/index.c > [3]PETSC ERROR: #2 MatFindOffBlockDiagonalEntries_MPIAIJ() line 229 in > /home/edaoktay/petsc-3.13.2/src/mat/impls/aij/mpi/mpiaij.c > [1]PETSC ERROR: #2 MatFindOffBlockDiagonalEntries_MPIAIJ() line 229 in > /home/edaoktay/petsc-3.13.2/src/mat/impls/aij/mpi/mpiaij.c > [1]PETSC ERROR: #3 MatFindOffBlockDiagonalEntries() line 9824 in > /home/edaoktay/petsc-3.13.2/src/mat/interface/matrix.c > #2 MatFindOffBlockDiagonalEntries_MPIAIJ() line 229 in > /home/edaoktay/petsc-3.13.2/src/mat/impls/aij/mpi/mpiaij.c > [2]PETSC ERROR: #3 MatFindOffBlockDiagonalEntries() line 9824 in > /home/edaoktay/petsc-3.13.2/src/mat/interface/matrix.c > [3]PETSC ERROR: #3 MatFindOffBlockDiagonalEntries() line 9824 in > /home/edaoktay/petsc-3.13.2/src/mat/interface/matrix.c > > I think this is again the same error when I used an older version. > Again, even though I am using a valid matrix and an index set, I don't > know why I am getting this error. > > About the small test example, I couldn't understand what Junchao ment > but this is the first time I am trying to use this function and my > current matrix is quite small 72*72. > > Mat Ais; > MatCreateSubMatrix(A,idUi,idUi,MAT_INITIAL_MATRIX,&Ais); > ISSetPermutation(idUi); > MatPermute(Ais,idUi,idUi,&PL); > IS offd; > MatFindOffBlockDiagonalEntries(PL,&offd); > > My matrix PL is working, I know that since I viewed it before the > line: IS offd. And although it says null object for parameter 1, I > believe that other lines below this error are about IS. Should I > create IS before using it in MatFindOffBlockDiagonalEntries? > > Mark Adams <mfad...@lbl.gov>, 19 Haz 2020 Cum, 20:37 tarihinde şunu yazdı: > > > > This is more useful and you are getting a new error. Maybe the old error > is fixed. > > It looks like you are trying to copy a matrix into another but they have > different communicators. For example, one is created with MPI_COMM_WOLD > and the other with PETSC_COMM_WORLD. > > > > On Fri, Jun 19, 2020 at 12:48 PM Eda Oktay <eda.ok...@metu.edu.tr> > wrote: > >> > >> Dear Mark, > >> > >> I updated Petsc and again, here's the error: > >> > >> [0]PETSC ERROR: --------------------- Error Message > >> -------------------------------------------------------------- > >> [0]PETSC ERROR: Arguments must have same communicators > >> [0]PETSC ERROR: [1]PETSC ERROR: --------------------- Error Message > >> -------------------------------------------------------------- > >> [1]PETSC ERROR: Arguments must have same communicators > >> [1]PETSC ERROR: Different communicators in the two objects: Argument # > >> 1 and 2 flag 3 > >> [1]PETSC ERROR: See > >> https://www.mcs.anl.gov/petsc/documentation/faq.html for trouble > >> shooting. > >> [1]PETSC ERROR: Petsc Release Version 3.13.2, Jun 02, 2020 > >> [1]PETSC ERROR: [2]PETSC ERROR: [3]PETSC ERROR: --------------------- > >> Error Message > -------------------------------------------------------------- > >> [3]PETSC ERROR: Arguments must have same communicators > >> [3]PETSC ERROR: Different communicators in the two objects: Argument # > >> 1 and 2 flag 3 > >> [3]PETSC ERROR: See > >> https://www.mcs.anl.gov/petsc/documentation/faq.html for trouble > >> shooting. > >> [3]PETSC ERROR: Petsc Release Version 3.13.2, Jun 02, 2020 > >> [3]PETSC ERROR: > >> ./approx_cut_deneme_clustering_son_final_edgecut_without_parmetis on a > >> arch-linux2-c-debug named b342.wls.metu.edu.tr by edaoktay Fri Jun 19 > >> 19:44:41 2020 > >> [3]PETSC ERROR: Configure options --with-cc=gcc --with-cxx=g++ > >> --download-mpich --download-openblas --download-slepc --download-metis > >> --download-parmetis --download-chaco --with-X=1 > >> Different communicators in the two objects: Argument # 1 and 2 flag 3 > >> [0]PETSC ERROR: See > >> https://www.mcs.anl.gov/petsc/documentation/faq.html for trouble > >> shooting. > >> [0]PETSC ERROR: Petsc Release Version 3.13.2, Jun 02, 2020 > >> [0]PETSC ERROR: > >> ./approx_cut_deneme_clustering_son_final_edgecut_without_parmetis on a > >> arch-linux2-c-debug named b342.wls.metu.edu.tr by edaoktay Fri Jun 19 > >> 19:44:41 2020 > >> [0]PETSC ERROR: Configure options --with-cc=gcc --with-cxx=g++ > >> --download-mpich --download-openblas --download-slepc --download-metis > >> --download-parmetis --download-chaco --with-X=1 > >> [0]PETSC ERROR: #1 MatCopy() line 4044 in > >> /home/edaoktay/petsc-3.13.2/src/mat/interface/matrix.c > >> ./approx_cut_deneme_clustering_son_final_edgecut_without_parmetis on a > >> arch-linux2-c-debug named b342.wls.metu.edu.tr by edaoktay Fri Jun 19 > >> 19:44:41 2020 > >> [1]PETSC ERROR: Configure options --with-cc=gcc --with-cxx=g++ > >> --download-mpich --download-openblas --download-slepc --download-metis > >> --download-parmetis --download-chaco --with-X=1 > >> [1]PETSC ERROR: #1 MatCopy() line 4044 in > >> /home/edaoktay/petsc-3.13.2/src/mat/interface/matrix.c > >> [3]PETSC ERROR: #1 MatCopy() line 4044 in > >> /home/edaoktay/petsc-3.13.2/src/mat/interface/matrix.c > >> [0]PETSC ERROR: --------------------- Error Message > >> -------------------------------------------------------------- > >> [0]PETSC ERROR: Nonconforming object sizes > >> [1]PETSC ERROR: --------------------- Error Message > >> -------------------------------------------------------------- > >> [1]PETSC ERROR: Nonconforming object sizes > >> [1]PETSC ERROR: Non conforming matrix add: global sizes 1 x 72, 1 x 72 > >> [3]PETSC ERROR: --------------------- Error Message > >> -------------------------------------------------------------- > >> [3]PETSC ERROR: Nonconforming object sizes > >> [3]PETSC ERROR: Non conforming matrix add: global sizes 16 x 72, 16 x 72 > >> [3]PETSC ERROR: [0]PETSC ERROR: Non conforming matrix add: global > >> sizes 52 x 72, 52 x 72 > >> [0]PETSC ERROR: See > >> https://www.mcs.anl.gov/petsc/documentation/faq.html for trouble > >> shooting. > >> [0]PETSC ERROR: Petsc Release Version 3.13.2, Jun 02, 2020 > >> [0]PETSC ERROR: [1]PETSC ERROR: See > >> https://www.mcs.anl.gov/petsc/documentation/faq.html for trouble > >> shooting. > >> [1]PETSC ERROR: Petsc Release Version 3.13.2, Jun 02, 2020 > >> [1]PETSC ERROR: > >> ./approx_cut_deneme_clustering_son_final_edgecut_without_parmetis on a > >> arch-linux2-c-debug named b342.wls.metu.edu.tr by edaoktay Fri Jun 19 > >> 19:44:41 2020 > >> See https://www.mcs.anl.gov/petsc/documentation/faq.html for trouble > shooting. > >> [3]PETSC ERROR: Petsc Release Version 3.13.2, Jun 02, 2020 > >> [3]PETSC ERROR: > >> ./approx_cut_deneme_clustering_son_final_edgecut_without_parmetis on a > >> arch-linux2-c-debug named b342.wls.metu.edu.tr by edaoktay Fri Jun 19 > >> 19:44:41 2020 > >> [3]PETSC ERROR: > >> ./approx_cut_deneme_clustering_son_final_edgecut_without_parmetis on a > >> arch-linux2-c-debug named b342.wls.metu.edu.tr by edaoktay Fri Jun 19 > >> 19:44:41 2020 > >> [0]PETSC ERROR: Configure options --with-cc=gcc --with-cxx=g++ > >> --download-mpich --download-openblas --download-slepc --download-metis > >> --download-parmetis --download-chaco --with-X=1 > >> [0]PETSC ERROR: [1]PETSC ERROR: Configure options --with-cc=gcc > >> --with-cxx=g++ --download-mpich --download-openblas --download-slepc > >> --download-metis --download-parmetis --download-chaco --with-X=1 > >> [1]PETSC ERROR: #2 MatAXPY() line 72 in > >> /home/edaoktay/petsc-3.13.2/src/mat/utils/axpy.c > >> Configure options --with-cc=gcc --with-cxx=g++ --download-mpich > >> --download-openblas --download-slepc --download-metis > >> --download-parmetis --download-chaco --with-X=1 > >> [3]PETSC ERROR: #2 MatAXPY() line 72 in > >> /home/edaoktay/petsc-3.13.2/src/mat/utils/axpy.c > >> #2 MatAXPY() line 72 in /home/edaoktay/petsc-3.13.2/src/mat/utils/axpy.c > >> > >> Elapsed time: 2.5e-01 > >> > >> --------------------- Error Message > >> -------------------------------------------------------------- > >> [2]PETSC ERROR: Arguments must have same communicators > >> [2]PETSC ERROR: Different communicators in the two objects: Argument # > >> 1 and 2 flag 3 > >> [2]PETSC ERROR: See > >> https://www.mcs.anl.gov/petsc/documentation/faq.html for trouble > >> shooting. > >> [2]PETSC ERROR: Petsc Release Version 3.13.2, Jun 02, 2020 > >> [2]PETSC ERROR: > >> ./approx_cut_deneme_clustering_son_final_edgecut_without_parmetis on a > >> arch-linux2-c-debug named b342.wls.metu.edu.tr by edaoktay Fri Jun 19 > >> 19:44:41 2020 > >> [2]PETSC ERROR: Configure options --with-cc=gcc --with-cxx=g++ > >> --download-mpich --download-openblas --download-slepc --download-metis > >> --download-parmetis --download-chaco --with-X=1 > >> [2]PETSC ERROR: #1 MatCopy() line 4044 in > >> /home/edaoktay/petsc-3.13.2/src/mat/interface/matrix.c > >> [2]PETSC ERROR: --------------------- Error Message > >> -------------------------------------------------------------- > >> [2]PETSC ERROR: Nonconforming object sizes > >> [2]PETSC ERROR: Non conforming matrix add: global sizes 3 x 72, 3 x 72 > >> [2]PETSC ERROR: See > >> https://www.mcs.anl.gov/petsc/documentation/faq.html for trouble > >> shooting. > >> [2]PETSC ERROR: Petsc Release Version 3.13.2, Jun 02, 2020 > >> [2]PETSC ERROR: > >> ./approx_cut_deneme_clustering_son_final_edgecut_without_parmetis on a > >> arch-linux2-c-debug named b342.wls.metu.edu.tr by edaoktay Fri Jun 19 > >> 19:44:41 2020 > >> [2]PETSC ERROR: Configure options --with-cc=gcc --with-cxx=g++ > >> --download-mpich --download-openblas --download-slepc --download-metis > >> --download-parmetis --download-chaco --with-X=1 > >> [2]PETSC ERROR: #2 MatAXPY() line 72 in > >> /home/edaoktay/petsc-3.13.2/src/mat/utils/axpy.c > >> > >> What can be the reason? > >> > >> Mark Adams <mfad...@lbl.gov>, 19 Haz 2020 Cum, 18:41 tarihinde şunu > yazdı: > >> > > >> > I don't know what is going on here. There was an update to this > function about a year ago, so that might fix your problem. > >> > > >> > We would need you to test with a current version. > >> > > >> > Mark > >> > > >> > On Fri, Jun 19, 2020 at 11:23 AM Eda Oktay <eda.ok...@metu.edu.tr> > wrote: > >> >> > >> >> Hi all, > >> >> > >> >> I am trying to find off block diagonal entries of a matrix and I am > >> >> trying to use MatFindOffBlockDiagonalEntries. However, although my > >> >> matrix is not NULL, I am getting an error message as follows: > >> >> > >> >> 1]PETSC ERROR: --------------------- Error Message > >> >> -------------------------------------------------------------- > >> >> [1]PETSC ERROR: [2]PETSC ERROR: --------------------- Error Message > >> >> -------------------------------------------------------------- > >> >> [2]PETSC ERROR: Null argument, when expecting valid pointer > >> >> [2]PETSC ERROR: Null Object: Parameter # 1 > >> >> [2]PETSC ERROR: See > >> >> http://www.mcs.anl.gov/petsc/documentation/faq.html for trouble > >> >> shooting. > >> >> [2]PETSC ERROR: Petsc Release Version 3.11.1, Apr, 12, 2019 > >> >> [2]PETSC ERROR: [3]PETSC ERROR: --------------------- Error Message > >> >> -------------------------------------------------------------- > >> >> [3]PETSC ERROR: Null argument, when expecting valid pointer > >> >> [3]PETSC ERROR: Null Object: Parameter # 1 > >> >> [3]PETSC ERROR: See > >> >> http://www.mcs.anl.gov/petsc/documentation/faq.html for trouble > >> >> shooting. > >> >> [3]PETSC ERROR: Null argument, when expecting valid pointer > >> >> [1]PETSC ERROR: Null Object: Parameter # 1 > >> >> [1]PETSC ERROR: See > >> >> http://www.mcs.anl.gov/petsc/documentation/faq.html for trouble > >> >> shooting. > >> >> [1]PETSC ERROR: Petsc Release Version 3.11.1, Apr, 12, 2019 > >> >> [1]PETSC ERROR: > >> >> ./approx_cut_deneme_clustering_son_final_edgecut_without_parmetis on > a > >> >> arch-linux2-c-debug named b342.wls.metu.edu.tr by edaoktay Fri Jun > 19 > >> >> 18:19:11 2020 > >> >> [1]PETSC ERROR: Configure options --download-mpich > --download-openblas > >> >> --download-slepc --download-metis --download-parmetis > --download-chaco > >> >> --with-X=1 > >> >> [1]PETSC ERROR: #1 ISGetSize() line 510 in > >> >> /home/edaoktay/petsc-3.11.1/src/vec/is/is/interface/index.c > >> >> [1]PETSC ERROR: > >> >> ./approx_cut_deneme_clustering_son_final_edgecut_without_parmetis on > a > >> >> arch-linux2-c-debug named b342.wls.metu.edu.tr by edaoktay Fri Jun > 19 > >> >> 18:19:11 2020 > >> >> [2]PETSC ERROR: Configure options --download-mpich > --download-openblas > >> >> --download-slepc --download-metis --download-parmetis > --download-chaco > >> >> --with-X=1 > >> >> [2]PETSC ERROR: #1 ISGetSize() line 510 in > >> >> /home/edaoktay/petsc-3.11.1/src/vec/is/is/interface/index.c > >> >> [2]PETSC ERROR: #2 MatFindOffBlockDiagonalEntries_MPIAIJ() line 211 > in > >> >> /home/edaoktay/petsc-3.11.1/src/mat/impls/aij/mpi/mpiaij.c > >> >> [2]PETSC ERROR: #2 MatFindOffBlockDiagonalEntries_MPIAIJ() line 211 > in > >> >> /home/edaoktay/petsc-3.11.1/src/mat/impls/aij/mpi/mpiaij.c > >> >> [1]PETSC ERROR: #3 MatFindOffBlockDiagonalEntries() line 10502 in > >> >> /home/edaoktay/petsc-3.11.1/src/mat/interface/matrix.c > >> >> #3 MatFindOffBlockDiagonalEntries() line 10502 in > >> >> /home/edaoktay/petsc-3.11.1/src/mat/interface/matrix.c > >> >> Petsc Release Version 3.11.1, Apr, 12, 2019 > >> >> [3]PETSC ERROR: > >> >> ./approx_cut_deneme_clustering_son_final_edgecut_without_parmetis on > a > >> >> arch-linux2-c-debug named b342.wls.metu.edu.tr by edaoktay Fri Jun > 19 > >> >> 18:19:11 2020 > >> >> [3]PETSC ERROR: Configure options --download-mpich > --download-openblas > >> >> --download-slepc --download-metis --download-parmetis > --download-chaco > >> >> --with-X=1 > >> >> [3]PETSC ERROR: #1 ISGetSize() line 510 in > >> >> /home/edaoktay/petsc-3.11.1/src/vec/is/is/interface/index.c > >> >> [3]PETSC ERROR: #2 MatFindOffBlockDiagonalEntries_MPIAIJ() line 211 > in > >> >> /home/edaoktay/petsc-3.11.1/src/mat/impls/aij/mpi/mpiaij.c > >> >> [3]PETSC ERROR: #3 MatFindOffBlockDiagonalEntries() line 10502 in > >> >> /home/edaoktay/petsc-3.11.1/src/mat/interface/matrix.c > >> >> > >> >> And the program is still running. > >> >> > >> >> The part of my program is : > >> >> > >> >> Mat Ais; > >> >> MatCreateSubMatrix(A,idUi,idUi,MAT_INITIAL_MATRIX,&Ais); > >> >> ISSetPermutation(idUi); > >> >> MatPermute(Ais,idUi,idUi,&PL); > >> >> IS offd; > >> >> MatFindOffBlockDiagonalEntries(PL,&offd); > >> >> > >> >> Since there isn't any examples for this function, I couldn't > >> >> understand what should I do. Since it says in the error that it is > >> >> about parameter 1, the matrix, I don't know how to solve the problem. > >> >> My matrix PL is parallel sparse. I know there shouldn't be any > problem > >> >> in PL since I can view PL by using MatView. > >> >> > >> >> Thanks! > >> >> > >> >> Eda > -- Stefano
