<div dir="ltr"><br><div class="gmail_extra"><br><div class="gmail_quote">On Sun, Sep 4, 2016 at 8:35 PM, Matthew Knepley <span dir="ltr"><<a href="mailto:knepley@gmail.com" target="_blank">knepley@gmail.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr"><div class="gmail_extra"><div class="gmail_quote"><span class="gmail-">On Sun, Sep 4, 2016 at 7:07 PM, 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:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr"><div class="gmail_extra"><div class="gmail_quote">On Sat, Sep 3, 2016 at 11:05 PM, Matthew Knepley <span dir="ltr"><<a href="mailto:knepley@gmail.com" target="_blank">knepley@gmail.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr"><div class="gmail_extra"><div class="gmail_quote"><span>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:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);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></span><div>The DM determines the block size based on the Section. There could be a bug, </div></div></div></div></blockquote><div><br></div><div>I suspect that I just found a hole in the combinatorics (finite!) that is Plex. I use -mat_block_size 3, so I could see how Plex might not pick it up.</div></div></div></div></blockquote><div><br></div></span><div>Perhaps. Give me the branch and exact command line to run.</div></div></div></div></blockquote><div><br></div><div>My pull request for snes/ex56 has:</div><div><br></div><div>ierr = VecNorm(xx,NORM_INFINITY,&mdisp[iter]);CHKERRQ(ierr);<br></div><div><br></div><div>When I tried to use VecStrideNorm I got the error that block size was not set. You can't set it manually so I was stuck (not a big deal, but there is either a bug in Plex or a bug in the way that I am using it).</div><div><br></div><div>Thanks,</div><div> </div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr"><div class="gmail_extra"><div class="gmail_quote"><div><br></div><div>  Matt</div><span class="gmail-"><div> </div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr"><div class="gmail_extra"><div class="gmail_quote"><div>All I do is call DMCreateMatrix(dm, &Amat);</div><div><br></div><div>And Amat get its w/o and explicit setfromoptions.</div><div><br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr"><div class="gmail_extra"><div class="gmail_quote"><div>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></div></div></blockquote><div><br></div><div>No</div><div> </div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr"><div class="gmail_extra"><div class="gmail_quote"><div><br></div><div>  Matt</div><span><div> </div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr"><div>snes/ex56 does use -mat_block_size 3.</div><span><font color="#888888"><div><br></div><div>Mark</div></font></span></div></blockquote></span></div><span><font color="#888888"><span><font color="#888888">-- <br><div>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>
</font></span></font></span></div></div>
</blockquote></div><br></div></div>
</blockquote></span></div><span class="gmail-"><br><br clear="all"><div><br></div>-- <br><div>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>
</span></div></div>
</blockquote></div><br></div></div>