1// Copyright 2009 The Go Authors. All rights reserved.
2// Use of this source code is governed by a BSD-style
3// license that can be found in the LICENSE file.
4
5#include "textflag.h"
6
7// func Cos(x float64) float64
8TEXT ·Cos(SB),NOSPLIT,$0
9	FMOVD   x+0(FP), F0  // F0=x
10	FCOS                 // F0=cos(x) if -2**63 < x < 2**63
11	FSTSW   AX           // AX=status word
12	ANDW    $0x0400, AX
13	JNE     3(PC)        // jump if x outside range
14	FMOVDP  F0, ret+8(FP)
15	RET
16	FLDPI                // F0=Pi, F1=x
17	FADDD   F0, F0       // F0=2*Pi, F1=x
18	FXCHD   F0, F1       // F0=x, F1=2*Pi
19	FPREM1               // F0=reduced_x, F1=2*Pi
20	FSTSW   AX           // AX=status word
21	ANDW    $0x0400, AX
22	JNE     -3(PC)       // jump if reduction incomplete
23	FMOVDP  F0, F1       // F0=reduced_x
24	FCOS                 // F0=cos(reduced_x)
25	FMOVDP  F0, ret+8(FP)
26	RET
27
28// func Sin(x float64) float64
29TEXT ·Sin(SB),NOSPLIT,$0
30	FMOVD   x+0(FP), F0  // F0=x
31	FSIN                 // F0=sin(x) if -2**63 < x < 2**63
32	FSTSW   AX           // AX=status word
33	ANDW    $0x0400, AX
34	JNE     3(PC)        // jump if x outside range
35	FMOVDP  F0, ret+8(FP)
36	RET
37	FLDPI                // F0=Pi, F1=x
38	FADDD   F0, F0       // F0=2*Pi, F1=x
39	FXCHD   F0, F1       // F0=x, F1=2*Pi
40	FPREM1               // F0=reduced_x, F1=2*Pi
41	FSTSW   AX           // AX=status word
42	ANDW    $0x0400, AX
43	JNE     -3(PC)       // jump if reduction incomplete
44	FMOVDP  F0, F1       // F0=reduced_x
45	FSIN                 // F0=sin(reduced_x)
46	FMOVDP  F0, ret+8(FP)
47	RET
48