<div dir="ltr"><div class="gmail_extra"><div class="gmail_quote">On Sat, Sep 3, 2016 at 11:12 AM, Mark Adams <span dir="ltr"><<a href="mailto:mfadams@lbl.gov" target="_blank">mfadams@lbl.gov</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><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></blockquote><div><br></div><div>The DM determines the block size based on the Section. There could be a bug, but first did you apply</div><div>a BC that would wreck the block size? Namely did you constrain just one or two unknowns at a point?</div><div><br></div><div>  Matt</div><div> </div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><div>snes/ex56 does use -mat_block_size 3.</div><span class="HOEnZb"><font color="#888888"><div><br></div><div>Mark</div></font></span></div></blockquote></div>-- <br><div class="gmail_signature" data-smartmail="gmail_signature">What most experimenters take for granted before they begin their experiments is infinitely more interesting than any results to which their experiments lead.<br>-- Norbert Wiener</div>
</div></div>