<p dir="ltr">Cool. Thank you all for your valuable suggestions.</p>
<div class="gmail_quote">On May 11, 2015 2:10 PM, "Jed Brown" <<a href="mailto:jed@jedbrown.org">jed@jedbrown.org</a>> wrote:<br type="attribution"><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Bikash Kanungo <<a href="mailto:bikash@umich.edu">bikash@umich.edu</a>> writes:<br>
<br>
> So I guess, in terms of performance creating M^{-1} through<br>
> MatGetLocalSubMatrix and then explicitly build M^{-1}H is a more efficient<br>
> way.<br>
<br>
Yes, unless this destroys sparsity that would otherwise be exploitable.<br>
</blockquote></div>