… | |
… | |
581 | || defined _M_IX86 || defined _M_AMD64 || defined _M_IA64 |
581 | || defined _M_IX86 || defined _M_AMD64 || defined _M_IA64 |
582 | #define ECB_STDFP 1 |
582 | #define ECB_STDFP 1 |
583 | #include <string.h> /* for memcpy */ |
583 | #include <string.h> /* for memcpy */ |
584 | #else |
584 | #else |
585 | #define ECB_STDFP 0 |
585 | #define ECB_STDFP 0 |
586 | #include <math.h> /* for frexp*, ldexp* */ |
|
|
587 | #endif |
586 | #endif |
588 | |
587 | |
589 | #ifndef ECB_NO_LIBM |
588 | #ifndef ECB_NO_LIBM |
|
|
589 | |
|
|
590 | #include <math.h> /* for frexp*, ldexp*, INFINITY, NAN */ |
|
|
591 | |
|
|
592 | /* only the oldest of old doesn't have this one. solaris. */ |
|
|
593 | #ifdef INFINITY |
|
|
594 | #define ECB_INFINITY INFINITY |
|
|
595 | #else |
|
|
596 | #define ECB_INFINITY HUGE_VAL |
|
|
597 | #endif |
|
|
598 | |
|
|
599 | #ifdef NAN |
|
|
600 | #define ECB_NAN NAN |
|
|
601 | #else |
|
|
602 | #define ECB_NAN ECB_INFINITY |
|
|
603 | #endif |
|
|
604 | |
|
|
605 | /* converts an ieee half/binary16 to a float */ |
|
|
606 | ecb_function_ float ecb_binary16_to_float (uint16_t x) ecb_const; |
|
|
607 | ecb_function_ float |
|
|
608 | ecb_binary16_to_float (uint16_t x) |
|
|
609 | { |
|
|
610 | int e = (x >> 10) & 0x1f; |
|
|
611 | int m = x & 0x3ff; |
|
|
612 | float r; |
|
|
613 | |
|
|
614 | if (!e ) r = ldexpf (m , -24); |
|
|
615 | else if (e != 31) r = ldexpf (m + 0x400, e - 25); |
|
|
616 | else if (m ) r = ECB_NAN; |
|
|
617 | else r = INFINITY; |
|
|
618 | |
|
|
619 | return x & 0x8000 ? -r : r; |
|
|
620 | } |
590 | |
621 | |
591 | /* convert a float to ieee single/binary32 */ |
622 | /* convert a float to ieee single/binary32 */ |
592 | ecb_function_ uint32_t ecb_float_to_binary32 (float x) ecb_const; |
623 | ecb_function_ uint32_t ecb_float_to_binary32 (float x) ecb_const; |
593 | ecb_function_ uint32_t |
624 | ecb_function_ uint32_t |
594 | ecb_float_to_binary32 (float x) |
625 | ecb_float_to_binary32 (float x) |