[petsc-dev] [petsc-users] MatSetBlockSize no longer works after MatCreateSeqAIJ
Matthew Knepley
knepley at gmail.com
Thu Jun 14 06:17:50 CDT 2012
On Thu, Jun 14, 2012 at 11:16 AM, Barry Smith <bsmith at mcs.anl.gov> wrote:
>
> On Jun 13, 2012, at 10:12 PM, Jed Brown wrote:
>
> > Okay, I thought switching to MatCreateAIJ was a good time to add it, bit
> it's kind of too late now.
>
> Good point. Do it, call it a bug when we didn't add the bs on that
> change so this is worthy of a patch for 3.3
Please no. This will break a TON of code, vs. the 1-3 messages we will get
now.
Matt
>
> Barry
>
> >
> > On Jun 13, 2012 10:06 PM, "Barry Smith" <bsmith at mcs.anl.gov> wrote:
> >
> > On Jun 13, 2012, at 8:03 PM, Jed Brown wrote:
> >
> > > What I do is
> > >
> > > MatCreate(comm,&A);
> > > MatSetSizes(A,...);
> > > MatSetBlockSize(A,bs);
> > > MatSetOptionsPrefix(A,"yourchoice_");
> > > MatSetFromOptions(A);
> > > MatXAIJSetPreallocation(A,...);
> > >
> > > This should work with any matrix type, parallel and sequential, and
> you can specify it without interfering with any other matrices using
> -yourchoice_mat_type aij {baij, sbaij}.
> > >
> > >
> > > Devs, is there a reason that block size wasn't added to the
> convenience constructor?
> >
> > To break less stuff and not make most people change their
> MatCreateMPIAIJ call()
> >
> > > Wasn't there a discussion about adding it?
> >
> > Yes
> >
> > | What should we do?
> >
> > Get rid of convenience constructors?
> >
> > Barry
> >
> >
> > >
> > > On Wed, Jun 13, 2012 at 7:55 PM, Geoffrey Irving <irving at naml.us>
> wrote:
> > > 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
> > >
> >
>
>
--
What most experimenters take for granted before they begin their
experiments is infinitely more interesting than any results to which their
experiments lead.
-- Norbert Wiener
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-dev/attachments/20120614/02b688ea/attachment.html>
More information about the petsc-dev
mailing list