<div dir="ltr">Thanks, Barry. <div><br></div><div>The reason I need to access these structures is that I am trying to interface PETSc with MKL through the PCSHELL interface. I need the access to the basic CSR structure in order to invoke the ILUT subroutine in the MKL package. <div>
<br></div><div style>BTW, should it be MatGetRowIJ() instead of MatGetIJ() ?</div><div style><br></div><div style>Thanks,</div><div style>Su</div><div style><br></div></div></div><div class="gmail_extra"><br><br><div class="gmail_quote">
On Tue, Jul 23, 2013 at 9:03 PM, Barry Smith <span dir="ltr"><<a href="mailto:bsmith@mcs.anl.gov" target="_blank">bsmith@mcs.anl.gov</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<div class="im"><br>
On Jul 23, 2013, at 8:59 PM, Su Yan <<a href="mailto:suyan0@gmail.com">suyan0@gmail.com</a>> wrote:<br>
<br>
> Hi, is there any way that I can extract the index information and element array from a AIJ matrix? Say, I have a Mat A, which is created using the CSR format with MatCreateSeqAIJ(). How can I get int *ia, int *ja, and double *array out of A?<br>

<br>
</div>   MatSeqAIJGetArray() for array and MatGetIJ() for ia and ja but we hope that generally users will never need to get access to these low-level underlying structures. Perhaps if you tell us why you want them we may suggest alternatives.<br>

<br>
    Barry<br>
<br>
><br>
> Thanks,<br>
> Su<br>
<br>
</blockquote></div><br></div>