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 <[email protected]> wrote: > Chung-Kan Huang <[email protected]> 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*
