| 1 | /* |
1 | /* |
| 2 | * libecb - http://software.schmorp.de/pkg/libecb |
2 | * libecb - http://software.schmorp.de/pkg/libecb |
| 3 | * |
3 | * |
| 4 | * Copyright (©) 2009-2012 Marc Alexander Lehmann <libecb@schmorp.de> |
4 | * Copyright (©) 2009-2013 Marc Alexander Lehmann <libecb@schmorp.de> |
| 5 | * Copyright (©) 2011 Emanuele Giaquinta |
5 | * Copyright (©) 2011 Emanuele Giaquinta |
| 6 | * All rights reserved. |
6 | * All rights reserved. |
| 7 | * |
7 | * |
| 8 | * Redistribution and use in source and binary forms, with or without modifica- |
8 | * Redistribution and use in source and binary forms, with or without modifica- |
| 9 | * tion, are permitted provided that the following conditions are met: |
9 | * tion, are permitted provided that the following conditions are met: |
| … | |
… | |
| 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 = ECB_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) |
