On Mon, Oct 29, 2012 at 11:15 AM, Lisandro Dalcin <span dir="ltr"><<a href="mailto:dalcinl@gmail.com" target="_blank">dalcinl@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">On 29 October 2012 09:46, Florian Rathgeber <<a href="mailto:florian.rathgeber@gmail.com">florian.rathgeber@gmail.com</a>> wrote:<br>
> -----BEGIN PGP SIGNED MESSAGE-----<br>
> Hash: SHA1<br>
><br>
> [Apologies if this should have gone to petsc-users instead]<br>
><br>
<br>
</div>I think this list is perfectly fine for the subject of your request.<br>
<div class="im"><br>
> I'd like to extend petsc4py's Mat extension type to expose<br>
> MatCreateSeqAIJWithArrays and MatGetArray/MatRestoreArray to python,<br>
> ideally via attributes/methods directly on the extension type.<br>
><br>
<br>
</div>A minor issue with MatGetArray/MatRestoreArray is to properly support<br>
AIJ and DENSE matrices.</blockquote><div><br></div><div>Don't do this. We have changed the interface so that it is specific to the</div><div>subtypes:</div><div><br></div><div>  <a href="http://www.mcs.anl.gov/petsc/documentation/changes/dev.html">http://www.mcs.anl.gov/petsc/documentation/changes/dev.html</a></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 class="im">
> I know these shouldn't normally be used and reading this thread [1]<br>
> I'm assuming Lisandro doesn't currently have any plans for adding<br>
> either of those to petsc4py?<br>
><br>
<br>
</div>I would love to add them, the only issue is that I'm not sure about<br>
how to support all the various sparse matrix formats.<br>
<br>
We could start by only supporting AIJ matrices. If you can help with<br>
an initial patch, I can continue to refine it.<br>
<br>
Are you using petsc-3.3? We can add the required support in branch<br>
release-3.3 of petsc4py, then make a new patch release.<br>
<span class="HOEnZb"><font color="#888888"><br>
<br>
--<br>
Lisandro Dalcin<br>
---------------<br>
CIMEC (INTEC/CONICET-UNL)<br>
Predio CONICET-Santa Fe<br>
Colectora RN 168 Km 472, Paraje El Pozo<br>
3000 Santa Fe, Argentina<br>
Tel: <a href="tel:%2B54-342-4511594" value="+543424511594">+54-342-4511594</a> (ext 1011)<br>
Tel/Fax: <a href="tel:%2B54-342-4511169" value="+543424511169">+54-342-4511169</a><br>
</font></span></blockquote></div><br><br clear="all"><div><br></div>-- <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>