[petsc-users] MatSetBlockSize no longer works after MatCreateSeqAIJ

[Geoffrey Irving]

Hello,

The following code used to work, but breaks after upgrading to 3.3:

    ....
    CHECK(MatCreateSeqAIJ(comm,3*n,3*n,0,scalar_view(nnz).data(),&mat));
    CHECK(MatSetBlockSize(mat,3));
    ...

The error follows.  Presumably MatSetBlockSize got more strict, and I
need to break the single MatCreateSeqAIJ into multiple lines (matrix
creation followed by preallocation) in order to insert MatSetBlockSize
in the middle.  What are the correct functions to call?

Thanks,
Geoffrey

[0]PETSC ERROR: --------------------- Error Message
------------------------------------
[0]PETSC ERROR: Arguments are incompatible!
[0]PETSC ERROR: Cannot change block size 1 to 3!
[0]PETSC ERROR:
------------------------------------------------------------------------
[0]PETSC ERROR: Petsc Release Version 3.3.0, Patch 0, Tue Jun  5
14:20:42 CDT 2012
[0]PETSC ERROR: See docs/changes/index.html for recent updates.
[0]PETSC ERROR: See docs/faq.html for hints about trouble shooting.
[0]PETSC ERROR: See docs/index.html for manual pages.
[0]PETSC ERROR:
------------------------------------------------------------------------
[0]PETSC ERROR: Unknown Name on a darwin named tile.local by irving
Wed Jun 13 17:51:34 2012
[0]PETSC ERROR: Libraries linked from
/opt/local/var/macports/build/_opt_local_var_macports_sources_rsync.macports.org_release_ports_math_petsc/petsc/work/destroot/opt/local/lib/petsc/lib
[0]PETSC ERROR: Configure run at Wed Jun 13 12:11:05 2012
[0]PETSC ERROR: Configure options --prefix=/opt/local --with-python
--with-debugging=0 --with-c-support=1 --with-c++-support=1
--with-pic=fPIC --with-shared-libraries=0 --with-mpi=1
--PETSC_ARCH=darwin
--prefix=/opt/local/var/macports/build/_opt_local_var_macports_sources_rsync.macports.org_release_ports_math_petsc/petsc/work/destroot/opt/local/lib/petsc
--with-cc=/opt/local/bin/openmpicc
--with-cxx=/opt/local/bin/openmpicxx
--with-mpiexec=/opt/local/bin/openmpiexec
--with-lapack-lib=/usr/lib/liblapack.dylib
--with-blas-lib=/usr/lib/libblas.dylib --download-umfpack=1
--download-ml --download-metis --download-parmetis --download-blacs
--download-scalapack --download-mumps --download-ptscotch
--with-fc=/opt/local/bin/openmpif90 --LIBS=-lstdc++
[0]PETSC ERROR:
------------------------------------------------------------------------
[0]PETSC ERROR: PetscLayoutSetBlockSize() line 460 in
/opt/local/var/macports/build/_opt_local_var_macports_sources_rsync.macports.org_release_ports_math_petsc/petsc/work/petsc-3.3-p0/src/vec/vec/impls/mpi/pmap.c
[0]PETSC ERROR: MatSetBlockSize() line 6691 in
/opt/local/var/macports/build/_opt_local_var_macports_sources_rsync.macports.org_release_ports_math_petsc/petsc/work/petsc-3.3-p0/src/mat/interface/matrix.c
[0]PETSC ERROR: --------------------- Error Message
------------------------------------
[0]PETSC ERROR: Arguments are incompatible!
[0]PETSC ERROR:  !
[0]PETSC ERROR:
------------------------------------------------------------------------
[0]PETSC ERROR: Petsc Release Version 3.3.0, Patch 0, Tue Jun  5
14:20:42 CDT 2012
[0]PETSC ERROR: See docs/changes/index.html for recent updates.
[0]PETSC ERROR: See docs/faq.html for hints about trouble shooting.
[0]PETSC ERROR: See docs/index.html for manual pages.
[0]PETSC ERROR:
------------------------------------------------------------------------
[0]PETSC ERROR: Unknown Name on a darwin named tile.local by irving
Wed Jun 13 17:51:34 2012
[0]PETSC ERROR: Libraries linked from
/opt/local/var/macports/build/_opt_local_var_macports_sources_rsync.macports.org_release_ports_math_petsc/petsc/work/destroot/opt/local/lib/petsc/lib
[0]PETSC ERROR: Configure run at Wed Jun 13 12:11:05 2012
[0]PETSC ERROR: Configure options --prefix=/opt/local --with-python
--with-debugging=0 --with-c-support=1 --with-c++-support=1
--with-pic=fPIC --with-shared-libraries=0 --with-mpi=1
--PETSC_ARCH=darwin
--prefix=/opt/local/var/macports/build/_opt_local_var_macports_sources_rsync.macports.org_release_ports_math_petsc/petsc/work/destroot/opt/local/lib/petsc
--with-cc=/opt/local/bin/openmpicc
--with-cxx=/opt/local/bin/openmpicxx
--with-mpiexec=/opt/local/bin/openmpiexec
--with-lapack-lib=/usr/lib/liblapack.dylib
--with-blas-lib=/usr/lib/libblas.dylib --download-umfpack=1
--download-ml --download-metis --download-parmetis --download-blacs
--download-scalapack --download-mumps --download-ptscotch
--with-fc=/opt/local/bin/openmpif90 --LIBS=-lstdc++
[0]PETSC ERROR:
------------------------------------------------------------------------
[0]PETSC ERROR: PetscSolidMat() line 114 in petsc/ksp.cpp