[petsc-users] Nonconforming object sizes in axpy

Wed Feb 27 18:19:05 CST 2019

I see the error now. You should pass SUBSET_NONZERO_PATTERN as MatStructure argument to MatAXPY
DIFFERENT_NONZERO_PATTERN now triggers preallocation, and MATELEMENTAL does not support MatGetRow.

> On Feb 27, 2019, at 3:59 PM, Marius Buerkle <mbuerkle at web.de> wrote:
> 
> Here is the error message and please find an example attached. 
> thanks for the efforts !
>  
> [0]PETSC ERROR: --------------------- Error Message --------------------------------------------------------------
> [0]PETSC ERROR: No support for this operation for this object type
> [0]PETSC ERROR: Mat type elemental
> [0]PETSC ERROR: See http://www.mcs.anl.gov/petsc/documentation/faq.html for trouble shooting.
> [0]PETSC ERROR: Petsc Development GIT revision: v3.10.3-1681-g73e1f9d014  GIT Date: 2019-02-24 19:52:00 +0100
> [0]PETSC ERROR: ./a.out on a  named tono-hpc1 by marius Thu Feb 28 08:56:48 2019
> [0]PETSC ERROR: Configure options --prefix=/home/marius/prog/petsc/intel_openmpi_openmp_mkl_dbg_with_slepc --with-fortran-kernels=0 --with-64-bit-indices=0 --CC=mpicc --COPTFLAGS="-g -O0 -qopenmp -std=c11" --CXX=mpicxx --CXXOPTFLAGS="-g -O0 -qopenmp -std=c++14" --FC=mpif90 --FOPTFLAGS="-g -O0 -qopenmp -traceback " --with-mpi=1 --with-x=0 --download-parmetis=1 --download-metis=1 --download-superlu_dist=1 --download-superlu_dist-commit=93a9aa53142253f2a56891e09c6f950d250ae251 --download-superlu_dist-cmake-arguments=-Denable_openmp=TRUE --download-hwloc=1 --download-sowing=1 --with-openmp=0 --with-pthread=0 --download-elemental=1 --download-elemental-commit=de7b5bea1abf5f626b91582f742cf99e2e551bff --with-cxx-dialect=c++11 --with-scalar-type=complex --with-debugging=1 --with-valgrind=0 --with-blaslapack-lib=" -L/home/marius/intel/compilers_and_libraries_2019.1.144/linux/mkl/lib/intel64 -lmkl_scalapack_lp64 -lmkl_intel_lp64 -lmkl_intel_thread -lmkl_core -lmkl_blacs_openmpi_lp64 -liomp5 -lpthread -lm -ldl" --with-scalapack-lib=" -L/home/marius/intel/compilers_and_libraries_2019.1.144/linux/mkl/lib/intel64 -lmkl_scalapack_lp64 -lmkl_intel_lp64 -lmkl_intel_thread -lmkl_core -lmkl_blacs_openmpi_lp64 -liomp5 -lpthread -lm -ldl" --with-mkl_pardiso-dir=/home/marius/intel/compilers_and_libraries_2019.1.144/linux/mkl --with-mkl_cpardiso-dir=/home/marius/intel/compilers_and_libraries_2019.1.144/linux/mkl --with-mkl_sparse-dir=/home/marius/intel/compilers_and_libraries_2019.1.144/linux/mkl --with-mkl_sparse_optimize-dir=/home/marius/intel/compilers_and_libraries_2019.1.144/linux/mkl --with-valgrind=1
> [0]PETSC ERROR: #1 MatGetRow() line 570 in /home/marius/prog/petsc/git/petsc/src/mat/interface/matrix.c
> [0]PETSC ERROR: #2 MatAXPY_Basic_Preallocate() line 95 in /home/marius/prog/petsc/git/petsc/src/mat/utils/axpy.c
> [0]PETSC ERROR: #3 MatAXPY() line 51 in /home/marius/prog/petsc/git/petsc/src/mat/utils/axpy.c
> [1]PETSC ERROR: --------------------- Error Message --------------------------------------------------------------
> [1]PETSC ERROR: No support for this operation for this object type
> [1]PETSC ERROR: Mat type elemental
> [1]PETSC ERROR: See http://www.mcs.anl.gov/petsc/documentation/faq.html for trouble shooting.
> [1]PETSC ERROR: Petsc Development GIT revision: v3.10.3-1681-g73e1f9d014  GIT Date: 2019-02-24 19:52:00 +0100
> [1]PETSC ERROR: ./a.out on a  named tono-hpc1 by marius Thu Feb 28 08:56:48 2019
> [1]PETSC ERROR: Configure options --prefix=/home/marius/prog/petsc/intel_openmpi_openmp_mkl_dbg_with_slepc --with-fortran-kernels=0 --with-64-bit-indices=0 --CC=mpicc --COPTFLAGS="-g -O0 -qopenmp -std=c11" --CXX=mpicxx --CXXOPTFLAGS="-g -O0 -qopenmp -std=c++14" --FC=mpif90 --FOPTFLAGS="-g -O0 -qopenmp -traceback " --with-mpi=1 --with-x=0 --download-parmetis=1 --download-metis=1 --download-superlu_dist=1 --download-superlu_dist-commit=93a9aa53142253f2a56891e09c6f950d250ae251 --download-superlu_dist-cmake-arguments=-Denable_openmp=TRUE --download-hwloc=1 --download-sowing=1 --with-openmp=0 --with-pthread=0 --download-elemental=1 --download-elemental-commit=de7b5bea1abf5f626b91582f742cf99e2e551bff --with-cxx-dialect=c++11 --with-scalar-type=complex --with-debugging=1 --with-valgrind=0 --with-blaslapack-lib=" -L/home/marius/intel/compilers_and_libraries_2019.1.144/linux/mkl/lib/intel64 -lmkl_scalapack_lp64 -lmkl_intel_lp64 -lmkl_intel_thread -lmkl_core -lmkl_blacs_openmpi_lp64 -liomp5 -lpthread -lm -ldl" --with-scalapack-lib=" -L/home/marius/intel/compilers_and_libraries_2019.1.144/linux/mkl/lib/intel64 -lmkl_scalapack_lp64 -lmkl_intel_lp64 -lmkl_intel_thread -lmkl_core -lmkl_blacs_openmpi_lp64 -liomp5 -lpthread -lm -ldl" --with-mkl_pardiso-dir=/home/marius/intel/compilers_and_libraries_2019.1.144/linux/mkl --with-mkl_cpardiso-dir=/home/marius/intel/compilers_and_libraries_2019.1.144/linux/mkl --with-mkl_sparse-dir=/home/marius/intel/compilers_and_libraries_2019.1.144/linux/mkl --with-mkl_sparse_optimize-dir=/home/marius/intel/compilers_and_libraries_2019.1.144/linux/mkl --with-valgrind=1
> [1]PETSC ERROR: #1 MatGetRow() line 570 in /home/marius/prog/petsc/git/petsc/src/mat/interface/matrix.c
> [1]PETSC ERROR: #2 MatAXPY_Basic_Preallocate() line 95 in /home/marius/prog/petsc/git/petsc/src/mat/utils/axpy.c
> [1]PETSC ERROR: #3 MatAXPY() line 51 in /home/marius/prog/petsc/git/petsc/src/mat/utils/axpy.c
> [2]PETSC ERROR: --------------------- Error Message --------------------------------------------------------------
> [2]PETSC ERROR: No support for this operation for this object type
> [2]PETSC ERROR: Mat type elemental
> [2]PETSC ERROR: See http://www.mcs.anl.gov/petsc/documentation/faq.html for trouble shooting.
> [2]PETSC ERROR: Petsc Development GIT revision: v3.10.3-1681-g73e1f9d014  GIT Date: 2019-02-24 19:52:00 +0100
> [2]PETSC ERROR: ./a.out on a  named tono-hpc1 by marius Thu Feb 28 08:56:48 2019
> [2]PETSC ERROR: Configure options --prefix=/home/marius/prog/petsc/intel_openmpi_openmp_mkl_dbg_with_slepc --with-fortran-kernels=0 --with-64-bit-indices=0 --CC=mpicc --COPTFLAGS="-g -O0 -qopenmp -std=c11" --CXX=mpicxx --CXXOPTFLAGS="-g -O0 -qopenmp -std=c++14" --FC=mpif90 --FOPTFLAGS="-g -O0 -qopenmp -traceback " --with-mpi=1 --with-x=0 --download-parmetis=1 --download-metis=1 --download-superlu_dist=1 --download-superlu_dist-commit=93a9aa53142253f2a56891e09c6f950d250ae251 --download-superlu_dist-cmake-arguments=-Denable_openmp=TRUE --download-hwloc=1 --download-sowing=1 --with-openmp=0 --with-pthread=0 --download-elemental=1 --download-elemental-commit=de7b5bea1abf5f626b91582f742cf99e2e551bff --with-cxx-dialect=c++11 --with-scalar-type=complex --with-debugging=1 --with-valgrind=0 --with-blaslapack-lib=" -L/home/marius/intel/compilers_and_libraries_2019.1.144/linux/mkl/lib/intel64 -lmkl_scalapack_lp64 -lmkl_intel_lp64 -lmkl_intel_thread -lmkl_core -lmkl_blacs_openmpi_lp64 -liomp5 -lpthread -lm -ldl" --with-scalapack-lib=" -L/home/marius/intel/compilers_and_libraries_2019.1.144/linux/mkl/lib/intel64 -lmkl_scalapack_lp64 -lmkl_intel_lp64 -lmkl_intel_thread -lmkl_core -lmkl_blacs_openmpi_lp64 -liomp5 -lpthread -lm -ldl" --with-mkl_pardiso-dir=/home/marius/intel/compilers_and_libraries_2019.1.144/linux/mkl --with-mkl_cpardiso-dir=/home/marius/intel/compilers_and_libraries_2019.1.144/linux/mkl --with-mkl_sparse-dir=/home/marius/intel/compilers_and_libraries_2019.1.144/linux/mkl --with-mkl_sparse_optimize-dir=/home/marius/intel/compilers_and_libraries_2019.1.144/linux/mkl --with-valgrind=1
> [2]PETSC ERROR: #1 MatGetRow() line 570 in /home/marius/prog/petsc/git/petsc/src/mat/interface/matrix.c
> [2]PETSC ERROR: #2 MatAXPY_Basic_Preallocate() line 95 in /home/marius/prog/petsc/git/petsc/src/mat/utils/axpy.c
> [2]PETSC ERROR: #3 MatAXPY() line 51 in /home/marius/prog/petsc/git/petsc/src/mat/utils/axpy.c
> [3]PETSC ERROR: --------------------- Error Message --------------------------------------------------------------
> [3]PETSC ERROR: No support for this operation for this object type
> [3]PETSC ERROR: Mat type elemental
> [3]PETSC ERROR: See http://www.mcs.anl.gov/petsc/documentation/faq.html for trouble shooting.
> [3]PETSC ERROR: Petsc Development GIT revision: v3.10.3-1681-g73e1f9d014  GIT Date: 2019-02-24 19:52:00 +0100
> [3]PETSC ERROR: ./a.out on a  named tono-hpc1 by marius Thu Feb 28 08:56:48 2019
> [3]PETSC ERROR: Configure options --prefix=/home/marius/prog/petsc/intel_openmpi_openmp_mkl_dbg_with_slepc --with-fortran-kernels=0 --with-64-bit-indices=0 --CC=mpicc --COPTFLAGS="-g -O0 -qopenmp -std=c11" --CXX=mpicxx --CXXOPTFLAGS="-g -O0 -qopenmp -std=c++14" --FC=mpif90 --FOPTFLAGS="-g -O0 -qopenmp -traceback " --with-mpi=1 --with-x=0 --download-parmetis=1 --download-metis=1 --download-superlu_dist=1 --download-superlu_dist-commit=93a9aa53142253f2a56891e09c6f950d250ae251 --download-superlu_dist-cmake-arguments=-Denable_openmp=TRUE --download-hwloc=1 --download-sowing=1 --with-openmp=0 --with-pthread=0 --download-elemental=1 --download-elemental-commit=de7b5bea1abf5f626b91582f742cf99e2e551bff --with-cxx-dialect=c++11 --with-scalar-type=complex --with-debugging=1 --with-valgrind=0 --with-blaslapack-lib=" -L/home/marius/intel/compilers_and_libraries_2019.1.144/linux/mkl/lib/intel64 -lmkl_scalapack_lp64 -lmkl_intel_lp64 -lmkl_intel_thread -lmkl_core -lmkl_blacs_openmpi_lp64 -liomp5 -lpthread -lm -ldl" --with-scalapack-lib=" -L/home/marius/intel/compilers_and_libraries_2019.1.144/linux/mkl/lib/intel64 -lmkl_scalapack_lp64 -lmkl_intel_lp64 -lmkl_intel_thread -lmkl_core -lmkl_blacs_openmpi_lp64 -liomp5 -lpthread -lm -ldl" --with-mkl_pardiso-dir=/home/marius/intel/compilers_and_libraries_2019.1.144/linux/mkl --with-mkl_cpardiso-dir=/home/marius/intel/compilers_and_libraries_2019.1.144/linux/mkl --with-mkl_sparse-dir=/home/marius/intel/compilers_and_libraries_2019.1.144/linux/mkl --with-mkl_sparse_optimize-dir=/home/marius/intel/compilers_and_libraries_2019.1.144/linux/mkl --with-valgrind=1
> [3]PETSC ERROR: #1 MatGetRow() line 570 in /home/marius/prog/petsc/git/petsc/src/mat/interface/matrix.c
> [3]PETSC ERROR: #2 MatAXPY_Basic_Preallocate() line 95 in /home/marius/prog/petsc/git/petsc/src/mat/utils/axpy.c
> [3]PETSC ERROR: #3 MatAXPY() line 51 in /home/marius/prog/petsc/git/petsc/src/mat/utils/axpy.c
>  
> 
> Von: "Stefano Zampini" <stefano.zampini at gmail.com>
> An: "Marius Buerkle" <mbuerkle at web.de>
> Cc: "Matthew Knepley" <knepley at gmail.com>, "PETSc users list" <petsc-users at mcs.anl.gov>
> Betreff: Re: [petsc-users] Nonconforming object sizes in axpy
> 
> 
> > On Feb 27, 2019, at 3:25 PM, Marius Buerkle <mbuerkle at web.de> wrote:
> >
> > Ok. I also tried without openmp which gave the same error.
> > The segfault left aside, is MatAXPY supposed to work when Y is of type ELEMENTAL and X of type MATMPIAIJ ? It worked before but now it stops with mat elemental not supported.
> >
> 
> I cannot find the error line that says Not supported for elemental.
> The error PETSc is raising in the log you sent is with respect to the nonconforming object sizes (which BTW are correct).
> Could you please send the backtrace (without openmp) for the error regarding “mat elemental not supported”?
> 
> 
> 
> >
> >> Von: "Stefano Zampini" <stefano.zampini at gmail.com>
> >> An: "Marius Buerkle" <mbuerkle at web.de>
> >> Cc: "Matthew Knepley" <knepley at gmail.com>, "PETSc users list" <petsc-users at mcs.anl.gov>
> >> Betreff: Re: [petsc-users] Nonconforming object sizes in axpy
> >>
> >> If you check ALL ERROR CODES by doing
> >> ierr = _AnyPetscCall_();CHKERRQ(ierr);
> >> The code will stop at the first error
> >>
> >> I see you are using openmp. Probably this is causing the problem, I don’t know what to say without a MWE.
> >>
> >>
> >>> On Feb 25, 2019, at 10:50 PM, Marius Buerkle <mbuerkle at web.de> wrote:
> >>>
> >>> <error_output.tar.gz>
> >>
> >>
>  
> <output.tar.gz>

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20190227/246a419a/attachment-0001.html>