! File TLATAN2.3.TEXT
!
! Tests of atan2(y,x)
!
! NANs
35	=d	Q	H	OK	Q
35	=d	Q	Hd1	OK	Q
35	=d	Q	1	OK	Q
35	=d	Q	E	OK	Q
35	=d	Q	Ed1	OK	Q
35	=d	Q	0i1	OK	Q
35	=d	Q	0	OK	Q
35	=d	Q	-H	OK	Q
35	=d	Q	-Hd1	OK	Q
35	=d	Q	-1	OK	Q
35	=d	Q	-E	OK	Q
35	=d	Q	-Ed1	OK	Q
35	=d	Q	-0i1	OK	Q
35	=d	Q	-0	OK	Q
35	=d	H	Q	OK	Q
35	=d	Hd1	Q	OK	Q
35	=d	1	Q	OK	Q
35	=d	E	Q	OK	Q
35	=d	Ed1	Q	OK	Q
35	=d	0i1	Q	OK	Q
35	=d	0	Q	OK	Q
35	=d	-H	Q	OK	Q
35	=d	-Hd1	Q	OK	Q
35	=d	-1	Q	OK	Q
35	=d	-E	Q	OK	Q
35	=d	-Ed1	Q	OK	Q
35	=d	-0i1	Q	OK	Q
35	=d	-0	Q	OK	Q
35	=d	Q	Q	OK	Q
!
! atan2(y,x) returns y for y = 0 and x > 0.
35	=d	0	0i1	OK	0
35	=d	-0	0i1	OK	-0
35	=d	0	E	OK	0
35	=d	-0	E	OK	-0
35	=d	0	1d1	OK	0
35	=d	-0	1d1	OK	-0
35	=d	0	Hd1	OK	0
35	=d	-0	Hd1	OK	-0
35	=d	0	H	OK	0
35	=d	-0	H	OK	-0
!
! atan2(y,+0) returns y for y = 0.
35	=d	0	0	OK	0
35	=d	-0	0	OK	-0
!
! atan2(y,x) returns copysign(pi,y) for y = 0 and x < 0.
35	=d	0	-0i1	x	P	!	flags??
35	=d	-0	-0i1	x	-P
35	=d	0	-E	x	P
35	=d	-0	-E	x	-P
35	=d	0	-1d1	x	P
35	=d	-0	-1d1	x	-P
35	=d	0	-Hd1	x	P
35	=d	-0	-Hd1	x	-P
35	=d	0	-H	x	P
35	=d	-0	-H	x	-P
!
! atan2(y,-0) returns copysign(pi,y) for y = 0.
35	=d	0	-0	x	P	!	flags??
35	=d	-0	-0	x	-P
!
! atan2(y,x) returns pi/2 for y > 0 and x = 0.
35	=d	0i1	0	x	Pm1	!	flags??
35	=d	0i1	-0	x	Pm1
35	=d	E	0	x	Pm1
35	=d	E	-0	x	Pm1
35	=d	1d1	0	x	Pm1
35	=d	1d1	-0	x	Pm1
35	=d	Hd1	0	x	Pm1
35	=d	Hd1	-0	x	Pm1
35	=d	H	0	x	Pm1
35	=d	H	-0	x	Pm1
!
! atan2(y,x) returns -pi/2 for y < 0 and x = 0.
35	=d	-0i1	0	x	-Pm1	!	flags??
35	=d	-0i1	-0	x	-Pm1
35	=d	-E	0	x	-Pm1
35	=d	-E	-0	x	-Pm1
35	=d	-1d1	0	x	-Pm1
35	=d	-1d1	-0	x	-Pm1
35	=d	-Hd1	0	x	-Pm1
35	=d	-Hd1	-0	x	-Pm1
35	=d	-H	0	x	-Pm1
35	=d	-H	-0	x	-Pm1
!
! atan2(y,+INFINITY) returns copysign(0,y) for finite y != 0.
35	=d	0i1	H	OK	0
35	=d	-0i1	H	OK	-0
35	=d	E	H	OK	0
35	=d	-E	H	OK	-0
35	=d	1d1	H	OK	0
35	=d	-1d1	H	OK	-0
35	=d	Hd1	H	OK	0
35	=d	-Hd1	H	OK	-0
!
! atan2(y,-INFINITY) returns copysign(pi,y) for finite y.
35	=d	0i1	-H	x	P	!	flags??
35	=d	-0i1	-H	x	-P
35	=d	E	-H	x	P
35	=d	-E	-H	x	-P
35	=d	1d1	-H	x	P
35	=d	-1d1	-H	x	-P
35	=d	Hd1	-H	x	P
35	=d	-Hd1	-H	x	-P
!
! atan2(+INFINITY,x) returns pi/2 for finite x.
35	=d	H	0i1	x	Pm1	!	flags??
35	=d	H	-0i1	x	Pm1
35	=d	H	E	x	Pm1
35	=d	H	-E	x	Pm1
35	=d	H	1d1	x	Pm1
35	=d	H	-1d1	x	Pm1
35	=d	H	Hd1	x	Pm1
35	=d	H	-Hd1	x	Pm1
!
! atan2(-INFINITY,x) returns -pi/2 for finite x.
35	=d	-H	0i1	x	-Pm1	!	flags??
35	=d	-H	-0i1	x	-Pm1
35	=d	-H	E	x	-Pm1
35	=d	-H	-E	x	-Pm1
35	=d	-H	1d1	x	-Pm1
35	=d	-H	-1d1	x	-Pm1
35	=d	-H	Hd1	x	-Pm1
35	=d	-H	-Hd1	x	-Pm1
!
! atan2(y,x) when both y and x are infinite.
35	=d	H	H	x	Pm2	!	pi/4	flags??
35	=d	-H	H	x	-Pm2	!	-pi/4	flags??
35	=d	H	-H	x	$400096CBE3F9990E91A8	!	3*pi/4	flags??
35	=d	-H	-H	x	$C00096CBE3F9990E91A8	!	-3*pi/4	flags??
!
! atan2(y,x) is y/x - noise when x,y > 0 and y/x is tiny (< 2^(-32).

! <scp> Disabled <0e and >0e variants, since function is computed round-nearest always.
35	=d	E	1	x	E
35	=d	Ei1	1	x	Ei1
35	=d	Ep1	1	x	Ep1
35	=d	Ep1i3	1	x	Ep1i3
35	=d	Ed2	1	ux	Ed2
35	=d	Ep1d5	1	x	Ep1d5
35	=d	0i1	1	ux	0i1
35	=d	0i9	1	ux	0i9
!
! atan2(y,x) is y/x + noise when x > 0, y < 0, and |y|/x is tiny.
35	=d	-E	1	x	-E
35	=d	-Ei1	1	x	-Ei1
35	=d	-Ep1	1	x	-Ep1
35	=d	-Ep1i3	1	x	-Ep1i3
35	=d	-Ed2	1	ux	-Ed2
35	=d	-Ep1d5	1	x	-Ep1d5
35	=d	-0i1	1	ux	-0i1
35	=d	-0i9	1	ux	-0i9

!
! atan2(y,x) is pi - (y/x - noise) when x < 0, y > 0 and y/|x| is tiny (< 2^(-32).
35	=d	E	-1	x	P
35	=d	Ei1	-1	x	P
35	=d	Ep1	-1	x	P
35	=d	Ep1i3	-1	x	P
35	=d	Ed2	-1	x	P
35	=d	Ep1d5	-1	x	P
35	=d	0i1	-1	x	P
35	=d	0i9	-1	x	P
!
! atan2(y,x) is -pi + (y/x - noise) when x < 0, y < 0 and y/x is tiny (< 2^(-32).
35	=d	-E	-1	x	-P
35	=d	-Ei1	-1	x	-P
35	=d	-Ep1	-1	x	-P
35	=d	-Ep1i3	-1	x	-P
35	=d	-Ed2	-1	x	-P
35	=d	-Ep1d5	-1	x	-P
35	=d	-0i1	-1	x	-P
35	=d	-0i9	-1	x	-P
!
! atan2(y,x) is a multiple of pi/3 when |y/x| = sqrt(3).

! <scp> harness appears to incorrectly convert the '$' constants below. 
! <scp> Use long double derived values instead.
!35	=e	  	$3FFFDDB3D742C265539E	1		x    $3FFF860A91C16B9B2C24	    ;pi/3
!35	=e	  	$BFFFDDB3D742C265539E	1		x    $BFFF860A91C16B9B2C24	    ;-pi/3
!35	=e	  	$3FFFDDB3D742C265539E	-1		x    $4000860A91C16B9B2C24-1    ;2*pi/3
!35	=e	  	$BFFFDDB3D742C265539E	-1		x    $C000860A91C16B9B2C24-1    ;-2*pi/3

35	=d	$3FFFDDB3D742C265539E	1	x	$3FFF860A91C16B9B2C23-1	;	pi/3
35	=d	$BFFFDDB3D742C265539E	1	x	$BFFF860A91C16B9B2C23-1	;	-pi/3
45	=d	0x3ffbb67a	0xe8584caa	0xbff00000	0x0	x	0x4000c152	0x382d7366	;	2*pi/3
45	=d	0xbffbb67a	0xe8584caa	0xbff00000	0x0	x	0xc000c152	0x382d7366	;	-2*pi/3