<div dir="ltr"><div class="gmail_extra"><div class="gmail_quote">On Wed, Jun 24, 2015 at 2:39 AM, domenico lahaye <span dir="ltr"><<a href="mailto:domenico_lahaye@yahoo.com" target="_blank">domenico_lahaye@yahoo.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div><div style="color:#000;background-color:#fff;font-family:times new roman,new york,times,serif;font-size:16px"><div>Would the same work for non-structured problems, for instance for power flow equations? </div></div></div></blockquote><div><br></div><div>Yes, PetscSection knows about components for fields, and this carries through to the DMNetwork setup.</div><div>It will, for instance, use block matrix formats if the block size is constant.</div><div><br></div><div>  Thanks,</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><div style="color:#000;background-color:#fff;font-family:times new roman,new york,times,serif;font-size:16px"><div>Thx. Domenico Lahaye. <br></div><div><span></span></div><br>  <div style="font-family:times new roman,new york,times,serif;font-size:16px"> <div style="font-family:HelveticaNeue,Helvetica Neue,Helvetica,Arial,Lucida Grande,sans-serif;font-size:16px"> <div dir="ltr"> <hr size="1">  <font face="Arial" size="2"> <b><span style="font-weight:bold">From:</span></b> Matthew Knepley <<a href="mailto:knepley@gmail.com" target="_blank">knepley@gmail.com</a>><br> <b><span style="font-weight:bold">To:</span></b> Gideon Simpson <<a href="mailto:gideon.simpson@gmail.com" target="_blank">gideon.simpson@gmail.com</a>> <br><b><span style="font-weight:bold">Cc:</span></b> petsc-users <<a href="mailto:petsc-users@mcs.anl.gov" target="_blank">petsc-users@mcs.anl.gov</a>> <br> <b><span style="font-weight:bold">Sent:</span></b> Wednesday, June 24, 2015 5:05 AM<br> <b><span style="font-weight:bold">Subject:</span></b> Re: [petsc-users] jacobians for multicomponent problems<br> </font> </div><span class=""> <div><br><div><div><div dir="ltr"><div><div>On Tue, Jun 23, 2015 at 7:52 PM, Gideon Simpson <span dir="ltr"><<a rel="nofollow" shape="rect" href="mailto:gideon.simpson@gmail.com" target="_blank">gideon.simpson@gmail.com</a>></span> wrote:<br clear="none"><blockquote style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div style="word-wrap:break-word">Suppose I have a problem for which the field variable has multiple components, in this case real and imaginary parts, and I have created a DM with two degrees of freedom.  When assembling a Jacobian associated with a nonlinear problem on this data, is there any data management accounting for the different components that makes this “easier?”</div></blockquote><div><br clear="none"></div><div>DMDA knows about components so that MatStencil has a .c slot for indices. Is that what you mean?</div><div><br clear="none"></div><div>  Thanks,</div><div><br clear="none"></div><div>    Matt</div><div><br><br></div><div><div> </div></div><blockquote style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div style="word-wrap:break-word"><div><span><font color="#888888"></font></span></div><div><div><div>
<span style="border-collapse:separate;color:rgb(0,0,0);font-family:Helvetica;font-style:normal;font-variant:normal;font-weight:normal;letter-spacing:normal;line-height:normal;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px">-gideon</span></div>

</div>
<br clear="none"></div></div></blockquote></div><br clear="none"><br clear="all"><div><br clear="none"></div>-- <br clear="none"><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 clear="none">-- Norbert Wiener</div><div>
</div></div></div></div></div><br><br></div> </span></div> </div>  </div></div></blockquote></div><br><br clear="all"><div><br></div>-- <br><div class="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>