Jak używać Gcc 4.6.0 libquadmath i __float128 na x86 i x86_64

Question

Mam program średniej wielkości C99, który używa typu long double (80bit) do obliczeń zmiennoprzecinkowych. Chcę poprawić precyzję z nowym rozszerzeniem GCC 4.6 __float128. Jak się okazuje, jest to 128-bitowa precyzyjna matematyka oparta na oprogramowaniu.

Jak skonwertować program z klasycznego long double 80-bitowego na quadowe 128-bitowe z emulacją oprogramowania z pełną precyzją? Co muszę zmienić? Flagi kompilatorów, źródła?

Mój program ma odczyt pełnej wartości dokładności z strtod, robi wiele różnych operacji na nich (tak jak + - */sin, cos, exp i drugiej z <math.h>) i printf -ing z nich.

PS: mimo że float128 jest zadeklarowany tylko dla Fortran (REAL * 16), libquadmath jest napisany w C i używa float128. Nie jestem pewien, czy GCC konwertuje operacje na float128 na bibliotekę wykonawczą, czy nie, i nie jestem pewien, jak przenieść z długich kopii do __float128 w moich źródłach.

PPS: Nie ma dokumentacji na "C" trybie język gcc: http://gcc.gnu.org/onlinedocs/gcc/Floating-Types.html

„kompilator GNU C wspomaga ... 128 bit (TFmode) pływających typów Wsparcie dla dodatkowych typów zawiera operatory arytmetyczne.: dodawanie, odejmowanie, mnożenie, dzielenie; unarne operatory arytmetyczne; operatory relacyjne; operatorów równości ... __float128 typy są obsługiwane na i386, x86_64"

Answer 1

Jak przekonwertować mój program z Classi c długa podwójna liczba bitów od 80 do 128 bitów z emulacją oprogramowania z pełną precyzją? Co muszę zmienić? Flagi kompilatorów, źródła?

Musisz najnowszą wersję oprogramowania, GCC przy wsparciu __float128 typu (4.6 i nowsze) oraz libquadmath (obsługiwanym tylko na x86 and x86_64 targets, w IA64 i HPPA z nowszymi GCC). Należy dodać łącznik flagę -lquadmath (The cannot find -lquadmath' pokaże, że masz zainstalowany bez libquadmath)

• Dodaj #include <quadmath.h> nagłówek mieć definicje makr i funkcji.
• Powinieneś zmodyfikować wszystkie definicje zmiennych long double na __float128.
• zmienne kompleks może być zmieniony na __complex128 typu (quadmath.h) lub bezpośrednio z typedef _Complex float __attribute__((mode(TC))) _Complex128;
• Wszystkie proste operacje arytmetyczne są automatycznie obsługiwane przez GCC (w przeliczeniu na wywołań funkcji pomocniczych, takich jak __*tf3()).
• W przypadku korzystania z makro jak LDBL_*, zastąpić je FLT128_* (pełna lista http://gcc.gnu.org/onlinedocs/libquadmath/Typedef-and-constants.html#Typedef-and-constants)
• Jeśli potrzebujesz konkretnych stałych jak PI (M_PI) lub E (M_E) z poczwórnej precyzji, należy użyć predefiniowanych stałych z q sufiksu (M_*q), jak M_PIq i M_Eq (pełna lista http://gcc.gnu.org/onlinedocs/libquadmath/Typedef-and-constants.html#Typedef-and-constants)
• stałe zdefiniowane przez użytkownika mogą być napisane z Q przyrostek, jak 1.3000011111111Q
• Wszystkie wywołania funkcji matematycznych należy wymienić dowcip h *q wersje, jak sqrtq(), sinq() (pełna lista http://gcc.gnu.org/onlinedocs/libquadmath/Math-Library-Routines.html#Math-Library-Routines)
• Reading czterordzeniowy pływak z łańcucha powinny być wykonane z __float128 strtoflt128 (const char *s, char **sp) - http://gcc.gnu.org/onlinedocs/libquadmath/strtoflt128.html#strtoflt128 (Uwaga, w starszych libquadmaths mogą wystąpić pewne błędy w strtoflt128, zrobić dokładnie sprawdzić czy)
• Drukowanie __float128 odbywa się za pomocą funkcji quadmath_snprintf. W dystrybucjach Linuksa z najnowszym glibc funkcja będzie automatycznie rejestrowana przez libquadmath, aby obsłużyć Q (może być również q) modyfikator długości ze specyfikatorów konwersji a, A, e, E, f, F, g, G we wszystkich printf s/sprintf s, tak jak zrobił to L dla długich debelów. Przykład: printf ("%Qe", 1.2Q), http://gcc.gnu.org/onlinedocs/libquadmath/quadmath_005fsnprintf.html#quadmath_005fsnprintf

Warto również wiedzieć, że od 4,6 gfortran użyje __float128 typ dla podwójnej precyzji, jeśli opcja -fdefault-real-8 dano i nie było opcji -fdefault-double-8. Może to być problem, ponieważ 128-bitowe podwójne jest znacznie wolniejsze niż standardowe podwójne podwójne na wielu platformach z powodu obliczeń programowych. (Dzięki postowi przez glennglockwood http://glennklockwood.blogspot.com/2014/02/linux-perf-libquadmath-and-gfortrans.html)