<html><head></head><body style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space; ">Dear Matt and Jed,<div><br></div><div>Thanks very much! That helps!</div><div><br></div><div>Best regards,</div><div><br></div><div>Rebecca</div><div><br></div><div><br></div><div><div>
</div>
<br><div><div>On Jul 14, 2011, at 11:36 AM, Matthew Knepley wrote:</div><br class="Apple-interchange-newline"><blockquote type="cite">On Thu, Jul 14, 2011 at 6:30 PM, Xuefei (Rebecca) Yuan <span dir="ltr">&lt;<a href="mailto:xyuan@lbl.gov">xyuan@lbl.gov</a>&gt;</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 style="word-wrap:break-word"><div><span style="border-collapse:separate;color:rgb(0, 0, 0);font-family:Helvetica;font-style:normal;font-variant:normal;font-weight:normal;letter-spacing:normal;line-height:normal;text-align:auto;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px;font-size:medium"><span style="border-collapse:separate;color:rgb(0, 0, 0);font-family:Helvetica;font-style:normal;font-variant:normal;font-weight:normal;letter-spacing:normal;line-height:normal;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px;font-size:medium"><div style="word-wrap:break-word">
Hi Jed,</div></span></span>
</div><div><br></div><div>I need to have this form for the simple input for PDSLin.</div><div><br></div><div>What do you mean by "COO"?</div>
<div><br></div>Any examples to show how to pull that new array for reading from the AIJ structure?</div></blockquote><div><br></div><div>First use MatGetSubmatrix() to pull the whole matrix into proc 0 if you are in parallel.</div>
<div>Then you MatGetRow() for each row, and write each entry in your form.</div><div><br></div><div>&nbsp; &nbsp;Matt</div><div>&nbsp;</div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex;">
<div style="word-wrap:break-word"><div>Thanks a lot!</div><div><br></div><font color="#888888"><div>R</div></font><div><div></div><div class="h5"><div><br><div><div>On Jul 14, 2011, at 11:28 AM, Jed Brown wrote:</div><br>
<blockquote type="cite"><div class="gmail_quote">On Thu, Jul 14, 2011 at 13:18, Xuefei (Rebecca) Yuan <span dir="ltr">&lt;<a href="mailto:xyuan@lbl.gov" target="_blank">xyuan@lbl.gov</a>&gt;</span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">

<div>I would like to store this sparse MXN matrix A in the following form:<br>
....................................................<br>
<br>
line1: &nbsp;M &nbsp; &nbsp; &nbsp; N &nbsp; &nbsp; &nbsp; numberofnonzeros<br>
line2: &nbsp;i &nbsp; &nbsp; &nbsp; j &nbsp; &nbsp; &nbsp; value<br>
line3: &nbsp;i &nbsp; &nbsp; &nbsp; j &nbsp; &nbsp; &nbsp; value</div></blockquote></div><br><div>This "COO" format is not an efficient representation. Why do you want to store it this way?</div><div><br></div><div>(You can produce this by getting the AIJ structure, allocating a new array for the row indices, and filling it redundantly from the AIJ structure.</div>

</blockquote></div><br></div></div></div></div></blockquote></div><br><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>
</blockquote></div><br></div></body></html>