[petsc-users] Questions about MatConvert BAIJ to AIJ

[Chung-Kan Huang]

In this case I expect the result should be byte-wise identical.  Should I
expect some roundoff? Why?

Thanks,

On Wed, Mar 18, 2015 at 12:17 PM, Jed Brown <jed at jedbrown.org> wrote:

> Chung-Kan Huang <ckhuangf at gmail.com> writes:
>
> > Hi,
> >
> > I created matrix with MATBAIJ and I use
> >     MatConvert(J, MATAIJ, MAT_REUSE_MATRIX, & J);
> > to convert it to MATAIJ.
> >
> > I ran some tests
> > MATBAIJ: create & use MATBAIJ for later assembly
> > MATBAIJ2MATAIJ: create as MATBAIJ then convert to MATAIJ & use MATAIJ for
> > later assembly
> > MATAIJ: create & use MATAIJ for later assembly
> >
> > and I found following results
> >
> >
> > MATBAIJ == MATBAIJ2MATAIJ != MATAIJ
> > as I was expected to see
> > MATBAIJ != MATBAIJ2MATAIJ == MATAIJ
>
> How are you defining equality?
>



-- 

*Cheers*
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20150318/e524f15f/attachment.html>