On Mon, Oct 29, 2012 at 11:29 AM, Florian Rathgeber <span dir="ltr"><<a href="mailto:florian.rathgeber@gmail.com" target="_blank">florian.rathgeber@gmail.com</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">> <a href="http://www.mcs.anl.gov/petsc/documentation/changes/dev.html" target="_blank">http://www.mcs.anl.gov/petsc/documentation/changes/dev.html</a><br>
<br>
</div>This change hasn't landed in the public repository yet, has it?<br></blockquote><div><br></div><div>Yes, there is no MatGetArray() in petsc-dev. It's MatSeqAIJGetArray(), etc.</div><div> </div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<div class="im">
<br>
</div>I'm using petsc-dev. In light of Matt's comment, we could initially<br>
support SeqAIJ matrices and when the interface change lands both<br>
SeqAIJ and dense. We could call the respective implementation<br>
depending on the matrix types and die for non-supported types. I'm<br>
happy to work on an initial patch.<br>
<br>
I'll also try to come up with a MatCreateAnyAIJWithArrays in the<br>
spirit of the existing MatCreateAnyAIJ that petsc4py supports.<br></blockquote><div><br></div><div>The upstream MatCreateMPIAIJ() was renamed MatCreateAIJ() because that's what it always did.</div><div><br></div><div>
You should be aware that MatCreateMPIAIJWithArrays() just does a loop over rows calling MatSetValues(). This copy is required because the internal format is different from the arrays. It's really much better to call MatSetValues() in your assembly so that there is no intermediate storage. The arrays that get passed to MatCreateMPIAIJWithArrays() are useless for computing with directly.</div>
<div>  </div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
Lisandro, do you prefer patches to this list, to your personal email<br>
or do you have some other platform for code review?</blockquote></div><br>