PETSc in double-double and quad-double...

2009-11-04 Thread Matthew Knepley
any results to which their >> experiments lead. >> -- Norbert Wiener >> > > -- What most experimenters take for granted before they begin their experiments is infinitely more interesting than any results to which their experiments lead. -- Norbert Wiener -- next part -- An HTML attachment was scrubbed... URL: <http://lists.mcs.anl.gov/pipermail/petsc-dev/attachments/20091104/eb82322e/attachment.html>

PETSc in double-double and quad-double...

2009-11-04 Thread Jed Brown
--- next part -- A non-text attachment was scrubbed... Name: signature.asc Type: application/pgp-signature Size: 261 bytes Desc: OpenPGP digital signature URL: <http://lists.mcs.anl.gov/pipermail/petsc-dev/attachments/20091104/ffa96917/attachment.pgp>

PETSc in double-double and quad-double...

2009-11-04 Thread Todd Munson
If we can do qd_real, we could also try using the arprec package and get arbitrary precision math. I might be interested in that as well... All are written as C++ types with operator overloading... Cheers, Todd. On Wed, 4 Nov 2009, Jed Brown wrote: > Matthew Knepley wrote: >> On Wed, Nov 4,

PETSc in double-double and quad-double...

2009-11-04 Thread Matthew Knepley
t now for several things in PETSc. Matt > Jed > > -- What most experimenters take for granted before they begin their experiments is infinitely more interesting than any results to which their experiments lead. -- Norbert Wiener -- next part -- An HTML attachment was scrubbed... URL: <http://lists.mcs.anl.gov/pipermail/petsc-dev/attachments/20091104/74a5998e/attachment.html>

PETSc in double-double and quad-double...

2009-11-04 Thread Jed Brown
text attachment was scrubbed... Name: signature.asc Type: application/pgp-signature Size: 261 bytes Desc: OpenPGP digital signature URL: <http://lists.mcs.anl.gov/pipermail/petsc-dev/attachments/20091104/daa035e4/attachment.pgp>

PETSc in double-double and quad-double...

2009-11-04 Thread Barry Smith
This is tricky because some C/C++ compilers long double uses 128 bits of space but actually (on Intel) only uses 80 of the bits (those are the bits that floating point unit handles). So you do not really get 128 bits. Also MPI may not properly handle the 128 bit doubles. I am invest