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 8 TEXT 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 29 TEXT 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
