1*25c28e83SPiotr Jasiukajtis /* 2*25c28e83SPiotr Jasiukajtis * CDDL HEADER START 3*25c28e83SPiotr Jasiukajtis * 4*25c28e83SPiotr Jasiukajtis * The contents of this file are subject to the terms of the 5*25c28e83SPiotr Jasiukajtis * Common Development and Distribution License (the "License"). 6*25c28e83SPiotr Jasiukajtis * You may not use this file except in compliance with the License. 7*25c28e83SPiotr Jasiukajtis * 8*25c28e83SPiotr Jasiukajtis * You can obtain a copy of the license at usr/src/OPENSOLARIS.LICENSE 9*25c28e83SPiotr Jasiukajtis * or http://www.opensolaris.org/os/licensing. 10*25c28e83SPiotr Jasiukajtis * See the License for the specific language governing permissions 11*25c28e83SPiotr Jasiukajtis * and limitations under the License. 12*25c28e83SPiotr Jasiukajtis * 13*25c28e83SPiotr Jasiukajtis * When distributing Covered Code, include this CDDL HEADER in each 14*25c28e83SPiotr Jasiukajtis * file and include the License file at usr/src/OPENSOLARIS.LICENSE. 15*25c28e83SPiotr Jasiukajtis * If applicable, add the following below this CDDL HEADER, with the 16*25c28e83SPiotr Jasiukajtis * fields enclosed by brackets "[]" replaced with your own identifying 17*25c28e83SPiotr Jasiukajtis * information: Portions Copyright [yyyy] [name of copyright owner] 18*25c28e83SPiotr Jasiukajtis * 19*25c28e83SPiotr Jasiukajtis * CDDL HEADER END 20*25c28e83SPiotr Jasiukajtis */ 21*25c28e83SPiotr Jasiukajtis /* 22*25c28e83SPiotr Jasiukajtis * Copyright 2011 Nexenta Systems, Inc. All rights reserved. 23*25c28e83SPiotr Jasiukajtis */ 24*25c28e83SPiotr Jasiukajtis /* 25*25c28e83SPiotr Jasiukajtis * Copyright 2005 Sun Microsystems, Inc. All rights reserved. 26*25c28e83SPiotr Jasiukajtis * Use is subject to license terms. 27*25c28e83SPiotr Jasiukajtis */ 28*25c28e83SPiotr Jasiukajtis 29*25c28e83SPiotr Jasiukajtis #pragma weak tanf = __tanf 30*25c28e83SPiotr Jasiukajtis 31*25c28e83SPiotr Jasiukajtis #include "libm.h" 32*25c28e83SPiotr Jasiukajtis 33*25c28e83SPiotr Jasiukajtis extern const int _TBL_ipio2_inf[]; 34*25c28e83SPiotr Jasiukajtis extern int __rem_pio2m(double *, double *, int, int, int, const int *); 35*25c28e83SPiotr Jasiukajtis #if defined(__i386) && !defined(__amd64) 36*25c28e83SPiotr Jasiukajtis extern int __swapRP(int); 37*25c28e83SPiotr Jasiukajtis #endif 38*25c28e83SPiotr Jasiukajtis 39*25c28e83SPiotr Jasiukajtis static const double C[] = { 40*25c28e83SPiotr Jasiukajtis 1.0, 41*25c28e83SPiotr Jasiukajtis 4.46066928428959230679140546271810308098793029785e-0003, 42*25c28e83SPiotr Jasiukajtis 4.92165316309189027066395283327437937259674072266e+0000, 43*25c28e83SPiotr Jasiukajtis -7.11410648161473480044492134766187518835067749023e-0001, 44*25c28e83SPiotr Jasiukajtis 4.08549808374053391446523164631798863410949707031e+0000, 45*25c28e83SPiotr Jasiukajtis 2.50411070398050927821032018982805311679840087891e+0000, 46*25c28e83SPiotr Jasiukajtis 1.11492064560251158411574579076841473579406738281e+0001, 47*25c28e83SPiotr Jasiukajtis -1.50565540968422650891511693771462887525558471680e+0000, 48*25c28e83SPiotr Jasiukajtis -1.81484378878349295050043110677506774663925170898e+0000, 49*25c28e83SPiotr Jasiukajtis 3.333335997532835641297409611782510896641e-0001, 50*25c28e83SPiotr Jasiukajtis 2.999997598248363761541668282006867229939e+00, 51*25c28e83SPiotr Jasiukajtis 0.636619772367581343075535, /* 2^ -1 * 1.45F306DC9C883 */ 52*25c28e83SPiotr Jasiukajtis 0.5, 53*25c28e83SPiotr Jasiukajtis 1.570796326734125614166, /* 2^ 0 * 1.921FB54400000 */ 54*25c28e83SPiotr Jasiukajtis 6.077100506506192601475e-11, /* 2^-34 * 1.0B4611A626331 */ 55*25c28e83SPiotr Jasiukajtis }; 56*25c28e83SPiotr Jasiukajtis 57*25c28e83SPiotr Jasiukajtis #define one C[0] 58*25c28e83SPiotr Jasiukajtis #define P0 C[1] 59*25c28e83SPiotr Jasiukajtis #define P1 C[2] 60*25c28e83SPiotr Jasiukajtis #define P2 C[3] 61*25c28e83SPiotr Jasiukajtis #define P3 C[4] 62*25c28e83SPiotr Jasiukajtis #define P4 C[5] 63*25c28e83SPiotr Jasiukajtis #define P5 C[6] 64*25c28e83SPiotr Jasiukajtis #define P6 C[7] 65*25c28e83SPiotr Jasiukajtis #define P7 C[8] 66*25c28e83SPiotr Jasiukajtis #define T0 C[9] 67*25c28e83SPiotr Jasiukajtis #define T1 C[10] 68*25c28e83SPiotr Jasiukajtis #define invpio2 C[11] 69*25c28e83SPiotr Jasiukajtis #define half C[12] 70*25c28e83SPiotr Jasiukajtis #define pio2_1 C[13] 71*25c28e83SPiotr Jasiukajtis #define pio2_t C[14] 72*25c28e83SPiotr Jasiukajtis 73*25c28e83SPiotr Jasiukajtis float 74*25c28e83SPiotr Jasiukajtis tanf(float x) 75*25c28e83SPiotr Jasiukajtis { 76*25c28e83SPiotr Jasiukajtis double y, z, w; 77*25c28e83SPiotr Jasiukajtis float f; 78*25c28e83SPiotr Jasiukajtis int n, ix, hx, hy; 79*25c28e83SPiotr Jasiukajtis volatile int i; 80*25c28e83SPiotr Jasiukajtis 81*25c28e83SPiotr Jasiukajtis hx = *((int *)&x); 82*25c28e83SPiotr Jasiukajtis ix = hx & 0x7fffffff; 83*25c28e83SPiotr Jasiukajtis 84*25c28e83SPiotr Jasiukajtis y = (double)x; 85*25c28e83SPiotr Jasiukajtis 86*25c28e83SPiotr Jasiukajtis if (ix <= 0x4016cbe4) { /* |x| < 3*pi/4 */ 87*25c28e83SPiotr Jasiukajtis if (ix <= 0x3f490fdb) { /* |x| < pi/4 */ 88*25c28e83SPiotr Jasiukajtis if (ix < 0x3c000000) { /* |x| < 2**-7 */ 89*25c28e83SPiotr Jasiukajtis if (ix <= 0x39800000) { /* |x| < 2**-12 */ 90*25c28e83SPiotr Jasiukajtis i = (int)y; 91*25c28e83SPiotr Jasiukajtis #ifdef lint 92*25c28e83SPiotr Jasiukajtis i = i; 93*25c28e83SPiotr Jasiukajtis #endif 94*25c28e83SPiotr Jasiukajtis return (x); 95*25c28e83SPiotr Jasiukajtis } 96*25c28e83SPiotr Jasiukajtis return ((float)((y * T0) * (T1 + y * y))); 97*25c28e83SPiotr Jasiukajtis } 98*25c28e83SPiotr Jasiukajtis z = y * y; 99*25c28e83SPiotr Jasiukajtis return ((float)(((P0 * y) * (P1 + z * (P2 + z)) * 100*25c28e83SPiotr Jasiukajtis (P3 + z * (P4 + z))) * 101*25c28e83SPiotr Jasiukajtis (P5 + z * (P6 + z * (P7 + z))))); 102*25c28e83SPiotr Jasiukajtis } 103*25c28e83SPiotr Jasiukajtis if (hx > 0) 104*25c28e83SPiotr Jasiukajtis y = (y - pio2_1) - pio2_t; 105*25c28e83SPiotr Jasiukajtis else 106*25c28e83SPiotr Jasiukajtis y = (y + pio2_1) + pio2_t; 107*25c28e83SPiotr Jasiukajtis hy = ((int *)&y)[HIWORD] & ~0x80000000; 108*25c28e83SPiotr Jasiukajtis if (hy < 0x3f800000) { /* |y| < 2**-7 */ 109*25c28e83SPiotr Jasiukajtis z = (y * T0) * (T1 + y * y); 110*25c28e83SPiotr Jasiukajtis return ((float)(-one / z)); 111*25c28e83SPiotr Jasiukajtis } 112*25c28e83SPiotr Jasiukajtis z = y * y; 113*25c28e83SPiotr Jasiukajtis w = ((P0 * y) * (P1 + z * (P2 + z)) * (P3 + z * (P4 + z))) * 114*25c28e83SPiotr Jasiukajtis (P5 + z * (P6 + z * (P7 + z))); 115*25c28e83SPiotr Jasiukajtis return ((float)(-one / w)); 116*25c28e83SPiotr Jasiukajtis } 117*25c28e83SPiotr Jasiukajtis 118*25c28e83SPiotr Jasiukajtis if (ix <= 0x49c90fdb) { /* |x| < 2^19*pi */ 119*25c28e83SPiotr Jasiukajtis #if defined(__i386) && !defined(__amd64) 120*25c28e83SPiotr Jasiukajtis int rp; 121*25c28e83SPiotr Jasiukajtis 122*25c28e83SPiotr Jasiukajtis rp = __swapRP(fp_extended); 123*25c28e83SPiotr Jasiukajtis #endif 124*25c28e83SPiotr Jasiukajtis w = y * invpio2; 125*25c28e83SPiotr Jasiukajtis if (hx < 0) 126*25c28e83SPiotr Jasiukajtis n = (int)(w - half); 127*25c28e83SPiotr Jasiukajtis else 128*25c28e83SPiotr Jasiukajtis n = (int)(w + half); 129*25c28e83SPiotr Jasiukajtis y = (y - n * pio2_1) - n * pio2_t; 130*25c28e83SPiotr Jasiukajtis #if defined(__i386) && !defined(__amd64) 131*25c28e83SPiotr Jasiukajtis if (rp != fp_extended) 132*25c28e83SPiotr Jasiukajtis (void) __swapRP(rp); 133*25c28e83SPiotr Jasiukajtis #endif 134*25c28e83SPiotr Jasiukajtis } else { 135*25c28e83SPiotr Jasiukajtis if (ix >= 0x7f800000) 136*25c28e83SPiotr Jasiukajtis return (x / x); /* sin(Inf or NaN) is NaN */ 137*25c28e83SPiotr Jasiukajtis hy = ((int *)&y)[HIWORD]; 138*25c28e83SPiotr Jasiukajtis n = ((hy >> 20) & 0x7ff) - 1046; 139*25c28e83SPiotr Jasiukajtis ((int *)&w)[HIWORD] = (hy & 0xfffff) | 0x41600000; 140*25c28e83SPiotr Jasiukajtis ((int *)&w)[LOWORD] = ((int *)&y)[LOWORD]; 141*25c28e83SPiotr Jasiukajtis n = __rem_pio2m(&w, &y, n, 1, 0, _TBL_ipio2_inf); 142*25c28e83SPiotr Jasiukajtis if (hy < 0) { 143*25c28e83SPiotr Jasiukajtis y = -y; 144*25c28e83SPiotr Jasiukajtis n = -n; 145*25c28e83SPiotr Jasiukajtis } 146*25c28e83SPiotr Jasiukajtis } 147*25c28e83SPiotr Jasiukajtis 148*25c28e83SPiotr Jasiukajtis hy = ((int *)&y)[HIWORD] & ~0x80000000; 149*25c28e83SPiotr Jasiukajtis if (hy < 0x3f800000) { /* |y| < 2**-7 */ 150*25c28e83SPiotr Jasiukajtis z = (y * T0) * (T1 + y * y); 151*25c28e83SPiotr Jasiukajtis f = ((n & 1) == 0)? (float)z : (float)(-one / z); 152*25c28e83SPiotr Jasiukajtis return (f); 153*25c28e83SPiotr Jasiukajtis } 154*25c28e83SPiotr Jasiukajtis z = y * y; 155*25c28e83SPiotr Jasiukajtis w = ((P0 * y) * (P1 + z * (P2 + z)) * (P3 + z * (P4 + z))) * 156*25c28e83SPiotr Jasiukajtis (P5 + z * (P6 + z * (P7 + z))); 157*25c28e83SPiotr Jasiukajtis f = ((n & 1) == 0)? (float)w : (float)(-one / w); 158*25c28e83SPiotr Jasiukajtis return (f); 159*25c28e83SPiotr Jasiukajtis } 160