libgcc-math/dbl-64/ulog.h

   1 /*
   2  * IBM Accurate Mathematical Library
   3  * Written by International Business Machines Corp.
   4  * Copyright (C) 2001 Free Software Foundation, Inc.
   5  *
   6  * This program is free software; you can redistribute it and/or modify
   7  * it under the terms of the GNU Lesser General Public License as published by
   8  * the Free Software Foundation; either version 2.1 of the License, or
   9  * (at your option) any later version.
  10  *
  11  * This program is distributed in the hope that it will be useful,
  12  * but WITHOUT ANY WARRANTY; without even the implied warranty of
  13  * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
  14  * GNU Lesser General Public License for more details.
  15  *
  16  * You should have received a copy of the GNU Lesser General Public License
  17  * along with this program; if not, write to the Free Software
  18  * Foundation, Inc., 59 Temple Place - Suite 330, Boston, MA 02111-1307, USA.
  19  */
  20
  21 /******************************************************************/
  22 /*                                                                */
  23 /* MODULE_NAME:ulog.h                                             */
  24 /*                                                                */
  25 /* common data and variables prototype and definition             */
  26 /******************************************************************/
  27
  28 #ifndef ULOG_H
  29 #define ULOG_H
  30
  31 #ifdef BIG_ENDI
  32   static const number
  33   /* polynomial I */
  34 /**/ a2             = {{0xbfe00000, 0x0001aa8f} }, /* -0.500... */
  35 /**/ a3             = {{0x3fd55555, 0x55588d2e} }, /*  0.333... */
  36   /* polynomial II */
  37 /**/ b0             = {{0x3fd55555, 0x55555555} }, /*  0.333... */
  38 /**/ b1             = {{0xbfcfffff, 0xffffffbb} }, /* -0.249... */
  39 /**/ b2             = {{0x3fc99999, 0x9999992f} }, /*  0.199... */
  40 /**/ b3             = {{0xbfc55555, 0x556503fd} }, /* -0.166... */
  41 /**/ b4             = {{0x3fc24924, 0x925b3d62} }, /*  0.142... */
  42 /**/ b5             = {{0xbfbffffe, 0x160472fc} }, /* -0.124... */
  43 /**/ b6             = {{0x3fbc71c5, 0x25db58ac} }, /*  0.111... */
  44 /**/ b7             = {{0xbfb9a4ac, 0x11a2a61c} }, /* -0.100... */
  45 /**/ b8             = {{0x3fb75077, 0x0df2b591} }, /*  0.091... */
  46   /* polynomial III */
  47 #if 0
  48 /**/ c1             = {{0x3ff00000, 0x00000000} }, /*  1        */
  49 #endif
  50 /**/ c2             = {{0xbfe00000, 0x00000000} }, /* -1/2      */
  51 /**/ c3             = {{0x3fd55555, 0x55555555} }, /*  1/3      */
  52 /**/ c4             = {{0xbfd00000, 0x00000000} }, /* -1/4      */
  53 /**/ c5             = {{0x3fc99999, 0x9999999a} }, /*  1/5      */
  54   /* polynomial IV */
  55 /**/ d2             = {{0xbfe00000, 0x00000000} }, /* -1/2      */
  56 /**/ dd2            = {{0x00000000, 0x00000000} }, /* -1/2-d2   */
  57 /**/ d3             = {{0x3fd55555, 0x55555555} }, /*  1/3      */
  58 /**/ dd3            = {{0x3c755555, 0x55555555} }, /*  1/3-d3   */
  59 /**/ d4             = {{0xbfd00000, 0x00000000} }, /* -1/4      */
  60 /**/ dd4            = {{0x00000000, 0x00000000} }, /* -1/4-d4   */
  61 /**/ d5             = {{0x3fc99999, 0x9999999a} }, /*  1/5      */
  62 /**/ dd5            = {{0xbc699999, 0x9999999a} }, /*  1/5-d5   */
  63 /**/ d6             = {{0xbfc55555, 0x55555555} }, /* -1/6      */
  64 /**/ dd6            = {{0xbc655555, 0x55555555} }, /* -1/6-d6   */
  65 /**/ d7             = {{0x3fc24924, 0x92492492} }, /*  1/7      */
  66 /**/ dd7            = {{0x3c624924, 0x92492492} }, /*  1/7-d7   */
  67 /**/ d8             = {{0xbfc00000, 0x00000000} }, /* -1/8      */
  68 /**/ dd8            = {{0x00000000, 0x00000000} }, /* -1/8-d8   */
  69 /**/ d9             = {{0x3fbc71c7, 0x1c71c71c} }, /*  1/9      */
  70 /**/ dd9            = {{0x3c5c71c7, 0x1c71c71c} }, /*  1/9-d9   */
  71 /**/ d10            = {{0xbfb99999, 0x9999999a} }, /* -1/10     */
  72 /**/ dd10           = {{0x3c599999, 0x9999999a} }, /* -1/10-d10 */
  73 /**/ d11            = {{0x3fb745d1, 0x745d1746} }, /*  1/11     */
  74 /**/ d12            = {{0xbfb55555, 0x55555555} }, /* -1/12     */
  75 /**/ d13            = {{0x3fb3b13b, 0x13b13b14} }, /*  1/13     */
  76 /**/ d14            = {{0xbfb24924, 0x92492492} }, /* -1/14     */
  77 /**/ d15            = {{0x3fb11111, 0x11111111} }, /*  1/15     */
  78 /**/ d16            = {{0xbfb00000, 0x00000000} }, /* -1/16     */
  79 /**/ d17            = {{0x3fae1e1e, 0x1e1e1e1e} }, /*  1/17     */
  80 /**/ d18            = {{0xbfac71c7, 0x1c71c71c} }, /* -1/18     */
  81 /**/ d19            = {{0x3faaf286, 0xbca1af28} }, /*  1/19     */
  82 /**/ d20            = {{0xbfa99999, 0x9999999a} }, /* -1/20     */
  83   /* constants    */
  84 /**/ zero           = {{0x00000000, 0x00000000} }, /* 0         */
  85 /**/ one            = {{0x3ff00000, 0x00000000} }, /* 1         */
  86 /**/ half           = {{0x3fe00000, 0x00000000} }, /* 1/2       */
  87 /**/ mhalf          = {{0xbfe00000, 0x00000000} }, /* -1/2      */
  88 /**/ sqrt_2         = {{0x3ff6a09e, 0x667f3bcc} }, /* sqrt(2)   */
  89 /**/ h1             = {{0x3fd2e000, 0x00000000} }, /* 151/2**9  */
  90 /**/ h2             = {{0x3f669000, 0x00000000} }, /* 361/2**17 */
  91 /**/ delu           = {{0x3f700000, 0x00000000} }, /* 1/2**8    */
  92 /**/ delv           = {{0x3ef00000, 0x00000000} }, /* 1/2**16   */
  93 /**/ ln2a           = {{0x3fe62e42, 0xfefa3800} }, /* ln(2) 43 bits */
  94 /**/ ln2b           = {{0x3d2ef357, 0x93c76730} }, /* ln(2)-ln2a    */
  95 /**/ e1             = {{0x3bbcc868, 0x00000000} }, /* 6.095e-21     */
  96 /**/ e2             = {{0x3c1138ce, 0x00000000} }, /* 2.334e-19     */
  97 /**/ e3             = {{0x3aa1565d, 0x00000000} }, /* 2.801e-26     */
  98 /**/ e4             = {{0x39809d88, 0x00000000} }, /* 1.024e-31     */
  99 /**/ e[M]           ={{{0x37da223a, 0x00000000} }, /* 1.2e-39       */
 100 /**/                  {{0x35c851c4, 0x00000000} }, /* 1.3e-49       */
 101 /**/                  {{0x2ab85e51, 0x00000000} }, /* 6.8e-103      */
 102 /**/                  {{0x17383827, 0x00000000} }},/* 8.1e-197      */
 103 /**/ two54          = {{0x43500000, 0x00000000} }, /* 2**54         */
 104 /**/ u03            = {{0x3f9eb851, 0xeb851eb8} }; /* 0.03          */
 105
 106 #else
 107 #ifdef LITTLE_ENDI
 108   static const number
 109   /* polynomial I */
 110 /**/ a2             = {{0x0001aa8f, 0xbfe00000} }, /* -0.500... */
 111 /**/ a3             = {{0x55588d2e, 0x3fd55555} }, /*  0.333... */
 112   /* polynomial II */
 113 /**/ b0             = {{0x55555555, 0x3fd55555} }, /*  0.333... */
 114 /**/ b1             = {{0xffffffbb, 0xbfcfffff} }, /* -0.249... */
 115 /**/ b2             = {{0x9999992f, 0x3fc99999} }, /*  0.199... */
 116 /**/ b3             = {{0x556503fd, 0xbfc55555} }, /* -0.166... */
 117 /**/ b4             = {{0x925b3d62, 0x3fc24924} }, /*  0.142... */
 118 /**/ b5             = {{0x160472fc, 0xbfbffffe} }, /* -0.124... */
 119 /**/ b6             = {{0x25db58ac, 0x3fbc71c5} }, /*  0.111... */
 120 /**/ b7             = {{0x11a2a61c, 0xbfb9a4ac} }, /* -0.100... */
 121 /**/ b8             = {{0x0df2b591, 0x3fb75077} }, /*  0.091... */
 122   /* polynomial III */
 123 #if 0
 124 /**/ c1             = {{0x00000000, 0x3ff00000} }, /*  1        */
 125 #endif
 126 /**/ c2             = {{0x00000000, 0xbfe00000} }, /* -1/2      */
 127 /**/ c3             = {{0x55555555, 0x3fd55555} }, /*  1/3      */
 128 /**/ c4             = {{0x00000000, 0xbfd00000} }, /* -1/4      */
 129 /**/ c5             = {{0x9999999a, 0x3fc99999} }, /*  1/5      */
 130   /* polynomial IV */
 131 /**/ d2             = {{0x00000000, 0xbfe00000} }, /* -1/2      */
 132 /**/ dd2            = {{0x00000000, 0x00000000} }, /* -1/2-d2   */
 133 /**/ d3             = {{0x55555555, 0x3fd55555} }, /*  1/3      */
 134 /**/ dd3            = {{0x55555555, 0x3c755555} }, /*  1/3-d3   */
 135 /**/ d4             = {{0x00000000, 0xbfd00000} }, /* -1/4      */
 136 /**/ dd4            = {{0x00000000, 0x00000000} }, /* -1/4-d4   */
 137 /**/ d5             = {{0x9999999a, 0x3fc99999} }, /*  1/5      */
 138 /**/ dd5            = {{0x9999999a, 0xbc699999} }, /*  1/5-d5   */
 139 /**/ d6             = {{0x55555555, 0xbfc55555} }, /* -1/6      */
 140 /**/ dd6            = {{0x55555555, 0xbc655555} }, /* -1/6-d6   */
 141 /**/ d7             = {{0x92492492, 0x3fc24924} }, /*  1/7      */
 142 /**/ dd7            = {{0x92492492, 0x3c624924} }, /*  1/7-d7   */
 143 /**/ d8             = {{0x00000000, 0xbfc00000} }, /* -1/8      */
 144 /**/ dd8            = {{0x00000000, 0x00000000} }, /* -1/8-d8   */
 145 /**/ d9             = {{0x1c71c71c, 0x3fbc71c7} }, /*  1/9      */
 146 /**/ dd9            = {{0x1c71c71c, 0x3c5c71c7} }, /*  1/9-d9   */
 147 /**/ d10            = {{0x9999999a, 0xbfb99999} }, /* -1/10     */
 148 /**/ dd10           = {{0x9999999a, 0x3c599999} }, /* -1/10-d10 */
 149 /**/ d11            = {{0x745d1746, 0x3fb745d1} }, /*  1/11     */
 150 /**/ d12            = {{0x55555555, 0xbfb55555} }, /* -1/12     */
 151 /**/ d13            = {{0x13b13b14, 0x3fb3b13b} }, /*  1/13     */
 152 /**/ d14            = {{0x92492492, 0xbfb24924} }, /* -1/14     */
 153 /**/ d15            = {{0x11111111, 0x3fb11111} }, /*  1/15     */
 154 /**/ d16            = {{0x00000000, 0xbfb00000} }, /* -1/16     */
 155 /**/ d17            = {{0x1e1e1e1e, 0x3fae1e1e} }, /*  1/17     */
 156 /**/ d18            = {{0x1c71c71c, 0xbfac71c7} }, /* -1/18     */
 157 /**/ d19            = {{0xbca1af28, 0x3faaf286} }, /*  1/19     */
 158 /**/ d20            = {{0x9999999a, 0xbfa99999} }, /* -1/20     */
 159   /* constants    */
 160 /**/ zero           = {{0x00000000, 0x00000000} }, /* 0         */
 161 /**/ one            = {{0x00000000, 0x3ff00000} }, /* 1         */
 162 /**/ half           = {{0x00000000, 0x3fe00000} }, /* 1/2       */
 163 /**/ mhalf          = {{0x00000000, 0xbfe00000} }, /* -1/2      */
 164 /**/ sqrt_2         = {{0x667f3bcc, 0x3ff6a09e} }, /* sqrt(2)   */
 165 /**/ h1             = {{0x00000000, 0x3fd2e000} }, /* 151/2**9  */
 166 /**/ h2             = {{0x00000000, 0x3f669000} }, /* 361/2**17 */
 167 /**/ delu           = {{0x00000000, 0x3f700000} }, /* 1/2**8    */
 168 /**/ delv           = {{0x00000000, 0x3ef00000} }, /* 1/2**16   */
 169 /**/ ln2a           = {{0xfefa3800, 0x3fe62e42} }, /* ln(2) 43 bits */
 170 /**/ ln2b           = {{0x93c76730, 0x3d2ef357} }, /* ln(2)-ln2a    */
 171 /**/ e1             = {{0x00000000, 0x3bbcc868} }, /* 6.095e-21     */
 172 /**/ e2             = {{0x00000000, 0x3c1138ce} }, /* 2.334e-19     */
 173 /**/ e3             = {{0x00000000, 0x3aa1565d} }, /* 2.801e-26     */
 174 /**/ e4             = {{0x00000000, 0x39809d88} }, /* 1.024e-31     */
 175 /**/ e[M]           ={{{0x00000000, 0x37da223a} }, /* 1.2e-39       */
 176 /**/                  {{0x00000000, 0x35c851c4} }, /* 1.3e-49       */
 177 /**/                  {{0x00000000, 0x2ab85e51} }, /* 6.8e-103      */
 178 /**/                  {{0x00000000, 0x17383827} }},/* 8.1e-197      */
 179 /**/ two54          = {{0x00000000, 0x43500000} }, /* 2**54         */
 180 /**/ u03            = {{0xeb851eb8, 0x3f9eb851} }; /* 0.03          */
 181
 182 #endif
 183 #endif
 184
 185 #define  ZERO      zero.d
 186 #define  ONE       one.d
 187 #define  HALF      half.d
 188 #define  MHALF     mhalf.d
 189 #define  SQRT_2    sqrt_2.d
 190 #define  DEL_U     delu.d
 191 #define  DEL_V     delv.d
 192 #define  LN2A      ln2a.d
 193 #define  LN2B      ln2b.d
 194 #define  E1        e1.d
 195 #define  E2        e2.d
 196 #define  E3        e3.d
 197 #define  E4        e4.d
 198 #define  U03       u03.d
 199
 200 #endif