On Wed, Nov 4, 2009 at 2:06 PM, Barry Smith <span dir="ltr"><<a href="mailto:bsmith@mcs.anl.gov">bsmith@mcs.anl.gov</a>></span> wrote:<br><div class="gmail_quote"><blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">
<br>
   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.<br>
</blockquote><div><br>I figured that his types were just C++ types once he said they were provided by a package.<br> </div><blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">

   I am investigating qd now. Is there any MPI implementation that uses qd? Do you need MPI?<font color="#888888"><br></font></blockquote><div><br>Can't we just make a type the correct number of bytes long?<br><br>  Matt<br>
 </div><blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;"><font color="#888888">
   Barry</font><div><div></div><div class="h5"><br>
<br>
<br>
On Nov 4, 2009, at 1:33 PM, Matthew Knepley wrote:<br>
<br>
<blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">
If you are alright with EVERYTHING running with that, it already works. Just set PetscScalar<br>
to that type. I would help you setup a configure switch for it if needed.<br>
<br>
If you only want it for some operations, we would have to hear more specifics.<br>
<br>
  Thanks,<br>
<br>
    Matt<br>
<br>
On Wed, Nov 4, 2009 at 10:45 AM, Todd Munson <<a href="mailto:tmunson@mcs.anl.gov" target="_blank">tmunson@mcs.anl.gov</a>> wrote:<br>
<br>
Barry and Matt,<br>
<br>
I was wondering how hard it would be to modify PETSc to support double-double and quad-double arithmetic using the QD package<br>
from LBL?<br>
<br>
Its looking more and more likely that I may need this functionality for some upcoming projects that have bad scaling issues and I don't want to have to reimplement KSP and PC methods...<br>
<br>
Cheers, Todd.<br>
<br>
<br>
<br>
<br>
-- <br>
What most experimenters take for granted before they begin their experiments is infinitely more interesting than any results to which their experiments lead.<br>
-- Norbert Wiener<br>
</blockquote>
<br>
</div></div></blockquote></div><br><br clear="all"><br>-- <br>What most experimenters take for granted before they begin their experiments is infinitely more interesting than any results to which their experiments lead.<br>
-- Norbert Wiener<br>