[petsc-users] Possible Bug ?

[Smith, Barry F.]

   I'll take a look at it. Compiling now. I note that this example is currently only run in the tests in complex, not sure why it is not tested in real

   Barry


> On Apr 22, 2019, at 5:00 PM, Sajid Ali via petsc-users <petsc-users at mcs.anl.gov> wrote:
> 
> Hi, 
> 
> I see that src/mat/examples/tests/ex112.c is failing for petsc at 3.11.1 configured without complex scalars. With complex scalars everything works fine. 
> 
> The error I see is : 
> ```
> [sajid at xrmlite bugfix]$  ./ex112                                                                                                                    
> No protocol specified                                                                                                                              
>                                                                                                                                                    
>  1-D: FFTW on vector of size 10                                                                                                                    
>   Error norm of |x - z| 5.37156                                                                                                                    
>   Error norm of |x - z| 5.90871                                                                                                                    
>   Error norm of |x - z| 5.96243                                                                                                                    
> [0]PETSC ERROR: --------------------- Error Message --------------------------------------------------------------                                 
> [0]PETSC ERROR: Nonconforming object sizes                                                                                                         
> [0]PETSC ERROR: Mat mat,Vec x: global dim 20 10                                                                                                    
> [0]PETSC ERROR: See http://www.mcs.anl.gov/petsc/documentation/faq.html for trouble shooting.                                                      
> [0]PETSC ERROR: Petsc Release Version 3.11.1, Apr, 12, 2019                                                                                        
> [0]PETSC ERROR: ./ex112 on a  named xrmlite.phys.northwestern.edu by sajid Mon Apr 22 16:58:41 2019                                                
> [0]PETSC ERROR: Configure options --prefix=/home/sajid/packages/spack/opt/spack/linux-centos7-x86_64/gcc-8.3.0/petsc-3.11.1-5bdbcozu3labtbbi7gtq4xa
> knay24lo6 --with-ssl=0 --download-c2html=0 --download-sowing=0 --download-hwloc=0 CFLAGS="-O2 -march=native" FFLAGS="-O2 -march=native" CXXFLAGS="-
> O2 -march=native" --with-cc=/home/sajid/packages/spack/opt/spack/linux-centos7-x86_64/gcc-8.3.0/mpich-3.3-ig4cr2xw2x63bqs5rnmhfshln4iv7av5/bin/mpic
> c --with-cxx=/home/sajid/packages/spack/opt/spack/linux-centos7-x86_64/gcc-8.3.0/mpich-3.3-ig4cr2xw2x63bqs5rnmhfshln4iv7av5/bin/mpic++ --with-fc=/h
> ome/sajid/packages/spack/opt/spack/linux-centos7-x86_64/gcc-8.3.0/mpich-3.3-ig4cr2xw2x63bqs5rnmhfshln4iv7av5/bin/mpif90 --with-precision=double --w
> ith-scalar-type=real --with-shared-libraries=1 --with-debugging=1 --with-64-bit-indices=0 --with-blaslapack-lib="/home/sajid/packages/spack/opt/spa
> ck/linux-centos7-x86_64/gcc-8.3.0/intel-mkl-2019.3.199-kzcly5rtcjbkwtnm6tri6kkexnwoat5m/compilers_and_libraries_2019.3.199/linux/mkl/lib/intel64/li
> bmkl_intel_lp64.so /home/sajid/packages/spack/opt/spack/linux-centos7-x86_64/gcc-8.3.0/intel-mkl-2019.3.199-kzcly5rtcjbkwtnm6tri6kkexnwoat5m/compil
> ers_and_libraries_2019.3.199/linux/mkl/lib/intel64/libmkl_sequential.so /home/sajid/packages/spack/opt/spack/linux-centos7-x86_64/gcc-8.3.0/intel-m
> kl-2019.3.199-kzcly5rtcjbkwtnm6tri6kkexnwoat5m/compilers_and_libraries_2019.3.199/linux/mkl/lib/intel64/libmkl_core.so /lib64/libpthread.so /lib64/
> libm.so /lib64/libdl.so" --with-x=0 --with-clanguage=C --with-scalapack=0 --with-metis=0 --with-hdf5=0 --with-hypre=0 --with-parmetis=0 --with-mump
> s=0 --with-trilinos=0 --with-fftw=1 --with-fftw-dir=/home/sajid/packages/spack/opt/spack/linux-centos7-x86_64/gcc-8.3.0/fftw-3.3.8-mkj4ho2jp6xfrnkl
> mrvhdfh73woer2s7 --with-superlu_dist=0 --with-suitesparse=0 --with-zlib-include=/home/sajid/packages/spack/opt/spack/linux-centos7-x86_64/gcc-8.3.0
> /zlib-1.2.11-jqrcjdjnrxvouufhjtxbfvfms23fsqpx/include --with-zlib-lib="-L/home/sajid/packages/spack/opt/spack/linux-centos7-x86_64/gcc-8.3.0/zlib-1
> .2.11-jqrcjdjnrxvouufhjtxbfvfms23fsqpx/lib -lz" --with-zlib=1                                                                                      
> [0]PETSC ERROR: #1 MatMult() line 2385 in /tmp/sajid/spack-stage/spack-stage-mTD0AV/petsc-3.11.1/src/mat/interface/matrix.c                        
> [0]PETSC ERROR: #2 main() line 94 in /home/sajid/packages/aux_xwp_petsc/bugfix/ex112.c                                                             
> [0]PETSC ERROR: No PETSc Option Table entries                                                                                                      
> [0]PETSC ERROR: ----------------End of Error Message -------send entire error message to petsc-maint at mcs.anl.gov----------                         
> application called MPI_Abort(MPI_COMM_WORLD, 60) - process 0                                                                                       
> [unset]: write_line error; fd=-1 buf=:cmd=abort exitcode=60                                                                                        
> :                                                                                                                                                  
> system msg for write_line failure : Bad file descriptor                                                                                            
> 
> ```
> 
> I came across this because I saw that MatMult was failing for a new test related to a PR I was working on. Is this a bug ? 
> 
> Thank You,
> Sajid Ali
> Applied Physics
> Northwestern University