diff options
| author | Keith Packard <keithp@keithp.com> | 2018-08-05 12:23:22 +0800 | 
|---|---|---|
| committer | Keith Packard <keithp@keithp.com> | 2018-10-13 08:21:58 -0700 | 
| commit | 7c04888cf9809e0c73f0813c74e8dd972facde3a (patch) | |
| tree | ebb43ee5aa61cc65ce28595c19b4dd3c19812ead /src/math/kf_sin.c | |
| parent | 56629222711ba3ef7853405c9b07ad614fb29b95 (diff) | |
altos: Switch to newlib-nano for libc on arm
Stop using pdclib
Signed-off-by: Keith Packard <keithp@keithp.com>
Diffstat (limited to 'src/math/kf_sin.c')
| -rw-r--r-- | src/math/kf_sin.c | 49 | 
1 files changed, 0 insertions, 49 deletions
diff --git a/src/math/kf_sin.c b/src/math/kf_sin.c deleted file mode 100644 index e81fa0bd..00000000 --- a/src/math/kf_sin.c +++ /dev/null @@ -1,49 +0,0 @@ -/* kf_sin.c -- float version of k_sin.c - * Conversion to float by Ian Lance Taylor, Cygnus Support, ian@cygnus.com. - */ - -/* - * ==================================================== - * Copyright (C) 1993 by Sun Microsystems, Inc. All rights reserved. - * - * Developed at SunPro, a Sun Microsystems, Inc. business. - * Permission to use, copy, modify, and distribute this - * software is freely granted, provided that this notice  - * is preserved. - * ==================================================== - */ - -#include "fdlibm.h" - -#ifdef __STDC__ -static const float  -#else -static float  -#endif -half =  5.0000000000e-01,/* 0x3f000000 */ -S1  = -1.6666667163e-01, /* 0xbe2aaaab */ -S2  =  8.3333337680e-03, /* 0x3c088889 */ -S3  = -1.9841270114e-04, /* 0xb9500d01 */ -S4  =  2.7557314297e-06, /* 0x3638ef1b */ -S5  = -2.5050759689e-08, /* 0xb2d72f34 */ -S6  =  1.5896910177e-10; /* 0x2f2ec9d3 */ - -#ifdef __STDC__ -	float __kernel_sinf(float x, float y, int iy) -#else -	float __kernel_sinf(x, y, iy) -	float x,y; int iy;		/* iy=0 if y is zero */ -#endif -{ -	float z,r,v; -	__int32_t ix; -	GET_FLOAT_WORD(ix,x); -	ix &= 0x7fffffff;			/* high word of x */ -	if(ix<0x32000000)			/* |x| < 2**-27 */ -	   {if((int)x==0) return x;}		/* generate inexact */ -	z	=  x*x; -	v	=  z*x; -	r	=  S2+z*(S3+z*(S4+z*(S5+z*S6))); -	if(iy==0) return x+v*(S1+z*r); -	else      return x-((z*(half*y-v*r)-y)-v*S1); -}  | 
