diff options
Diffstat (limited to 'src/pkg/math/exp2_386.s')
-rw-r--r-- | src/pkg/math/exp2_386.s | 38 |
1 files changed, 38 insertions, 0 deletions
diff --git a/src/pkg/math/exp2_386.s b/src/pkg/math/exp2_386.s new file mode 100644 index 000000000..ed82a4dd3 --- /dev/null +++ b/src/pkg/math/exp2_386.s @@ -0,0 +1,38 @@ +// 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. + +// func Exp2(x float64) float64 +TEXT ·Exp2(SB),7,$0 +// test bits for not-finite + MOVL x+4(FP), AX + ANDL $0x7ff00000, AX + CMPL AX, $0x7ff00000 + JEQ not_finite + FMOVD x+0(FP), F0 // F0=x + FMOVD F0, F1 // F0=x, F1=x + FRNDINT // F0=int(x), F1=x + FSUBD F0, F1 // F0=int(x), F1=x-int(x) + FXCHD F0, F1 // F0=x-int(x), F1=int(x) + F2XM1 // F0=2**(x-int(x))-1, F1=int(x) + FLD1 // F0=1, F1=2**(x-int(x))-1, F2=int(x) + FADDDP F0, F1 // F0=2**(x-int(x)), F1=int(x) + FSCALE // F0=2**x, F1=int(x) + FMOVDP F0, F1 // F0=2**x + FMOVDP F0, r+8(FP) + RET +not_finite: +// test bits for -Inf + MOVL x+4(FP), BX + MOVL x+0(FP), CX + CMPL BX, $0xfff00000 + JNE not_neginf + CMPL CX, $0 + JNE not_neginf + MOVL $0, r+8(FP) + MOVL $0, r+12(FP) + RET +not_neginf: + MOVL CX, r+8(FP) + MOVL BX, r+12(FP) + RET |