2010-02-19 00:33:15 -07:00
|
|
|
// 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.
|
|
|
|
|
2013-08-12 11:25:18 -06:00
|
|
|
#include "../../cmd/ld/textflag.h"
|
|
|
|
|
2010-02-19 00:33:15 -07:00
|
|
|
// func Sincos(x float64) (sin, cos float64)
|
2013-08-12 11:25:18 -06:00
|
|
|
TEXT ·Sincos(SB),NOSPLIT,$0
|
2010-02-19 00:33:15 -07:00
|
|
|
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
|
2012-09-20 22:35:56 -06:00
|
|
|
FMOVDP F0, cos+16(FP) // F0=sin(x)
|
|
|
|
FMOVDP F0, sin+8(FP)
|
2010-02-19 00:33:15 -07:00
|
|
|
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)
|
2012-09-20 22:35:56 -06:00
|
|
|
FMOVDP F0, cos+16(FP) // F0=sin(reduced_x)
|
|
|
|
FMOVDP F0, sin+8(FP)
|
2010-02-19 00:33:15 -07:00
|
|
|
RET
|