On Tue, Feb 7, 2012 at 3:34 PM, Jed Brown <span dir="ltr"><<a href="mailto:jedbrown@mcs.anl.gov">jedbrown@mcs.anl.gov</a>></span> wrote:<br><div class="gmail_quote"><blockquote class="gmail_quote" style="margin:0pt 0pt 0pt 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">
<div dir="ltr"><div class="im"><div class="gmail_quote">On Tue, Feb 7, 2012 at 23:18, Barry Smith <span dir="ltr"><<a href="mailto:bsmith@mcs.anl.gov" target="_blank">bsmith@mcs.anl.gov</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0pt 0pt 0pt 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">
I have not tried complex quad precision. If you are right than it may not be possible right out of the box with std:complex. We are not in position to provide a proper complex class for this case (for the same reasons I listed before) so if someone wants PETSc with complex quad they may have their work cut out for them.</blockquote>
</div><br></div><div>It's spelled __complex128. <a href="http://gcc.gnu.org/onlinedocs/gcc-4.6.2/libquadmath.pdf" target="_blank">http://gcc.gnu.org/onlinedocs/gcc-4.6.2/libquadmath.pdf</a></div><div><br></div><div>Perhaps it's a shame they didn't use C99 complex __float128, but such is life.</div>
</div></blockquote><div><br>Good to know; I didn't know libquadmath provided that. The catch is that it does not allow one to easily extend templates that previously handled 32-bit and 64-bit base types. On the other hand, a custom complex class could easily just call the built-in functions for __complex128 when instantiated with a base type of __float128. <br>
<br>Jack<br></div></div><br>