So I guess I misunderstood the meaning of " inplace conversion". and what do MAT_INITIAL_MATRIX and MAT_REUSE_MATRIX do anyway?
On Wed, Mar 18, 2015 at 2:48 PM, Jed Brown <[email protected]> wrote: > Chung-Kan Huang <[email protected]> writes: > > > Yes, I was suprised that I can still use MatSetValuesBlocked after I did > > MatConvert MATBAIJ to MATAIJ. > > > > Then I just found the mat type is still seqbaij which means the > MatConvert > > didn't work. > > > > This is how I did MatConvert > > > > MatConvert(J, MATAIJ, MAT_REUSE_MATRIX, & J); > > > > where J is my Jacobian matrix. > > Seems to me this should fail because you didn't MatConvert this matrix > yet. Anyway, you should use MAT_INITIAL_MATRIX and put the result in a > new Mat. > -- *Cheers*
