Yes I am doing: call MatMPIBAIJSetPreallocation(this%A, 4-bdim, flubioSolvers%d_nz, mesh%d_nnz, flubioSolvers%o_nz, mesh%o_nnz, ierr)
with d_nnz the number of diagonal blocks and o_nnz the number of off-diagonal blocks. However I am getting this: [0]PETSC ERROR: --------------------- Error Message -------------------------------------------------------------- [0]PETSC ERROR: Argument out of range [0]PETSC ERROR: New nonzero at (3,3) caused a malloc Use MatSetOption(A, MAT_NEW_NONZERO_ALLOCATION_ERR, PETSC_FALSE) to turn off this check [0]PETSC ERROR: See https://petsc.org/release/faq/ for trouble shooting. [0]PETSC ERROR: Petsc Release Version 3.18.0, Sep 30, 2022 [0]PETSC ERROR: flubio_coupled on a gnu named alienware by edo Thu Nov 3 18:19:30 2022 [0]PETSC ERROR: Configure options PETSC_ARCH=gnu FOPTFLAGS=-O3 COPTFLAGS=-O3 CXXOPTFLAGS=-O3 -with-debugging=no -download-fblaslapack=1 -download-superlu_dist -download-mumps -download-hypre -download-metis -download-parmetis -download-scalapack -download-ml -download-slepc -download-hpddm -download-cmake -with-mpi-dir=/home/edo/software/openmpi-4.1.1/build/ [0]PETSC ERROR: #1 MatSetValuesBlocked_SeqBAIJ_Inlined() at /home/edo/software/petsc-3.18.0/src/mat/impls/baij/mpi/mpibaij.c:318 [0]PETSC ERROR: #2 MatSetValuesBlocked_MPIBAIJ() at /home/edo/software/petsc-3.18.0/src/mat/impls/baij/mpi/mpibaij.c:389 [1]PETSC ERROR: --------------------- Error Message -------------------------------------------------------------- [1]PETSC ERROR: Argument out of range [1]PETSC ERROR: New nonzero at (0,0) caused a malloc Use MatSetOption(A, MAT_NEW_NONZERO_ALLOCATION_ERR, PETSC_FALSE) to turn off this check [1]PETSC ERROR: See https://petsc.org/release/faq/ for trouble shooting. [1]PETSC ERROR: Petsc Release Version 3.18.0, Sep 30, 2022 [1]PETSC ERROR: flubio_coupled on a gnu named alienware by edo Thu Nov 3 18:19:30 2022 [1]PETSC ERROR: Configure options PETSC_ARCH=gnu FOPTFLAGS=-O3 COPTFLAGS=-O3 CXXOPTFLAGS=-O3 -with-debugging=no -download-fblaslapack=1 -download-superlu_dist -download-mumps -download-hypre -download-metis -download-parmetis -download-scalapack -download-ml -download-slepc -download-hpddm -download-cmake -with-mpi-dir=/home/edo/software/openmpi-4.1.1/build/ [1]PETSC ERROR: #1 MatSetValuesBlocked_SeqBAIJ_Inlined() at /home/edo/software/petsc-3.18.0/src/mat/impls/baij/mpi/mpibaij.c:318 [1]PETSC ERROR: #2 MatSetValuesBlocked_MPIBAIJ() at /home/edo/software/petsc-3.18.0/src/mat/impls/baij/mpi/mpibaij.c:419 [1]PETSC ERROR: #3 MatSetValuesBlocked() at /home/edo/software/petsc-3.18.0/src/mat/interface/matrix.c:1978 [1]PETSC ERROR: #4 MatSetValuesBlocked_SeqBAIJ_Inlined() at /home/edo/software/petsc-3.18.0/src/mat/impls/baij/mpi/mpibaij.c:318 [1]PETSC ERROR: #5 MatSetValuesBlocked_MPIBAIJ() at /home/edo/software/petsc-3.18.0/src/mat/impls/baij/mpi/mpibaij.c:419 [1]PETSC ERROR: #6 MatSetValuesBlocked() at /home/edo/software/petsc-3.18.0/src/mat/interface/matrix.c:1978 [0]PETSC ERROR: #3 MatSetValuesBlocked() at /home/edo/software/petsc-3.18.0/src/mat/interface/matrix.c:1978 [0]PETSC ERROR: #4 MatSetValuesBlocked_SeqBAIJ_Inlined() at /home/edo/software/petsc-3.18.0/src/mat/impls/baij/mpi/mpibaij.c:318 [0]PETSC ERROR: #5 MatSetValuesBlocked_MPIBAIJ() at /home/edo/software/petsc-3.18.0/src/mat/impls/baij/mpi/mpibaij.c:419 [0]PETSC ERROR: #6 MatAssemblyEnd_MPIBAIJ() at /home/edo/software/petsc-3.18.0/src/mat/impls/baij/mpi/mpibaij.c:906 [0]PETSC ERROR: #7 MatAssemblyEnd() at /home/edo/software/petsc-3.18.0/src/mat/interface/matrix.c:5696
