diff options
Diffstat (limited to 'src/math/sincos_386.s')
| -rw-r--r-- | src/math/sincos_386.s | 28 |
1 files changed, 28 insertions, 0 deletions
diff --git a/src/math/sincos_386.s b/src/math/sincos_386.s new file mode 100644 index 000000000..83af5016e --- /dev/null +++ b/src/math/sincos_386.s @@ -0,0 +1,28 @@ +// Copyright 2010 The Go Authors. All rights reserved. +// Use of this source code is governed by a BSD-style +// license that can be found in the LICENSE file. + +#include "textflag.h" + +// func Sincos(x float64) (sin, cos float64) +TEXT ·Sincos(SB),NOSPLIT,$0 + FMOVD x+0(FP), F0 // F0=x + FSINCOS // F0=cos(x), F1=sin(x) if -2**63 < x < 2**63 + FSTSW AX // AX=status word + ANDW $0x0400, AX + JNE 4(PC) // jump if x outside range + FMOVDP F0, cos+16(FP) // F0=sin(x) + FMOVDP F0, sin+8(FP) + RET + FLDPI // F0=Pi, F1=x + FADDD F0, F0 // F0=2*Pi, F1=x + FXCHD F0, F1 // F0=x, F1=2*Pi + FPREM1 // F0=reduced_x, F1=2*Pi + FSTSW AX // AX=status word + ANDW $0x0400, AX + JNE -3(PC) // jump if reduction incomplete + FMOVDP F0, F1 // F0=reduced_x + FSINCOS // F0=cos(reduced_x), F1=sin(reduced_x) + FMOVDP F0, cos+16(FP) // F0=sin(reduced_x) + FMOVDP F0, sin+8(FP) + RET |