<div dir="ltr"><div>On Thu, Sep 3, 2015 at 8:11 PM, Barry Smith <span dir="ltr"><<a href="mailto:bsmith@mcs.anl.gov" target="_blank">bsmith@mcs.anl.gov</a>></span> wrote:<br></div><div class="gmail_extra"><div class="gmail_quote"><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left-width:1px;border-left-color:rgb(204,204,204);border-left-style:solid;padding-left:1ex"><br>
  You might also consider the MatNest object. It's job is to help construct matrices that have a natural set of sub matrices such as yours. Also look at MatGetLocalSubMat()<br></blockquote><div><br></div><div><br>OK, thanks for the pointers.<div><br></div><div>Best regards,</div><div>David</div><div><br></div><div></div></div><div><br></div><div> </div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left-width:1px;border-left-color:rgb(204,204,204);border-left-style:solid;padding-left:1ex"><div><div>
<br>
> On Sep 3, 2015, at 6:52 PM, David Knezevic <<a href="mailto:david.knezevic@akselos.com" target="_blank">david.knezevic@akselos.com</a>> wrote:<br>
><br>
> I need to add Lagrange multipliers to an existing system of equations, which requires me to expand the matrix.<br>
><br>
> I figured I would allocate a new larger matrix with the appropriate non-zero structure, and then copy the entries from the smaller matrix over to it. But I was wondering if there might be a better way to do this?<br>
><br>
> Thanks,<br>
> David<br>
><br>
><br>
<br>
</div></div></blockquote></div><br></div></div>