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

源码开发语言/平台
当前位置：首页 > 源码/资料 > 嵌入式/单片机编程 > VxWorks > 查看源码
log__L.c
资源名称：vxwork_src.rar　[点击查看]
上传用户：nvosite88
上传日期：2007-01-17
资源大小：4983k
文件大小：5k
源码类别：

VxWorks

开发平台：

C/C++

log__L.c：源码内容
  1. /* log__L.c - math routine */
  3. /* Copyright 1992 Wind River Systems, Inc. */
  5. /*
  6. modification history
  7. --------------------
  8. 01a,08jul92,smb  documentation.
  9. */
  11. /*
  12. DESCRIPTION
  13. * Copyright (c) 1985 Regents of the University of California.
  14. * All rights reserved.
  15. *
  16. * Redistribution and use in source and binary forms are permitted
  17. * provided that the above copyright notice and this paragraph are
  18. * duplicated in all such forms and that any documentation,
  19. * advertising materials, and other materials related to such
  20. * distribution and use acknowledge that the software was developed
  21. * by the University of California, Berkeley.  The name of the
  22. * University may not be used to endorse or promote products derived
  23. * from this software without specific prior written permission.
  24. * THIS SOFTWARE IS PROVIDED ``AS IS'' AND WITHOUT ANY EXPRESS OR
  25. * IMPLIED WARRANTIES, INCLUDING, WITHOUT LIMITATION, THE IMPLIED
  26. * WARRANTIES OF MERCHANTIBILITY AND FITNESS FOR A PARTICULAR PURPOSE.
  27. *
  28. * All recipients should regard themselves as participants in an ongoing
  29. * research project and hence should feel obligated to report their
  30. * experiences (good or bad) with these elementary function codes, using
  31. * the sendbug(8) program, to the authors.
  32. *
  33. SEE ALSO: American National Standard X3.159-1989
  34. NOMANUAL
  35. */
  37. #include "vxWorks.h"
  38. #include "math.h"
  40. #if defined(vax)||defined(tahoe) /* VAX D format (56 bits) */
  41. #ifdef vax
  42. #define _0x(A,B) 0x/**/A/**/B
  43. #else /* vax */
  44. #define _0x(A,B) 0x/**/B/**/A
  45. #endif /* vax */
  46. /* static double */
  47. /* L1     =  6.6666666666666703212E-1    , Hex  2^  0   *  .AAAAAAAAAAAAC5 */
  48. /* L2     =  3.9999999999970461961E-1    , Hex  2^ -1   *  .CCCCCCCCCC2684 */
  49. /* L3     =  2.8571428579395698188E-1    , Hex  2^ -1   *  .92492492F85782 */
  50. /* L4     =  2.2222221233634724402E-1    , Hex  2^ -2   *  .E38E3839B7AF2C */
  51. /* L5     =  1.8181879517064680057E-1    , Hex  2^ -2   *  .BA2EB4CC39655E */
  52. /* L6     =  1.5382888777946145467E-1    , Hex  2^ -2   *  .9D8551E8C5781D */
  53. /* L7     =  1.3338356561139403517E-1    , Hex  2^ -2   *  .8895B3907FCD92 */
  54. /* L8     =  1.2500000000000000000E-1    , Hex  2^ -2   *  .80000000000000 */
  55. static long        L1x[] = { _0x(aaaa,402a), _0x(aac5,aaaa)};
  56. static long        L2x[] = { _0x(cccc,3fcc), _0x(2684,cccc)};
  57. static long        L3x[] = { _0x(4924,3f92), _0x(5782,92f8)};
  58. static long        L4x[] = { _0x(8e38,3f63), _0x(af2c,39b7)};
  59. static long        L5x[] = { _0x(2eb4,3f3a), _0x(655e,cc39)};
  60. static long        L6x[] = { _0x(8551,3f1d), _0x(781d,e8c5)};
  61. static long        L7x[] = { _0x(95b3,3f08), _0x(cd92,907f)};
  62. static long        L8x[] = { _0x(0000,3f00), _0x(0000,0000)};
  63. #define       L1    (*(double*)L1x)
  64. #define       L2    (*(double*)L2x)
  65. #define       L3    (*(double*)L3x)
  66. #define       L4    (*(double*)L4x)
  67. #define       L5    (*(double*)L5x)
  68. #define       L6    (*(double*)L6x)
  69. #define       L7    (*(double*)L7x)
  70. #define       L8    (*(double*)L8x)
  71. #else /* defined(vax)||defined(tahoe) */
  72. static double
  73. L1     =  6.6666666666667340202E-1    , /*Hex  2^ -1   *  1.5555555555592 */
  74. L2     =  3.9999999999416702146E-1    , /*Hex  2^ -2   *  1.999999997FF24 */
  75. L3     =  2.8571428742008753154E-1    , /*Hex  2^ -2   *  1.24924941E07B4 */
  76. L4     =  2.2222198607186277597E-1    , /*Hex  2^ -3   *  1.C71C52150BEA6 */
  77. L5     =  1.8183562745289935658E-1    , /*Hex  2^ -3   *  1.74663CC94342F */
  78. L6     =  1.5314087275331442206E-1    , /*Hex  2^ -3   *  1.39A1EC014045B */
  79. L7     =  1.4795612545334174692E-1    ; /*Hex  2^ -3   *  1.2F039F0085122 */
  80. #endif /* defined(vax)||defined(tahoe) */
  82. /*****************************************************************************
  83. * log__l -
  84. *
  85. * log__L(Z)
  86. * LOG(1+X) - 2S        X
  87. * RETURN      ---------------  WHERE Z = S*S,  S = ------- , 0 <= Z <= .0294...
  88. *       S      2 + X
  89. *
  90. * DOUBLE PRECISION (VAX D FORMAT 56 bits or IEEE DOUBLE 53 BITS)
  91. * KERNEL FUNCTION FOR LOG; TO BE USED IN LOG1P, LOG, AND POW FUNCTIONS
  92. * CODED IN C BY K.C. NG, 1/19/85;
  93. * REVISED BY K.C. Ng, 2/3/85, 4/16/85.
  94. *
  95. * Method :
  96. * 1. Polynomial approximation: let s = x/(2+x).
  97. *    Based on log(1+x) = log(1+s) - log(1-s)
  98. *  = 2s + 2/3 s**3 + 2/5 s**5 + .....,
  99. *
  100. *    (log(1+x) - 2s)/s is computed by
  101. *
  102. *        z*(L1 + z*(L2 + z*(... (L7 + z*L8)...)))
  103. *
  104. *    where z=s*s. (See the listing below for Lk's values.) The
  105. *    coefficients are obtained by a special Remez algorithm.
  106. *
  107. * Accuracy:
  108. * Assuming no rounding error, the maximum magnitude of the approximation
  109. * error (absolute) is 2**(-58.49) for IEEE double, and 2**(-63.63)
  110. * for VAX D format.
  111. *
  112. * Constants:
  113. * The hexadecimal values are the intended ones for the following constants.
  114. * The decimal values may be used, provided that the compiler will convert
  115. * from decimal to binary accurately enough to produce the hexadecimal values
  116. * shown.
  117. NOMANUAL
  118. */
  120. double log__L(z)
  121. double z;
  122. {
  123. #if defined(vax)||defined(tahoe)
  124.     return(z*(L1+z*(L2+z*(L3+z*(L4+z*(L5+z*(L6+z*(L7+z*L8))))))));
  125. #else /* defined(vax)||defined(tahoe) */
  126.     return(z*(L1+z*(L2+z*(L3+z*(L4+z*(L5+z*(L6+z*L7)))))));
  127. #endif /* defined(vax)||defined(tahoe) */
  128. }