n_sincos.S revision 1.3 
1/*	$NetBSD: n_sincos.S,v 1.3 1999/07/02 15:37:35 simonb Exp $	*/
2/*
3 * Copyright (c) 1985, 1993
4 *	The Regents of the University of California.  All rights reserved.
5 *
6 * Redistribution and use in source and binary forms, with or without
7 * modification, are permitted provided that the following conditions
8 * are met:
9 * 1. Redistributions of source code must retain the above copyright
10 *    notice, this list of conditions and the following disclaimer.
11 * 2. Redistributions in binary form must reproduce the above copyright
12 *    notice, this list of conditions and the following disclaimer in the
13 *    documentation and/or other materials provided with the distribution.
14 * 3. All advertising materials mentioning features or use of this software
15 *    must display the following acknowledgement:
16 *	This product includes software developed by the University of
17 *	California, Berkeley and its contributors.
18 * 4. Neither the name of the University nor the names of its contributors
19 *    may be used to endorse or promote products derived from this software
20 *    without specific prior written permission.
21 *
22 * THIS SOFTWARE IS PROVIDED BY THE REGENTS AND CONTRIBUTORS ``AS IS'' AND
23 * ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE
24 * IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE
25 * ARE DISCLAIMED.  IN NO EVENT SHALL THE REGENTS OR CONTRIBUTORS BE LIABLE
26 * FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR CONSEQUENTIAL
27 * DAMAGES (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS
28 * OR SERVICES; LOSS OF USE, DATA, OR PROFITS; OR BUSINESS INTERRUPTION)
29 * HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT
30 * LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY
31 * OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF
32 * SUCH DAMAGE.
33 *
34 *	@(#)sincos.s	8.1 (Berkeley) 6/4/93
35 */
36
37/*
38 *  This is the implementation of Peter Tang's double precision
39 *  sine and cosine for the VAX using Bob Corbett's argument reduction.
40 *
41 *  Notes:
42 *       under 1,024,000 random arguments testing on [0,2*pi]
43 *       sin() observed maximum error = 0.814 ulps
44 *       cos() observed maximum error = 0.792 ulps
45 *
46 * double sin(arg)
47 * double arg;
48 * method: true range reduction to [-pi/4,pi/4], P. Tang  &  B. Corbett
49 * S. McDonald, April 4,  1985
50 */
51#include <machine/asm.h>
52	.text
53	.align	1
54	.globl	_sin
55	.type	_sin@function
56_sin:
57	.word	0xffc		# save r2-r11
58	movq	4(ap),r0
59	bicw3	$0x807f,r0,r2
60	beql	1f		# if x is zero or reserved operand then return x
61/*
62 * Save the PSL's IV & FU bits on the stack.
63 */
64	movpsl	r2
65	bicw3	$0xff9f,r2,-(sp)
66/*
67 * Clear the IV & FU bits.
68 */
69	bicpsw	$0x0060
70/*
71 *  Entered by  sine    ; save  0  in  r4 .
72 */
73	jsb	libm$argred
74	movl	$0,r4
75	jsb	libm$sincos
76	bispsw	(sp)+
771:	ret
78
79/*
80 * double cos(arg)
81 * double arg;
82 * method: true range reduction to [-pi/4,pi/4], P. Tang  &  B. Corbett
83 * S. McDonald, April 4,  1985
84 */
85	.text
86	.align	1
87	.globl	_cos
88	.type	_cos,@function
89_cos:
90	.word	0xffc		# save r2-r11
91	movq	4(ap),r0
92	bicw3	$0x7f,r0,r2
93	cmpw	$0x8000,r2
94	beql	1f		# if x is reserved operand then return x
95/*
96 * Save the PSL's IV & FU bits on the stack.
97 */
98	movpsl	r2
99	bicw3	$0xff9f,r2,-(sp)
100/*
101 * Clear the IV & FU bits.
102 */
103	bicpsw	$0x0060
104/*
105 *  Entered by  cosine  ; save  1  in  r4 .
106 */
107	jsb	libm$argred
108	movl	$1,r4
109	jsb	libm$sincos
110	bispsw	(sp)+
1111:	ret
112