开通博客赚积分
发布资源赚积分
分类

源码开发语言/平台
当前位置：首页 > 源码/资料 > 多媒体 > DVD > 查看源码
fdlibm.c
资源名称：8202s.rar　[点击查看]
上传用户：poi891205
上传日期：2013-07-15
资源大小：9745k
文件大小：18k
源码类别：

DVD

开发平台：

C/C++

fdlibm.c：源码内容
  1. # 1 "fdlibm.h"
  3. #ifndef __IEEE_BIG_ENDIAN
  4. #ifndef __IEEE_LITTLE_ENDIAN
  6. #if defined(__arm__) || defined(__thumb__)
  7. /* ARM always has big-endian words.  Within those words the byte ordering
  8.    appears to be big or little endian.  Newlib doesn't seem to care about
  9.    the byte ordering within words.  */
  10. #define __IEEE_BIG_ENDIAN
  11. #endif
  13. #ifdef __hppa__
  14. #define __IEEE_BIG_ENDIAN
  15. #endif
  17. #ifdef __sparc__
  18. #define __IEEE_BIG_ENDIAN
  19. #endif
  21. //#if defined(__m68k__) || defined(__mc68000__)
  22. //#define __IEEE_BIG_ENDIAN
  23. //#endif
  25. #if defined (__H8300__) || defined (__H8300H__) || defined (__H8300S__)
  26. #define __IEEE_BIG_ENDIAN
  27. #define __SMALL_BITFIELDS
  28. #define _DOUBLE_IS_32BITS
  29. #endif
  31. #ifdef __H8500__
  32. #define __IEEE_BIG_ENDIAN
  33. #define __SMALL_BITFIELDS
  34. #define _DOUBLE_IS_32BITS
  35. #endif
  37. #ifdef __sh__
  38. #ifdef __LITTLE_ENDIAN__
  39. #define __IEEE_LITTLE_ENDIAN
  40. #else
  41. #define __IEEE_BIG_ENDIAN
  42. #endif
  43. #if defined(__SH3E__) || defined(__SH4_SINGLE_ONLY__)
  44. #define _DOUBLE_IS_32BITS
  45. #endif
  46. #endif
  48. #ifdef _AM29K
  49. #define __IEEE_BIG_ENDIAN
  50. #endif
  52. #ifdef __i386__
  53. #define __IEEE_LITTLE_ENDIAN
  54. #endif
  56. #ifdef __i960__
  57. #define __IEEE_LITTLE_ENDIAN
  58. #endif
  60. #ifdef __M32R__
  61. #define __IEEE_BIG_ENDIAN
  62. #endif
  65. #ifdef __MIPSEL__
  66. #define __IEEE_LITTLE_ENDIAN
  67. #endif
  68. #ifdef __MIPSEB__
  69. #define __IEEE_BIG_ENDIAN
  70. #endif
  72. /* necv70 was __IEEE_LITTLE_ENDIAN. */
  74. #ifdef __W65__
  75. #define __IEEE_LITTLE_ENDIAN
  76. #define __SMALL_BITFIELDS
  77. #define _DOUBLE_IS_32BITS
  78. #endif
  80. #if defined(__Z8001__) || defined(__Z8002__)
  81. #define __IEEE_BIG_ENDIAN
  82. #endif
  84. #ifdef __m88k__
  85. #define __IEEE_BIG_ENDIAN
  86. #endif
  88. #ifdef __mn10300__
  89. #define __IEEE_LITTLE_ENDIAN
  90. #endif
  92. #ifdef __mn10200__
  93. #define __IEEE_LITTLE_ENDIAN
  94. #define __SMALL_BITFIELDS
  95. #define _DOUBLE_IS_32BITS
  96. #endif
  98. #ifdef __v800
  99. #define __IEEE_LITTLE_ENDIAN
  100. #endif
  102. #ifdef __v850
  103. #define __IEEE_LITTLE_ENDIAN
  104. #endif
  106. #ifdef __D10V__
  107. #define __IEEE_BIG_ENDIAN
  108. #define _DOUBLE_IS_32BITS
  109. #define __SMALL_BITFIELDS
  110. #endif
  112. #ifdef __PPC__
  113. #if (defined(_BIG_ENDIAN) && _BIG_ENDIAN) || (defined(_AIX) && _AIX)
  114. #define __IEEE_BIG_ENDIAN
  115. #else
  116. #if (defined(_LITTLE_ENDIAN) && _LITTLE_ENDIAN) || (defined(__sun__) && __sun__) || (defined(_WIN32) && _WIN32)
  117. #define __IEEE_LITTLE_ENDIAN
  118. #endif
  119. #endif
  120. #endif
  122. #ifdef __arc__
  123. #ifdef __big_endian__
  124. #define __IEEE_BIG_ENDIAN
  125. #else
  126. #define __IEEE_LITTLE_ENDIAN
  127. #endif
  128. #endif
  130. #ifndef __IEEE_BIG_ENDIAN
  131. #ifndef __IEEE_LITTLE_ENDIAN
  132. #error Endianess not declared!!
  133. #endif /* not __IEEE_LITTLE_ENDIAN */
  134. #endif /* not __IEEE_BIG_ENDIAN */
  136. #endif /* not __IEEE_LITTLE_ENDIAN */
  137. #endif /* not __IEEE_BIG_ENDIAN */
  140.  
  141. # 1 "/home/potatooo/fcc/lib/gcc-lib/mips-ecoff/2.8.1/include/math.h" 1 3
  142.  
  150. # 1 "/home/potatooo/fcc/mips-ecoff/include/sys/reent.h" 1 3
  151.  
  153.  
  163. # 1 "/home/potatooo/fcc/lib/gcc-lib/mips-ecoff/2.8.1/include/_ansi.h" 1 3
  164.  
  166.  
  178. # 1 "/home/potatooo/fcc/mips-ecoff/include/sys/config.h" 1 3
  182.  
  183.  
  184. # 14 "/home/potatooo/fcc/mips-ecoff/include/sys/config.h" 3
  187. # 25 "/home/potatooo/fcc/mips-ecoff/include/sys/config.h" 3
  190.  
  199. # 44 "/home/potatooo/fcc/mips-ecoff/include/sys/config.h" 3
  239.  
  248. typedef int __int32_t;
  249. typedef unsigned int __uint32_t;
  257. # 15 "/home/potatooo/fcc/lib/gcc-lib/mips-ecoff/2.8.1/include/_ansi.h" 2 3
  260.  
  261.  
  287. # 61 "/home/potatooo/fcc/lib/gcc-lib/mips-ecoff/2.8.1/include/_ansi.h" 3
  290.  
  299. # 13 "/home/potatooo/fcc/mips-ecoff/include/sys/reent.h" 2 3
  302. struct _glue 
  303. {
  304.   struct _glue *_next;
  305.   int _niobs;
  306.   struct __sFILE *_iobs;
  307. };
  309. struct _Bigint 
  310. {
  311.   struct _Bigint *_next;
  312.   int _k, _maxwds, _sign, _wds;
  313.   unsigned long _x[1];
  314. };
  316.  
  322. struct _atexit {
  323. struct _atexit *_next;  
  324. int _ind;  
  325. void (*_fns[32 ])();  
  326. };
  328.  
  335. struct __sbuf {
  336. unsigned char *_base;
  337. int _size;
  338. };
  340.  
  345. typedef long _fpos_t;  
  346.  
  348.  
  373. struct __sFILE {
  374.   unsigned char *_p;  
  375.   int _r;  
  376.   int _w;  
  377.   short _flags;  
  378.   short _file;  
  379.   struct __sbuf _bf;  
  380.   int _lbfsize;  
  382.    
  383.   void *  _cookie;  
  385.   int  (*_read)   (void *  _cookie, char *_buf, int _n)  ;
  386.   int  (*_write)   (void *  _cookie, const char *_buf, int _n)  ;
  387.   _fpos_t  (*_seek)   (void *  _cookie, _fpos_t _offset, int _whence)  ;
  388.   int  (*_close)   (void *  _cookie)  ;
  390.    
  391.   struct __sbuf _ub;  
  392.   unsigned char *_up;  
  393.   int _ur;  
  395.    
  396.   unsigned char _ubuf[3];  
  397.   unsigned char _nbuf[1];  
  399.    
  400.   struct __sbuf _lb;  
  402.    
  403.   int _blksize;  
  404.   int _offset;  
  406.   struct _reent *_data;
  407. };
  409.  
  417. struct _reent
  418. {
  419.    
  420.   int _errno;
  422.    
  425.   struct __sFILE *_stdin, *_stdout, *_stderr;
  427.   int  _inc;  
  428.   char _emergency[25];
  429.  
  430.   int _current_category;  
  431.   const  char *_current_locale;
  433.   int __sdidinit;  
  435.   void  (*__cleanup)   (struct _reent *)  ;
  437.    
  438.   struct _Bigint *_result;
  439.   int _result_k;
  440.   struct _Bigint *_p5s;
  441.   struct _Bigint **_freelist;
  443.    
  444.   int _cvtlen;  
  445.   char *_cvtbuf;
  447.    
  451.   unsigned char * _nextf[30 ];
  452.   unsigned int _nmalloc[30 ];
  454.    
  455.   struct _atexit *_atexit;  
  456.   struct _atexit _atexit0;  
  458.    
  459.   void (**(_sig_func))();
  461.    
  464.   struct _glue __sglue;  
  465.   struct __sFILE __sf[3];  
  466. };
  471.  
  480. extern struct _reent *_impure_ptr  ;
  482. void _reclaim_reent  (struct _reent *)  ;
  484.  
  494. # 9 "/home/potatooo/fcc/lib/gcc-lib/mips-ecoff/2.8.1/include/math.h" 2 3
  496. # 1 "/home/potatooo/fcc/mips-ecoff/include/machine/ieeefp.h" 1 3
  531. # 44 "/home/potatooo/fcc/mips-ecoff/include/machine/ieeefp.h" 3
  558.  
  598. # 118 "/home/potatooo/fcc/mips-ecoff/include/machine/ieeefp.h" 3
  618. # 10 "/home/potatooo/fcc/lib/gcc-lib/mips-ecoff/2.8.1/include/math.h" 2 3
  624.  
  627. union __dmath
  628. {
  629.   __uint32_t i[2];
  630.   double d;
  631. };
  633. extern const union __dmath __infinity;
  641.  
  644. extern double atan  (double)  ;
  645. extern double cos  (double)  ;
  646. extern double sin  (double)  ;
  647. extern double tan  (double)  ;
  648. extern double tanh  (double)  ;
  649. extern double frexp  (double, int *)  ;
  650. extern double modf  (double, double *)  ;
  651. extern double ceil  (double)  ;
  652. extern double fabs  (double)  ;
  653. extern double floor  (double)  ;
  656.  
  660. extern double acos  (double)  ;
  661. extern double asin  (double)  ;
  662. extern double atan2  (double, double)  ;
  663. extern double cosh  (double)  ;
  664. extern double sinh  (double)  ;
  665. extern double exp  (double)  ;
  666. extern double ldexp  (double, int)  ;
  667. extern double log  (double)  ;
  668. extern double log10  (double)  ;
  669. extern double pow  (double, double)  ;
  670. extern double sqrt  (double)  ;
  671. extern double fmod  (double, double)  ;
  677.  
  679. extern double infinity  (void)  ;
  680. extern double nan  (void)  ;
  681. extern int isnan  (double)  ;
  682. extern int isinf  (double)  ;
  683. extern int finite  (double)  ;
  684. extern double copysign  (double, double)  ;
  685. extern int ilogb  (double)  ;
  687. extern double asinh  (double)  ;
  688. extern double cbrt  (double)  ;
  689. extern double nextafter  (double, double)  ;
  690. extern double rint  (double)  ;
  691. extern double scalbn  (double, int)  ;
  694. extern double log1p  (double)  ;
  695. extern double expm1  (double)  ;
  699. extern double acosh  (double)  ;
  700. extern double atanh  (double)  ;
  701. extern double remainder  (double, double)  ;
  702. extern double gamma  (double)  ;
  703. extern double gamma_r  (double, int *)  ;
  704. extern double lgamma  (double)  ;
  705. extern double lgamma_r  (double, int *)  ;
  706. extern double erf  (double)  ;
  707. extern double erfc  (double)  ;
  708. extern double y0  (double)  ;
  709. extern double y1  (double)  ;
  710. extern double yn  (int, double)  ;
  711. extern double j0  (double)  ;
  712. extern double j1  (double)  ;
  713. extern double jn  (int, double)  ;
  717. extern double hypot  (double, double)  ;
  720. extern double cabs();
  721. extern double drem  (double, double)  ;
  725.  
  727. extern float atanf  (float)  ;
  728. extern float cosf  (float)  ;
  729. extern float sinf  (float)  ;
  730. extern float tanf  (float)  ;
  731. extern float tanhf  (float)  ;
  732. extern float frexpf  (float, int *)  ;
  733. extern float modff  (float, float *)  ;
  734. extern float ceilf  (float)  ;
  735. extern float fabsf  (float)  ;
  736. extern float floorf  (float)  ;
  739. extern float acosf  (float)  ;
  740. extern float asinf  (float)  ;
  741. extern float atan2f  (float, float)  ;
  742. extern float coshf  (float)  ;
  743. extern float sinhf  (float)  ;
  744. extern float expf  (float)  ;
  745. extern float ldexpf  (float, int)  ;
  746. extern float logf  (float)  ;
  747. extern float log10f  (float)  ;
  748. extern float powf  (float, float)  ;
  749. extern float sqrtf  (float)  ;
  750. extern float fmodf  (float, float)  ;
  753.  
  755. extern float infinityf  (void)  ;
  756. extern float nanf  (void)  ;
  757. extern int isnanf  (float)  ;
  758. extern int isinff  (float)  ;
  759. extern int finitef  (float)  ;
  760. extern float copysignf  (float, float)  ;
  761. extern int ilogbf  (float)  ;
  763. extern float asinhf  (float)  ;
  764. extern float cbrtf  (float)  ;
  765. extern float nextafterf  (float, float)  ;
  766. extern float rintf  (float)  ;
  767. extern float scalbnf  (float, int)  ;
  768. extern float log1pf  (float)  ;
  769. extern float expm1f  (float)  ;
  772. extern float acoshf  (float)  ;
  773. extern float atanhf  (float)  ;
  774. extern float remainderf  (float, float)  ;
  775. extern float gammaf  (float)  ;
  776. extern float gammaf_r  (float, int *)  ;
  777. extern float lgammaf  (float)  ;
  778. extern float lgammaf_r  (float, int *)  ;
  779. extern float erff  (float)  ;
  780. extern float erfcf  (float)  ;
  781. extern float y0f  (float)  ;
  782. extern float y1f  (float)  ;
  783. extern float ynf  (int, float)  ;
  784. extern float j0f  (float)  ;
  785. extern float j1f  (float)  ;
  786. extern float jnf  (int, float)  ;
  788. extern float hypotf  (float, float)  ;
  790. extern float cabsf();
  791. extern float dremf  (float, float)  ;
  795.  
  797. extern int signgam;
  802. struct exception;
  809.  
  820. struct exception 
  825. {
  826.   int type;
  827.   char *name;
  828.   double arg1;
  829.   double arg2;
  830.   double retval;
  831.   int err;
  832. };
  834. # 233 "/home/potatooo/fcc/lib/gcc-lib/mips-ecoff/2.8.1/include/math.h" 3
  842. extern int matherr  (struct exception *e)  ;
  854.  
  866.  
  891.  
  893. enum __fdlibm_version
  894. {
  895.   __fdlibm_ieee = -1,
  896.   __fdlibm_svid,
  897.   __fdlibm_xopen,
  898.   __fdlibm_posix
  899. };
  904. extern enum __fdlibm_version  __fdlib_version ;
  917. # 15 "fdlibm.h" 2
  919. # 1 "/home/potatooo/fcc/mips-ecoff/include/machine/ieeefp.h" 1 3
  954. # 44 "/home/potatooo/fcc/mips-ecoff/include/machine/ieeefp.h" 3
  981.  
  1021. # 118 "/home/potatooo/fcc/mips-ecoff/include/machine/ieeefp.h" 3
  1041. # 16 "fdlibm.h" 2
  1044.  
  1055.  
  1062.  
  1064. extern double logb  (double)  ;
  1068. extern double scalb  (double, double)  ;
  1070. extern double significand  (double)  ;
  1072.  
  1073. extern double __ieee754_sqrt  (double)  ;
  1074. extern double __ieee754_acos  (double)  ;
  1075. extern double __ieee754_acosh  (double)  ;
  1076. extern double __ieee754_log  (double)  ;
  1077. extern double __ieee754_atanh  (double)  ;
  1078. extern double __ieee754_asin  (double)  ;
  1079. extern double __ieee754_atan2  (double,double)  ;
  1080. extern double __ieee754_exp  (double)  ;
  1081. extern double __ieee754_cosh  (double)  ;
  1082. extern double __ieee754_fmod  (double,double)  ;
  1083. extern double __ieee754_pow  (double,double)  ;
  1084. extern double __ieee754_lgamma_r  (double,int *)  ;
  1085. extern double __ieee754_gamma_r  (double,int *)  ;
  1086. extern double __ieee754_log10  (double)  ;
  1087. extern double __ieee754_sinh  (double)  ;
  1088. extern double __ieee754_hypot  (double,double)  ;
  1089. extern double __ieee754_j0  (double)  ;
  1090. extern double __ieee754_j1  (double)  ;
  1091. extern double __ieee754_y0  (double)  ;
  1092. extern double __ieee754_y1  (double)  ;
  1093. extern double __ieee754_jn  (int,double)  ;
  1094. extern double __ieee754_yn  (int,double)  ;
  1095. extern double __ieee754_remainder  (double,double)  ;
  1096. extern __int32_t __ieee754_rem_pio2  (double,double*)  ;
  1100. extern double __ieee754_scalb  (double,double)  ;
  1103.  
  1104. extern double __kernel_standard  (double,double,int)  ;
  1105. extern double __kernel_sin  (double,double,int)  ;
  1106. extern double __kernel_cos  (double,double)  ;
  1107. extern double __kernel_tan  (double,double,int)  ;
  1108. extern int    __kernel_rem_pio2  (double*,double*,int,int,int,const __int32_t*)  ;
  1110.  
  1111. extern float logbf  (float)  ;
  1115. extern float scalbf  (float, float)  ;
  1117. extern float significandf  (float)  ;
  1119.  
  1120. extern float __ieee754_sqrtf  (float)  ;
  1121. extern float __ieee754_acosf  (float)  ;
  1122. extern float __ieee754_acoshf  (float)  ;
  1123. extern float __ieee754_logf  (float)  ;
  1124. extern float __ieee754_atanhf  (float)  ;
  1125. extern float __ieee754_asinf  (float)  ;
  1126. extern float __ieee754_atan2f  (float,float)  ;
  1127. extern float __ieee754_expf  (float)  ;
  1128. extern float __ieee754_coshf  (float)  ;
  1129. extern float __ieee754_fmodf  (float,float)  ;
  1130. extern float __ieee754_powf  (float,float)  ;
  1131. extern float __ieee754_lgammaf_r  (float,int *)  ;
  1132. extern float __ieee754_gammaf_r  (float,int *)  ;
  1133. extern float __ieee754_log10f  (float)  ;
  1134. extern float __ieee754_sinhf  (float)  ;
  1135. extern float __ieee754_hypotf  (float,float)  ;
  1136. extern float __ieee754_j0f  (float)  ;
  1137. extern float __ieee754_j1f  (float)  ;
  1138. extern float __ieee754_y0f  (float)  ;
  1139. extern float __ieee754_y1f  (float)  ;
  1140. extern float __ieee754_jnf  (int,float)  ;
  1141. extern float __ieee754_ynf  (int,float)  ;
  1142. extern float __ieee754_remainderf  (float,float)  ;
  1143. extern __int32_t __ieee754_rem_pio2f  (float,float*)  ;
  1147. extern float __ieee754_scalbf  (float,float)  ;
  1150.  
  1151. extern float __kernel_sinf  (float,float,int)  ;
  1152. extern float __kernel_cosf  (float,float)  ;
  1153. extern float __kernel_tanf  (float,float,int)  ;
  1154. extern int   __kernel_rem_pio2f  (float*,float*,int,int,int,const __int32_t*)  ;
  1156.  
  1160. #ifdef __IEEE_BIG_ENDIAN
  1162. typedef union
  1163. {
  1164.   double value;
  1165.   struct
  1166.   {
  1167.     __uint32_t msw;
  1168.     __uint32_t lsw;
  1169.   } parts;
  1170. } ieee_double_shape_type;
  1172. #endif
  1174. #ifdef __IEEE_LITTLE_ENDIAN
  1176. typedef union
  1177. {
  1178.   double value;
  1179.   struct
  1180.   {
  1181.     __uint32_t lsw;
  1182.     __uint32_t msw;
  1183.   } parts;
  1184. } ieee_double_shape_type;
  1186. #endif
  1188. /* Get two 32 bit ints from a double.  */
  1190. #define EXTRACT_WORDS(ix0,ix1,d)                                
  1191. do {                                                            
  1192.   ieee_double_shape_type ew_u;                                  
  1193.   ew_u.value = (d);                                             
  1194.   (ix0) = ew_u.parts.msw;                                       
  1195.   (ix1) = ew_u.parts.lsw;                                       
  1196. } while (0)
  1198. /* Get the more significant 32 bit int from a double.  */
  1200. #define GET_HIGH_WORD(i,d)                                      
  1201. do {                                                            
  1202.   ieee_double_shape_type gh_u;                                  
  1203.   gh_u.value = (d);                                             
  1204.   (i) = gh_u.parts.msw;                                         
  1205. } while (0)
  1208. /* Get the less significant 32 bit int from a double.  */
  1210. #define GET_LOW_WORD(i,d)                                       
  1211. do {                                                            
  1212.   ieee_double_shape_type gl_u;                                  
  1213.   gl_u.value = (d);                                             
  1214.   (i) = gl_u.parts.lsw;                                         
  1215. } while (0)
  1217. /* Set a double from two 32 bit ints.  */
  1219. #define INSERT_WORDS(d,ix0,ix1)                                 
  1220. do {                                                            
  1221.   ieee_double_shape_type iw_u;                                  
  1222.   iw_u.parts.msw = (ix0);                                       
  1223.   iw_u.parts.lsw = (ix1);                                       
  1224.   (d) = iw_u.value;                                             
  1225. } while (0)
  1227. /* Set the more significant 32 bits of a double from an int.  */
  1229. #define SET_HIGH_WORD(d,v)                                      
  1230. do {                                                            
  1231.   ieee_double_shape_type sh_u;                                  
  1232.   sh_u.value = (d);                                             
  1233.   sh_u.parts.msw = (v);                                         
  1234.   (d) = sh_u.value;                                             
  1235. } while (0)
  1237. /* Set the less significant 32 bits of a double from an int.  */
  1239. #define SET_LOW_WORD(d,v)                                       
  1240. do {                                                            
  1241.   ieee_double_shape_type sl_u;                                  
  1242.   sl_u.value = (d);                                             
  1243.   sl_u.parts.lsw = (v);                                         
  1244.   (d) = sl_u.value;                                             
  1245. } while (0)
  1246.  
  1249. typedef union
  1250. {
  1251.   float value;
  1252.   __uint32_t word;
  1253. } ieee_float_shape_type;
  1255.  
  1256. #define GET_FLOAT_WORD(i,d)                                     
  1257. do {                                                            
  1258.   ieee_float_shape_type gf_u;                                   
  1259.   gf_u.value = (d);                                             
  1260.   (i) = gf_u.word;                                              
  1261. } while (0)
  1263. /* Set a float from a 32 bit int.  */
  1265. #define SET_FLOAT_WORD(d,i)                                     
  1266. do {                                                            
  1267.   ieee_float_shape_type sf_u;                                   
  1268.   sf_u.word = (i);                                              
  1269.   (d) = sf_u.value;                                             
  1270. } while (0)