[petsc-dev] How to get block size into a DM

[Mark Adams]

The DM in snes/ex56 does elasticity but the block size does not get into
vectors. It does something like this:

PetscFECreateDefault(dm, dim, dim, PETSC_FALSE, NULL, 3,
&fe);CHKERRQ(ierr); /* elasticity */
DMGetDS(dm, &prob);CHKERRQ(ierr);
PetscDSSetDiscretization(prob, 0, (PetscObject) fe);CHKERRQ(ierr);
DMCreateGlobalVector(dm,&xx);

The block size of xx is not 3. (dim==3). It is not clear to me how to make
this work.

snes/ex56 does use -mat_block_size 3.

Mark
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-dev/attachments/20160903/976dd48e/attachment.html>