9  * or http://www.opensolaris.org/os/licensing.
  10  * See the License for the specific language governing permissions
  11  * and limitations under the License.
  12  *
  13  * When distributing Covered Code, include this CDDL HEADER in each
  14  * file and include the License file at usr/src/OPENSOLARIS.LICENSE.
  15  * If applicable, add the following below this CDDL HEADER, with the
  16  * fields enclosed by brackets "[]" replaced with your own identifying
  17  * information: Portions Copyright [yyyy] [name of copyright owner]
  18  *
  19  * CDDL HEADER END
  20  */
  21 /*
  22  * Copyright 2011 Nexenta Systems, Inc.  All rights reserved.
  23  */
  24 /*
  25  * Copyright 2005 Sun Microsystems, Inc.  All rights reserved.
  26  * Use is subject to license terms.
  27  */
  28 
  29 #pragma weak sincos = __sincos
  30 
  31 /* INDENT OFF */
  32 /*
  33  * sincos(x,s,c)
  34  * Accurate Table look-up algorithm by K.C. Ng, 2000.
  35  *
  36  * 1. Reduce x to x>0 by cos(-x)=cos(x), sin(-x)=-sin(x).
  37  * 2. For 0<= x < 8, let i = (64*x chopped)-10. Let d = x - a[i], where
  38  *    a[i] is a double that is close to (i+10.5)/64 (and hence |d|< 10.5/64)
  39  *    and such that sin(a[i]) and cos(a[i]) is close to a double (with error
  40  *    less than 2**-8 ulp). Then
  41  *
  42  *      cos(x) = cos(a[i]+d) = cos(a[i])cos(d) - sin(a[i])*sin(d)
  43  *             = TBL_cos_a[i]*(1+QQ1*d^2+QQ2*d^4) -
  44  *                      TBL_sin_a[i]*(d+PP1*d^3+PP2*d^5)
  45  *             = TBL_cos_a[i] + (TBL_cos_a[i]*d^2*(QQ1+QQ2*d^2) -
  46  *                      TBL_sin_a[i]*(d+PP1*d^3+PP2*d^5))
  47  *
  48  *      sin(x) = sin(a[i]+d) = sin(a[i])cos(d) + cos(a[i])*sin(d)
  49  *             = TBL_sin_a[i]*(1+QQ1*d^2+QQ2*d^4) +

   9  * or http://www.opensolaris.org/os/licensing.
  10  * See the License for the specific language governing permissions
  11  * and limitations under the License.
  12  *
  13  * When distributing Covered Code, include this CDDL HEADER in each
  14  * file and include the License file at usr/src/OPENSOLARIS.LICENSE.
  15  * If applicable, add the following below this CDDL HEADER, with the
  16  * fields enclosed by brackets "[]" replaced with your own identifying
  17  * information: Portions Copyright [yyyy] [name of copyright owner]
  18  *
  19  * CDDL HEADER END
  20  */
  21 /*
  22  * Copyright 2011 Nexenta Systems, Inc.  All rights reserved.
  23  */
  24 /*
  25  * Copyright 2005 Sun Microsystems, Inc.  All rights reserved.
  26  * Use is subject to license terms.
  27  */
  28 
  29 #pragma weak __sincos = sincos
  30 
  31 /* INDENT OFF */
  32 /*
  33  * sincos(x,s,c)
  34  * Accurate Table look-up algorithm by K.C. Ng, 2000.
  35  *
  36  * 1. Reduce x to x>0 by cos(-x)=cos(x), sin(-x)=-sin(x).
  37  * 2. For 0<= x < 8, let i = (64*x chopped)-10. Let d = x - a[i], where
  38  *    a[i] is a double that is close to (i+10.5)/64 (and hence |d|< 10.5/64)
  39  *    and such that sin(a[i]) and cos(a[i]) is close to a double (with error
  40  *    less than 2**-8 ulp). Then
  41  *
  42  *      cos(x) = cos(a[i]+d) = cos(a[i])cos(d) - sin(a[i])*sin(d)
  43  *             = TBL_cos_a[i]*(1+QQ1*d^2+QQ2*d^4) -
  44  *                      TBL_sin_a[i]*(d+PP1*d^3+PP2*d^5)
  45  *             = TBL_cos_a[i] + (TBL_cos_a[i]*d^2*(QQ1+QQ2*d^2) -
  46  *                      TBL_sin_a[i]*(d+PP1*d^3+PP2*d^5))
  47  *
  48  *      sin(x) = sin(a[i]+d) = sin(a[i])cos(d) + cos(a[i])*sin(d)
  49  *             = TBL_sin_a[i]*(1+QQ1*d^2+QQ2*d^4) +