[petsc-users] MatMult() returning different values depending on # of processors?

Jed Brown jed at jedbrown.org
Thu Jun 26 22:48:09 CDT 2014


Barry Smith <bsmith at mcs.anl.gov> writes:
>    Note also that the Intel floating point registers are 80 bits wide
>    while double precision numbers in memory are 64 bit. Thus moving a
>    partially computed sum to memory and then back to register wipes
>    out 16 bits of the number.

This comment refers to the x87 unit which fell out of favor when SSE
provided better performance and vectorization.  It is only used now if
you compile for x86 (32-bit) without SSE support.
-------------- next part --------------
A non-text attachment was scrubbed...
Name: not available
Type: application/pgp-signature
Size: 818 bytes
Desc: not available
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20140626/7d99bc89/attachment-0001.pgp>


More information about the petsc-users mailing list