<div dir="ltr">Yes I am doing: <div><br></div><div>call MatMPIBAIJSetPreallocation(this%A, 4-bdim, flubioSolvers%d_nz, mesh%d_nnz, flubioSolvers%o_nz, mesh%o_nnz, ierr)<br></div><div><br></div><div>with d_nnz the number of diagonal blocks and o_nnz the number of off-diagonal blocks. However I am getting this: </div><div><br></div><div>[0]PETSC ERROR: --------------------- Error Message --------------------------------------------------------------<br>[0]PETSC ERROR: Argument out of range<br>[0]PETSC ERROR: New nonzero at (3,3) caused a malloc<br>Use MatSetOption(A, MAT_NEW_NONZERO_ALLOCATION_ERR, PETSC_FALSE) to turn off this check<br>[0]PETSC ERROR: See <a href="https://petsc.org/release/faq/">https://petsc.org/release/faq/</a> for trouble shooting.<br>[0]PETSC ERROR: Petsc Release Version 3.18.0, Sep 30, 2022 <br>[0]PETSC ERROR: flubio_coupled on a gnu named alienware by edo Thu Nov  3 18:19:30 2022<br>[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/<br>[0]PETSC ERROR: #1 MatSetValuesBlocked_SeqBAIJ_Inlined() at /home/edo/software/petsc-3.18.0/src/mat/impls/baij/mpi/mpibaij.c:318<br>[0]PETSC ERROR: #2 MatSetValuesBlocked_MPIBAIJ() at /home/edo/software/petsc-3.18.0/src/mat/impls/baij/mpi/mpibaij.c:389<br>[1]PETSC ERROR: --------------------- Error Message --------------------------------------------------------------<br>[1]PETSC ERROR: Argument out of range<br>[1]PETSC ERROR: New nonzero at (0,0) caused a malloc<br>Use MatSetOption(A, MAT_NEW_NONZERO_ALLOCATION_ERR, PETSC_FALSE) to turn off this check<br>[1]PETSC ERROR: See <a href="https://petsc.org/release/faq/">https://petsc.org/release/faq/</a> for trouble shooting.<br>[1]PETSC ERROR: Petsc Release Version 3.18.0, Sep 30, 2022 <br>[1]PETSC ERROR: flubio_coupled on a gnu named alienware by edo Thu Nov  3 18:19:30 2022<br>[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/<br>[1]PETSC ERROR: #1 MatSetValuesBlocked_SeqBAIJ_Inlined() at /home/edo/software/petsc-3.18.0/src/mat/impls/baij/mpi/mpibaij.c:318<br>[1]PETSC ERROR: #2 MatSetValuesBlocked_MPIBAIJ() at /home/edo/software/petsc-3.18.0/src/mat/impls/baij/mpi/mpibaij.c:419<br>[1]PETSC ERROR: #3 MatSetValuesBlocked() at /home/edo/software/petsc-3.18.0/src/mat/interface/matrix.c:1978<br>[1]PETSC ERROR: #4 MatSetValuesBlocked_SeqBAIJ_Inlined() at /home/edo/software/petsc-3.18.0/src/mat/impls/baij/mpi/mpibaij.c:318<br>[1]PETSC ERROR: #5 MatSetValuesBlocked_MPIBAIJ() at /home/edo/software/petsc-3.18.0/src/mat/impls/baij/mpi/mpibaij.c:419<br>[1]PETSC ERROR: #6 MatSetValuesBlocked() at /home/edo/software/petsc-3.18.0/src/mat/interface/matrix.c:1978<br>[0]PETSC ERROR: #3 MatSetValuesBlocked() at /home/edo/software/petsc-3.18.0/src/mat/interface/matrix.c:1978<br>[0]PETSC ERROR: #4 MatSetValuesBlocked_SeqBAIJ_Inlined() at /home/edo/software/petsc-3.18.0/src/mat/impls/baij/mpi/mpibaij.c:318<br>[0]PETSC ERROR: #5 MatSetValuesBlocked_MPIBAIJ() at /home/edo/software/petsc-3.18.0/src/mat/impls/baij/mpi/mpibaij.c:419<br>[0]PETSC ERROR: #6 MatAssemblyEnd_MPIBAIJ() at /home/edo/software/petsc-3.18.0/src/mat/impls/baij/mpi/mpibaij.c:906<br>[0]PETSC ERROR: #7 MatAssemblyEnd() at /home/edo/software/petsc-3.18.0/src/mat/interface/matrix.c:5696</div></div>