On Wed, Nov 4, 2009 at 2:21 PM, Jed Brown <span dir="ltr"><<a href="mailto:jed@59a2.org">jed@59a2.org</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;">
<div class="im">Matthew Knepley wrote:<br>
> On Wed, Nov 4, 2009 at 2:06 PM, Barry Smith <<a href="mailto:bsmith@mcs.anl.gov">bsmith@mcs.anl.gov</a><br>
</div><div class="im">> <mailto:<a href="mailto:bsmith@mcs.anl.gov">bsmith@mcs.anl.gov</a>>> wrote:<br>
><br>
><br>
>       This is tricky because some C/C++ compilers long double uses 128<br>
>     bits of space but actually (on Intel) only uses 80 of the bits<br>
>     (those are the bits that floating point unit handles). So you do not<br>
>     really get 128 bits. Also MPI may not properly handle the 128 bit<br>
>     doubles.<br>
><br>
><br>
> I figured that his types were just C++ types once he said they were<br>
> provided by a package.<br>
<br>
</div>QD is honest quad-double (32 bytes, all of them significant)<br>
<br>
  <a href="http://crd.lbl.gov/%7Edhbailey/mpdist/" target="_blank">http://crd.lbl.gov/~dhbailey/mpdist/</a><br>
<br>
They are C++ types with appropriate overloading, but the operations are<br>
implemented in terms of conventional double-precision where as "long<br>
double" is a different beast.<br><div class="im"></div></blockquote><div><br>Okay, then this should be no harder than complex.<br> </div><blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">
<div class="im">
>       I am investigating qd now. Is there any MPI implementation that<br>
>     uses qd? Do you need MPI?<br>
><br>
><br>
> Can't we just make a type the correct number of bytes long?<br>
<br>
</div>I think so (on homogeneous hardware), but you have to create your own<br>
version of each MPI_Op that does the right thing with this longer type.<font color="#888888"><br></font></blockquote><div><br>Yep, we do that now for several things in PETSc.<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">
Jed<br>
<br>
</font></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>