On Fri, May 25, 2012 at 9:49 AM, Jed Brown <span dir="ltr"><<a href="mailto:jedbrown@mcs.anl.gov" target="_blank">jedbrown@mcs.anl.gov</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">
<p>No, its not okay.</p>
<p>I think MatGetArray needs to be removed. Pretty much any use of it is wrong. It is not even the right thing for dense matrices because they should be distributed differently than PETSc does now.</p></blockquote><div>Without it, people cannot get directly to the AIJ data structure. No matter how wrong it is,</div>
<div>it is a constant request. I would keep it in just to satisfy that large number of petsc-maints.</div><div><br></div><div>   Matt </div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<p>On May 25, 2012 8:44 AM, "Hui Zhang" <<a href="mailto:mike.hui.zhang@hotmail.com" target="_blank">mike.hui.zhang@hotmail.com</a>> wrote:</p><div class="HOEnZb"><div class="h5"><div class="gmail_quote"><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">

A new question: can I use MatSetValues in the following way?<br>
<br>
MatGetArray(A,&a);<br>
a[0]= 1;<br>
MatSetValues(A,...);    /* in particular, is this allowed inside GetArray/RestoreArray? */<br>
MatRestoreArray(A,&a);<br>
<br>
Thanks!<br>
Hui<br>
<br>
On May 14, 2012, at 2:29 PM, Barry Smith wrote:<br>
<br>
><br>
>   If B is 50% dense then store it in a dense format. It will be much faster and the extra memory is minimal.  Even 30% dense.<br>
><br>
>   Barry<br>
><br>
> On May 14, 2012, at 5:13 AM, Hui Zhang wrote:<br>
><br>
>> I have two matrices A and B stored in sparse format.  While A is really sparse, B is relatively<br>
>> dense (say non-zeros entries about 50%).  Now to multiply the two matrices, among<br>
>><br>
>>     A*B<br>
>> or<br>
>>    (B^T * A^T)^T<br>
>><br>
>> which is better, or no big difference ?<br>
>><br>
>> Thanks!<br>
>><br>
><br>
><br>
<br>
</blockquote></div>
</div></div></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>