On Tue, Mar 15, 2011 at 4:10 PM, Jed Brown <span dir="ltr"><<a href="mailto:jed@59a2.org">jed@59a2.org</a>></span> wrote:<br><div class="gmail_quote"><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex;">
<div class="im"><div class="gmail_quote">On Tue, Mar 15, 2011 at 22:07, 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:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">

<div>I guess I have no problem with MatGetSubMatrix() being the interface for this, since all solvers get the Mat object. However, my point was</div><div>really that the interface for creating such things is horrible. There is no easy way for me to do what I want, and the hard way is too hard for</div>


<div>me (a developer) without looking at a lot of code, trying things and seeing where they break, and probably some crying. This will have to</div><div>change or no one will use this.</div></blockquote></div><br></div><div>
I'd love to hear suggestions for something with similar or better modularity.</div>
</blockquote></div><br>What I have now is a MatShell that implements MatGetSubMatrix(), and it asks<div><br></div><div>  Is this the field I am replacing?</div><div><br></div><div>    Yes - give back my stored matrix</div>
<div><br></div><div>     No  - Dispatch to the original matrix</div><div><br></div><div>For this I need IS retrieval from FS which I am adding. I think you can do this in</div><div>general very easily. Basically you have a function SwapOutBlockOfMatrix().</div>
<div><br></div><div>   Matt<br clear="all"><br>-- <br>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<br>

</div>