Book a Demo!
CoCalc Logo Icon
StoreFeaturesDocsShareSupportNewsAboutPoliciesSign UpSign In
wine-mirror
GitHub Repository: wine-mirror/wine
 Path: blob/master/libs/musl/src/math/__sindf.c
4397 views
1
/* origin: FreeBSD /usr/src/lib/msun/src/k_sinf.c */

2
/*

3
 * Conversion to float by Ian Lance Taylor, Cygnus Support, [email protected].

4
 * Optimized by Bruce D. Evans.

5
 */

6
/*

7
 * ====================================================

8
 * Copyright (C) 1993 by Sun Microsystems, Inc. All rights reserved.

9
 *

10
 * Developed at SunPro, a Sun Microsystems, Inc. business.

11
 * Permission to use, copy, modify, and distribute this

12
 * software is freely granted, provided that this notice

13
 * is preserved.

14
 * ====================================================

15
 */

16


17
#include "libm.h"

18


19
/* |sin(x)/x - s(x)| < 2**-37.5 (~[-4.89e-12, 4.824e-12]). */

20
static const double

21
S1 = -0x1.5555555555555p-3,

22
S2 =  0x1.1111111111111p-7,

23
S3 = -0x1.a01a01a01a01ap-13,

24
S4 =  0x1.71de3a556c734p-19;

25


26
float __sindf(double x)

27
{

28
	double_t r, s, w, z;

29


30
	/* Try to optimize for parallel evaluation as in __tandf.c. */

31
	z = x*x;

32
	if (x > -7.8175831586122513e-03 && x < 7.8175831586122513e-03)

33
		return x * (1 + S1 * z);

34
	w = z*z;

35
	r = S3 + z*S4;

36
	s = z*x;

37
	return (x + s*(S1 + z*S2)) + s*w*r;

38
}

39


40