<div dir="ltr"><br><div class="gmail_extra"><br><div class="gmail_quote">On Wed, Mar 18, 2015 at 3:09 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:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><span class=""><br>
> On Mar 18, 2015, at 3:07 PM, Chung-Kan Huang <<a href="mailto:ckhuangf@gmail.com">ckhuangf@gmail.com</a>> wrote:<br>
><br>
</span><span class="">> So I guess I misunderstood the meaning of " inplace conversion".<br>
> and what do MAT_INITIAL_MATRIX and MAT_REUSE_MATRIX do anyway?<br>
<br>
</span>  No you have not misunderstood it.<br>
<br>
Currently<br>
   MAT_REUSE_MATRIX is only supported for inplace conversion, otherwise use<br>
   MAT_INITIAL_MATRIX.<br>
<br>
  For testing purposes you can use MAT_INITIAL_MATRIX.<br>
<div class=""><div class="h5"><br>
</div><div></div></div></blockquote><div><br>I will try to convert the matrix in a new matrix using MAT_INITIAL_MATRIX. <br></div><div>while <br>MatConvert(J, MATAIJ, MAT_INITIAL_MATRIX, & A); <br></div><div>will convert the matrix and port the results into a new<br></div><div>what does <br>MatConvert(J, MATAIJ, MAT_REUSE_MATRIX, & J);<br></div><div>actually do?<br><br></div><div>Thanks,<br></div><div> </div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div class=""><div><br></div></div></blockquote><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div class=""><div class="h5">
><br>
><br>
><br>
><br>
> On Wed, Mar 18, 2015 at 2:48 PM, Jed Brown <<a href="mailto:jed@jedbrown.org">jed@jedbrown.org</a>> wrote:<br>
> 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>
> 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>
><br>
><br>
><br>
> --<br>
> Cheers<br>
><br>
<br>
</div></div></blockquote></div><br><br clear="all"><br>-- <br><div class="gmail_signature"><p><b>Cheers</b></p></div>
</div></div>