sysdeps/ieee754/dbl-64/atnat2.h

   1
   2 /*
   3  * IBM Accurate Mathematical Library
   4  * Written by International Business Machines Corp.
   5  * Copyright (C) 2001 Free Software Foundation, Inc.
   6  *
   7  * This program is free software; you can redistribute it and/or modify
   8  * it under the terms of the GNU Lesser General Public License as published by
   9  * the Free Software Foundation; either version 2.1 of the License, or
  10  * (at your option) any later version.
  11  *
  12  * This program is distributed in the hope that it will be useful,
  13  * but WITHOUT ANY WARRANTY; without even the implied warranty of
  14  * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
  15  * GNU Lesser General Public License for more details.
  16  *
  17  * You should have received a copy of the GNU Lesser General Public License
  18  * along with this program; if not, see <http://www.gnu.org/licenses/>.
  19  */
  20
  21 /************************************************************************/
  22 /*  MODULE_NAME: atnat2.h                                                */
  23 /*                                                                      */
  24 /*                                                                      */
  25 /*      common data and variables definition for BIG or LITTLE ENDIAN   */
  26 /************************************************************************/
  27
  28
  29
  30 #ifndef ATNAT2_H
  31 #define ATNAT2_H
  32
  33
  34 #define MM 5
  35 #ifdef BIG_ENDI
  36
  37   static const number
  38   /* polynomial I */
  39 /**/ d3             = {{0xbfd55555, 0x55555555} }, /* -0.333... */
  40 /**/ d5             = {{0x3fc99999, 0x999997fd} }, /*  0.199... */
  41 /**/ d7             = {{0xbfc24924, 0x923f7603} }, /* -0.142... */
  42 /**/ d9             = {{0x3fbc71c6, 0xe5129a3b} }, /*  0.111... */
  43 /**/ d11            = {{0xbfb74580, 0x22b13c25} }, /* -0.090... */
  44 /**/ d13            = {{0x3fb375f0, 0x8b31cbce} }, /*  0.076... */
  45   /* polynomial II */
  46 /**/ f3             = {{0xbfd55555, 0x55555555} }, /* -1/3      */
  47 /**/ ff3            = {{0xbc755555, 0x55555555} }, /* -1/3-f3   */
  48 /**/ f5             = {{0x3fc99999, 0x9999999a} }, /*  1/5      */
  49 /**/ ff5            = {{0xbc699999, 0x9999999a} }, /*  1/5-f5   */
  50 /**/ f7             = {{0xbfc24924, 0x92492492} }, /* -1/7      */
  51 /**/ ff7            = {{0xbc624924, 0x92492492} }, /* -1/7-f7   */
  52 /**/ f9             = {{0x3fbc71c7, 0x1c71c71c} }, /*  1/9      */
  53 /**/ ff9            = {{0x3c5c71c7, 0x1c71c71c} }, /*  1/9-f9   */
  54 /**/ f11            = {{0xbfb745d1, 0x745d1746} }, /* -1/11     */
  55 /**/ f13            = {{0x3fb3b13b, 0x13b13b14} }, /*  1/13     */
  56 /**/ f15            = {{0xbfb11111, 0x11111111} }, /* -1/15     */
  57 /**/ f17            = {{0x3fae1e1e, 0x1e1e1e1e} }, /*  1/17     */
  58 /**/ f19            = {{0xbfaaf286, 0xbca1af28} }, /* -1/19     */
  59   /* constants    */
  60 /**/ zero           = {{0x00000000, 0x00000000} }, /*  0            */
  61 /**/ mzero          = {{0x80000000, 0x00000000} }, /* -0            */
  62 /**/ one            = {{0x3ff00000, 0x00000000} }, /*  1            */
  63 /**/ inv16          = {{0x3fb00000, 0x00000000} }, /*  1/16         */
  64 /**/ opi            = {{0x400921fb, 0x54442d18} }, /*  pi           */
  65 /**/ opi1           = {{0x3ca1a626, 0x33145c07} }, /*  pi-opi       */
  66 /**/ mopi           = {{0xc00921fb, 0x54442d18} }, /* -pi           */
  67 /**/ hpi            = {{0x3ff921fb, 0x54442d18} }, /*  pi/2         */
  68 /**/ hpi1           = {{0x3c91a626, 0x33145c07} }, /*  pi/2-hpi     */
  69 /**/ mhpi           = {{0xbff921fb, 0x54442d18} }, /* -pi/2         */
  70 /**/ qpi            = {{0x3fe921fb, 0x54442d18} }, /*  pi/4         */
  71 /**/ qpi1           = {{0x3c81a626, 0x33145c07} }, /*  pi/4-qpi     */
  72 /**/ mqpi           = {{0xbfe921fb, 0x54442d18} }, /* -pi/4         */
  73 /**/ tqpi           = {{0x4002d97c, 0x7f3321d2} }, /*  3pi/4        */
  74 /**/ tqpi1          = {{0x3c9a7939, 0x4c9e8a0a} }, /*  3pi/4-tqpi   */
  75 /**/ mtqpi          = {{0xc002d97c, 0x7f3321d2} }, /* -3pi/4        */
  76 /**/ u1             = {{0x3c314c2a, 0x00000000} }, /*  9.377e-19    */
  77 /**/ u2             = {{0x3bf955e4, 0x00000000} }, /*  8.584e-20    */
  78 /**/ u3             = {{0x3bf955e4, 0x00000000} }, /*  8.584e-20    */
  79 /**/ u4             = {{0x3bf955e4, 0x00000000} }, /*  8.584e-20    */
  80 /**/ u5             = {{0x3aaef2d1, 0x00000000} }, /*  5e-26        */
  81 /**/ u6             = {{0x3a6eeb36, 0x00000000} }, /*  3.122e-27    */
  82 /**/ u7             = {{0x3a6eeb36, 0x00000000} }, /*  3.122e-27    */
  83 /**/ u8             = {{0x3a6eeb36, 0x00000000} }, /*  3.122e-27    */
  84 /**/ u91            = {{0x3c6dffc0, 0x00000000} }, /*  1.301e-17    */
  85 /**/ u92            = {{0x3c527bd0, 0x00000000} }, /*  4.008e-18    */
  86 /**/ u93            = {{0x3c3cd057, 0x00000000} }, /*  1.562e-18    */
  87 /**/ u94            = {{0x3c329cdf, 0x00000000} }, /*  1.009e-18    */
  88 /**/ ua1            = {{0x3c3a1edf, 0x00000000} }, /*  1.416e-18    */
  89 /**/ ua2            = {{0x3c33f0e1, 0x00000000} }, /*  1.081e-18    */
  90 /**/ ub             = {{0x3a98c56d, 0x00000000} }, /*  2.001e-26    */
  91 /**/ uc             = {{0x3a9375de, 0x00000000} }, /*  1.572e-26    */
  92 /**/ ud[MM]         ={{{0x38c6eddf, 0x00000000} }, /*  3.450e-35    */
  93 /**/                  {{0x35c6ef60, 0x00000000} }, /*  1.226e-49    */
  94 /**/                  {{0x32c6ed2f, 0x00000000} }, /*  4.354e-64    */
  95 /**/                  {{0x23c6eee8, 0x00000000} }, /*  2.465e-136   */
  96 /**/                  {{0x11c6ed16, 0x00000000} }},/*  4.955e-223   */
  97 /**/ ue             = {{0x38900e9d, 0x00000000} }, /*  3.02e-36     */
  98 /**/ two8           = {{0x40700000, 0x00000000} }, /*  2**8=256     */
  99 /**/ two52          = {{0x43300000, 0x00000000} }, /*  2**52        */
 100 /**/ two500         = {{0x5f300000, 0x00000000} }, /*  2**500       */
 101 /**/ twom500        = {{0x20b00000, 0x00000000} }, /*  2**(-500)    */
 102 /**/ twom1022       = {{0x00100000, 0x00000000} }; /*  2**(-1022)   */
 103
 104 #else
 105 #ifdef LITTLE_ENDI
 106
 107   static const number
 108   /* polynomial I */
 109 /**/ d3             = {{0x55555555, 0xbfd55555} }, /* -0.333... */
 110 /**/ d5             = {{0x999997fd, 0x3fc99999} }, /*  0.199... */
 111 /**/ d7             = {{0x923f7603, 0xbfc24924} }, /* -0.142... */
 112 /**/ d9             = {{0xe5129a3b, 0x3fbc71c6} }, /*  0.111... */
 113 /**/ d11            = {{0x22b13c25, 0xbfb74580} }, /* -0.090... */
 114 /**/ d13            = {{0x8b31cbce, 0x3fb375f0} }, /*  0.076... */
 115   /* polynomial II */
 116 /**/ f3             = {{0x55555555, 0xbfd55555} }, /* -1/3      */
 117 /**/ ff3            = {{0x55555555, 0xbc755555} }, /* -1/3-f3   */
 118 /**/ f5             = {{0x9999999a, 0x3fc99999} }, /*  1/5      */
 119 /**/ ff5            = {{0x9999999a, 0xbc699999} }, /*  1/5-f5   */
 120 /**/ f7             = {{0x92492492, 0xbfc24924} }, /* -1/7      */
 121 /**/ ff7            = {{0x92492492, 0xbc624924} }, /* -1/7-f7   */
 122 /**/ f9             = {{0x1c71c71c, 0x3fbc71c7} }, /*  1/9      */
 123 /**/ ff9            = {{0x1c71c71c, 0x3c5c71c7} }, /*  1/9-f9   */
 124 /**/ f11            = {{0x745d1746, 0xbfb745d1} }, /* -1/11     */
 125 /**/ f13            = {{0x13b13b14, 0x3fb3b13b} }, /*  1/13     */
 126 /**/ f15            = {{0x11111111, 0xbfb11111} }, /* -1/15     */
 127 /**/ f17            = {{0x1e1e1e1e, 0x3fae1e1e} }, /*  1/17     */
 128 /**/ f19            = {{0xbca1af28, 0xbfaaf286} }, /* -1/19     */
 129   /* constants    */
 130 /**/ zero           = {{0x00000000, 0x00000000} }, /*  0            */
 131 /**/ mzero          = {{0x00000000, 0x80000000} }, /* -0            */
 132 /**/ one            = {{0x00000000, 0x3ff00000} }, /*  1            */
 133 /**/ inv16          = {{0x00000000, 0x3fb00000} }, /*  1/16         */
 134 /**/ opi            = {{0x54442d18, 0x400921fb} }, /*  pi           */
 135 /**/ opi1           = {{0x33145c07, 0x3ca1a626} }, /*  pi-opi       */
 136 /**/ mopi           = {{0x54442d18, 0xc00921fb} }, /* -pi           */
 137 /**/ hpi            = {{0x54442d18, 0x3ff921fb} }, /*  pi/2         */
 138 /**/ hpi1           = {{0x33145c07, 0x3c91a626} }, /*  pi/2-hpi     */
 139 /**/ mhpi           = {{0x54442d18, 0xbff921fb} }, /* -pi/2         */
 140 /**/ qpi            = {{0x54442d18, 0x3fe921fb} }, /*  pi/4         */
 141 /**/ qpi1           = {{0x33145c07, 0x3c81a626} }, /*  pi/4-qpi     */
 142 /**/ mqpi           = {{0x54442d18, 0xbfe921fb} }, /* -pi/4         */
 143 /**/ tqpi           = {{0x7f3321d2, 0x4002d97c} }, /*  3pi/4        */
 144 /**/ tqpi1          = {{0x4c9e8a0a, 0x3c9a7939} }, /*  3pi/4-tqpi   */
 145 /**/ mtqpi          = {{0x7f3321d2, 0xc002d97c} }, /* -3pi/4        */
 146 /**/ u1             = {{0x00000000, 0x3c314c2a} }, /*  9.377e-19    */
 147 /**/ u2             = {{0x00000000, 0x3bf955e4} }, /*  8.584e-20    */
 148 /**/ u3             = {{0x00000000, 0x3bf955e4} }, /*  8.584e-20    */
 149 /**/ u4             = {{0x00000000, 0x3bf955e4} }, /*  8.584e-20    */
 150 /**/ u5             = {{0x00000000, 0x3aaef2d1} }, /*  5e-26        */
 151 /**/ u6             = {{0x00000000, 0x3a6eeb36} }, /*  3.122e-27    */
 152 /**/ u7             = {{0x00000000, 0x3a6eeb36} }, /*  3.122e-27    */
 153 /**/ u8             = {{0x00000000, 0x3a6eeb36} }, /*  3.122e-27    */
 154 /**/ u91            = {{0x00000000, 0x3c6dffc0} }, /*  1.301e-17    */
 155 /**/ u92            = {{0x00000000, 0x3c527bd0} }, /*  4.008e-18    */
 156 /**/ u93            = {{0x00000000, 0x3c3cd057} }, /*  1.562e-18    */
 157 /**/ u94            = {{0x00000000, 0x3c329cdf} }, /*  1.009e-18    */
 158 /**/ ua1            = {{0x00000000, 0x3c3a1edf} }, /*  1.416e-18    */
 159 /**/ ua2            = {{0x00000000, 0x3c33f0e1} }, /*  1.081e-18    */
 160 /**/ ub             = {{0x00000000, 0x3a98c56d} }, /*  2.001e-26    */
 161 /**/ uc             = {{0x00000000, 0x3a9375de} }, /*  1.572e-26    */
 162 /**/ ud[MM]         ={{{0x00000000, 0x38c6eddf} }, /*  3.450e-35    */
 163 /**/                  {{0x00000000, 0x35c6ef60} }, /*  1.226e-49    */
 164 /**/                  {{0x00000000, 0x32c6ed2f} }, /*  4.354e-64    */
 165 /**/                  {{0x00000000, 0x23c6eee8} }, /*  2.465e-136   */
 166 /**/                  {{0x00000000, 0x11c6ed16} }},/*  4.955e-223   */
 167 /**/ ue             = {{0x00000000, 0x38900e9d} }, /*  3.02e-36     */
 168 /**/ two8           = {{0x00000000, 0x40700000} }, /*  2**8=256     */
 169 /**/ two52          = {{0x00000000, 0x43300000} }, /*  2**52        */
 170 /**/ two500         = {{0x00000000, 0x5f300000} }, /*  2**500       */
 171 /**/ twom500        = {{0x00000000, 0x20b00000} }, /*  2**(-500)    */
 172 /**/ twom1022       = {{0x00000000, 0x00100000} }; /*  2**(-1022)   */
 173
 174 #endif
 175 #endif
 176
 177 #define  ZERO      zero.d
 178 #define  MZERO     mzero.d
 179 #define  ONE       one.d
 180 #define  TWO8      two8.d
 181 #define  TWO52     two52.d
 182 #define  TWOM1022  twom1022.d
 183
 184 #endif