<div dir="ltr"><div class="gmail_extra"><div class="gmail_quote">On Thu, Oct 15, 2015 at 12:10 PM, K. N. Ramachandran <span dir="ltr"><<a href="mailto:knram06@gmail.com" target="_blank">knram06@gmail.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr">Hello Matt,<br><div class="gmail_extra"><br><div class="gmail_quote">On Thu, Oct 15, 2015 at 12:18 PM, Matthew Knepley <span dir="ltr"><<a href="mailto:knepley@gmail.com" target="_blank">knepley@gmail.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><br><div class="gmail_extra"><div class="gmail_quote"><span></span><div>The bad design decision is using a storage format designed for high speed MatMult() to use for data interchange. You should be</div><div>using the API for value insertion rather than some storage format.</div></div></div></div></blockquote><div><br></div><div>Ok, noted. Instead of using CreateMatSeqAIJWithArrays, I have also tried setting up the Mat passed into the ComputeMatrix function, using MatSetValue where the cols and values are referred to values in the CSR format, but I still get the same error. It looks like the Mat is not preallocated properly in the DM, even when I use DMDACreate3d.<br></div></div></div></div></blockquote><div><br></div><div>Another explanation is that the matrix is preallocated correctly, but you are mapping the row and columns indices differently than the</div><div>DM. This what I was refering to when I said that getting the parallel partition correct would be challenging.</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 dir="ltr"><div class="gmail_extra"><div class="gmail_quote"><div></div><div>Could you suggest what I can do to resolve this?<br></div></div><br></div><div class="gmail_extra">Thanks,<br></div><div class="gmail_extra"><div><div dir="ltr">K.N.Ramachandran<br><div>Ph: <a href="tel:814-441-4279" value="+18144414279" target="_blank">814-441-4279</a></div></div></div>
</div></div>
</blockquote></div><br><br clear="all"><div><br></div>-- <br><div class="gmail_signature">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</div>
</div></div>