<div dir="ltr"><div>The DM in snes/ex56 does elasticity but the block size does not get into vectors. It does something like this:</div><div><br></div><div>PetscFECreateDefault(dm, dim, dim, PETSC_FALSE, NULL, 3, &fe);CHKERRQ(ierr); /* elasticity */ </div><div>DMGetDS(dm, &prob);CHKERRQ(ierr);</div><div>PetscDSSetDiscretization(prob, 0, (PetscObject) fe);CHKERRQ(ierr);</div><div>DMCreateGlobalVector(dm,&xx);<br></div><div><br></div><div>The block size of xx is not 3. (dim==3). It is not clear to me how to make this work.</div><div><br></div><div>snes/ex56 does use -mat_block_size 3.</div><div><br></div><div>Mark</div></div>