[petsc-users] Nonconforming object sizes in axpy

[Stefano Zampini]

We need a stacktrace or an MWE to understand what is going on

> On Feb 25, 2019, at 9:51 PM, Marius Buerkle <mbuerkle at web.de> 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 elemental
>  
>  
> Prior 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 <petsc-users at mcs.anl.gov <mailto:petsc-users at mcs.anl.gov>> 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 <petsc-users at mcs.anl.gov <mailto:petsc-users at mcs.anl.gov>> 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 54
>  
> This seems impossible since it reports above that the sizes are the same, but the error says they are different. Is it possible you source
> got accidentally changed?
>  
>   Matt
>  
> [6]PETSC ERROR: See http://www.mcs.anl.gov/petsc/documentation/faq.html <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.c
>  
> best,
> 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
>  
> 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/20190225/88fbb724/attachment-0001.html>