Please find attached the error output from petsc and the valgrind output. I could not reporduce the actual error in an MWE as it stops before calling MatLUFactorNumeric with the error at MatAXPY complaining that Mat type elemental is not supported, why it continues in my original code despite a nonzero error I don't know.
We need a stacktrace or an MWE to understand what is going on
On Feb 25, 2019, at 9:51 PM, Marius Buerkle <[email protected]> wrote:I checked the problem further. I am not sure if this is the cause but what happens is that I use MatAXPY with Y = MATELEMENTAL and X= MATMPIAIJ before calling MatLUFactorNumeric, it turns out that MatAXPY fails but did not raise an exception and I did not check the error code returned from MatAXPY. In my code it give error code 60 in a MWE I actually get error code 56, with[1]PETSC ERROR: --------------------- Error Message --------------------------------------------------------------
[1]PETSC ERROR: No support for this operation for this object type
[1]PETSC ERROR: Mat type elementalPrior to 646531bb84 MatAXPY worked with Y = MATELEMENTAL and X= MATMPIAIJ not raising a nonzero error code.I don't see any issue with the commit.Did you recompile petsc ? Can you send us a MWE to reproduce the error? Have you checked under valgrind?Il Lun 25 Feb 2019, 06:58 Marius Buerkle via petsc-users <[email protected]> ha scritto:No excatly same source on my side. It works with commits prior to 646531bb84 and fails starting from 646531bb84.
On Mon, Feb 25, 2019 at 9:00 AM Marius Buerkle via petsc-users <[email protected]> wrote:Hi,After commit 646531bb84 I get the following error message when calling MatLUFactorNumeric. It worked without a problem before.[6]PETSC ERROR: --------------------- Error Message --------------------------------------------------------------
[6]PETSC ERROR: Nonconforming object sizes
[6]PETSC ERROR: Non conforming matrix add: local sizes 51 x 54, 51 x 54This seems impossible since it reports above that the sizes are the same, but the error says they are different. Is it possible you sourcegot accidentally changed?Matt[6]PETSC ERROR: See http://www.mcs.anl.gov/petsc/documentation/faq.html for trouble shooting.
[6]PETSC ERROR: Petsc Development GIT revision: v3.10.3-1412-g646531bb84 GIT Date: 2019-02-07 21:47:41 +0300
[6]PETSC ERROR: /home/marius/prog/ownstuff/fortran/programs/transomat_dev/transomat/transomat on a named tono-hpc1 by marius Mon Feb 25 22:43:23 2019
[6]PETSC ERROR: Configure options --prefix=/home/marius/prog/petsc/intel_openmpi_openmp_mkl_opt_with_slepc --with-fortran-kernels=1 --with-64-bit-indices=0 --CC=mpi
cc --COPTFLAGS="-g -O0 -std=c11" --CXX=mpicxx --CXXOPTFLAGS="-g -O0 -std=c++14" --FC=mpif90 --FOPTFLAGS="-g -O0 -traceback " --with-mpi=1 --with-x=0 --download-p
armetis=1 --download-metis=1 --download-superlu_dist=1 --download-superlu_dist-commit=93a9aa53142253f2a56891e09c6f950d250ae251 --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-scala
r-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_scal
apack_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_an
d_libraries_2019.1.144/linux/mkl
[6]PETSC ERROR: #1 MatAXPY() line 39 in /home/marius/prog/petsc/git/petsc/src/mat/utils/axpy.cbest,Marus--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
error_output.tar.gz
Description: Binary data
