Cómo usar Gcc 4.6.0 libquadmath y __float128 en x86 y x86_64

Question

Tengo un programa de tamaño medio C99 que usa el tipo long double (80 bits) para el cálculo en coma flotante. Quiero mejorar la precisión con la nueva extensión GCC 4.6 __float128. Según entiendo, es una matemática de precisión de 128 bits emulada por software.

¿Cómo debo convertir mi programa de doble larga clásica de 80 bits a cuádruple flotadores de 128 bits con emulación de software de precisión total? ¿Qué necesito cambiar? Banderas del compilador, fuentes?

Mi programa tienen la lectura de los valores de plena precisión con strtod, haciendo un montón de diferentes operaciones sobre ellos (como + - */sen, cos, exp y otra de <math.h>) y printf -ing de ellos.

PD: a pesar de que float128 está declarado solo para Fortran (REAL * 16), el libquadmath está escrito en C y usa float128. No estoy seguro si GCC convertirá las operaciones en float128 a la biblioteca de tiempo de ejecución o no y no estoy seguro de cómo migrar de doble largo a __float128 en mis fuentes.

PPS: Hay una documentación sobre "C" modo de idioma gcc: http://gcc.gnu.org/onlinedocs/gcc/Floating-Types.html

"compilador GNU C soporta ... 128 bits (TFmode) tipos flotantes soporte para tipos adicionales incluye los operadores aritméticos.: sumar, restar, multiplicar, dividir; operadores aritméticos unitarios; operadores relacionales; operadores de igualdad ... __float128 tipos son compatibles en i386, x86_64"

Answer 1

¿Cómo debo convertir mi programa de clasificación ¿c largo doble de 80 bit a quad flota de 128 bit con emulación de software de precisión total? ¿Qué necesito cambiar? Banderas del compilador, fuentes?

se necesita un software reciente, versión de GCC con el apoyo de __float128 tipo (4.6 y posteriores) y libquadmath (compatible sólo con x86 and x86_64 targets; en IA64 y HPPA con nuevo GCC). Usted debe agregar el indicador de enlazador -lquadmath (el cannot find -lquadmath' mostrará que ha instalado ningún libquadmath)

• Añadir #include <quadmath.h> encabezado tener definiciones de macros y de función.
• Debe modificar todas las definiciones de las variables long double a __float128.
• las variables complejas pueden ser cambiados para __complex128 tipo (quadmath.h) o directamente con typedef _Complex float __attribute__((mode(TC))) _Complex128;
• Todas las operaciones aritméticas sencillas son manejados automáticamente por GCC (convertidos a las llamadas de funciones auxiliares como __*tf3()).
• Si usa cualquier macro como LDBL_*, reemplazarlos con FLT128_* (la lista completa http://gcc.gnu.org/onlinedocs/libquadmath/Typedef-and-constants.html#Typedef-and-constants)
• Si necesita algunas constantes específicas como pi (M_PI) o E (M_E) con una precisión cuádruple, utilice constantes predefinidas con q sufijo (M_*q ingenio), al igual que M_PIq y M_Eq (la lista completa http://gcc.gnu.org/onlinedocs/libquadmath/Typedef-and-constants.html#Typedef-and-constants)
• constantes definidas por el usuario pueden ser escritos con Q sufijo, como 1.3000011111111Q
• Todas las llamadas a funciones matemáticas deben ser reemplazados h *q versiones, como sqrtq(), sinq() (la lista completa http://gcc.gnu.org/onlinedocs/libquadmath/Math-Library-Routines.html#Math-Library-Routines)
• lectura quad-float de cadena se debe hacer con __float128 strtoflt128 (const char *s, char **sp)-http://gcc.gnu.org/onlinedocs/libquadmath/strtoflt128.html#strtoflt128 (Advertencia, en libquadmaths de edad avanzada que puede haber algunos errores en strtoflt128, hacer un doble control)
• La impresión del __float128 se realiza con la ayuda de la función quadmath_snprintf. En las distribuciones de Linux con glibc reciente, la función se registrará automágicamente por libquadmath para manejar Q (también puede ser q) modificador de longitud de a, A, e, E, f, F, g, G especificadores de conversión en todos printf s/sprintf s, como lo hizo L para dobles largos. Ejemplo: printf ("%Qe", 1.2Q), http://gcc.gnu.org/onlinedocs/libquadmath/quadmath_005fsnprintf.html#quadmath_005fsnprintf

También debe saber, que desde 4.6 gfortran utilizará __float128 tipo para doble precisión, si la opción -fdefault-real-8 se le dio y no había ninguna opción -fdefault-double-8. Esto puede ser un problema, ya que 128 long double es mucho más lento que el double estándar doble en muchas plataformas debido al cálculo del software. (Gracias a la publicación de glennglockwood http://glennklockwood.blogspot.com/2014/02/linux-perf-libquadmath-and-gfortrans.html)