[petsc-users] Questions about MatConvert BAIJ to AIJ

[Chung-Kan Huang]

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 <jed at jedbrown.org> wrote:

> Chung-Kan Huang <ckhuangf at gmail.com> 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*
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20150318/4272db46/attachment.html>