sysdeps/ieee754/dbl-64/atnat.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 /*  MODULE_NAME: atnat.h                                                */
  23 /*                                                                      */
  24 /*                                                                      */
  25 /*      common data and variables definition for BIG or LITTLE ENDIAN   */
  26 /************************************************************************/
  27 #ifndef ATNAT_H
  28 #define ATNAT_H
  29
  30 #define M 4
  31
  32 #ifdef BIG_ENDI
  33   static const number
  34   /* polynomial I */
  35 /**/ d3             = {{0xbfd55555, 0x55555555} }, /* -0.333... */
  36 /**/ d5             = {{0x3fc99999, 0x999997fd} }, /*  0.199... */
  37 /**/ d7             = {{0xbfc24924, 0x923f7603} }, /* -0.142... */
  38 /**/ d9             = {{0x3fbc71c6, 0xe5129a3b} }, /*  0.111... */
  39 /**/ d11            = {{0xbfb74580, 0x22b13c25} }, /* -0.090... */
  40 /**/ d13            = {{0x3fb375f0, 0x8b31cbce} }, /*  0.076... */
  41   /* polynomial II */
  42 /**/ f3             = {{0xbfd55555, 0x55555555} }, /* -1/3      */
  43 /**/ ff3            = {{0xbc755555, 0x55555555} }, /* -1/3-f3   */
  44 /**/ f5             = {{0x3fc99999, 0x9999999a} }, /*  1/5      */
  45 /**/ ff5            = {{0xbc699999, 0x9999999a} }, /*  1/5-f5   */
  46 /**/ f7             = {{0xbfc24924, 0x92492492} }, /* -1/7      */
  47 /**/ ff7            = {{0xbc624924, 0x92492492} }, /* -1/7-f7   */
  48 /**/ f9             = {{0x3fbc71c7, 0x1c71c71c} }, /*  1/9      */
  49 /**/ ff9            = {{0x3c5c71c7, 0x1c71c71c} }, /*  1/9-f9   */
  50 /**/ f11            = {{0xbfb745d1, 0x745d1746} }, /* -1/11     */
  51 /**/ f13            = {{0x3fb3b13b, 0x13b13b14} }, /*  1/13     */
  52 /**/ f15            = {{0xbfb11111, 0x11111111} }, /* -1/15     */
  53 /**/ f17            = {{0x3fae1e1e, 0x1e1e1e1e} }, /*  1/17     */
  54 /**/ f19            = {{0xbfaaf286, 0xbca1af28} }, /* -1/19     */
  55   /* constants    */
  56 /**/ zero           = {{0x00000000, 0x00000000} }, /*  0            */
  57 /**/ one            = {{0x3ff00000, 0x00000000} }, /*  1            */
  58 /**/ a              = {{0x3e4bb67a, 0x00000000} }, /*  1.290e-8     */
  59 /**/ b              = {{0x3fb00000, 0x00000000} }, /*  1/16         */
  60 /**/ c              = {{0x3ff00000, 0x00000000} }, /*  1            */
  61 /**/ d              = {{0x40300000, 0x00000000} }, /*  16           */
  62 /**/ e              = {{0x43349ff2, 0x00000000} }, /*  5.805e15     */
  63 /**/ hpi            = {{0x3ff921fb, 0x54442d18} }, /*  pi/2         */
  64 /**/ mhpi           = {{0xbff921fb, 0x54442d18} }, /* -pi/2         */
  65 /**/ hpi1           = {{0x3c91a626, 0x33145c07} }, /*  pi/2-hpi     */
  66 /**/ u1             = {{0x3c2d3382, 0x00000000} }, /*  7.915e-19    */
  67 /**/ u21            = {{0x3c6dffc0, 0x00000000} }, /*  1.301e-17    */
  68 /**/ u22            = {{0x3c527bd0, 0x00000000} }, /*  4.008e-18    */
  69 /**/ u23            = {{0x3c3cd057, 0x00000000} }, /*  1.562e-18    */
  70 /**/ u24            = {{0x3c329cdf, 0x00000000} }, /*  1.009e-18    */
  71 /**/ u31            = {{0x3c3a1edf, 0x00000000} }, /*  1.416e-18    */
  72 /**/ u32            = {{0x3c33f0e1, 0x00000000} }, /*  1.081e-18    */
  73 /**/ u4             = {{0x3bf955e4, 0x00000000} }, /*  8.584e-20    */
  74 /**/ u5             = {{0x3aaef2d1, 0x00000000} }, /*  5e-26        */
  75 /**/ u6             = {{0x3a98c56d, 0x00000000} }, /*  2.001e-26    */
  76 /**/ u7             = {{0x3a9375de, 0x00000000} }, /*  1.572e-26    */
  77 /**/ u8             = {{0x3a6eeb36, 0x00000000} }, /*  3.122e-27    */
  78 /**/ u9[M]          ={{{0x38c1aa5b, 0x00000000} }, /* 2.658e-35     */
  79 /**/                  {{0x35c1aa4d, 0x00000000} }, /* 9.443e-50     */
  80 /**/                  {{0x32c1aa88, 0x00000000} }, /* 3.355e-64     */
  81 /**/                  {{0x11c1aa56, 0x00000000} }},/* 3.818e-223    */
  82 /**/ two8           = {{0x40700000, 0x00000000} }, /*  2**8=256     */
  83 /**/ two52          = {{0x43300000, 0x00000000} }; /*  2**52        */
  84
  85 #else
  86 #ifdef LITTLE_ENDI
  87   static const number
  88   /* polynomial I */
  89 /**/ d3             = {{0x55555555, 0xbfd55555} }, /* -0.333... */
  90 /**/ d5             = {{0x999997fd, 0x3fc99999} }, /*  0.199... */
  91 /**/ d7             = {{0x923f7603, 0xbfc24924} }, /* -0.142... */
  92 /**/ d9             = {{0xe5129a3b, 0x3fbc71c6} }, /*  0.111... */
  93 /**/ d11            = {{0x22b13c25, 0xbfb74580} }, /* -0.090... */
  94 /**/ d13            = {{0x8b31cbce, 0x3fb375f0} }, /*  0.076... */
  95   /* polynomial II */
  96 /**/ f3             = {{0x55555555, 0xbfd55555} }, /* -1/3      */
  97 /**/ ff3            = {{0x55555555, 0xbc755555} }, /* -1/3-f3   */
  98 /**/ f5             = {{0x9999999a, 0x3fc99999} }, /*  1/5      */
  99 /**/ ff5            = {{0x9999999a, 0xbc699999} }, /*  1/5-f5   */
 100 /**/ f7             = {{0x92492492, 0xbfc24924} }, /* -1/7      */
 101 /**/ ff7            = {{0x92492492, 0xbc624924} }, /* -1/7-f7   */
 102 /**/ f9             = {{0x1c71c71c, 0x3fbc71c7} }, /*  1/9      */
 103 /**/ ff9            = {{0x1c71c71c, 0x3c5c71c7} }, /*  1/9-f9   */
 104 /**/ f11            = {{0x745d1746, 0xbfb745d1} }, /* -1/11     */
 105 /**/ f13            = {{0x13b13b14, 0x3fb3b13b} }, /*  1/13     */
 106 /**/ f15            = {{0x11111111, 0xbfb11111} }, /* -1/15     */
 107 /**/ f17            = {{0x1e1e1e1e, 0x3fae1e1e} }, /*  1/17     */
 108 /**/ f19            = {{0xbca1af28, 0xbfaaf286} }, /* -1/19     */
 109   /* constants    */
 110 /**/ zero           = {{0x00000000, 0x00000000} }, /*  0            */
 111 /**/ one            = {{0x00000000, 0x3ff00000} }, /*  1            */
 112 /**/ a              = {{0x00000000, 0x3e4bb67a} }, /*  1.290e-8     */
 113 /**/ b              = {{0x00000000, 0x3fb00000} }, /*  1/16         */
 114 /**/ c              = {{0x00000000, 0x3ff00000} }, /*  1            */
 115 /**/ d              = {{0x00000000, 0x40300000} }, /*  16           */
 116 /**/ e              = {{0x00000000, 0x43349ff2} }, /*  5.805e15     */
 117 /**/ hpi            = {{0x54442d18, 0x3ff921fb} }, /*  pi/2         */
 118 /**/ mhpi           = {{0x54442d18, 0xbff921fb} }, /* -pi/2         */
 119 /**/ hpi1           = {{0x33145c07, 0x3c91a626} }, /*  pi/2-hpi     */
 120 /**/ u1             = {{0x00000000, 0x3c2d3382} }, /*  7.915e-19    */
 121 /**/ u21            = {{0x00000000, 0x3c6dffc0} }, /*  1.301e-17    */
 122 /**/ u22            = {{0x00000000, 0x3c527bd0} }, /*  4.008e-18    */
 123 /**/ u23            = {{0x00000000, 0x3c3cd057} }, /*  1.562e-18    */
 124 /**/ u24            = {{0x00000000, 0x3c329cdf} }, /*  1.009e-18    */
 125 /**/ u31            = {{0x00000000, 0x3c3a1edf} }, /*  1.416e-18    */
 126 /**/ u32            = {{0x00000000, 0x3c33f0e1} }, /*  1.081e-18    */
 127 /**/ u4             = {{0x00000000, 0x3bf955e4} }, /*  8.584e-20    */
 128 /**/ u5             = {{0x00000000, 0x3aaef2d1} }, /*  5e-26        */
 129 /**/ u6             = {{0x00000000, 0x3a98c56d} }, /*  2.001e-26    */
 130 /**/ u7             = {{0x00000000, 0x3a9375de} }, /*  1.572e-26    */
 131 /**/ u8             = {{0x00000000, 0x3a6eeb36} }, /*  3.122e-27    */
 132 /**/ u9[M]          ={{{0x00000000, 0x38c1aa5b} }, /* 2.658e-35     */
 133 /**/                  {{0x00000000, 0x35c1aa4d} }, /* 9.443e-50     */
 134 /**/                  {{0x00000000, 0x32c1aa88} }, /* 3.355e-64     */
 135 /**/                  {{0x00000000, 0x11c1aa56} }},/* 3.818e-223    */
 136 /**/ two8           = {{0x00000000, 0x40700000} }, /*  2**8=256     */
 137 /**/ two52          = {{0x00000000, 0x43300000} }; /*  2**52        */
 138
 139 #endif
 140 #endif
 141
 142 #define  ZERO      zero.d
 143 #define  ONE       one.d
 144 #define  A         a.d
 145 #define  B         b.d
 146 #define  C         c.d
 147 #define  D         d.d
 148 #define  E         e.d
 149 #define  HPI       hpi.d
 150 #define  MHPI      mhpi.d
 151 #define  HPI1      hpi1.d
 152 #define  U1        u1.d
 153 #define  U21       u21.d
 154 #define  U22       u22.d
 155 #define  U23       u23.d
 156 #define  U24       u24.d
 157 #define  U31       u31.d
 158 #define  U32       u32.d
 159 #define  U4        u4.d
 160 #define  U5        u5.d
 161 #define  U6        u6.d
 162 #define  U7        u7.d
 163 #define  U8        u8.d
 164 #define  TWO8      two8.d
 165 #define  TWO52     two52.d
 166
 167 #endif