ViewVC Help
View File | Revision Log | Show Annotations | Download File
/cvs/libecb/ecb.h
(Generate patch)

Comparing libecb/ecb.h (file contents):
Revision 1.161 by sf-exg, Thu Mar 19 20:55:22 2015 UTC vs.
Revision 1.177 by root, Mon Nov 19 00:27:38 2018 UTC

… …
40 40
41#ifndef ECB_H 41#ifndef ECB_H
42#define ECB_H 42#define ECB_H
43 43
44/* 16 bits major, 16 bits minor */ 44/* 16 bits major, 16 bits minor */
45#define ECB_VERSION 0x00010004 45#define ECB_VERSION 0x00010005
46 46
47#ifdef _WIN32 47#ifdef _WIN32
48 typedef signed char int8_t; 48 typedef signed char int8_t;
49 typedef unsigned char uint8_t; 49 typedef unsigned char uint8_t;
50 typedef signed short int16_t; 50 typedef signed short int16_t;
… …
67 typedef uint32_t uintptr_t; 67 typedef uint32_t uintptr_t;
68 typedef int32_t intptr_t; 68 typedef int32_t intptr_t;
69 #endif 69 #endif
70#else 70#else
71 #include <inttypes.h> 71 #include <inttypes.h>
72 #if UINTMAX_MAX > 0xffffffffU 72 #if (defined INTPTR_MAX ? INTPTR_MAX : ULONG_MAX) > 0xffffffffU
73 #define ECB_PTRSIZE 8 73 #define ECB_PTRSIZE 8
74 #else 74 #else
75 #define ECB_PTRSIZE 4 75 #define ECB_PTRSIZE 4
76 #endif 76 #endif
77#endif 77#endif
… …
115 #define ECB_CLANG_EXTENSION(x) 0 115 #define ECB_CLANG_EXTENSION(x) 0
116#endif 116#endif
117 117
118#define ECB_CPP (__cplusplus+0) 118#define ECB_CPP (__cplusplus+0)
119#define ECB_CPP11 (__cplusplus >= 201103L) 119#define ECB_CPP11 (__cplusplus >= 201103L)
120#define ECB_CPP14 (__cplusplus >= 201402L)
121#define ECB_CPP17 (__cplusplus >= 201703L)
120 122
121#if ECB_CPP 123#if ECB_CPP
122 #define ECB_C 0 124 #define ECB_C 0
123 #define ECB_STDC_VERSION 0 125 #define ECB_STDC_VERSION 0
124#else 126#else
… …
126 #define ECB_STDC_VERSION __STDC_VERSION__ 128 #define ECB_STDC_VERSION __STDC_VERSION__
127#endif 129#endif
128 130
129#define ECB_C99 (ECB_STDC_VERSION >= 199901L) 131#define ECB_C99 (ECB_STDC_VERSION >= 199901L)
130#define ECB_C11 (ECB_STDC_VERSION >= 201112L) 132#define ECB_C11 (ECB_STDC_VERSION >= 201112L)
133#define ECB_C17 (ECB_STDC_VERSION >= 201710L)
131 134
132#if ECB_CPP 135#if ECB_CPP
133 #define ECB_EXTERN_C extern "C" 136 #define ECB_EXTERN_C extern "C"
134 #define ECB_EXTERN_C_BEG ECB_EXTERN_C { 137 #define ECB_EXTERN_C_BEG ECB_EXTERN_C {
135 #define ECB_EXTERN_C_END } 138 #define ECB_EXTERN_C_END }
… …
148 #define ECB_NO_SMP 1 151 #define ECB_NO_SMP 1
149#endif 152#endif
150 153
151#if ECB_NO_SMP 154#if ECB_NO_SMP
152 #define ECB_MEMORY_FENCE do { } while (0) 155 #define ECB_MEMORY_FENCE do { } while (0)
156#endif
157
158/* http://www-01.ibm.com/support/knowledgecenter/SSGH3R_13.1.0/com.ibm.xlcpp131.aix.doc/compiler_ref/compiler_builtins.html */
159#if __xlC__ && ECB_CPP
160 #include <builtins.h>
161#endif
162
163#if 1400 <= _MSC_VER
164 #include <intrin.h> /* fence functions _ReadBarrier, also bit search functions _BitScanReverse */
153#endif 165#endif
154 166
155#ifndef ECB_MEMORY_FENCE 167#ifndef ECB_MEMORY_FENCE
156 #if ECB_GCC_VERSION(2,5) || defined __INTEL_COMPILER || (__llvm__ && __GNUC__) || __SUNPRO_C >= 0x5110 || __SUNPRO_CC >= 0x5110 168 #if ECB_GCC_VERSION(2,5) || defined __INTEL_COMPILER || (__llvm__ && __GNUC__) || __SUNPRO_C >= 0x5110 || __SUNPRO_CC >= 0x5110
157 #if __i386 || __i386__ 169 #if __i386 || __i386__
158 #define ECB_MEMORY_FENCE __asm__ __volatile__ ("lock; orb $0, -1(%%esp)" : : : "memory") 170 #define ECB_MEMORY_FENCE __asm__ __volatile__ ("lock; orb $0, -1(%%esp)" : : : "memory")
159 #define ECB_MEMORY_FENCE_ACQUIRE __asm__ __volatile__ ("" : : : "memory") 171 #define ECB_MEMORY_FENCE_ACQUIRE __asm__ __volatile__ ("" : : : "memory")
160 #define ECB_MEMORY_FENCE_RELEASE __asm__ __volatile__ ("") 172 #define ECB_MEMORY_FENCE_RELEASE __asm__ __volatile__ ("" : : : "memory")
161 #elif ECB_GCC_AMD64 173 #elif ECB_GCC_AMD64
162 #define ECB_MEMORY_FENCE __asm__ __volatile__ ("mfence" : : : "memory") 174 #define ECB_MEMORY_FENCE __asm__ __volatile__ ("mfence" : : : "memory")
163 #define ECB_MEMORY_FENCE_ACQUIRE __asm__ __volatile__ ("" : : : "memory") 175 #define ECB_MEMORY_FENCE_ACQUIRE __asm__ __volatile__ ("" : : : "memory")
164 #define ECB_MEMORY_FENCE_RELEASE __asm__ __volatile__ ("") 176 #define ECB_MEMORY_FENCE_RELEASE __asm__ __volatile__ ("" : : : "memory")
165 #elif __powerpc__ || __ppc__ || __powerpc64__ || __ppc64__ 177 #elif __powerpc__ || __ppc__ || __powerpc64__ || __ppc64__
166 #define ECB_MEMORY_FENCE __asm__ __volatile__ ("sync" : : : "memory") 178 #define ECB_MEMORY_FENCE __asm__ __volatile__ ("sync" : : : "memory")
179 #elif defined __ARM_ARCH_2__ \
180 || defined __ARM_ARCH_3__ || defined __ARM_ARCH_3M__ \
181 || defined __ARM_ARCH_4__ || defined __ARM_ARCH_4T__ \
182 || defined __ARM_ARCH_5__ || defined __ARM_ARCH_5E__ \
183 || defined __ARM_ARCH_5T__ || defined __ARM_ARCH_5TE__ \
184 || defined __ARM_ARCH_5TEJ__
185 /* should not need any, unless running old code on newer cpu - arm doesn't support that */
167 #elif defined __ARM_ARCH_6__ || defined __ARM_ARCH_6J__ \ 186 #elif defined __ARM_ARCH_6__ || defined __ARM_ARCH_6J__ \
168 || defined __ARM_ARCH_6K__ || defined __ARM_ARCH_6ZK__ 187 || defined __ARM_ARCH_6K__ || defined __ARM_ARCH_6ZK__ \
188 || defined __ARM_ARCH_6T2__
169 #define ECB_MEMORY_FENCE __asm__ __volatile__ ("mcr p15,0,%0,c7,c10,5" : : "r" (0) : "memory") 189 #define ECB_MEMORY_FENCE __asm__ __volatile__ ("mcr p15,0,%0,c7,c10,5" : : "r" (0) : "memory")
170 #elif defined __ARM_ARCH_7__ || defined __ARM_ARCH_7A__ \ 190 #elif defined __ARM_ARCH_7__ || defined __ARM_ARCH_7A__ \
171 || defined __ARM_ARCH_7M__ || defined __ARM_ARCH_7R__ 191 || defined __ARM_ARCH_7R__ || defined __ARM_ARCH_7M__
172 #define ECB_MEMORY_FENCE __asm__ __volatile__ ("dmb" : : : "memory") 192 #define ECB_MEMORY_FENCE __asm__ __volatile__ ("dmb" : : : "memory")
173 #elif __aarch64__ 193 #elif __aarch64__
174 #define ECB_MEMORY_FENCE __asm__ __volatile__ ("dmb ish" : : : "memory") 194 #define ECB_MEMORY_FENCE __asm__ __volatile__ ("dmb ish" : : : "memory")
175 #elif (__sparc || __sparc__) && !__sparcv8 195 #elif (__sparc || __sparc__) && !(__sparc_v8__ || defined __sparcv8)
176 #define ECB_MEMORY_FENCE __asm__ __volatile__ ("membar #LoadStore | #LoadLoad | #StoreStore | #StoreLoad" : : : "memory") 196 #define ECB_MEMORY_FENCE __asm__ __volatile__ ("membar #LoadStore | #LoadLoad | #StoreStore | #StoreLoad" : : : "memory")
177 #define ECB_MEMORY_FENCE_ACQUIRE __asm__ __volatile__ ("membar #LoadStore | #LoadLoad" : : : "memory") 197 #define ECB_MEMORY_FENCE_ACQUIRE __asm__ __volatile__ ("membar #LoadStore | #LoadLoad" : : : "memory")
178 #define ECB_MEMORY_FENCE_RELEASE __asm__ __volatile__ ("membar #LoadStore | #StoreStore") 198 #define ECB_MEMORY_FENCE_RELEASE __asm__ __volatile__ ("membar #LoadStore | #StoreStore")
179 #elif defined __s390__ || defined __s390x__ 199 #elif defined __s390__ || defined __s390x__
180 #define ECB_MEMORY_FENCE __asm__ __volatile__ ("bcr 15,0" : : : "memory") 200 #define ECB_MEMORY_FENCE __asm__ __volatile__ ("bcr 15,0" : : : "memory")
… …
352 #define ecb_deprecated __declspec (deprecated) 372 #define ecb_deprecated __declspec (deprecated)
353#else 373#else
354 #define ecb_deprecated ecb_attribute ((__deprecated__)) 374 #define ecb_deprecated ecb_attribute ((__deprecated__))
355#endif 375#endif
356 376
357#if __MSC_VER >= 1500 377#if _MSC_VER >= 1500
358 #define ecb_deprecated_message(msg) __declspec (deprecated (msg)) 378 #define ecb_deprecated_message(msg) __declspec (deprecated (msg))
359#elif ECB_GCC_VERSION(4,5) 379#elif ECB_GCC_VERSION(4,5)
360 #define ecb_deprecated_message(msg) ecb_attribute ((__deprecated__ (msg)) 380 #define ecb_deprecated_message(msg) ecb_attribute ((__deprecated__ (msg))
361#else 381#else
362 #define ecb_deprecated_message(msg) ecb_deprecated 382 #define ecb_deprecated_message(msg) ecb_deprecated
… …
371#define ecb_unused ecb_attribute ((__unused__)) 391#define ecb_unused ecb_attribute ((__unused__))
372#define ecb_const ecb_attribute ((__const__)) 392#define ecb_const ecb_attribute ((__const__))
373#define ecb_pure ecb_attribute ((__pure__)) 393#define ecb_pure ecb_attribute ((__pure__))
374 394
375#if ECB_C11 || __IBMC_NORETURN 395#if ECB_C11 || __IBMC_NORETURN
376 /* http://pic.dhe.ibm.com/infocenter/compbg/v121v141/topic/com.ibm.xlcpp121.bg.doc/language_ref/noreturn.html */ 396 /* http://www-01.ibm.com/support/knowledgecenter/SSGH3R_13.1.0/com.ibm.xlcpp131.aix.doc/language_ref/noreturn.html */
377 #define ecb_noreturn _Noreturn 397 #define ecb_noreturn _Noreturn
378#elif ECB_CPP11 398#elif ECB_CPP11
379 #define ecb_noreturn [[noreturn]] 399 #define ecb_noreturn [[noreturn]]
380#elif _MSC_VER >= 1200 400#elif _MSC_VER >= 1200
381 /* http://msdn.microsoft.com/en-us/library/k6ktzx3s.aspx */ 401 /* http://msdn.microsoft.com/en-us/library/k6ktzx3s.aspx */
… …
418#else 438#else
419 ecb_function_ ecb_const int ecb_ctz32 (uint32_t x); 439 ecb_function_ ecb_const int ecb_ctz32 (uint32_t x);
420 ecb_function_ ecb_const int 440 ecb_function_ ecb_const int
421 ecb_ctz32 (uint32_t x) 441 ecb_ctz32 (uint32_t x)
422 { 442 {
443#if 1400 <= _MSC_VER && (_M_IX86 || _M_X64 || _M_IA64 || _M_ARM)
444 unsigned long r;
445 _BitScanForward (&r, x);
446 return (int)r;
447#else
423 int r = 0; 448 int r = 0;
424 449
425 x &= ~x + 1; /* this isolates the lowest bit */ 450 x &= ~x + 1; /* this isolates the lowest bit */
426 451
427#if ECB_branchless_on_i386 452#if ECB_branchless_on_i386
… …
437 if (x & 0xff00ff00) r += 8; 462 if (x & 0xff00ff00) r += 8;
438 if (x & 0xffff0000) r += 16; 463 if (x & 0xffff0000) r += 16;
439#endif 464#endif
440 465
441 return r; 466 return r;
467#endif
442 } 468 }
443 469
444 ecb_function_ ecb_const int ecb_ctz64 (uint64_t x); 470 ecb_function_ ecb_const int ecb_ctz64 (uint64_t x);
445 ecb_function_ ecb_const int 471 ecb_function_ ecb_const int
446 ecb_ctz64 (uint64_t x) 472 ecb_ctz64 (uint64_t x)
447 { 473 {
474#if 1400 <= _MSC_VER && (_M_X64 || _M_IA64 || _M_ARM)
475 unsigned long r;
476 _BitScanForward64 (&r, x);
477 return (int)r;
478#else
448 int shift = x & 0xffffffffU ? 0 : 32; 479 int shift = x & 0xffffffff ? 0 : 32;
449 return ecb_ctz32 (x >> shift) + shift; 480 return ecb_ctz32 (x >> shift) + shift;
481#endif
450 } 482 }
451 483
452 ecb_function_ ecb_const int ecb_popcount32 (uint32_t x); 484 ecb_function_ ecb_const int ecb_popcount32 (uint32_t x);
453 ecb_function_ ecb_const int 485 ecb_function_ ecb_const int
454 ecb_popcount32 (uint32_t x) 486 ecb_popcount32 (uint32_t x)
… …
462 } 494 }
463 495
464 ecb_function_ ecb_const int ecb_ld32 (uint32_t x); 496 ecb_function_ ecb_const int ecb_ld32 (uint32_t x);
465 ecb_function_ ecb_const int ecb_ld32 (uint32_t x) 497 ecb_function_ ecb_const int ecb_ld32 (uint32_t x)
466 { 498 {
499#if 1400 <= _MSC_VER && (_M_IX86 || _M_X64 || _M_IA64 || _M_ARM)
500 unsigned long r;
501 _BitScanReverse (&r, x);
502 return (int)r;
503#else
467 int r = 0; 504 int r = 0;
468 505
469 if (x >> 16) { x >>= 16; r += 16; } 506 if (x >> 16) { x >>= 16; r += 16; }
470 if (x >> 8) { x >>= 8; r += 8; } 507 if (x >> 8) { x >>= 8; r += 8; }
471 if (x >> 4) { x >>= 4; r += 4; } 508 if (x >> 4) { x >>= 4; r += 4; }
472 if (x >> 2) { x >>= 2; r += 2; } 509 if (x >> 2) { x >>= 2; r += 2; }
473 if (x >> 1) { r += 1; } 510 if (x >> 1) { r += 1; }
474 511
475 return r; 512 return r;
513#endif
476 } 514 }
477 515
478 ecb_function_ ecb_const int ecb_ld64 (uint64_t x); 516 ecb_function_ ecb_const int ecb_ld64 (uint64_t x);
479 ecb_function_ ecb_const int ecb_ld64 (uint64_t x) 517 ecb_function_ ecb_const int ecb_ld64 (uint64_t x)
480 { 518 {
519#if 1400 <= _MSC_VER && (_M_X64 || _M_IA64 || _M_ARM)
520 unsigned long r;
521 _BitScanReverse64 (&r, x);
522 return (int)r;
523#else
481 int r = 0; 524 int r = 0;
482 525
483 if (x >> 32) { x >>= 32; r += 32; } 526 if (x >> 32) { x >>= 32; r += 32; }
484 527
485 return r + ecb_ld32 (x); 528 return r + ecb_ld32 (x);
529#endif
486 } 530 }
487#endif 531#endif
488 532
489ecb_function_ ecb_const ecb_bool ecb_is_pot32 (uint32_t x); 533ecb_function_ ecb_const ecb_bool ecb_is_pot32 (uint32_t x);
490ecb_function_ ecb_const ecb_bool ecb_is_pot32 (uint32_t x) { return !(x & (x - 1)); } 534ecb_function_ ecb_const ecb_bool ecb_is_pot32 (uint32_t x) { return !(x & (x - 1)); }
… …
547ecb_inline ecb_const uint32_t ecb_rotr32 (uint32_t x, unsigned int count) { return (x << (32 - count)) | (x >> count); } 591ecb_inline ecb_const uint32_t ecb_rotr32 (uint32_t x, unsigned int count) { return (x << (32 - count)) | (x >> count); }
548ecb_inline ecb_const uint64_t ecb_rotl64 (uint64_t x, unsigned int count) { return (x >> (64 - count)) | (x << count); } 592ecb_inline ecb_const uint64_t ecb_rotl64 (uint64_t x, unsigned int count) { return (x >> (64 - count)) | (x << count); }
549ecb_inline ecb_const uint64_t ecb_rotr64 (uint64_t x, unsigned int count) { return (x << (64 - count)) | (x >> count); } 593ecb_inline ecb_const uint64_t ecb_rotr64 (uint64_t x, unsigned int count) { return (x << (64 - count)) | (x >> count); }
550 594
551#if ECB_GCC_VERSION(4,3) || (ECB_CLANG_BUILTIN(__builtin_bswap32) && ECB_CLANG_BUILTIN(__builtin_bswap64)) 595#if ECB_GCC_VERSION(4,3) || (ECB_CLANG_BUILTIN(__builtin_bswap32) && ECB_CLANG_BUILTIN(__builtin_bswap64))
596 #if ECB_GCC_VERSION(4,8) || ECB_CLANG_BUILTIN(__builtin_bswap16)
597 #define ecb_bswap16(x) __builtin_bswap16 (x)
598 #else
552 #define ecb_bswap16(x) (__builtin_bswap32 (x) >> 16) 599 #define ecb_bswap16(x) (__builtin_bswap32 (x) >> 16)
600 #endif
553 #define ecb_bswap32(x) __builtin_bswap32 (x) 601 #define ecb_bswap32(x) __builtin_bswap32 (x)
554 #define ecb_bswap64(x) __builtin_bswap64 (x) 602 #define ecb_bswap64(x) __builtin_bswap64 (x)
603#elif _MSC_VER
604 #include <stdlib.h>
605 #define ecb_bswap16(x) ((uint16_t)_byteswap_ushort ((uint16_t)(x)))
606 #define ecb_bswap32(x) ((uint32_t)_byteswap_ulong ((uint32_t)(x)))
607 #define ecb_bswap64(x) ((uint64_t)_byteswap_uint64 ((uint64_t)(x)))
555#else 608#else
556 ecb_function_ ecb_const uint16_t ecb_bswap16 (uint16_t x); 609 ecb_function_ ecb_const uint16_t ecb_bswap16 (uint16_t x);
557 ecb_function_ ecb_const uint16_t 610 ecb_function_ ecb_const uint16_t
558 ecb_bswap16 (uint16_t x) 611 ecb_bswap16 (uint16_t x)
559 { 612 {
… …
584#endif 637#endif
585 638
586/* try to tell the compiler that some condition is definitely true */ 639/* try to tell the compiler that some condition is definitely true */
587#define ecb_assume(cond) if (!(cond)) ecb_unreachable (); else 0 640#define ecb_assume(cond) if (!(cond)) ecb_unreachable (); else 0
588 641
589ecb_inline ecb_const unsigned char ecb_byteorder_helper (void); 642ecb_inline ecb_const uint32_t ecb_byteorder_helper (void);
590ecb_inline ecb_const unsigned char 643ecb_inline ecb_const uint32_t
591ecb_byteorder_helper (void) 644ecb_byteorder_helper (void)
592{ 645{
593 /* the union code still generates code under pressure in gcc, */ 646 /* the union code still generates code under pressure in gcc, */
594 /* but less than using pointers, and always seems to */ 647 /* but less than using pointers, and always seems to */
595 /* successfully return a constant. */ 648 /* successfully return a constant. */
596 /* the reason why we have this horrible preprocessor mess */ 649 /* the reason why we have this horrible preprocessor mess */
597 /* is to avoid it in all cases, at least on common architectures */ 650 /* is to avoid it in all cases, at least on common architectures */
598 /* or when using a recent enough gcc version (>= 4.6) */ 651 /* or when using a recent enough gcc version (>= 4.6) */
599#if ((__i386 || __i386__) && !__VOS__) || _M_IX86 || ECB_GCC_AMD64 || ECB_MSVC_AMD64
600 return 0x44;
601#elif __BYTE_ORDER__ && __BYTE_ORDER__ == __ORDER_LITTLE_ENDIAN__ 652#if (defined __BYTE_ORDER__ && __BYTE_ORDER__ == __ORDER_LITTLE_ENDIAN__) \
653 || ((__i386 || __i386__ || _M_IX86 || ECB_GCC_AMD64 || ECB_MSVC_AMD64) && !__VOS__)
654 #define ECB_LITTLE_ENDIAN 1
602 return 0x44; 655 return 0x44332211;
603#elif __BYTE_ORDER__ && __BYTE_ORDER__ == __ORDER_BIG_ENDIAN__ 656#elif (defined __BYTE_ORDER__ && __BYTE_ORDER__ == __ORDER_BIG_ENDIAN__) \
657 || ((__AARCH64EB__ || __MIPSEB__ || __ARMEB__) && !__VOS__)
658 #define ECB_BIG_ENDIAN 1
604 return 0x11; 659 return 0x11223344;
605#else 660#else
606 union 661 union
607 { 662 {
663 uint8_t c[4];
608 uint32_t i; 664 uint32_t u;
609 uint8_t c;
610 } u = { 0x11223344 }; 665 } u = { 0x11, 0x22, 0x33, 0x44 };
611 return u.c; 666 return u.u;
612#endif 667#endif
613} 668}
614 669
615ecb_inline ecb_const ecb_bool ecb_big_endian (void); 670ecb_inline ecb_const ecb_bool ecb_big_endian (void);
616ecb_inline ecb_const ecb_bool ecb_big_endian (void) { return ecb_byteorder_helper () == 0x11; } 671ecb_inline ecb_const ecb_bool ecb_big_endian (void) { return ecb_byteorder_helper () == 0x11223344; }
617ecb_inline ecb_const ecb_bool ecb_little_endian (void); 672ecb_inline ecb_const ecb_bool ecb_little_endian (void);
618ecb_inline ecb_const ecb_bool ecb_little_endian (void) { return ecb_byteorder_helper () == 0x44; } 673ecb_inline ecb_const ecb_bool ecb_little_endian (void) { return ecb_byteorder_helper () == 0x44332211; }
619 674
620#if ECB_GCC_VERSION(3,0) || ECB_C99 675#if ECB_GCC_VERSION(3,0) || ECB_C99
621 #define ecb_mod(m,n) ((m) % (n) + ((m) % (n) < 0 ? (n) : 0)) 676 #define ecb_mod(m,n) ((m) % (n) + ((m) % (n) < 0 ? (n) : 0))
622#else 677#else
623 #define ecb_mod(m,n) ((m) < 0 ? ((n) - 1 - ((-1 - (m)) % (n))) : ((m) % (n))) 678 #define ecb_mod(m,n) ((m) < 0 ? ((n) - 1 - ((-1 - (m)) % (n))) : ((m) % (n)))
… …
647 return N; 702 return N;
648 } 703 }
649#else 704#else
650 #define ecb_array_length(name) (sizeof (name) / sizeof (name [0])) 705 #define ecb_array_length(name) (sizeof (name) / sizeof (name [0]))
651#endif 706#endif
707
708ecb_function_ ecb_const uint32_t ecb_binary16_to_binary32 (uint32_t x);
709ecb_function_ ecb_const uint32_t
710ecb_binary16_to_binary32 (uint32_t x)
711{
712 unsigned int s = (x & 0x8000) << (31 - 15);
713 int e = (x >> 10) & 0x001f;
714 unsigned int m = x & 0x03ff;
715
716 if (ecb_expect_false (e == 31))
717 /* infinity or NaN */
718 e = 255 - (127 - 15);
719 else if (ecb_expect_false (!e))
720 {
721 if (ecb_expect_true (!m))
722 /* zero, handled by code below by forcing e to 0 */
723 e = 0 - (127 - 15);
724 else
725 {
726 /* subnormal, renormalise */
727 unsigned int s = 10 - ecb_ld32 (m);
728
729 m = (m << s) & 0x3ff; /* mask implicit bit */
730 e -= s - 1;
731 }
732 }
733
734 /* e and m now are normalised, or zero, (or inf or nan) */
735 e += 127 - 15;
736
737 return s | (e << 23) | (m << (23 - 10));
738}
739
740ecb_function_ ecb_const uint16_t ecb_binary32_to_binary16 (uint32_t x);
741ecb_function_ ecb_const uint16_t
742ecb_binary32_to_binary16 (uint32_t x)
743{
744 unsigned int s = (x >> 16) & 0x00008000; /* sign bit, the easy part */
745 unsigned int e = ((x >> 23) & 0x000000ff) - (127 - 15); /* the desired exponent */
746 unsigned int m = x & 0x007fffff;
747
748 x &= 0x7fffffff;
749
750 /* if it's within range of binary16 normals, use fast path */
751 if (ecb_expect_true (0x38800000 <= x && x <= 0x477fefff))
752 {
753 /* mantissa round-to-even */
754 m += 0x00000fff + ((m >> (23 - 10)) & 1);
755
756 /* handle overflow */
757 if (ecb_expect_false (m >= 0x00800000))
758 {
759 m >>= 1;
760 e += 1;
761 }
762
763 return s | (e << 10) | (m >> (23 - 10));
764 }
765
766 /* handle large numbers and infinity */
767 if (ecb_expect_true (0x477fefff < x && x <= 0x7f800000))
768 return s | 0x7c00;
769
770 /* handle zero, subnormals and small numbers */
771 if (ecb_expect_true (x < 0x38800000))
772 {
773 /* zero */
774 if (ecb_expect_true (!x))
775 return s;
776
777 /* handle subnormals */
778
779 /* too small, will be zero */
780 if (e < (14 - 24)) /* might not be sharp, but is good enough */
781 return s;
782
783 m |= 0x00800000; /* make implicit bit explicit */
784
785 /* very tricky - we need to round to the nearest e (+10) bit value */
786 {
787 unsigned int bits = 14 - e;
788 unsigned int half = (1 << (bits - 1)) - 1;
789 unsigned int even = (m >> bits) & 1;
790
791 /* if this overflows, we will end up with a normalised number */
792 m = (m + half + even) >> bits;
793 }
794
795 return s | m;
796 }
797
798 /* handle NaNs, preserve leftmost nan bits, but make sure we don't turn them into infinities */
799 m >>= 13;
800
801 return s | 0x7c00 | m | !m;
802}
652 803
653/*******************************************************************************/ 804/*******************************************************************************/
654/* floating point stuff, can be disabled by defining ECB_NO_LIBM */ 805/* floating point stuff, can be disabled by defining ECB_NO_LIBM */
655 806
656/* basically, everything uses "ieee pure-endian" floating point numbers */ 807/* basically, everything uses "ieee pure-endian" floating point numbers */
… …
693 #define ECB_NAN ECB_INFINITY 844 #define ECB_NAN ECB_INFINITY
694 #endif 845 #endif
695 846
696 #if ECB_C99 || _XOPEN_VERSION >= 600 || _POSIX_VERSION >= 200112L 847 #if ECB_C99 || _XOPEN_VERSION >= 600 || _POSIX_VERSION >= 200112L
697 #define ecb_ldexpf(x,e) ldexpf ((x), (e)) 848 #define ecb_ldexpf(x,e) ldexpf ((x), (e))
849 #define ecb_frexpf(x,e) frexpf ((x), (e))
698 #else 850 #else
699 #define ecb_ldexpf(x,e) (float) ldexp ((double) (x), (e)) 851 #define ecb_ldexpf(x,e) (float) ldexp ((double) (x), (e))
852 #define ecb_frexpf(x,e) (float) frexp ((double) (x), (e))
700 #endif 853 #endif
701
702 /* converts an ieee half/binary16 to a float */
703 ecb_function_ ecb_const float ecb_binary16_to_float (uint16_t x);
704 ecb_function_ ecb_const float
705 ecb_binary16_to_float (uint16_t x)
706 {
707 int e = (x >> 10) & 0x1f;
708 int m = x & 0x3ff;
709 float r;
710
711 if (!e ) r = ecb_ldexpf (m , -24);
712 else if (e != 31) r = ecb_ldexpf (m + 0x400, e - 25);
713 else if (m ) r = ECB_NAN;
714 else r = ECB_INFINITY;
715
716 return x & 0x8000 ? -r : r;
717 }
718 854
719 /* convert a float to ieee single/binary32 */ 855 /* convert a float to ieee single/binary32 */
720 ecb_function_ ecb_const uint32_t ecb_float_to_binary32 (float x); 856 ecb_function_ ecb_const uint32_t ecb_float_to_binary32 (float x);
721 ecb_function_ ecb_const uint32_t 857 ecb_function_ ecb_const uint32_t
722 ecb_float_to_binary32 (float x) 858 ecb_float_to_binary32 (float x)
… …
733 if (x == 0e0f ) return 0x00000000U; 869 if (x == 0e0f ) return 0x00000000U;
734 if (x > +3.40282346638528860e+38f) return 0x7f800000U; 870 if (x > +3.40282346638528860e+38f) return 0x7f800000U;
735 if (x < -3.40282346638528860e+38f) return 0xff800000U; 871 if (x < -3.40282346638528860e+38f) return 0xff800000U;
736 if (x != x ) return 0x7fbfffffU; 872 if (x != x ) return 0x7fbfffffU;
737 873
738 m = frexpf (x, &e) * 0x1000000U; 874 m = ecb_frexpf (x, &e) * 0x1000000U;
739 875
740 r = m & 0x80000000U; 876 r = m & 0x80000000U;
741 877
742 if (r) 878 if (r)
743 m = -m; 879 m = -m;
… …
854 #endif 990 #endif
855 991
856 return r; 992 return r;
857 } 993 }
858 994
859#endif 995 /* convert a float to ieee half/binary16 */
996 ecb_function_ ecb_const uint16_t ecb_float_to_binary16 (float x);
997 ecb_function_ ecb_const uint16_t
998 ecb_float_to_binary16 (float x)
999 {
1000 return ecb_binary32_to_binary16 (ecb_float_to_binary32 (x));
1001 }
860 1002
861#endif 1003 /* convert an ieee half/binary16 to float */
1004 ecb_function_ ecb_const float ecb_binary16_to_float (uint16_t x);
1005 ecb_function_ ecb_const float
1006 ecb_binary16_to_float (uint16_t x)
1007 {
1008 return ecb_binary32_to_float (ecb_binary16_to_binary32 (x));
1009 }
862 1010
1011#endif
1012
1013#endif
1014

Diff Legend

– Removed lines
+ Added lines
< Changed lines
> Changed lines