diff options
| author | Ondřej Surý <ondrej@sury.org> | 2011-01-17 12:40:45 +0100 |
|---|---|---|
| committer | Ondřej Surý <ondrej@sury.org> | 2011-01-17 12:40:45 +0100 |
| commit | 3e45412327a2654a77944249962b3652e6142299 (patch) | |
| tree | bc3bf69452afa055423cbe0c5cfa8ca357df6ccf /src/pkg/math/hypot_amd64.s | |
| parent | c533680039762cacbc37db8dc7eed074c3e497be (diff) | |
| download | golang-3e45412327a2654a77944249962b3652e6142299.tar.gz | |
Imported Upstream version 2011.01.12upstream/2011.01.12
Diffstat (limited to 'src/pkg/math/hypot_amd64.s')
| -rw-r--r-- | src/pkg/math/hypot_amd64.s | 50 |
1 files changed, 50 insertions, 0 deletions
diff --git a/src/pkg/math/hypot_amd64.s b/src/pkg/math/hypot_amd64.s new file mode 100644 index 000000000..1f691e70e --- /dev/null +++ b/src/pkg/math/hypot_amd64.s @@ -0,0 +1,50 @@ +// 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. + +#define PosInf 0x7ff0000000000000 +#define NaN 0x7FF0000000000001 + +// func Hypot(x, y float64) float64 +TEXT ·Hypot(SB),7,$0 + // test bits for special cases + MOVQ x+0(FP), BX + MOVQ $~(1<<63), AX + ANDQ AX, BX // x = |x| + MOVQ y+8(FP), CX + ANDQ AX, CX // y = |y| + MOVQ $PosInf, AX + CMPQ AX, BX + JLE isInfOrNaN + CMPQ AX, CX + JLE isInfOrNaN + // hypot = max * sqrt(1 + (min/max)**2) + MOVQ BX, X0 + MOVQ CX, X1 + ORQ CX, BX + JEQ isZero + MOVAPD X0, X2 + MAXSD X1, X0 + MINSD X2, X1 + DIVSD X0, X1 + MULSD X1, X1 + ADDSD $1.0, X1 + SQRTSD X1, X1 + MULSD X1, X0 + MOVSD X0, r+16(FP) + RET +isInfOrNaN: + CMPQ AX, BX + JEQ isInf + CMPQ AX, CX + JEQ isInf + MOVQ $NaN, AX + MOVQ AX, r+16(FP) // return NaN + RET +isInf: + MOVQ AX, r+16(FP) // return +Inf + RET +isZero: + MOVQ $0, AX + MOVQ AX, r+16(FP) // return 0 + RET |