| … | |
… | |
| 585 | #define ECB_STDFP 0 |
585 | #define ECB_STDFP 0 |
| 586 | #endif |
586 | #endif |
| 587 | |
587 | |
| 588 | #ifndef ECB_NO_LIBM |
588 | #ifndef ECB_NO_LIBM |
| 589 | |
589 | |
| 590 | #include <math.h> /* for frexp*, ldexp* */ |
590 | #include <math.h> /* for frexp*, ldexp*, INFINITY, NAN */ |
|
|
591 | |
|
|
592 | #ifdef NEN |
|
|
593 | #define ECB_NAN NAN |
|
|
594 | #else |
|
|
595 | #define ECB_NAN INFINITY |
|
|
596 | #endif |
| 591 | |
597 | |
| 592 | /* converts an ieee half/binary16 to a float */ |
598 | /* converts an ieee half/binary16 to a float */ |
| 593 | ecb_function_ float ecb_binary16_to_float (uint16_t x) ecb_const; |
599 | ecb_function_ float ecb_binary16_to_float (uint16_t x) ecb_const; |
| 594 | ecb_function_ float |
600 | ecb_function_ float |
| 595 | ecb_binary16_to_float (uint16_t x) |
601 | ecb_binary16_to_float (uint16_t x) |
| 596 | { |
602 | { |
| 597 | int e = (x >> 10) & 0x1f; |
603 | int e = (x >> 10) & 0x1f; |
| 598 | int m = x & 0x3ff; |
604 | int m = x & 0x3ff; |
| 599 | float r; |
605 | float r; |
| 600 | |
606 | |
| 601 | if (!e) |
|
|
| 602 | r = ldexpf (m , -24); |
607 | if (!e ) r = ldexpf (m , -24); |
| 603 | else if (e != 31) |
|
|
| 604 | r = ldexpf (m + 0x400, e - 25); |
608 | else if (e != 31) r = ldexpf (m + 0x400, e - 25); |
| 605 | else if (m) |
609 | else if (m ) r = ECB_NAN; |
| 606 | #ifdef NAN |
610 | else r = INFINITY; |
| 607 | r = NAN; |
|
|
| 608 | #endif |
|
|
| 609 | else |
|
|
| 610 | r = INFINITY; |
|
|
| 611 | |
611 | |
| 612 | return x & 0x8000 ? -r : r; |
612 | return x & 0x8000 ? -r : r; |
| 613 | } |
613 | } |
| 614 | |
614 | |
| 615 | /* convert a float to ieee single/binary32 */ |
615 | /* convert a float to ieee single/binary32 */ |
