[petsc-users] On the usage of MatSetValuesBlocked

[Edoardo alinovi]

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
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20221103/48997d81/attachment-0001.html>