From 50944bca4bfeefc5de2ec205ad549eb669ed1008 Mon Sep 17 00:00:00 2001 From: Ulrich Drepper Date: Mon, 12 Mar 2001 07:57:09 +0000 Subject: Fix warnings. --- sysdeps/ieee754/dbl-64/asincos.tbl | 8 +- sysdeps/ieee754/dbl-64/atnat.h | 194 +- sysdeps/ieee754/dbl-64/atnat2.h | 254 +- sysdeps/ieee754/dbl-64/branred.c | 53 +- sysdeps/ieee754/dbl-64/branred.h | 72 +- sysdeps/ieee754/dbl-64/doasin.h | 46 +- sysdeps/ieee754/dbl-64/dosincos.c | 84 +- sysdeps/ieee754/dbl-64/dosincos.h | 74 +- sysdeps/ieee754/dbl-64/e_asin.c | 42 +- sysdeps/ieee754/dbl-64/e_atan2.c | 16 +- sysdeps/ieee754/dbl-64/e_exp.c | 14 +- sysdeps/ieee754/dbl-64/e_log.c | 7 +- sysdeps/ieee754/dbl-64/e_pow.c | 30 +- sysdeps/ieee754/dbl-64/e_remainder.c | 12 +- sysdeps/ieee754/dbl-64/endian.h | 4 + sysdeps/ieee754/dbl-64/mpa.c | 43 +- sysdeps/ieee754/dbl-64/mpa.h | 9 +- sysdeps/ieee754/dbl-64/mpa2.h | 74 +- sysdeps/ieee754/dbl-64/mpatan.c | 36 +- sysdeps/ieee754/dbl-64/mpatan.h | 212 +- sysdeps/ieee754/dbl-64/mpatan2.c | 12 +- sysdeps/ieee754/dbl-64/mpexp.c | 15 +- sysdeps/ieee754/dbl-64/mpexp.h | 198 +- sysdeps/ieee754/dbl-64/mplog.c | 15 +- sysdeps/ieee754/dbl-64/mplog.h | 10 +- sysdeps/ieee754/dbl-64/mpsqrt.c | 18 +- sysdeps/ieee754/dbl-64/mpsqrt.h | 14 +- sysdeps/ieee754/dbl-64/mptan.c | 10 +- sysdeps/ieee754/dbl-64/s_atan.c | 15 +- sysdeps/ieee754/dbl-64/s_sin.c | 49 +- sysdeps/ieee754/dbl-64/s_tan.c | 9 +- sysdeps/ieee754/dbl-64/sincos.tbl | 104 +- sysdeps/ieee754/dbl-64/sincos32.c | 116 +- sysdeps/ieee754/dbl-64/sincos32.h | 26 +- sysdeps/ieee754/dbl-64/slowexp.c | 27 +- sysdeps/ieee754/dbl-64/slowpow.c | 26 +- sysdeps/ieee754/dbl-64/uasncs.h | 22 +- sysdeps/ieee754/dbl-64/uatan.tbl | 22173 ++++++++++++++++----------------- sysdeps/ieee754/dbl-64/uexp.h | 44 +- sysdeps/ieee754/dbl-64/uexp.tbl | 16 +- sysdeps/ieee754/dbl-64/ulog.h | 270 +- sysdeps/ieee754/dbl-64/ulog.tbl | 6528 +++++----- sysdeps/ieee754/dbl-64/upow.h | 46 +- sysdeps/ieee754/dbl-64/upow.tbl | 16 +- sysdeps/ieee754/dbl-64/urem.h | 36 +- sysdeps/ieee754/dbl-64/usncs.h | 78 +- sysdeps/ieee754/dbl-64/utan.h | 422 +- sysdeps/ieee754/dbl-64/utan.tbl | 2977 +++-- 48 files changed, 17358 insertions(+), 17218 deletions(-) diff --git a/sysdeps/ieee754/dbl-64/asincos.tbl b/sysdeps/ieee754/dbl-64/asincos.tbl index ffae0bc288..5f41e54d2b 100644 --- a/sysdeps/ieee754/dbl-64/asincos.tbl +++ b/sysdeps/ieee754/dbl-64/asincos.tbl @@ -3,7 +3,7 @@ /***************************************************************************/ #ifdef BIG_ENDI -static const union {int4 i[5136];double x[2568];} asncs = { +static const union {int4 i[5136];double x[2568];} asncs = { .i = { /**/ 0x3FC04000, 0x00000000, /**/ 0x3FF02169, 0x88994424, /**/ 0x3FB0A6A2, 0xB799B115, @@ -2571,11 +2571,11 @@ static const union {int4 i[5136];double x[2568];} asncs = { /**/ 0x3C935690, 0xE395EEA6, /**/ 0x3FF4FF86, 0x2E5965A2, /**/ 0x400F4000, 0x00000000, -/**/ 0x3F7C6B82, 0xD36A5E70 }; +/**/ 0x3F7C6B82, 0xD36A5E70 } }; #else #ifdef LITTLE_ENDI -static const union {int4 i[5136];double x[2568];} asncs = { +static const union {int4 i[5136];double x[2568];} asncs = { .i = { /**/ 0x00000000, 0x3FC04000, /**/ 0x88994424, 0x3FF02169, /**/ 0xB799B115, 0x3FB0A6A2, @@ -5143,7 +5143,7 @@ static const union {int4 i[5136];double x[2568];} asncs = { /**/ 0xE395EEA6, 0x3C935690, /**/ 0x2E5965A2, 0x3FF4FF86, /**/ 0x00000000, 0x400F4000, -/**/ 0xD36A5E70, 0x3F7C6B82 }; +/**/ 0xD36A5E70, 0x3F7C6B82 } }; #endif #endif diff --git a/sysdeps/ieee754/dbl-64/atnat.h b/sysdeps/ieee754/dbl-64/atnat.h index 7ddf235f03..2a3d70d363 100644 --- a/sysdeps/ieee754/dbl-64/atnat.h +++ b/sysdeps/ieee754/dbl-64/atnat.h @@ -5,9 +5,9 @@ * * This program is free software; you can redistribute it and/or modify * it under the terms of the GNU Lesser General Public License as published by - * the Free Software Foundation; either version 2 of the License, or + * the Free Software Foundation; either version 2 of the License, or * (at your option) any later version. - * + * * This program is distributed in the hope that it will be useful, * but WITHOUT ANY WARRANTY; without even the implied warranty of * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the @@ -15,7 +15,7 @@ * * You should have received a copy of the GNU Lesser General Public License * along with this program; if not, write to the Free Software - * Foundation, Inc., 59 Temple Place - Suite 330, Boston, MA 02111-1307, USA. + * Foundation, Inc., 59 Temple Place - Suite 330, Boston, MA 02111-1307, USA. */ /************************************************************************/ /* MODULE_NAME: atnat.h */ @@ -31,109 +31,109 @@ #ifdef BIG_ENDI static const number /* polynomial I */ -/**/ d3 = {0xbfd55555, 0x55555555, }, /* -0.333... */ -/**/ d5 = {0x3fc99999, 0x999997fd, }, /* 0.199... */ -/**/ d7 = {0xbfc24924, 0x923f7603, }, /* -0.142... */ -/**/ d9 = {0x3fbc71c6, 0xe5129a3b, }, /* 0.111... */ -/**/ d11 = {0xbfb74580, 0x22b13c25, }, /* -0.090... */ -/**/ d13 = {0x3fb375f0, 0x8b31cbce, }, /* 0.076... */ +/**/ d3 = {{0xbfd55555, 0x55555555} }, /* -0.333... */ +/**/ d5 = {{0x3fc99999, 0x999997fd} }, /* 0.199... */ +/**/ d7 = {{0xbfc24924, 0x923f7603} }, /* -0.142... */ +/**/ d9 = {{0x3fbc71c6, 0xe5129a3b} }, /* 0.111... */ +/**/ d11 = {{0xbfb74580, 0x22b13c25} }, /* -0.090... */ +/**/ d13 = {{0x3fb375f0, 0x8b31cbce} }, /* 0.076... */ /* polynomial II */ -/**/ f3 = {0xbfd55555, 0x55555555, }, /* -1/3 */ -/**/ ff3 = {0xbc755555, 0x55555555, }, /* -1/3-f3 */ -/**/ f5 = {0x3fc99999, 0x9999999a, }, /* 1/5 */ -/**/ ff5 = {0xbc699999, 0x9999999a, }, /* 1/5-f5 */ -/**/ f7 = {0xbfc24924, 0x92492492, }, /* -1/7 */ -/**/ ff7 = {0xbc624924, 0x92492492, }, /* -1/7-f7 */ -/**/ f9 = {0x3fbc71c7, 0x1c71c71c, }, /* 1/9 */ -/**/ ff9 = {0x3c5c71c7, 0x1c71c71c, }, /* 1/9-f9 */ -/**/ f11 = {0xbfb745d1, 0x745d1746, }, /* -1/11 */ -/**/ f13 = {0x3fb3b13b, 0x13b13b14, }, /* 1/13 */ -/**/ f15 = {0xbfb11111, 0x11111111, }, /* -1/15 */ -/**/ f17 = {0x3fae1e1e, 0x1e1e1e1e, }, /* 1/17 */ -/**/ f19 = {0xbfaaf286, 0xbca1af28, }, /* -1/19 */ +/**/ f3 = {{0xbfd55555, 0x55555555} }, /* -1/3 */ +/**/ ff3 = {{0xbc755555, 0x55555555} }, /* -1/3-f3 */ +/**/ f5 = {{0x3fc99999, 0x9999999a} }, /* 1/5 */ +/**/ ff5 = {{0xbc699999, 0x9999999a} }, /* 1/5-f5 */ +/**/ f7 = {{0xbfc24924, 0x92492492} }, /* -1/7 */ +/**/ ff7 = {{0xbc624924, 0x92492492} }, /* -1/7-f7 */ +/**/ f9 = {{0x3fbc71c7, 0x1c71c71c} }, /* 1/9 */ +/**/ ff9 = {{0x3c5c71c7, 0x1c71c71c} }, /* 1/9-f9 */ +/**/ f11 = {{0xbfb745d1, 0x745d1746} }, /* -1/11 */ +/**/ f13 = {{0x3fb3b13b, 0x13b13b14} }, /* 1/13 */ +/**/ f15 = {{0xbfb11111, 0x11111111} }, /* -1/15 */ +/**/ f17 = {{0x3fae1e1e, 0x1e1e1e1e} }, /* 1/17 */ +/**/ f19 = {{0xbfaaf286, 0xbca1af28} }, /* -1/19 */ /* constants */ -/**/ zero = {0x00000000, 0x00000000, }, /* 0 */ -/**/ one = {0x3ff00000, 0x00000000, }, /* 1 */ -/**/ a = {0x3e4bb67a, 0x00000000, }, /* 1.290e-8 */ -/**/ b = {0x3fb00000, 0x00000000, }, /* 1/16 */ -/**/ c = {0x3ff00000, 0x00000000, }, /* 1 */ -/**/ d = {0x40300000, 0x00000000, }, /* 16 */ -/**/ e = {0x43349ff2, 0x00000000, }, /* 5.805e15 */ -/**/ hpi = {0x3ff921fb, 0x54442d18, }, /* pi/2 */ -/**/ mhpi = {0xbff921fb, 0x54442d18, }, /* -pi/2 */ -/**/ hpi1 = {0x3c91a626, 0x33145c07, }, /* pi/2-hpi */ -/**/ u1 = {0x3c2d3382, 0x00000000, }, /* 7.915e-19 */ -/**/ u21 = {0x3c6dffc0, 0x00000000, }, /* 1.301e-17 */ -/**/ u22 = {0x3c527bd0, 0x00000000, }, /* 4.008e-18 */ -/**/ u23 = {0x3c3cd057, 0x00000000, }, /* 1.562e-18 */ -/**/ u24 = {0x3c329cdf, 0x00000000, }, /* 1.009e-18 */ -/**/ u31 = {0x3c3a1edf, 0x00000000, }, /* 1.416e-18 */ -/**/ u32 = {0x3c33f0e1, 0x00000000, }, /* 1.081e-18 */ -/**/ u4 = {0x3bf955e4, 0x00000000, }, /* 8.584e-20 */ -/**/ u5 = {0x3aaef2d1, 0x00000000, }, /* 5e-26 */ -/**/ u6 = {0x3a98c56d, 0x00000000, }, /* 2.001e-26 */ -/**/ u7 = {0x3a9375de, 0x00000000, }, /* 1.572e-26 */ -/**/ u8 = {0x3a6eeb36, 0x00000000, }, /* 3.122e-27 */ -/**/ u9[M] ={{0x38c1aa5b, 0x00000000, }, /* 2.658e-35 */ -/**/ {0x35c1aa4d, 0x00000000, }, /* 9.443e-50 */ -/**/ {0x32c1aa88, 0x00000000, }, /* 3.355e-64 */ -/**/ {0x11c1aa56, 0x00000000, }},/* 3.818e-223 */ -/**/ two8 = {0x40700000, 0x00000000, }, /* 2**8=256 */ -/**/ two52 = {0x43300000, 0x00000000, }; /* 2**52 */ +/**/ zero = {{0x00000000, 0x00000000} }, /* 0 */ +/**/ one = {{0x3ff00000, 0x00000000} }, /* 1 */ +/**/ a = {{0x3e4bb67a, 0x00000000} }, /* 1.290e-8 */ +/**/ b = {{0x3fb00000, 0x00000000} }, /* 1/16 */ +/**/ c = {{0x3ff00000, 0x00000000} }, /* 1 */ +/**/ d = {{0x40300000, 0x00000000} }, /* 16 */ +/**/ e = {{0x43349ff2, 0x00000000} }, /* 5.805e15 */ +/**/ hpi = {{0x3ff921fb, 0x54442d18} }, /* pi/2 */ +/**/ mhpi = {{0xbff921fb, 0x54442d18} }, /* -pi/2 */ +/**/ hpi1 = {{0x3c91a626, 0x33145c07} }, /* pi/2-hpi */ +/**/ u1 = {{0x3c2d3382, 0x00000000} }, /* 7.915e-19 */ +/**/ u21 = {{0x3c6dffc0, 0x00000000} }, /* 1.301e-17 */ +/**/ u22 = {{0x3c527bd0, 0x00000000} }, /* 4.008e-18 */ +/**/ u23 = {{0x3c3cd057, 0x00000000} }, /* 1.562e-18 */ +/**/ u24 = {{0x3c329cdf, 0x00000000} }, /* 1.009e-18 */ +/**/ u31 = {{0x3c3a1edf, 0x00000000} }, /* 1.416e-18 */ +/**/ u32 = {{0x3c33f0e1, 0x00000000} }, /* 1.081e-18 */ +/**/ u4 = {{0x3bf955e4, 0x00000000} }, /* 8.584e-20 */ +/**/ u5 = {{0x3aaef2d1, 0x00000000} }, /* 5e-26 */ +/**/ u6 = {{0x3a98c56d, 0x00000000} }, /* 2.001e-26 */ +/**/ u7 = {{0x3a9375de, 0x00000000} }, /* 1.572e-26 */ +/**/ u8 = {{0x3a6eeb36, 0x00000000} }, /* 3.122e-27 */ +/**/ u9[M] ={{{0x38c1aa5b, 0x00000000} }, /* 2.658e-35 */ +/**/ {{0x35c1aa4d, 0x00000000} }, /* 9.443e-50 */ +/**/ {{0x32c1aa88, 0x00000000} }, /* 3.355e-64 */ +/**/ {{0x11c1aa56, 0x00000000} }},/* 3.818e-223 */ +/**/ two8 = {{0x40700000, 0x00000000} }, /* 2**8=256 */ +/**/ two52 = {{0x43300000, 0x00000000} }; /* 2**52 */ #else #ifdef LITTLE_ENDI static const number /* polynomial I */ -/**/ d3 = {0x55555555, 0xbfd55555, }, /* -0.333... */ -/**/ d5 = {0x999997fd, 0x3fc99999, }, /* 0.199... */ -/**/ d7 = {0x923f7603, 0xbfc24924, }, /* -0.142... */ -/**/ d9 = {0xe5129a3b, 0x3fbc71c6, }, /* 0.111... */ -/**/ d11 = {0x22b13c25, 0xbfb74580, }, /* -0.090... */ -/**/ d13 = {0x8b31cbce, 0x3fb375f0, }, /* 0.076... */ +/**/ d3 = {{0x55555555, 0xbfd55555} }, /* -0.333... */ +/**/ d5 = {{0x999997fd, 0x3fc99999} }, /* 0.199... */ +/**/ d7 = {{0x923f7603, 0xbfc24924} }, /* -0.142... */ +/**/ d9 = {{0xe5129a3b, 0x3fbc71c6} }, /* 0.111... */ +/**/ d11 = {{0x22b13c25, 0xbfb74580} }, /* -0.090... */ +/**/ d13 = {{0x8b31cbce, 0x3fb375f0} }, /* 0.076... */ /* polynomial II */ -/**/ f3 = {0x55555555, 0xbfd55555, }, /* -1/3 */ -/**/ ff3 = {0x55555555, 0xbc755555, }, /* -1/3-f3 */ -/**/ f5 = {0x9999999a, 0x3fc99999, }, /* 1/5 */ -/**/ ff5 = {0x9999999a, 0xbc699999, }, /* 1/5-f5 */ -/**/ f7 = {0x92492492, 0xbfc24924, }, /* -1/7 */ -/**/ ff7 = {0x92492492, 0xbc624924, }, /* -1/7-f7 */ -/**/ f9 = {0x1c71c71c, 0x3fbc71c7, }, /* 1/9 */ -/**/ ff9 = {0x1c71c71c, 0x3c5c71c7, }, /* 1/9-f9 */ -/**/ f11 = {0x745d1746, 0xbfb745d1, }, /* -1/11 */ -/**/ f13 = {0x13b13b14, 0x3fb3b13b, }, /* 1/13 */ -/**/ f15 = {0x11111111, 0xbfb11111, }, /* -1/15 */ -/**/ f17 = {0x1e1e1e1e, 0x3fae1e1e, }, /* 1/17 */ -/**/ f19 = {0xbca1af28, 0xbfaaf286, }, /* -1/19 */ +/**/ f3 = {{0x55555555, 0xbfd55555} }, /* -1/3 */ +/**/ ff3 = {{0x55555555, 0xbc755555} }, /* -1/3-f3 */ +/**/ f5 = {{0x9999999a, 0x3fc99999} }, /* 1/5 */ +/**/ ff5 = {{0x9999999a, 0xbc699999} }, /* 1/5-f5 */ +/**/ f7 = {{0x92492492, 0xbfc24924} }, /* -1/7 */ +/**/ ff7 = {{0x92492492, 0xbc624924} }, /* -1/7-f7 */ +/**/ f9 = {{0x1c71c71c, 0x3fbc71c7} }, /* 1/9 */ +/**/ ff9 = {{0x1c71c71c, 0x3c5c71c7} }, /* 1/9-f9 */ +/**/ f11 = {{0x745d1746, 0xbfb745d1} }, /* -1/11 */ +/**/ f13 = {{0x13b13b14, 0x3fb3b13b} }, /* 1/13 */ +/**/ f15 = {{0x11111111, 0xbfb11111} }, /* -1/15 */ +/**/ f17 = {{0x1e1e1e1e, 0x3fae1e1e} }, /* 1/17 */ +/**/ f19 = {{0xbca1af28, 0xbfaaf286} }, /* -1/19 */ /* constants */ -/**/ zero = {0x00000000, 0x00000000, }, /* 0 */ -/**/ one = {0x00000000, 0x3ff00000, }, /* 1 */ -/**/ a = {0x00000000, 0x3e4bb67a, }, /* 1.290e-8 */ -/**/ b = {0x00000000, 0x3fb00000, }, /* 1/16 */ -/**/ c = {0x00000000, 0x3ff00000, }, /* 1 */ -/**/ d = {0x00000000, 0x40300000, }, /* 16 */ -/**/ e = {0x00000000, 0x43349ff2, }, /* 5.805e15 */ -/**/ hpi = {0x54442d18, 0x3ff921fb, }, /* pi/2 */ -/**/ mhpi = {0x54442d18, 0xbff921fb, }, /* -pi/2 */ -/**/ hpi1 = {0x33145c07, 0x3c91a626, }, /* pi/2-hpi */ -/**/ u1 = {0x00000000, 0x3c2d3382, }, /* 7.915e-19 */ -/**/ u21 = {0x00000000, 0x3c6dffc0, }, /* 1.301e-17 */ -/**/ u22 = {0x00000000, 0x3c527bd0, }, /* 4.008e-18 */ -/**/ u23 = {0x00000000, 0x3c3cd057, }, /* 1.562e-18 */ -/**/ u24 = {0x00000000, 0x3c329cdf, }, /* 1.009e-18 */ -/**/ u31 = {0x00000000, 0x3c3a1edf, }, /* 1.416e-18 */ -/**/ u32 = {0x00000000, 0x3c33f0e1, }, /* 1.081e-18 */ -/**/ u4 = {0x00000000, 0x3bf955e4, }, /* 8.584e-20 */ -/**/ u5 = {0x00000000, 0x3aaef2d1, }, /* 5e-26 */ -/**/ u6 = {0x00000000, 0x3a98c56d, }, /* 2.001e-26 */ -/**/ u7 = {0x00000000, 0x3a9375de, }, /* 1.572e-26 */ -/**/ u8 = {0x00000000, 0x3a6eeb36, }, /* 3.122e-27 */ -/**/ u9[M] ={{0x00000000, 0x38c1aa5b, }, /* 2.658e-35 */ -/**/ {0x00000000, 0x35c1aa4d, }, /* 9.443e-50 */ -/**/ {0x00000000, 0x32c1aa88, }, /* 3.355e-64 */ -/**/ {0x00000000, 0x11c1aa56, }},/* 3.818e-223 */ -/**/ two8 = {0x00000000, 0x40700000, }, /* 2**8=256 */ -/**/ two52 = {0x00000000, 0x43300000, }; /* 2**52 */ +/**/ zero = {{0x00000000, 0x00000000} }, /* 0 */ +/**/ one = {{0x00000000, 0x3ff00000} }, /* 1 */ +/**/ a = {{0x00000000, 0x3e4bb67a} }, /* 1.290e-8 */ +/**/ b = {{0x00000000, 0x3fb00000} }, /* 1/16 */ +/**/ c = {{0x00000000, 0x3ff00000} }, /* 1 */ +/**/ d = {{0x00000000, 0x40300000} }, /* 16 */ +/**/ e = {{0x00000000, 0x43349ff2} }, /* 5.805e15 */ +/**/ hpi = {{0x54442d18, 0x3ff921fb} }, /* pi/2 */ +/**/ mhpi = {{0x54442d18, 0xbff921fb} }, /* -pi/2 */ +/**/ hpi1 = {{0x33145c07, 0x3c91a626} }, /* pi/2-hpi */ +/**/ u1 = {{0x00000000, 0x3c2d3382} }, /* 7.915e-19 */ +/**/ u21 = {{0x00000000, 0x3c6dffc0} }, /* 1.301e-17 */ +/**/ u22 = {{0x00000000, 0x3c527bd0} }, /* 4.008e-18 */ +/**/ u23 = {{0x00000000, 0x3c3cd057} }, /* 1.562e-18 */ +/**/ u24 = {{0x00000000, 0x3c329cdf} }, /* 1.009e-18 */ +/**/ u31 = {{0x00000000, 0x3c3a1edf} }, /* 1.416e-18 */ +/**/ u32 = {{0x00000000, 0x3c33f0e1} }, /* 1.081e-18 */ +/**/ u4 = {{0x00000000, 0x3bf955e4} }, /* 8.584e-20 */ +/**/ u5 = {{0x00000000, 0x3aaef2d1} }, /* 5e-26 */ +/**/ u6 = {{0x00000000, 0x3a98c56d} }, /* 2.001e-26 */ +/**/ u7 = {{0x00000000, 0x3a9375de} }, /* 1.572e-26 */ +/**/ u8 = {{0x00000000, 0x3a6eeb36} }, /* 3.122e-27 */ +/**/ u9[M] ={{{0x00000000, 0x38c1aa5b} }, /* 2.658e-35 */ +/**/ {{0x00000000, 0x35c1aa4d} }, /* 9.443e-50 */ +/**/ {{0x00000000, 0x32c1aa88} }, /* 3.355e-64 */ +/**/ {{0x00000000, 0x11c1aa56} }},/* 3.818e-223 */ +/**/ two8 = {{0x00000000, 0x40700000} }, /* 2**8=256 */ +/**/ two52 = {{0x00000000, 0x43300000} }; /* 2**52 */ #endif #endif diff --git a/sysdeps/ieee754/dbl-64/atnat2.h b/sysdeps/ieee754/dbl-64/atnat2.h index 967bbc30b9..6972bac71e 100644 --- a/sysdeps/ieee754/dbl-64/atnat2.h +++ b/sysdeps/ieee754/dbl-64/atnat2.h @@ -5,9 +5,9 @@ * * This program is free software; you can redistribute it and/or modify * it under the terms of the GNU Lesser General Public License as published by - * the Free Software Foundation; either version 2 of the License, or + * the Free Software Foundation; either version 2 of the License, or * (at your option) any later version. - * + * * This program is distributed in the hope that it will be useful, * but WITHOUT ANY WARRANTY; without even the implied warranty of * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the @@ -15,7 +15,7 @@ * * You should have received a copy of the GNU Lesser General Public License * along with this program; if not, write to the Free Software - * Foundation, Inc., 59 Temple Place - Suite 330, Boston, MA 02111-1307, USA. + * Foundation, Inc., 59 Temple Place - Suite 330, Boston, MA 02111-1307, USA. */ /************************************************************************/ /* MODULE_NAME: atnat2.h */ @@ -35,140 +35,140 @@ static const number /* polynomial I */ -/**/ d3 = {0xbfd55555, 0x55555555, }, /* -0.333... */ -/**/ d5 = {0x3fc99999, 0x999997fd, }, /* 0.199... */ -/**/ d7 = {0xbfc24924, 0x923f7603, }, /* -0.142... */ -/**/ d9 = {0x3fbc71c6, 0xe5129a3b, }, /* 0.111... */ -/**/ d11 = {0xbfb74580, 0x22b13c25, }, /* -0.090... */ -/**/ d13 = {0x3fb375f0, 0x8b31cbce, }, /* 0.076... */ +/**/ d3 = {{0xbfd55555, 0x55555555} }, /* -0.333... */ +/**/ d5 = {{0x3fc99999, 0x999997fd} }, /* 0.199... */ +/**/ d7 = {{0xbfc24924, 0x923f7603} }, /* -0.142... */ +/**/ d9 = {{0x3fbc71c6, 0xe5129a3b} }, /* 0.111... */ +/**/ d11 = {{0xbfb74580, 0x22b13c25} }, /* -0.090... */ +/**/ d13 = {{0x3fb375f0, 0x8b31cbce} }, /* 0.076... */ /* polynomial II */ -/**/ f3 = {0xbfd55555, 0x55555555, }, /* -1/3 */ -/**/ ff3 = {0xbc755555, 0x55555555, }, /* -1/3-f3 */ -/**/ f5 = {0x3fc99999, 0x9999999a, }, /* 1/5 */ -/**/ ff5 = {0xbc699999, 0x9999999a, }, /* 1/5-f5 */ -/**/ f7 = {0xbfc24924, 0x92492492, }, /* -1/7 */ -/**/ ff7 = {0xbc624924, 0x92492492, }, /* -1/7-f7 */ -/**/ f9 = {0x3fbc71c7, 0x1c71c71c, }, /* 1/9 */ -/**/ ff9 = {0x3c5c71c7, 0x1c71c71c, }, /* 1/9-f9 */ -/**/ f11 = {0xbfb745d1, 0x745d1746, }, /* -1/11 */ -/**/ f13 = {0x3fb3b13b, 0x13b13b14, }, /* 1/13 */ -/**/ f15 = {0xbfb11111, 0x11111111, }, /* -1/15 */ -/**/ f17 = {0x3fae1e1e, 0x1e1e1e1e, }, /* 1/17 */ -/**/ f19 = {0xbfaaf286, 0xbca1af28, }, /* -1/19 */ +/**/ f3 = {{0xbfd55555, 0x55555555} }, /* -1/3 */ +/**/ ff3 = {{0xbc755555, 0x55555555} }, /* -1/3-f3 */ +/**/ f5 = {{0x3fc99999, 0x9999999a} }, /* 1/5 */ +/**/ ff5 = {{0xbc699999, 0x9999999a} }, /* 1/5-f5 */ +/**/ f7 = {{0xbfc24924, 0x92492492} }, /* -1/7 */ +/**/ ff7 = {{0xbc624924, 0x92492492} }, /* -1/7-f7 */ +/**/ f9 = {{0x3fbc71c7, 0x1c71c71c} }, /* 1/9 */ +/**/ ff9 = {{0x3c5c71c7, 0x1c71c71c} }, /* 1/9-f9 */ +/**/ f11 = {{0xbfb745d1, 0x745d1746} }, /* -1/11 */ +/**/ f13 = {{0x3fb3b13b, 0x13b13b14} }, /* 1/13 */ +/**/ f15 = {{0xbfb11111, 0x11111111} }, /* -1/15 */ +/**/ f17 = {{0x3fae1e1e, 0x1e1e1e1e} }, /* 1/17 */ +/**/ f19 = {{0xbfaaf286, 0xbca1af28} }, /* -1/19 */ /* constants */ -/**/ zero = {0x00000000, 0x00000000, }, /* 0 */ -/**/ mzero = {0x80000000, 0x00000000, }, /* -0 */ -/**/ one = {0x3ff00000, 0x00000000, }, /* 1 */ -/**/ inv16 = {0x3fb00000, 0x00000000, }, /* 1/16 */ -/**/ opi = {0x400921fb, 0x54442d18, }, /* pi */ -/**/ opi1 = {0x3ca1a626, 0x33145c07, }, /* pi-opi */ -/**/ mopi = {0xc00921fb, 0x54442d18, }, /* -pi */ -/**/ hpi = {0x3ff921fb, 0x54442d18, }, /* pi/2 */ -/**/ hpi1 = {0x3c91a626, 0x33145c07, }, /* pi/2-hpi */ -/**/ mhpi = {0xbff921fb, 0x54442d18, }, /* -pi/2 */ -/**/ qpi = {0x3fe921fb, 0x54442d18, }, /* pi/4 */ -/**/ qpi1 = {0x3c81a626, 0x33145c07, }, /* pi/4-qpi */ -/**/ mqpi = {0xbfe921fb, 0x54442d18, }, /* -pi/4 */ -/**/ tqpi = {0x4002d97c, 0x7f3321d2, }, /* 3pi/4 */ -/**/ tqpi1 = {0x3c9a7939, 0x4c9e8a0a, }, /* 3pi/4-tqpi */ -/**/ mtqpi = {0xc002d97c, 0x7f3321d2, }, /* -3pi/4 */ -/**/ u1 = {0x3c314c2a, 0x00000000, }, /* 9.377e-19 */ -/**/ u2 = {0x3bf955e4, 0x00000000, }, /* 8.584e-20 */ -/**/ u3 = {0x3bf955e4, 0x00000000, }, /* 8.584e-20 */ -/**/ u4 = {0x3bf955e4, 0x00000000, }, /* 8.584e-20 */ -/**/ u5 = {0x3aaef2d1, 0x00000000, }, /* 5e-26 */ -/**/ u6 = {0x3a6eeb36, 0x00000000, }, /* 3.122e-27 */ -/**/ u7 = {0x3a6eeb36, 0x00000000, }, /* 3.122e-27 */ -/**/ u8 = {0x3a6eeb36, 0x00000000, }, /* 3.122e-27 */ -/**/ u91 = {0x3c6dffc0, 0x00000000, }, /* 1.301e-17 */ -/**/ u92 = {0x3c527bd0, 0x00000000, }, /* 4.008e-18 */ -/**/ u93 = {0x3c3cd057, 0x00000000, }, /* 1.562e-18 */ -/**/ u94 = {0x3c329cdf, 0x00000000, }, /* 1.009e-18 */ -/**/ ua1 = {0x3c3a1edf, 0x00000000, }, /* 1.416e-18 */ -/**/ ua2 = {0x3c33f0e1, 0x00000000, }, /* 1.081e-18 */ -/**/ ub = {0x3a98c56d, 0x00000000, }, /* 2.001e-26 */ -/**/ uc = {0x3a9375de, 0x00000000, }, /* 1.572e-26 */ -/**/ ud[MM] ={{0x38c6eddf, 0x00000000, }, /* 3.450e-35 */ -/**/ {0x35c6ef60, 0x00000000, }, /* 1.226e-49 */ -/**/ {0x32c6ed2f, 0x00000000, }, /* 4.354e-64 */ -/**/ {0x23c6eee8, 0x00000000, }, /* 2.465e-136 */ -/**/ {0x11c6ed16, 0x00000000, }},/* 4.955e-223 */ -/**/ ue = {0x38900e9d, 0x00000000, }, /* 3.02e-36 */ -/**/ two8 = {0x40700000, 0x00000000, }, /* 2**8=256 */ -/**/ two52 = {0x43300000, 0x00000000, }, /* 2**52 */ -/**/ two500 = {0x5f300000, 0x00000000, }, /* 2**500 */ -/**/ twom500 = {0x20b00000, 0x00000000, }, /* 2**(-500) */ -/**/ twom1022 = {0x00100000, 0x00000000, }; /* 2**(-1022) */ +/**/ zero = {{0x00000000, 0x00000000} }, /* 0 */ +/**/ mzero = {{0x80000000, 0x00000000} }, /* -0 */ +/**/ one = {{0x3ff00000, 0x00000000} }, /* 1 */ +/**/ inv16 = {{0x3fb00000, 0x00000000} }, /* 1/16 */ +/**/ opi = {{0x400921fb, 0x54442d18} }, /* pi */ +/**/ opi1 = {{0x3ca1a626, 0x33145c07} }, /* pi-opi */ +/**/ mopi = {{0xc00921fb, 0x54442d18} }, /* -pi */ +/**/ hpi = {{0x3ff921fb, 0x54442d18} }, /* pi/2 */ +/**/ hpi1 = {{0x3c91a626, 0x33145c07} }, /* pi/2-hpi */ +/**/ mhpi = {{0xbff921fb, 0x54442d18} }, /* -pi/2 */ +/**/ qpi = {{0x3fe921fb, 0x54442d18} }, /* pi/4 */ +/**/ qpi1 = {{0x3c81a626, 0x33145c07} }, /* pi/4-qpi */ +/**/ mqpi = {{0xbfe921fb, 0x54442d18} }, /* -pi/4 */ +/**/ tqpi = {{0x4002d97c, 0x7f3321d2} }, /* 3pi/4 */ +/**/ tqpi1 = {{0x3c9a7939, 0x4c9e8a0a} }, /* 3pi/4-tqpi */ +/**/ mtqpi = {{0xc002d97c, 0x7f3321d2} }, /* -3pi/4 */ +/**/ u1 = {{0x3c314c2a, 0x00000000} }, /* 9.377e-19 */ +/**/ u2 = {{0x3bf955e4, 0x00000000} }, /* 8.584e-20 */ +/**/ u3 = {{0x3bf955e4, 0x00000000} }, /* 8.584e-20 */ +/**/ u4 = {{0x3bf955e4, 0x00000000} }, /* 8.584e-20 */ +/**/ u5 = {{0x3aaef2d1, 0x00000000} }, /* 5e-26 */ +/**/ u6 = {{0x3a6eeb36, 0x00000000} }, /* 3.122e-27 */ +/**/ u7 = {{0x3a6eeb36, 0x00000000} }, /* 3.122e-27 */ +/**/ u8 = {{0x3a6eeb36, 0x00000000} }, /* 3.122e-27 */ +/**/ u91 = {{0x3c6dffc0, 0x00000000} }, /* 1.301e-17 */ +/**/ u92 = {{0x3c527bd0, 0x00000000} }, /* 4.008e-18 */ +/**/ u93 = {{0x3c3cd057, 0x00000000} }, /* 1.562e-18 */ +/**/ u94 = {{0x3c329cdf, 0x00000000} }, /* 1.009e-18 */ +/**/ ua1 = {{0x3c3a1edf, 0x00000000} }, /* 1.416e-18 */ +/**/ ua2 = {{0x3c33f0e1, 0x00000000} }, /* 1.081e-18 */ +/**/ ub = {{0x3a98c56d, 0x00000000} }, /* 2.001e-26 */ +/**/ uc = {{0x3a9375de, 0x00000000} }, /* 1.572e-26 */ +/**/ ud[MM] ={{{0x38c6eddf, 0x00000000} }, /* 3.450e-35 */ +/**/ {{0x35c6ef60, 0x00000000} }, /* 1.226e-49 */ +/**/ {{0x32c6ed2f, 0x00000000} }, /* 4.354e-64 */ +/**/ {{0x23c6eee8, 0x00000000} }, /* 2.465e-136 */ +/**/ {{0x11c6ed16, 0x00000000} }},/* 4.955e-223 */ +/**/ ue = {{0x38900e9d, 0x00000000} }, /* 3.02e-36 */ +/**/ two8 = {{0x40700000, 0x00000000} }, /* 2**8=256 */ +/**/ two52 = {{0x43300000, 0x00000000} }, /* 2**52 */ +/**/ two500 = {{0x5f300000, 0x00000000} }, /* 2**500 */ +/**/ twom500 = {{0x20b00000, 0x00000000} }, /* 2**(-500) */ +/**/ twom1022 = {{0x00100000, 0x00000000} }; /* 2**(-1022) */ #else #ifdef LITTLE_ENDI static const number /* polynomial I */ -/**/ d3 = {0x55555555, 0xbfd55555, }, /* -0.333... */ -/**/ d5 = {0x999997fd, 0x3fc99999, }, /* 0.199... */ -/**/ d7 = {0x923f7603, 0xbfc24924, }, /* -0.142... */ -/**/ d9 = {0xe5129a3b, 0x3fbc71c6, }, /* 0.111... */ -/**/ d11 = {0x22b13c25, 0xbfb74580, }, /* -0.090... */ -/**/ d13 = {0x8b31cbce, 0x3fb375f0, }, /* 0.076... */ +/**/ d3 = {{0x55555555, 0xbfd55555} }, /* -0.333... */ +/**/ d5 = {{0x999997fd, 0x3fc99999} }, /* 0.199... */ +/**/ d7 = {{0x923f7603, 0xbfc24924} }, /* -0.142... */ +/**/ d9 = {{0xe5129a3b, 0x3fbc71c6} }, /* 0.111... */ +/**/ d11 = {{0x22b13c25, 0xbfb74580} }, /* -0.090... */ +/**/ d13 = {{0x8b31cbce, 0x3fb375f0} }, /* 0.076... */ /* polynomial II */ -/**/ f3 = {0x55555555, 0xbfd55555, }, /* -1/3 */ -/**/ ff3 = {0x55555555, 0xbc755555, }, /* -1/3-f3 */ -/**/ f5 = {0x9999999a, 0x3fc99999, }, /* 1/5 */ -/**/ ff5 = {0x9999999a, 0xbc699999, }, /* 1/5-f5 */ -/**/ f7 = {0x92492492, 0xbfc24924, }, /* -1/7 */ -/**/ ff7 = {0x92492492, 0xbc624924, }, /* -1/7-f7 */ -/**/ f9 = {0x1c71c71c, 0x3fbc71c7, }, /* 1/9 */ -/**/ ff9 = {0x1c71c71c, 0x3c5c71c7, }, /* 1/9-f9 */ -/**/ f11 = {0x745d1746, 0xbfb745d1, }, /* -1/11 */ -/**/ f13 = {0x13b13b14, 0x3fb3b13b, }, /* 1/13 */ -/**/ f15 = {0x11111111, 0xbfb11111, }, /* -1/15 */ -/**/ f17 = {0x1e1e1e1e, 0x3fae1e1e, }, /* 1/17 */ -/**/ f19 = {0xbca1af28, 0xbfaaf286, }, /* -1/19 */ +/**/ f3 = {{0x55555555, 0xbfd55555} }, /* -1/3 */ +/**/ ff3 = {{0x55555555, 0xbc755555} }, /* -1/3-f3 */ +/**/ f5 = {{0x9999999a, 0x3fc99999} }, /* 1/5 */ +/**/ ff5 = {{0x9999999a, 0xbc699999} }, /* 1/5-f5 */ +/**/ f7 = {{0x92492492, 0xbfc24924} }, /* -1/7 */ +/**/ ff7 = {{0x92492492, 0xbc624924} }, /* -1/7-f7 */ +/**/ f9 = {{0x1c71c71c, 0x3fbc71c7} }, /* 1/9 */ +/**/ ff9 = {{0x1c71c71c, 0x3c5c71c7} }, /* 1/9-f9 */ +/**/ f11 = {{0x745d1746, 0xbfb745d1} }, /* -1/11 */ +/**/ f13 = {{0x13b13b14, 0x3fb3b13b} }, /* 1/13 */ +/**/ f15 = {{0x11111111, 0xbfb11111} }, /* -1/15 */ +/**/ f17 = {{0x1e1e1e1e, 0x3fae1e1e} }, /* 1/17 */ +/**/ f19 = {{0xbca1af28, 0xbfaaf286} }, /* -1/19 */ /* constants */ -/**/ zero = {0x00000000, 0x00000000, }, /* 0 */ -/**/ mzero = {0x00000000, 0x80000000, }, /* -0 */ -/**/ one = {0x00000000, 0x3ff00000, }, /* 1 */ -/**/ inv16 = {0x00000000, 0x3fb00000, }, /* 1/16 */ -/**/ opi = {0x54442d18, 0x400921fb, }, /* pi */ -/**/ opi1 = {0x33145c07, 0x3ca1a626, }, /* pi-opi */ -/**/ mopi = {0x54442d18, 0xc00921fb, }, /* -pi */ -/**/ hpi = {0x54442d18, 0x3ff921fb, }, /* pi/2 */ -/**/ hpi1 = {0x33145c07, 0x3c91a626, }, /* pi/2-hpi */ -/**/ mhpi = {0x54442d18, 0xbff921fb, }, /* -pi/2 */ -/**/ qpi = {0x54442d18, 0x3fe921fb, }, /* pi/4 */ -/**/ qpi1 = {0x33145c07, 0x3c81a626, }, /* pi/4-qpi */ -/**/ mqpi = {0x54442d18, 0xbfe921fb, }, /* -pi/4 */ -/**/ tqpi = {0x7f3321d2, 0x4002d97c, }, /* 3pi/4 */ -/**/ tqpi1 = {0x4c9e8a0a, 0x3c9a7939, }, /* 3pi/4-tqpi */ -/**/ mtqpi = {0x7f3321d2, 0xc002d97c, }, /* -3pi/4 */ -/**/ u1 = {0x00000000, 0x3c314c2a, }, /* 9.377e-19 */ -/**/ u2 = {0x00000000, 0x3bf955e4, }, /* 8.584e-20 */ -/**/ u3 = {0x00000000, 0x3bf955e4, }, /* 8.584e-20 */ -/**/ u4 = {0x00000000, 0x3bf955e4, }, /* 8.584e-20 */ -/**/ u5 = {0x00000000, 0x3aaef2d1, }, /* 5e-26 */ -/**/ u6 = {0x00000000, 0x3a6eeb36, }, /* 3.122e-27 */ -/**/ u7 = {0x00000000, 0x3a6eeb36, }, /* 3.122e-27 */ -/**/ u8 = {0x00000000, 0x3a6eeb36, }, /* 3.122e-27 */ -/**/ u91 = {0x00000000, 0x3c6dffc0, }, /* 1.301e-17 */ -/**/ u92 = {0x00000000, 0x3c527bd0, }, /* 4.008e-18 */ -/**/ u93 = {0x00000000, 0x3c3cd057, }, /* 1.562e-18 */ -/**/ u94 = {0x00000000, 0x3c329cdf, }, /* 1.009e-18 */ -/**/ ua1 = {0x00000000, 0x3c3a1edf, }, /* 1.416e-18 */ -/**/ ua2 = {0x00000000, 0x3c33f0e1, }, /* 1.081e-18 */ -/**/ ub = {0x00000000, 0x3a98c56d, }, /* 2.001e-26 */ -/**/ uc = {0x00000000, 0x3a9375de, }, /* 1.572e-26 */ -/**/ ud[MM] ={{0x00000000, 0x38c6eddf, }, /* 3.450e-35 */ -/**/ {0x00000000, 0x35c6ef60, }, /* 1.226e-49 */ -/**/ {0x00000000, 0x32c6ed2f, }, /* 4.354e-64 */ -/**/ {0x00000000, 0x23c6eee8, }, /* 2.465e-136 */ -/**/ {0x00000000, 0x11c6ed16, }},/* 4.955e-223 */ -/**/ ue = {0x00000000, 0x38900e9d, }, /* 3.02e-36 */ -/**/ two8 = {0x00000000, 0x40700000, }, /* 2**8=256 */ -/**/ two52 = {0x00000000, 0x43300000, }, /* 2**52 */ -/**/ two500 = {0x00000000, 0x5f300000, }, /* 2**500 */ -/**/ twom500 = {0x00000000, 0x20b00000, }, /* 2**(-500) */ -/**/ twom1022 = {0x00000000, 0x00100000, }; /* 2**(-1022) */ +/**/ zero = {{0x00000000, 0x00000000} }, /* 0 */ +/**/ mzero = {{0x00000000, 0x80000000} }, /* -0 */ +/**/ one = {{0x00000000, 0x3ff00000} }, /* 1 */ +/**/ inv16 = {{0x00000000, 0x3fb00000} }, /* 1/16 */ +/**/ opi = {{0x54442d18, 0x400921fb} }, /* pi */ +/**/ opi1 = {{0x33145c07, 0x3ca1a626} }, /* pi-opi */ +/**/ mopi = {{0x54442d18, 0xc00921fb} }, /* -pi */ +/**/ hpi = {{0x54442d18, 0x3ff921fb} }, /* pi/2 */ +/**/ hpi1 = {{0x33145c07, 0x3c91a626} }, /* pi/2-hpi */ +/**/ mhpi = {{0x54442d18, 0xbff921fb} }, /* -pi/2 */ +/**/ qpi = {{0x54442d18, 0x3fe921fb} }, /* pi/4 */ +/**/ qpi1 = {{0x33145c07, 0x3c81a626} }, /* pi/4-qpi */ +/**/ mqpi = {{0x54442d18, 0xbfe921fb} }, /* -pi/4 */ +/**/ tqpi = {{0x7f3321d2, 0x4002d97c} }, /* 3pi/4 */ +/**/ tqpi1 = {{0x4c9e8a0a, 0x3c9a7939} }, /* 3pi/4-tqpi */ +/**/ mtqpi = {{0x7f3321d2, 0xc002d97c} }, /* -3pi/4 */ +/**/ u1 = {{0x00000000, 0x3c314c2a} }, /* 9.377e-19 */ +/**/ u2 = {{0x00000000, 0x3bf955e4} }, /* 8.584e-20 */ +/**/ u3 = {{0x00000000, 0x3bf955e4} }, /* 8.584e-20 */ +/**/ u4 = {{0x00000000, 0x3bf955e4} }, /* 8.584e-20 */ +/**/ u5 = {{0x00000000, 0x3aaef2d1} }, /* 5e-26 */ +/**/ u6 = {{0x00000000, 0x3a6eeb36} }, /* 3.122e-27 */ +/**/ u7 = {{0x00000000, 0x3a6eeb36} }, /* 3.122e-27 */ +/**/ u8 = {{0x00000000, 0x3a6eeb36} }, /* 3.122e-27 */ +/**/ u91 = {{0x00000000, 0x3c6dffc0} }, /* 1.301e-17 */ +/**/ u92 = {{0x00000000, 0x3c527bd0} }, /* 4.008e-18 */ +/**/ u93 = {{0x00000000, 0x3c3cd057} }, /* 1.562e-18 */ +/**/ u94 = {{0x00000000, 0x3c329cdf} }, /* 1.009e-18 */ +/**/ ua1 = {{0x00000000, 0x3c3a1edf} }, /* 1.416e-18 */ +/**/ ua2 = {{0x00000000, 0x3c33f0e1} }, /* 1.081e-18 */ +/**/ ub = {{0x00000000, 0x3a98c56d} }, /* 2.001e-26 */ +/**/ uc = {{0x00000000, 0x3a9375de} }, /* 1.572e-26 */ +/**/ ud[MM] ={{{0x00000000, 0x38c6eddf} }, /* 3.450e-35 */ +/**/ {{0x00000000, 0x35c6ef60} }, /* 1.226e-49 */ +/**/ {{0x00000000, 0x32c6ed2f} }, /* 4.354e-64 */ +/**/ {{0x00000000, 0x23c6eee8} }, /* 2.465e-136 */ +/**/ {{0x00000000, 0x11c6ed16} }},/* 4.955e-223 */ +/**/ ue = {{0x00000000, 0x38900e9d} }, /* 3.02e-36 */ +/**/ two8 = {{0x00000000, 0x40700000} }, /* 2**8=256 */ +/**/ two52 = {{0x00000000, 0x43300000} }, /* 2**52 */ +/**/ two500 = {{0x00000000, 0x5f300000} }, /* 2**500 */ +/**/ twom500 = {{0x00000000, 0x20b00000} }, /* 2**(-500) */ +/**/ twom1022 = {{0x00000000, 0x00100000} }; /* 2**(-1022) */ #endif #endif diff --git a/sysdeps/ieee754/dbl-64/branred.c b/sysdeps/ieee754/dbl-64/branred.c index 09b55602bf..2bceb947c0 100644 --- a/sysdeps/ieee754/dbl-64/branred.c +++ b/sysdeps/ieee754/dbl-64/branred.c @@ -5,9 +5,9 @@ * * This program is free software; you can redistribute it and/or modify * it under the terms of the GNU Lesser General Public License as published by - * the Free Software Foundation; either version 2 of the License, or + * the Free Software Foundation; either version 2 of the License, or * (at your option) any later version. - * + * * This program is distributed in the hope that it will be useful, * but WITHOUT ANY WARRANTY; without even the implied warranty of * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the @@ -15,14 +15,14 @@ * * You should have received a copy of the GNU Lesser General Public License * along with this program; if not, write to the Free Software - * Foundation, Inc., 59 Temple Place - Suite 330, Boston, MA 02111-1307, USA. + * Foundation, Inc., 59 Temple Place - Suite 330, Boston, MA 02111-1307, USA. */ /*******************************************************************/ /* */ -/* MODULE_NAME: branred.c */ +/* MODULE_NAME: branred.c */ /* */ /* FUNCTIONS: branred */ -/* */ +/* */ /* FILES NEEDED: branred.h mydefs.h endian.h mpa.h */ /* mha.c */ /* */ @@ -45,35 +45,41 @@ /* x=n*pi/2+(a+aa), abs(a+aa)>20)&2047; + k = (u.i[HIGH_HALF]>>20)&2047; k = (k-450)/24; if (k<0) k=0; gor.x = t576.x; gor.i[HIGH_HALF] -= ((k*24)<<20); - for (i=0;i<6;i++) + for (i=0;i<6;i++) { r[i] = x1*toverp[k+i]*gor.x; gor.x *= tm24.x; } for (i=0;i<3;i++) { - s=(r[i]+big.x)-big.x; + s=(r[i]+big.x)-big.x; sum+=s; r[i]-=s; } t=0; - for (i=0;i<6;i++) + for (i=0;i<6;i++) t+=r[5-i]; - bb=(((((r[0]-t)+r[1])+r[2])+r[3])+r[4])+r[5]; + bb=(((((r[0]-t)+r[1])+r[2])+r[3])+r[4])+r[5]; s=(t+big.x)-big.x; sum+=s; t-=s; @@ -85,7 +91,7 @@ int branred(double x, double *a, double *aa) bb1=bb; sum1=sum; sum=0; - + u.x = x2; k = (u.i[HIGH_HALF]>>20)&2047; k = (k-450)/24; @@ -96,14 +102,14 @@ int branred(double x, double *a, double *aa) for (i=0;i<6;i++) { r[i] = x2*toverp[k+i]*gor.x; gor.x *= tm24.x; } for (i=0;i<3;i++) { - s=(r[i]+big.x)-big.x; + s=(r[i]+big.x)-big.x; sum+=s; r[i]-=s; } t=0; for (i=0;i<6;i++) t+=r[5-i]; - bb=(((((r[0]-t)+r[1])+r[2])+r[3])+r[4])+r[5]; + bb=(((((r[0]-t)+r[1])+r[2])+r[3])+r[4])+r[5]; s=(t+big.x)-big.x; sum+=s; t-=s; @@ -111,17 +117,17 @@ int branred(double x, double *a, double *aa) bb=(t-b)+bb; s=(sum+big1.x)-big1.x; sum-=s; - + b2=b; bb2=bb; sum2=sum; - + sum=sum1+sum2; b=b1+b2; - bb = (ABS(b1)>ABS(b2))? (b1-b)+b2 : (b2-b)+b1; - if (b > 0.5) + bb = (ABS(b1)>ABS(b2))? (b1-b)+b2 : (b2-b)+b1; + if (b > 0.5) {b-=1.0; sum+=1.0;} - else if (b < -0.5) + else if (b < -0.5) {b+=1.0; sum-=1.0;} s=b+(bb+bb1+bb2); t=((b-s)+bb)+(bb1+bb2); @@ -136,4 +142,3 @@ int branred(double x, double *a, double *aa) *aa=t; return ((int) sum)&3; /* return quater of unit circle */ } - diff --git a/sysdeps/ieee754/dbl-64/branred.h b/sysdeps/ieee754/dbl-64/branred.h index 2258e58eb1..f5e9439e5a 100644 --- a/sysdeps/ieee754/dbl-64/branred.h +++ b/sysdeps/ieee754/dbl-64/branred.h @@ -5,9 +5,9 @@ * * This program is free software; you can redistribute it and/or modify * it under the terms of the GNU Lesser General Public License as published by - * the Free Software Foundation; either version 2 of the License, or + * the Free Software Foundation; either version 2 of the License, or * (at your option) any later version. - * + * * This program is distributed in the hope that it will be useful, * but WITHOUT ANY WARRANTY; without even the implied warranty of * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the @@ -15,7 +15,7 @@ * * You should have received a copy of the GNU Lesser General Public License * along with this program; if not, write to the Free Software - * Foundation, Inc., 59 Temple Place - Suite 330, Boston, MA 02111-1307, USA. + * Foundation, Inc., 59 Temple Place - Suite 330, Boston, MA 02111-1307, USA. */ /************************************************************************/ /* MODULE_NAME: branred.h */ @@ -31,49 +31,49 @@ #ifdef BIG_ENDI static const mynumber -/**/ t576 = {0x63f00000, 0x00000000 }, /* 2 ^ 576 */ -/**/ tm600 = {0x1a700000, 0x00000000 }, /* 2 ^- 600 */ -/**/ tm24 = {0x3e700000, 0x00000000 }, /* 2 ^- 24 */ -/**/ big = {0x43380000, 0x00000000 }, /* 6755399441055744 */ -/**/ big1 = {0x43580000, 0x00000000 }, /* 27021597764222976 */ -/**/ hp0 = {0x3FF921FB, 0x54442D18 } ,/* 1.5707963267948966 */ -/**/ hp1 = {0x3C91A626, 0x33145C07 } ,/* 6.123233995736766e-17 */ -/**/ mp1 = {0x3FF921FB, 0x58000000 }, /* 1.5707963407039642 */ -/**/ mp2 = {0xBE4DDE97, 0x40000000 }; /*-1.3909067675399456e-08 */ +/**/ t576 = {{0x63f00000, 0x00000000}}, /* 2 ^ 576 */ +/**/ tm600 = {{0x1a700000, 0x00000000}}, /* 2 ^- 600 */ +/**/ tm24 = {{0x3e700000, 0x00000000}}, /* 2 ^- 24 */ +/**/ big = {{0x43380000, 0x00000000}}, /* 6755399441055744 */ +/**/ big1 = {{0x43580000, 0x00000000}}, /* 27021597764222976 */ +/**/ hp0 = {{0x3FF921FB, 0x54442D18}} ,/* 1.5707963267948966 */ +/**/ hp1 = {{0x3C91A626, 0x33145C07}} ,/* 6.123233995736766e-17 */ +/**/ mp1 = {{0x3FF921FB, 0x58000000}}, /* 1.5707963407039642 */ +/**/ mp2 = {{0xBE4DDE97, 0x40000000}}; /*-1.3909067675399456e-08 */ #else #ifdef LITTLE_ENDI static const mynumber -/**/ t576 = {0x00000000, 0x63f00000 }, /* 2 ^ 576 */ -/**/ tm600 = {0x00000000, 0x1a700000 }, /* 2 ^- 600 */ -/**/ tm24 = {0x00000000, 0x3e700000 }, /* 2 ^- 24 */ -/**/ big = {0x00000000, 0x43380000 }, /* 6755399441055744 */ -/**/ big1 = {0x00000000, 0x43580000 }, /* 27021597764222976 */ -/**/ hp0 = {0x54442D18, 0x3FF921FB }, /* 1.5707963267948966 */ -/**/ hp1 = {0x33145C07, 0x3C91A626 }, /* 6.123233995736766e-17 */ -/**/ mp1 = {0x58000000, 0x3FF921FB }, /* 1.5707963407039642 */ -/**/ mp2 = {0x40000000, 0xBE4DDE97 }; /*-1.3909067675399456e-08 */ +/**/ t576 = {{0x00000000, 0x63f00000}}, /* 2 ^ 576 */ +/**/ tm600 = {{0x00000000, 0x1a700000}}, /* 2 ^- 600 */ +/**/ tm24 = {{0x00000000, 0x3e700000}}, /* 2 ^- 24 */ +/**/ big = {{0x00000000, 0x43380000}}, /* 6755399441055744 */ +/**/ big1 = {{0x00000000, 0x43580000}}, /* 27021597764222976 */ +/**/ hp0 = {{0x54442D18, 0x3FF921FB}}, /* 1.5707963267948966 */ +/**/ hp1 = {{0x33145C07, 0x3C91A626}}, /* 6.123233995736766e-17 */ +/**/ mp1 = {{0x58000000, 0x3FF921FB}}, /* 1.5707963407039642 */ +/**/ mp2 = {{0x40000000, 0xBE4DDE97}}; /*-1.3909067675399456e-08 */ #endif #endif static const double toverp[75] = { /* 2/ PI base 24*/ - 10680707.0, 7228996.0, 1387004.0, 2578385.0, 16069853.0, - 12639074.0, 9804092.0, 4427841.0, 16666979.0, 11263675.0, - 12935607.0, 2387514.0, 4345298.0, 14681673.0, 3074569.0, - 13734428.0, 16653803.0, 1880361.0, 10960616.0, 8533493.0, - 3062596.0, 8710556.0, 7349940.0, 6258241.0, 3772886.0, - 3769171.0, 3798172.0, 8675211.0, 12450088.0, 3874808.0, - 9961438.0, 366607.0, 15675153.0, 9132554.0, 7151469.0, - 3571407.0, 2607881.0, 12013382.0, 4155038.0, 6285869.0, - 7677882.0, 13102053.0, 15825725.0, 473591.0, 9065106.0, - 15363067.0, 6271263.0, 9264392.0, 5636912.0, 4652155.0, - 7056368.0, 13614112.0, 10155062.0, 1944035.0, 9527646.0, - 15080200.0, 6658437.0, 6231200.0, 6832269.0, 16767104.0, - 5075751.0, 3212806.0, 1398474.0, 7579849.0, 6349435.0, - 12618859.0, 4703257.0, 12806093.0, 14477321.0, 2786137.0, - 12875403.0, 9837734.0, 14528324.0, 13719321.0, 343717.0 }; + 10680707.0, 7228996.0, 1387004.0, 2578385.0, 16069853.0, + 12639074.0, 9804092.0, 4427841.0, 16666979.0, 11263675.0, + 12935607.0, 2387514.0, 4345298.0, 14681673.0, 3074569.0, + 13734428.0, 16653803.0, 1880361.0, 10960616.0, 8533493.0, + 3062596.0, 8710556.0, 7349940.0, 6258241.0, 3772886.0, + 3769171.0, 3798172.0, 8675211.0, 12450088.0, 3874808.0, + 9961438.0, 366607.0, 15675153.0, 9132554.0, 7151469.0, + 3571407.0, 2607881.0, 12013382.0, 4155038.0, 6285869.0, + 7677882.0, 13102053.0, 15825725.0, 473591.0, 9065106.0, + 15363067.0, 6271263.0, 9264392.0, 5636912.0, 4652155.0, + 7056368.0, 13614112.0, 10155062.0, 1944035.0, 9527646.0, + 15080200.0, 6658437.0, 6231200.0, 6832269.0, 16767104.0, + 5075751.0, 3212806.0, 1398474.0, 7579849.0, 6349435.0, + 12618859.0, 4703257.0, 12806093.0, 14477321.0, 2786137.0, + 12875403.0, 9837734.0, 14528324.0, 13719321.0, 343717.0 }; static const double split = 134217729.0; diff --git a/sysdeps/ieee754/dbl-64/doasin.h b/sysdeps/ieee754/dbl-64/doasin.h index 01f44e0adb..078716dc60 100644 --- a/sysdeps/ieee754/dbl-64/doasin.h +++ b/sysdeps/ieee754/dbl-64/doasin.h @@ -5,9 +5,9 @@ * * This program is free software; you can redistribute it and/or modify * it under the terms of the GNU Lesser General Public License as published by - * the Free Software Foundation; either version 2 of the License, or + * the Free Software Foundation; either version 2 of the License, or * (at your option) any later version. - * + * * This program is distributed in the hope that it will be useful, * but WITHOUT ANY WARRANTY; without even the implied warranty of * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the @@ -15,7 +15,7 @@ * * You should have received a copy of the GNU Lesser General Public License * along with this program; if not, write to the Free Software - * Foundation, Inc., 59 Temple Place - Suite 330, Boston, MA 02111-1307, USA. + * Foundation, Inc., 59 Temple Place - Suite 330, Boston, MA 02111-1307, USA. */ /************************************************************************/ /* MODULE_NAME: doasin.h */ @@ -31,30 +31,30 @@ #ifdef BIG_ENDI - static const mynumber -/**/ c1 = {0x3FC55555, 0x55555555}, /* 0.16666666666666666 */ -/**/ cc1 = {0x3C655555, 0x55775389}, /* 9.2518585419753846e-18 */ -/**/ c2 = {0x3FB33333, 0x33333333}, /* 0.074999999999999997 */ -/**/ cc2 = {0x3C499993, 0x63F1A115}, /* 2.7755472886508899e-18 */ -/**/ c3 = {0x3FA6DB6D, 0xB6DB6DB7}, /* 0.044642857142857144 */ -/**/ cc3 = {0xBC320FC0, 0x3D5CF0C5}, /* -9.7911734574147224e-19 */ -/**/ c4 = {0x3F9F1C71, 0xC71C71C5}, /* 0.030381944444444437 */ -/**/ cc4 = {0xBC02B240, 0xFF23ED1E}; /* -1.2669108566898312e-19 */ + static const mynumber +/**/ c1 = {{0x3FC55555, 0x55555555}}, /* 0.16666666666666666 */ +/**/ cc1 = {{0x3C655555, 0x55775389}}, /* 9.2518585419753846e-18 */ +/**/ c2 = {{0x3FB33333, 0x33333333}}, /* 0.074999999999999997 */ +/**/ cc2 = {{0x3C499993, 0x63F1A115}}, /* 2.7755472886508899e-18 */ +/**/ c3 = {{0x3FA6DB6D, 0xB6DB6DB7}}, /* 0.044642857142857144 */ +/**/ cc3 = {{0xBC320FC0, 0x3D5CF0C5}}, /* -9.7911734574147224e-19 */ +/**/ c4 = {{0x3F9F1C71, 0xC71C71C5}}, /* 0.030381944444444437 */ +/**/ cc4 = {{0xBC02B240, 0xFF23ED1E}}; /* -1.2669108566898312e-19 */ #else #ifdef LITTLE_ENDI - static const mynumber -/**/ c1 = {0x55555555, 0x3FC55555}, /* 0.16666666666666666 */ -/**/ cc1 = {0x55775389, 0x3C655555}, /* 9.2518585419753846e-18 */ -/**/ c2 = {0x33333333, 0x3FB33333}, /* 0.074999999999999997 */ -/**/ cc2 = {0x63F1A115, 0x3C499993}, /* 2.7755472886508899e-18 */ -/**/ c3 = {0xB6DB6DB7, 0x3FA6DB6D}, /* 0.044642857142857144 */ -/**/ cc3 = {0x3D5CF0C5, 0xBC320FC0}, /* -9.7911734574147224e-19 */ -/**/ c4 = {0xC71C71C5, 0x3F9F1C71}, /* 0.030381944444444437 */ -/**/ cc4 = {0xFF23ED1E, 0xBC02B240}; /* -1.2669108566898312e-19 */ - - + static const mynumber +/**/ c1 = {{0x55555555, 0x3FC55555}}, /* 0.16666666666666666 */ +/**/ cc1 = {{0x55775389, 0x3C655555}}, /* 9.2518585419753846e-18 */ +/**/ c2 = {{0x33333333, 0x3FB33333}}, /* 0.074999999999999997 */ +/**/ cc2 = {{0x63F1A115, 0x3C499993}}, /* 2.7755472886508899e-18 */ +/**/ c3 = {{0xB6DB6DB7, 0x3FA6DB6D}}, /* 0.044642857142857144 */ +/**/ cc3 = {{0x3D5CF0C5, 0xBC320FC0}}, /* -9.7911734574147224e-19 */ +/**/ c4 = {{0xC71C71C5, 0x3F9F1C71}}, /* 0.030381944444444437 */ +/**/ cc4 = {{0xFF23ED1E, 0xBC02B240}}; /* -1.2669108566898312e-19 */ + + #endif #endif diff --git a/sysdeps/ieee754/dbl-64/dosincos.c b/sysdeps/ieee754/dbl-64/dosincos.c index 2a5f8fcc34..692fc85889 100644 --- a/sysdeps/ieee754/dbl-64/dosincos.c +++ b/sysdeps/ieee754/dbl-64/dosincos.c @@ -5,9 +5,9 @@ * * This program is free software; you can redistribute it and/or modify * it under the terms of the GNU Lesser General Public License as published by - * the Free Software Foundation; either version 2 of the License, or + * the Free Software Foundation; either version 2 of the License, or * (at your option) any later version. - * + * * This program is distributed in the hope that it will be useful, * but WITHOUT ANY WARRANTY; without even the implied warranty of * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the @@ -15,10 +15,10 @@ * * You should have received a copy of the GNU Lesser General Public License * along with this program; if not, write to the Free Software - * Foundation, Inc., 59 Temple Place - Suite 330, Boston, MA 02111-1307, USA. + * Foundation, Inc., 59 Temple Place - Suite 330, Boston, MA 02111-1307, USA. */ /********************************************************************/ -/* */ +/* */ /* MODULE_NAME: dosincos.c */ /* */ /* */ @@ -34,7 +34,7 @@ #include "endian.h" -#include "mydefs.h" +#include "mydefs.h" #include "sincos.tbl" #include "dla.h" #include "dosincos.h" @@ -42,15 +42,18 @@ /* Routine receive Double-Length number (x+dx) and computing sin(x+dx) */ /* as Double-Length number and store it at array v .It computes it by */ /* arithmetic action on Double-Length numbers */ -/*(x+dx) between 0 and PI/4 */ +/*(x+dx) between 0 and PI/4 */ /***********************************************************************/ void dubsin(double x, double dx, double v[]) { - double xx,y,yy,z,zz,r,s,p,hx,tx,hy,ty,q,c,cc,d,dd,d2,dd2,e,ee, + double r,s,p,hx,tx,hy,ty,q,c,cc,d,dd,d2,dd2,e,ee, sn,ssn,cs,ccs,ds,dss,dc,dcc; +#if 0 + double xx,y,yy,z,zz; +#endif mynumber u; int4 k; - + u.x=x+big.x; k = u.i[LOW_HALF]<<2; x=x-(u.x-big.x); @@ -63,26 +66,26 @@ void dubsin(double x, double dx, double v[]) { cs=sincos.x[k+2]; /* */ ccs=sincos.x[k+3]; /* */ MUL2(d2,dd2,s7.x,ss7.x,ds,dss,p,hx,tx,hy,ty,q,c,cc); /* Taylor */ - ADD2(ds,dss,s5.x,ss5.x,ds,dss,r,s); + ADD2(ds,dss,s5.x,ss5.x,ds,dss,r,s); MUL2(d2,dd2,ds,dss,ds,dss,p,hx,tx,hy,ty,q,c,cc); /* series */ - ADD2(ds,dss,s3.x,ss3.x,ds,dss,r,s); + ADD2(ds,dss,s3.x,ss3.x,ds,dss,r,s); MUL2(d2,dd2,ds,dss,ds,dss,p,hx,tx,hy,ty,q,c,cc); /* for sin */ MUL2(d,dd,ds,dss,ds,dss,p,hx,tx,hy,ty,q,c,cc); ADD2(ds,dss,d,dd,ds,dss,r,s); /* ds=sin(t) */ - + MUL2(d2,dd2,c8.x,cc8.x,dc,dcc,p,hx,tx,hy,ty,q,c,cc); ;/* Taylor */ - ADD2(dc,dcc,c6.x,cc6.x,dc,dcc,r,s); + ADD2(dc,dcc,c6.x,cc6.x,dc,dcc,r,s); MUL2(d2,dd2,dc,dcc,dc,dcc,p,hx,tx,hy,ty,q,c,cc); /* series */ - ADD2(dc,dcc,c4.x,cc4.x,dc,dcc,r,s); + ADD2(dc,dcc,c4.x,cc4.x,dc,dcc,r,s); MUL2(d2,dd2,dc,dcc,dc,dcc,p,hx,tx,hy,ty,q,c,cc); /* for cos */ - ADD2(dc,dcc,c2.x,cc2.x,dc,dcc,r,s); + ADD2(dc,dcc,c2.x,cc2.x,dc,dcc,r,s); MUL2(d2,dd2,dc,dcc,dc,dcc,p,hx,tx,hy,ty,q,c,cc); /* dc=cos(t) */ - + MUL2(cs,ccs,ds,dss,e,ee,p,hx,tx,hy,ty,q,c,cc); MUL2(dc,dcc,sn,ssn,dc,dcc,p,hx,tx,hy,ty,q,c,cc); - SUB2(e,ee,dc,dcc,e,ee,r,s); + SUB2(e,ee,dc,dcc,e,ee,r,s); ADD2(e,ee,sn,ssn,e,ee,r,s); /* e+ee=sin(x+dx) */ - + v[0]=e; v[1]=ee; } @@ -94,46 +97,49 @@ void dubsin(double x, double dx, double v[]) { /**********************************************************************/ void dubcos(double x, double dx, double v[]) { - double xx,y,yy,z,zz,r,s,p,hx,tx,hy,ty,q,c,cc,d,dd,d2,dd2,e,ee, + double r,s,p,hx,tx,hy,ty,q,c,cc,d,dd,d2,dd2,e,ee, sn,ssn,cs,ccs,ds,dss,dc,dcc; +#if 0 + double xx,y,yy,z,zz; +#endif mynumber u; - int4 k; - u.x=x+big.x; + int4 k; + u.x=x+big.x; k = u.i[LOW_HALF]<<2; x=x-(u.x-big.x); - d=x+dx; + d=x+dx; dd=(x-d)+dx; /* cos(x+dx)=cos(Xi+t)=cos(Xi)cos(t) - sin(Xi)sin(t) */ - MUL2(d,dd,d,dd,d2,dd2,p,hx,tx,hy,ty,q,c,cc); + MUL2(d,dd,d,dd,d2,dd2,p,hx,tx,hy,ty,q,c,cc); sn=sincos.x[k]; /* */ ssn=sincos.x[k+1]; /* sin(Xi) and cos(Xi) */ cs=sincos.x[k+2]; /* */ ccs=sincos.x[k+3]; /* */ MUL2(d2,dd2,s7.x,ss7.x,ds,dss,p,hx,tx,hy,ty,q,c,cc); - ADD2(ds,dss,s5.x,ss5.x,ds,dss,r,s); + ADD2(ds,dss,s5.x,ss5.x,ds,dss,r,s); MUL2(d2,dd2,ds,dss,ds,dss,p,hx,tx,hy,ty,q,c,cc); - ADD2(ds,dss,s3.x,ss3.x,ds,dss,r,s); + ADD2(ds,dss,s3.x,ss3.x,ds,dss,r,s); MUL2(d2,dd2,ds,dss,ds,dss,p,hx,tx,hy,ty,q,c,cc); MUL2(d,dd,ds,dss,ds,dss,p,hx,tx,hy,ty,q,c,cc); - ADD2(ds,dss,d,dd,ds,dss,r,s); - + ADD2(ds,dss,d,dd,ds,dss,r,s); + MUL2(d2,dd2,c8.x,cc8.x,dc,dcc,p,hx,tx,hy,ty,q,c,cc); - ADD2(dc,dcc,c6.x,cc6.x,dc,dcc,r,s); + ADD2(dc,dcc,c6.x,cc6.x,dc,dcc,r,s); MUL2(d2,dd2,dc,dcc,dc,dcc,p,hx,tx,hy,ty,q,c,cc); - ADD2(dc,dcc,c4.x,cc4.x,dc,dcc,r,s); + ADD2(dc,dcc,c4.x,cc4.x,dc,dcc,r,s); MUL2(d2,dd2,dc,dcc,dc,dcc,p,hx,tx,hy,ty,q,c,cc); - ADD2(dc,dcc,c2.x,cc2.x,dc,dcc,r,s); + ADD2(dc,dcc,c2.x,cc2.x,dc,dcc,r,s); MUL2(d2,dd2,dc,dcc,dc,dcc,p,hx,tx,hy,ty,q,c,cc); - + MUL2(cs,ccs,ds,dss,e,ee,p,hx,tx,hy,ty,q,c,cc); MUL2(dc,dcc,sn,ssn,dc,dcc,p,hx,tx,hy,ty,q,c,cc); MUL2(d2,dd2,s7.x,ss7.x,ds,dss,p,hx,tx,hy,ty,q,c,cc); ADD2(ds,dss,s5.x,ss5.x,ds,dss,r,s); MUL2(d2,dd2,ds,dss,ds,dss,p,hx,tx,hy,ty,q,c,cc); - ADD2(ds,dss,s3.x,ss3.x,ds,dss,r,s); + ADD2(ds,dss,s3.x,ss3.x,ds,dss,r,s); MUL2(d2,dd2,ds,dss,ds,dss,p,hx,tx,hy,ty,q,c,cc); - MUL2(d,dd,ds,dss,ds,dss,p,hx,tx,hy,ty,q,c,cc); - ADD2(ds,dss,d,dd,ds,dss,r,s); + MUL2(d,dd,ds,dss,ds,dss,p,hx,tx,hy,ty,q,c,cc); + ADD2(ds,dss,d,dd,ds,dss,r,s); MUL2(d2,dd2,c8.x,cc8.x,dc,dcc,p,hx,tx,hy,ty,q,c,cc); ADD2(dc,dcc,c6.x,cc6.x,dc,dcc,r,s); MUL2(d2,dd2,dc,dcc,dc,dcc,p,hx,tx,hy,ty,q,c,cc); @@ -143,9 +149,9 @@ void dubcos(double x, double dx, double v[]) { MUL2(d2,dd2,dc,dcc,dc,dcc,p,hx,tx,hy,ty,q,c,cc); MUL2(sn,ssn,ds,dss,e,ee,p,hx,tx,hy,ty,q,c,cc); MUL2(dc,dcc,cs,ccs,dc,dcc,p,hx,tx,hy,ty,q,c,cc); - ADD2(e,ee,dc,dcc,e,ee,r,s); - SUB2(cs,ccs,e,ee,e,ee,r,s); - + ADD2(e,ee,dc,dcc,e,ee,r,s); + SUB2(cs,ccs,e,ee,e,ee,r,s); + v[0]=e; v[1]=ee; } @@ -164,14 +170,14 @@ void docos(double x, double dx, double v[]) { yy=hp1.x-yy; y=p+yy; yy=(p-y)+yy; - if (y>0) {dubsin(y,yy,w); v[0]=w[0]; v[1]=w[1];} + if (y>0) {dubsin(y,yy,w); v[0]=w[0]; v[1]=w[1];} /* cos(x) = sin ( 90 - x ) */ - else {dubsin(-y,-yy,w); v[0]=-w[0]; v[1]=-w[1]; + else {dubsin(-y,-yy,w); v[0]=-w[0]; v[1]=-w[1]; } } else { /* y>= 3/4 * PI */ p=2.0*hp0.x-y; /* p = PI- y */ - yy=2.0*hp1.x-yy; + yy=2.0*hp1.x-yy; y=p+yy; yy=(p-y)+yy; dubcos(y,yy,w); diff --git a/sysdeps/ieee754/dbl-64/dosincos.h b/sysdeps/ieee754/dbl-64/dosincos.h index 20f4ede78c..ce7df09dfa 100644 --- a/sysdeps/ieee754/dbl-64/dosincos.h +++ b/sysdeps/ieee754/dbl-64/dosincos.h @@ -5,9 +5,9 @@ * * This program is free software; you can redistribute it and/or modify * it under the terms of the GNU Lesser General Public License as published by - * the Free Software Foundation; either version 2 of the License, or + * the Free Software Foundation; either version 2 of the License, or * (at your option) any later version. - * + * * This program is distributed in the hope that it will be useful, * but WITHOUT ANY WARRANTY; without even the implied warranty of * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the @@ -15,7 +15,7 @@ * * You should have received a copy of the GNU Lesser General Public License * along with this program; if not, write to the Free Software - * Foundation, Inc., 59 Temple Place - Suite 330, Boston, MA 02111-1307, USA. + * Foundation, Inc., 59 Temple Place - Suite 330, Boston, MA 02111-1307, USA. */ /************************************************************************/ /* MODULE_NAME: dosincos.h */ @@ -32,47 +32,47 @@ #ifdef BIG_ENDI static const mynumber -/**/ s3 = {0xBFC55555, 0x55555555}, /* -0.16666666666666666 */ -/**/ ss3 = {0xBC6553AA, 0xE77EE482}, /* -9.2490366677784492e-18 */ -/**/ s5 = {0x3F811111, 0x11110F15}, /* 0.008333333333332452 */ -/**/ ss5 = {0xBC21AC06, 0xDA488820}, /* -4.7899996586987931e-19 */ -/**/ s7 = {0xBF2A019F, 0x5816C78D}, /* -0.00019841261022928957 */ -/**/ ss7 = {0x3BCDCEC9, 0x6A18BF2A}, /* 1.2624077757871259e-20 */ -/**/ c2 = {0x3FE00000, 0x00000000}, /* 0.5 */ -/**/ cc2 = {0xBA282FD8, 0x00000000}, /* -1.5264073330037701e-28 */ -/**/ c4 = {0xBFA55555, 0x55555555}, /* -0.041666666666666664 */ -/**/ cc4 = {0xBC4554BC, 0x2FFF257E}, /* -2.312711276085743e-18 */ -/**/ c6 = {0x3F56C16C, 0x16C16A96}, /* 0.0013888888888888055 */ -/**/ cc6 = {0xBBD2E846, 0xE6346F14}, /* -1.6015133010194884e-20 */ -/**/ c8 = {0xBEFA019F, 0x821D5987}, /* -2.480157866754367e-05 */ -/**/ cc8 = {0x3B7AB71E, 0x72FFE5CC}, /* 3.5357416224857556e-22 */ +/**/ s3 = {{0xBFC55555, 0x55555555}},/* -0.16666666666666666 */ +/**/ ss3 = {{0xBC6553AA, 0xE77EE482}},/* -9.2490366677784492e-18 */ +/**/ s5 = {{0x3F811111, 0x11110F15}},/* 0.008333333333332452 */ +/**/ ss5 = {{0xBC21AC06, 0xDA488820}},/* -4.7899996586987931e-19 */ +/**/ s7 = {{0xBF2A019F, 0x5816C78D}},/* -0.00019841261022928957 */ +/**/ ss7 = {{0x3BCDCEC9, 0x6A18BF2A}},/* 1.2624077757871259e-20 */ +/**/ c2 = {{0x3FE00000, 0x00000000}},/* 0.5 */ +/**/ cc2 = {{0xBA282FD8, 0x00000000}},/* -1.5264073330037701e-28 */ +/**/ c4 = {{0xBFA55555, 0x55555555}},/* -0.041666666666666664 */ +/**/ cc4 = {{0xBC4554BC, 0x2FFF257E}},/* -2.312711276085743e-18 */ +/**/ c6 = {{0x3F56C16C, 0x16C16A96}},/* 0.0013888888888888055 */ +/**/ cc6 = {{0xBBD2E846, 0xE6346F14}},/* -1.6015133010194884e-20 */ +/**/ c8 = {{0xBEFA019F, 0x821D5987}},/* -2.480157866754367e-05 */ +/**/ cc8 = {{0x3B7AB71E, 0x72FFE5CC}},/* 3.5357416224857556e-22 */ -/**/ big = {0x42c80000, 0x00000000}, /* 52776558133248 */ +/**/ big = {{0x42c80000, 0x00000000}}, /* 52776558133248 */ -/**/ hp0 = {0x3FF921FB, 0x54442D18 }, /* PI / 2 */ -/**/ hp1 = {0x3C91A626, 0x33145C07 }; /* 6.123233995736766e-17 */ +/**/ hp0 = {{0x3FF921FB, 0x54442D18}}, /* PI / 2 */ +/**/ hp1 = {{0x3C91A626, 0x33145C07}}; /* 6.123233995736766e-17 */ #else #ifdef LITTLE_ENDI static const mynumber -/**/ s3 = {0x55555555, 0xBFC55555}, /* -0.16666666666666666 */ -/**/ ss3 = {0xE77EE482, 0xBC6553AA}, /* -9.2490366677784492e-18 */ -/**/ s5 = {0x11110F15, 0x3F811111}, /* 0.008333333333332452 */ -/**/ ss5 = {0xDA488820, 0xBC21AC06}, /* -4.7899996586987931e-19 */ -/**/ s7 = {0x5816C78D, 0xBF2A019F}, /* -0.00019841261022928957 */ -/**/ ss7 = {0x6A18BF2A, 0x3BCDCEC9}, /* 1.2624077757871259e-20 */ -/**/ c2 = {0x00000000, 0x3FE00000}, /* 0.5 */ -/**/ cc2 = {0x00000000, 0xBA282FD8}, /* -1.5264073330037701e-28 */ -/**/ c4 = {0x55555555, 0xBFA55555}, /* -0.041666666666666664 */ -/**/ cc4 = {0x2FFF257E, 0xBC4554BC}, /* -2.312711276085743e-18 */ -/**/ c6 = {0x16C16A96, 0x3F56C16C}, /* 0.0013888888888888055 */ -/**/ cc6 = {0xE6346F14, 0xBBD2E846}, /* -1.6015133010194884e-20 */ -/**/ c8 = {0x821D5987, 0xBEFA019F}, /* -2.480157866754367e-05 */ -/**/ cc8 = {0x72FFE5CC, 0x3B7AB71E}, /* 3.5357416224857556e-22 */ +/**/ s3 = {{0x55555555, 0xBFC55555}},/* -0.16666666666666666 */ +/**/ ss3 = {{0xE77EE482, 0xBC6553AA}},/* -9.2490366677784492e-18 */ +/**/ s5 = {{0x11110F15, 0x3F811111}},/* 0.008333333333332452 */ +/**/ ss5 = {{0xDA488820, 0xBC21AC06}},/* -4.7899996586987931e-19 */ +/**/