<div dir="ltr"><div>So I guess I misunderstood the meaning of " inplace conversion".<br></div>and what do MAT_INITIAL_MATRIX and <span class="im">MAT_REUSE_MATRIX do </span>anyway?<br><br><br><br></div><div class="gmail_extra"><br><div class="gmail_quote">On Wed, Mar 18, 2015 at 2:48 PM, Jed Brown <span dir="ltr"><<a href="mailto:jed@jedbrown.org" target="_blank">jed@jedbrown.org</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><span class="">Chung-Kan Huang <<a href="mailto:ckhuangf@gmail.com">ckhuangf@gmail.com</a>> writes:<br>
<br>
> Yes, I was suprised that I can still use MatSetValuesBlocked after I did<br>
> MatConvert MATBAIJ to MATAIJ.<br>
><br>
> Then I just found the mat type is still seqbaij which means the MatConvert<br>
> didn't work.<br>
><br>
> This is how I did MatConvert<br>
><br>
> MatConvert(J, MATAIJ, MAT_REUSE_MATRIX, & J);<br>
><br>
> where J is my Jacobian matrix.<br>
<br>
</span>Seems to me this should fail because you didn't MatConvert this matrix<br>
yet.  Anyway, you should use MAT_INITIAL_MATRIX and put the result in a<br>
new Mat.<br>
</blockquote></div><br><br clear="all"><br>-- <br><div class="gmail_signature"><p><strong>Cheers</strong></p></div>
</div>