Blame - src/math/hypot_amd64.s - go

blob: a68eebc8ca619810107fac75a91c0afa89e89ebb [file] [log] [blame]

Charles L. Dorian	c8c2bdb	2010-08-06 16:50:48 -0700	[diff] [blame]	1	// Copyright 2010 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
Rob Pike	8bca148	2014-08-12 17:04:45 -0700	[diff] [blame]	5	#include "textflag.h"
Keith Randall	1f79663	2013-08-12 10:25:18 -0700	[diff] [blame]	6
Charles L. Dorian	a5c4e0f	2012-08-09 09:40:05 +1000	[diff] [blame]	7	#define PosInf 0x7FF0000000000000
				8	#define NaN 0x7FF8000000000001
Charles L. Dorian	c8c2bdb	2010-08-06 16:50:48 -0700	[diff] [blame]	9
Russ Cox	8572950	2012-09-21 00:35:56 -0400	[diff] [blame]	10	// func Hypot(p, q float64) float64
Keith Randall	1f79663	2013-08-12 10:25:18 -0700	[diff] [blame]	11	TEXT ·Hypot(SB),NOSPLIT,$0
Charles L. Dorian	c8c2bdb	2010-08-06 16:50:48 -0700	[diff] [blame]	12	// test bits for special cases
Russ Cox	8572950	2012-09-21 00:35:56 -0400	[diff] [blame]	13	MOVQ p+0(FP), BX
Charles L. Dorian	c8c2bdb	2010-08-06 16:50:48 -0700	[diff] [blame]	14	MOVQ $~(1<<63), AX
Russ Cox	8572950	2012-09-21 00:35:56 -0400	[diff] [blame]	15	ANDQ AX, BX // p = \|p\|
				16	MOVQ q+8(FP), CX
				17	ANDQ AX, CX // q = \|q\|
Charles L. Dorian	c8c2bdb	2010-08-06 16:50:48 -0700	[diff] [blame]	18	MOVQ $PosInf, AX
				19	CMPQ AX, BX
				20	JLE isInfOrNaN
				21	CMPQ AX, CX
				22	JLE isInfOrNaN
				23	// hypot = max * sqrt(1 + (min/max)**2)
				24	MOVQ BX, X0
				25	MOVQ CX, X1
				26	ORQ CX, BX
				27	JEQ isZero
				28	MOVAPD X0, X2
				29	MAXSD X1, X0
				30	MINSD X2, X1
				31	DIVSD X0, X1
				32	MULSD X1, X1
				33	ADDSD $1.0, X1
				34	SQRTSD X1, X1
				35	MULSD X1, X0
Russ Cox	07720b6	2013-03-22 12:57:55 -0400	[diff] [blame]	36	MOVSD X0, ret+16(FP)
Charles L. Dorian	c8c2bdb	2010-08-06 16:50:48 -0700	[diff] [blame]	37	RET
				38	isInfOrNaN:
				39	CMPQ AX, BX
				40	JEQ isInf
				41	CMPQ AX, CX
				42	JEQ isInf
				43	MOVQ $NaN, AX
Russ Cox	07720b6	2013-03-22 12:57:55 -0400	[diff] [blame]	44	MOVQ AX, ret+16(FP) // return NaN
Charles L. Dorian	c8c2bdb	2010-08-06 16:50:48 -0700	[diff] [blame]	45	RET
				46	isInf:
Russ Cox	07720b6	2013-03-22 12:57:55 -0400	[diff] [blame]	47	MOVQ AX, ret+16(FP) // return +Inf
Charles L. Dorian	c8c2bdb	2010-08-06 16:50:48 -0700	[diff] [blame]	48	RET
				49	isZero:
				50	MOVQ $0, AX
Russ Cox	07720b6	2013-03-22 12:57:55 -0400	[diff] [blame]	51	MOVQ AX, ret+16(FP) // return 0
Charles L. Dorian	c8c2bdb	2010-08-06 16:50:48 -0700	[diff] [blame]	52	RET