fpmath.h 3.2 KB

123456789101112131415161718192021222324252627282930313233343536373839404142434445464748495051525354555657585960616263646566676869707172737475767778798081828384858687888990919293949596979899100101102103104105106107108109110111
  1. /*-
  2. * Copyright (c) 2003 Mike Barcroft <mike@FreeBSD.org>
  3. * Copyright (c) 2002 David Schultz <das@FreeBSD.ORG>
  4. * All rights reserved.
  5. *
  6. * Redistribution and use in source and binary forms, with or without
  7. * modification, are permitted provided that the following conditions
  8. * are met:
  9. * 1. Redistributions of source code must retain the above copyright
  10. * notice, this list of conditions and the following disclaimer.
  11. * 2. Redistributions in binary form must reproduce the above copyright
  12. * notice, this list of conditions and the following disclaimer in the
  13. * documentation and/or other materials provided with the distribution.
  14. *
  15. * THIS SOFTWARE IS PROVIDED BY THE AUTHOR AND CONTRIBUTORS ``AS IS'' AND
  16. * ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE
  17. * IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE
  18. * ARE DISCLAIMED. IN NO EVENT SHALL THE AUTHOR OR CONTRIBUTORS BE LIABLE
  19. * FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR CONSEQUENTIAL
  20. * DAMAGES (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS
  21. * OR SERVICES; LOSS OF USE, DATA, OR PROFITS; OR BUSINESS INTERRUPTION)
  22. * HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT
  23. * LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY
  24. * OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF
  25. * SUCH DAMAGE.
  26. *
  27. * $FreeBSD: src/lib/libc/include/fpmath.h,v 1.4 2008/12/23 22:20:59 marcel Exp $
  28. */
  29. // Currently assumes Intel platform
  30. #if defined (__i386__) || defined(__x86_64__)
  31. #ifdef __LP64__
  32. #include "amd64_fpmath.h"
  33. #else
  34. #include "i386_fpmath.h"
  35. #endif
  36. #endif
  37. #ifdef __linux
  38. #include <features.h>
  39. #include <endian.h>
  40. #define _LITTLE_ENDIAN __LITTLE_ENDIAN
  41. #define _BIG_ENDIAN __BIG_ENDIAN
  42. #define _PDP_ENDIAN __PDP_ENDIAN
  43. #define _BYTE_ORDER __BYTE_ORDER
  44. #endif
  45. #ifdef __APPLE__
  46. #include <machine/endian.h>
  47. #define _LITTLE_ENDIAN LITTLE_ENDIAN
  48. #define _BIG_ENDIAN BIG_ENDIAN
  49. #define _PDP_ENDIAN PDP_ENDIAN
  50. #define _BYTE_ORDER BYTE_ORDER
  51. #endif
  52. #ifdef WIN32
  53. #define _LITTLE_ENDIAN 1234
  54. #define _BIG_ENDIAN 4321
  55. #define _PDP_ENDIAN 3412
  56. #define _BYTE_ORDER __LITTLE_ENDIAN
  57. #define _FLOAT_WORD_ORDER __LITTLE_ENDIAN
  58. #define LITTLE_ENDIAN __LITTLE_ENDIAN
  59. #define BIG_ENDIAN __BIG_ENDIAN
  60. #define PDP_ENDIAN __PDP_ENDIAN
  61. #define BYTE_ORDER __BYTE_ORDER
  62. #endif
  63. #ifndef _IEEE_WORD_ORDER
  64. #define _IEEE_WORD_ORDER _BYTE_ORDER
  65. #endif
  66. union IEEEf2bits {
  67. float f;
  68. struct {
  69. #if _BYTE_ORDER == _LITTLE_ENDIAN
  70. unsigned int man :23;
  71. unsigned int exp :8;
  72. unsigned int sign :1;
  73. #else /* _BIG_ENDIAN */
  74. unsigned int sign :1;
  75. unsigned int exp :8;
  76. unsigned int man :23;
  77. #endif
  78. } bits;
  79. };
  80. #define DBL_MANH_SIZE 20
  81. #define DBL_MANL_SIZE 32
  82. union IEEEd2bits {
  83. double d;
  84. struct {
  85. #if _BYTE_ORDER == _LITTLE_ENDIAN
  86. #if _IEEE_WORD_ORDER == _LITTLE_ENDIAN
  87. unsigned int manl :32;
  88. #endif
  89. unsigned int manh :20;
  90. unsigned int exp :11;
  91. unsigned int sign :1;
  92. #if _IEEE_WORD_ORDER == _BIG_ENDIAN
  93. unsigned int manl :32;
  94. #endif
  95. #else /* _BIG_ENDIAN */
  96. unsigned int sign :1;
  97. unsigned int exp :11;
  98. unsigned int manh :20;
  99. unsigned int manl :32;
  100. #endif
  101. } bits;
  102. };