Confusing error message about index in matrix.c : PETSC ERROR: Column block index 5 (index 730) great than column length 750
Hi,
Using Petsc 3.19.6, I got an error that confused me:
[1]PETSC ERROR: Column block index 5 (index 730) great than column length 750
telling "730 is greater than 750"...
Here is the full log:
[1]PETSC ERROR: --------------------- Error Message --------------------------------------------------------------
[1]PETSC ERROR: Argument out of range
[1]PETSC ERROR: Column block index 5 (index 730) great than column length 750
[1]PETSC ERROR: WARNING! There are option(s) set that were not used! Could be the program crashed before they were used or a spelling mistake, etc!
[1]PETSC ERROR: Option left: name:-bddcericksp_atol value: 1e-20 source: code
[1]PETSC ERROR: Option left: name:-bddcericksp_divtol value: 1e30 source: code
[1]PETSC ERROR: Option left: name:-bddcericksp_max_it value: 4000 source: code
[1]PETSC ERROR: Option left: name:-bddcericksp_monitor (no value) source: code
[1]PETSC ERROR: Option left: name:-bddcericksp_rtol value: 1e-12 source: code
[1]PETSC ERROR: Option left: name:-bddcericksp_view (no value) source: code
[1]PETSC ERROR: Option left: name:-gmres_hpddmksp_atol value: 1e-20 source: code
[1]PETSC ERROR: Option left: name:-gmres_hpddmksp_divtol value: 1e30 source: code
[1]PETSC ERROR: Option left: name:-gmres_hpddmksp_gmres_modifiedgramschmidt (no value) source: code
[1]PETSC ERROR: Option left: name:-gmres_hpddmksp_max_it value: 4000 source: code
[1]PETSC ERROR: Option left: name:-gmres_hpddmksp_monitor (no value) source: code
[1]PETSC ERROR: Option left: name:-gmres_hpddmksp_monitor_singular_value (no value) source: code
[1]PETSC ERROR: Option left: name:-gmres_hpddmksp_rtol value: 1e-9 source: code
[1]PETSC ERROR: Option left: name:-gmres_hpddmksp_view (no value) source: code
[1]PETSC ERROR: Option left: name:-gmres_hpddmmat_block_size value: 3 source: code
[1]PETSC ERROR: Option left: name:-gmres_hpddmmat_type value: is source: code
[1]PETSC ERROR: Option left: name:-gmres_hpddmpc_hpddm_coarse_correction value: balanced source: code
[1]PETSC ERROR: Option left: name:-gmres_hpddmpc_hpddm_coarse_p value: 2 source: code
[1]PETSC ERROR: Option left: name:-gmres_hpddmpc_hpddm_coarse_sub_pc_type value: lu source: code
[1]PETSC ERROR: Option left: name:-gmres_hpddmpc_hpddm_define_subdomains (no value) source: code
[1]PETSC ERROR: Option left: name:-gmres_hpddmpc_hpddm_has_neumann (no value) source: code
[1]PETSC ERROR: Option left: name:-gmres_hpddmpc_hpddm_levels_1_eps_gen_non_hermitian (no value) source: code
[1]PETSC ERROR: Option left: name:-gmres_hpddmpc_hpddm_levels_1_st_share_sub_ksp (no value) source: code
[1]PETSC ERROR: Option left: name:-gmres_hpddmpc_hpddm_levels_1_sub_pc_type value: lu source: code
[1]PETSC ERROR: Option left: name:-gmres_hpddmpc_hpddm_schur_precondition value: geneo source: code
[1]PETSC ERROR: Option left: name:-gmres_hpddmpc_type value: hpddm source: code
[1]PETSC ERROR: Option left: name:-grosLUksp_atol value: 1e-20 source: code
[1]PETSC ERROR: Option left: name:-grosLUksp_divtol value: 1e30 source: code
[1]PETSC ERROR: Option left: name:-grosLUksp_max_it value: 4000 source: code
[1]PETSC ERROR: Option left: name:-grosLUksp_monitor (no value) source: code
[1]PETSC ERROR: Option left: name:-grosLUksp_monitor_singular_value (no value) source: code
[1]PETSC ERROR: Option left: name:-grosLUksp_rtol value: 1e-9 source: code
[1]PETSC ERROR: Option left: name:-grosLUksp_view value: ::ascii_info_detail source: code
[1]PETSC ERROR: Option left: name:-grosLUpc_type value: lu source: code
[1]PETSC ERROR: See https://petsc.org/release/faq/ for trouble shooting.
[1]PETSC ERROR: Petsc Release Version 3.19.6, Sep 28, 2023
[1]PETSC ERROR: MEF++.dev on a named rohan by ericc Wed Mar 27 16:24:20 2024
[1]PETSC ERROR: Configure options --prefix=/opt/petsc-3.19.6-debug_openmpi-4.1.6 --with-mpi-compilers=1 --with-mpi-dir=/opt/openmpi-4.1.6 --with-make-np=12 --with-shared-libraries=1 --with-debugging=yes --with-memalign=64 --with-visibility=0 --with-64-bit-indices=0 --download-chaco=yes --download-ml=yes --download-party=yes --download-spai=yes --download-superlu=yes --download-cmake=yes --download-hpddm=yes --download-hypre=yes --download-metis=yes --download-mumps=yes --download-parmetis=yes --download-ptscotch=yes --download-slepc=yes --download-strumpack=yes --download-suitesparse=yes --download-superlu_dist=yes --with-blaslapack-dir=/opt/intel/oneapi/mkl/2024.0/lib/intel64 --with-mkl_pardiso-dir=/opt/intel/oneapi/mkl/2024.0 --with-scalapack=1 --with-scalapack-include=/opt/intel/oneapi/mkl/2024.0/include --with-mkl_cpardiso-dir=/opt/intel/oneapi/mkl/2024.0 --with-scalapack-lib="-L/opt/intel/oneapi/mkl/2024.0/lib/intel64 -lmkl_scalapack_lp64 -lmkl_blacs_openmpi_lp64"
[1]PETSC ERROR: #1 MatSetValuesBlocked() at /home/mefpp_ericc/petsc/ompi-opt/petsc-3.19.6-debug/src/mat/interface/matrix.c:2005
[1]PETSC ERROR: #2 MatSetValuesBlocked_IS() at /home/mefpp_ericc/petsc/ompi-opt/petsc-3.19.6-debug/src/mat/impls/is/matis.c:2589
[1]PETSC ERROR: #3 MatSetValuesBlocked() at /home/mefpp_ericc/petsc/ompi-opt/petsc-3.19.6-debug/src/mat/interface/matrix.c:2013
Maybe the message is fixed in recent PETSc?
Anyway, I know it is all about block sizes, but fixing the message should not be a big challenge...
Thanks,
Eric