Back to index

glibc  2.9
atnat2.h
Go to the documentation of this file.
00001 
00002 /*
00003  * IBM Accurate Mathematical Library
00004  * Written by International Business Machines Corp.
00005  * Copyright (C) 2001 Free Software Foundation, Inc.
00006  *
00007  * This program is free software; you can redistribute it and/or modify
00008  * it under the terms of the GNU Lesser General Public License as published by
00009  * the Free Software Foundation; either version 2.1 of the License, or
00010  * (at your option) any later version.
00011  *
00012  * This program is distributed in the hope that it will be useful,
00013  * but WITHOUT ANY WARRANTY; without even the implied warranty of
00014  * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
00015  * GNU Lesser General Public License for more details.
00016  *
00017  * You should have received a copy of the GNU Lesser General Public License
00018  * along with this program; if not, write to the Free Software
00019  * Foundation, Inc., 59 Temple Place - Suite 330, Boston, MA 02111-1307, USA.
00020  */
00021 
00022 /************************************************************************/
00023 /*  MODULE_NAME: atnat2.h                                                */
00024 /*                                                                      */
00025 /*                                                                      */
00026 /*     common data and variables definition for BIG or LITTLE ENDIAN   */
00027 /************************************************************************/
00028 
00029 
00030 
00031 #ifndef ATNAT2_H
00032 #define ATNAT2_H
00033 
00034 
00035 #define MM 5
00036 #ifdef BIG_ENDI
00037 
00038   static const number
00039   /* polynomial I */
00040  d3             = {{0xbfd55555, 0x55555555} }, /* -0.333... */
00041  d5             = {{0x3fc99999, 0x999997fd} }, /*  0.199... */
00042  d7             = {{0xbfc24924, 0x923f7603} }, /* -0.142... */
00043  d9             = {{0x3fbc71c6, 0xe5129a3b} }, /*  0.111... */
00044  d11            = {{0xbfb74580, 0x22b13c25} }, /* -0.090... */
00045  d13            = {{0x3fb375f0, 0x8b31cbce} }, /*  0.076... */
00046   /* polynomial II */
00047  f3             = {{0xbfd55555, 0x55555555} }, /* -1/3      */
00048  ff3            = {{0xbc755555, 0x55555555} }, /* -1/3-f3   */
00049  f5             = {{0x3fc99999, 0x9999999a} }, /*  1/5      */
00050  ff5            = {{0xbc699999, 0x9999999a} }, /*  1/5-f5   */
00051  f7             = {{0xbfc24924, 0x92492492} }, /* -1/7      */
00052  ff7            = {{0xbc624924, 0x92492492} }, /* -1/7-f7   */
00053  f9             = {{0x3fbc71c7, 0x1c71c71c} }, /*  1/9      */
00054  ff9            = {{0x3c5c71c7, 0x1c71c71c} }, /*  1/9-f9   */
00055  f11            = {{0xbfb745d1, 0x745d1746} }, /* -1/11     */
00056  f13            = {{0x3fb3b13b, 0x13b13b14} }, /*  1/13     */
00057  f15            = {{0xbfb11111, 0x11111111} }, /* -1/15     */
00058  f17            = {{0x3fae1e1e, 0x1e1e1e1e} }, /*  1/17     */
00059  f19            = {{0xbfaaf286, 0xbca1af28} }, /* -1/19     */
00060   /* constants    */
00061  zero           = {{0x00000000, 0x00000000} }, /*  0            */
00062  mzero          = {{0x80000000, 0x00000000} }, /* -0            */
00063  one            = {{0x3ff00000, 0x00000000} }, /*  1            */
00064  inv16          = {{0x3fb00000, 0x00000000} }, /*  1/16         */
00065  opi            = {{0x400921fb, 0x54442d18} }, /*  pi           */
00066  opi1           = {{0x3ca1a626, 0x33145c07} }, /*  pi-opi       */
00067  mopi           = {{0xc00921fb, 0x54442d18} }, /* -pi           */
00068  hpi            = {{0x3ff921fb, 0x54442d18} }, /*  pi/2         */
00069  hpi1           = {{0x3c91a626, 0x33145c07} }, /*  pi/2-hpi     */
00070  mhpi           = {{0xbff921fb, 0x54442d18} }, /* -pi/2         */
00071  qpi            = {{0x3fe921fb, 0x54442d18} }, /*  pi/4         */
00072  qpi1           = {{0x3c81a626, 0x33145c07} }, /*  pi/4-qpi     */
00073  mqpi           = {{0xbfe921fb, 0x54442d18} }, /* -pi/4         */
00074  tqpi           = {{0x4002d97c, 0x7f3321d2} }, /*  3pi/4        */
00075  tqpi1          = {{0x3c9a7939, 0x4c9e8a0a} }, /*  3pi/4-tqpi   */
00076  mtqpi          = {{0xc002d97c, 0x7f3321d2} }, /* -3pi/4        */
00077  u1             = {{0x3c314c2a, 0x00000000} }, /*  9.377e-19    */
00078  u2             = {{0x3bf955e4, 0x00000000} }, /*  8.584e-20    */
00079  u3             = {{0x3bf955e4, 0x00000000} }, /*  8.584e-20    */
00080  u4             = {{0x3bf955e4, 0x00000000} }, /*  8.584e-20    */
00081  u5             = {{0x3aaef2d1, 0x00000000} }, /*  5e-26        */
00082  u6             = {{0x3a6eeb36, 0x00000000} }, /*  3.122e-27    */
00083  u7             = {{0x3a6eeb36, 0x00000000} }, /*  3.122e-27    */
00084  u8             = {{0x3a6eeb36, 0x00000000} }, /*  3.122e-27    */
00085  u91            = {{0x3c6dffc0, 0x00000000} }, /*  1.301e-17    */
00086  u92            = {{0x3c527bd0, 0x00000000} }, /*  4.008e-18    */
00087  u93            = {{0x3c3cd057, 0x00000000} }, /*  1.562e-18    */
00088  u94            = {{0x3c329cdf, 0x00000000} }, /*  1.009e-18    */
00089  ua1            = {{0x3c3a1edf, 0x00000000} }, /*  1.416e-18    */
00090  ua2            = {{0x3c33f0e1, 0x00000000} }, /*  1.081e-18    */
00091  ub             = {{0x3a98c56d, 0x00000000} }, /*  2.001e-26    */
00092  uc             = {{0x3a9375de, 0x00000000} }, /*  1.572e-26    */
00093  ud[MM]         ={{{0x38c6eddf, 0x00000000} }, /*  3.450e-35    */
00094                   {{0x35c6ef60, 0x00000000} }, /*  1.226e-49    */
00095                   {{0x32c6ed2f, 0x00000000} }, /*  4.354e-64    */
00096                   {{0x23c6eee8, 0x00000000} }, /*  2.465e-136   */
00097                   {{0x11c6ed16, 0x00000000} }},/*  4.955e-223   */
00098  ue             = {{0x38900e9d, 0x00000000} }, /*  3.02e-36     */
00099  two8           = {{0x40700000, 0x00000000} }, /*  2**8=256     */
00100  two52          = {{0x43300000, 0x00000000} }, /*  2**52        */
00101  two500         = {{0x5f300000, 0x00000000} }, /*  2**500       */
00102  twom500        = {{0x20b00000, 0x00000000} }, /*  2**(-500)    */
00103  twom1022       = {{0x00100000, 0x00000000} }; /*  2**(-1022)   */
00104 
00105 #else
00106 #ifdef LITTLE_ENDI
00107 
00108   static const number
00109   /* polynomial I */
00110  d3             = {{0x55555555, 0xbfd55555} }, /* -0.333... */
00111  d5             = {{0x999997fd, 0x3fc99999} }, /*  0.199... */
00112  d7             = {{0x923f7603, 0xbfc24924} }, /* -0.142... */
00113  d9             = {{0xe5129a3b, 0x3fbc71c6} }, /*  0.111... */
00114  d11            = {{0x22b13c25, 0xbfb74580} }, /* -0.090... */
00115  d13            = {{0x8b31cbce, 0x3fb375f0} }, /*  0.076... */
00116   /* polynomial II */
00117  f3             = {{0x55555555, 0xbfd55555} }, /* -1/3      */
00118  ff3            = {{0x55555555, 0xbc755555} }, /* -1/3-f3   */
00119  f5             = {{0x9999999a, 0x3fc99999} }, /*  1/5      */
00120  ff5            = {{0x9999999a, 0xbc699999} }, /*  1/5-f5   */
00121  f7             = {{0x92492492, 0xbfc24924} }, /* -1/7      */
00122  ff7            = {{0x92492492, 0xbc624924} }, /* -1/7-f7   */
00123  f9             = {{0x1c71c71c, 0x3fbc71c7} }, /*  1/9      */
00124  ff9            = {{0x1c71c71c, 0x3c5c71c7} }, /*  1/9-f9   */
00125  f11            = {{0x745d1746, 0xbfb745d1} }, /* -1/11     */
00126  f13            = {{0x13b13b14, 0x3fb3b13b} }, /*  1/13     */
00127  f15            = {{0x11111111, 0xbfb11111} }, /* -1/15     */
00128  f17            = {{0x1e1e1e1e, 0x3fae1e1e} }, /*  1/17     */
00129  f19            = {{0xbca1af28, 0xbfaaf286} }, /* -1/19     */
00130   /* constants    */
00131  zero           = {{0x00000000, 0x00000000} }, /*  0            */
00132  mzero          = {{0x00000000, 0x80000000} }, /* -0            */
00133  one            = {{0x00000000, 0x3ff00000} }, /*  1            */
00134  inv16          = {{0x00000000, 0x3fb00000} }, /*  1/16         */
00135  opi            = {{0x54442d18, 0x400921fb} }, /*  pi           */
00136  opi1           = {{0x33145c07, 0x3ca1a626} }, /*  pi-opi       */
00137  mopi           = {{0x54442d18, 0xc00921fb} }, /* -pi           */
00138  hpi            = {{0x54442d18, 0x3ff921fb} }, /*  pi/2         */
00139  hpi1           = {{0x33145c07, 0x3c91a626} }, /*  pi/2-hpi     */
00140  mhpi           = {{0x54442d18, 0xbff921fb} }, /* -pi/2         */
00141  qpi            = {{0x54442d18, 0x3fe921fb} }, /*  pi/4         */
00142  qpi1           = {{0x33145c07, 0x3c81a626} }, /*  pi/4-qpi     */
00143  mqpi           = {{0x54442d18, 0xbfe921fb} }, /* -pi/4         */
00144  tqpi           = {{0x7f3321d2, 0x4002d97c} }, /*  3pi/4        */
00145  tqpi1          = {{0x4c9e8a0a, 0x3c9a7939} }, /*  3pi/4-tqpi   */
00146  mtqpi          = {{0x7f3321d2, 0xc002d97c} }, /* -3pi/4        */
00147  u1             = {{0x00000000, 0x3c314c2a} }, /*  9.377e-19    */
00148  u2             = {{0x00000000, 0x3bf955e4} }, /*  8.584e-20    */
00149  u3             = {{0x00000000, 0x3bf955e4} }, /*  8.584e-20    */
00150  u4             = {{0x00000000, 0x3bf955e4} }, /*  8.584e-20    */
00151  u5             = {{0x00000000, 0x3aaef2d1} }, /*  5e-26        */
00152  u6             = {{0x00000000, 0x3a6eeb36} }, /*  3.122e-27    */
00153  u7             = {{0x00000000, 0x3a6eeb36} }, /*  3.122e-27    */
00154  u8             = {{0x00000000, 0x3a6eeb36} }, /*  3.122e-27    */
00155  u91            = {{0x00000000, 0x3c6dffc0} }, /*  1.301e-17    */
00156  u92            = {{0x00000000, 0x3c527bd0} }, /*  4.008e-18    */
00157  u93            = {{0x00000000, 0x3c3cd057} }, /*  1.562e-18    */
00158  u94            = {{0x00000000, 0x3c329cdf} }, /*  1.009e-18    */
00159  ua1            = {{0x00000000, 0x3c3a1edf} }, /*  1.416e-18    */
00160  ua2            = {{0x00000000, 0x3c33f0e1} }, /*  1.081e-18    */
00161  ub             = {{0x00000000, 0x3a98c56d} }, /*  2.001e-26    */
00162  uc             = {{0x00000000, 0x3a9375de} }, /*  1.572e-26    */
00163  ud[MM]         ={{{0x00000000, 0x38c6eddf} }, /*  3.450e-35    */
00164                   {{0x00000000, 0x35c6ef60} }, /*  1.226e-49    */
00165                   {{0x00000000, 0x32c6ed2f} }, /*  4.354e-64    */
00166                   {{0x00000000, 0x23c6eee8} }, /*  2.465e-136   */
00167                   {{0x00000000, 0x11c6ed16} }},/*  4.955e-223   */
00168  ue             = {{0x00000000, 0x38900e9d} }, /*  3.02e-36     */
00169  two8           = {{0x00000000, 0x40700000} }, /*  2**8=256     */
00170  two52          = {{0x00000000, 0x43300000} }, /*  2**52        */
00171  two500         = {{0x00000000, 0x5f300000} }, /*  2**500       */
00172  twom500        = {{0x00000000, 0x20b00000} }, /*  2**(-500)    */
00173  twom1022       = {{0x00000000, 0x00100000} }; /*  2**(-1022)   */
00174 
00175 #endif
00176 #endif
00177 
00178 #define  ZERO      zero.d
00179 #define  MZERO     mzero.d
00180 #define  ONE       one.d
00181 #define  TWO8      two8.d
00182 #define  TWO52     two52.d
00183 #define  TWOM1022  twom1022.d
00184 
00185 #endif