1/*-
2 * Copyright (c) 2011 The NetBSD Foundation, Inc.
3 * All rights reserved.
4 *
5 * This code is derived from software contributed to The NetBSD Foundation
6 * by Matt Thomas of 3am Software Foundry.
7 *
8 * Redistribution and use in source and binary forms, with or without
9 * modification, are permitted provided that the following conditions
10 * are met:
11 * 1. Redistributions of source code must retain the above copyright
12 *    notice, this list of conditions and the following disclaimer.
13 * 2. Redistributions in binary form must reproduce the above copyright
14 *    notice, this list of conditions and the following disclaimer in the
15 *    documentation and/or other materials provided with the distribution.
16 *
17 * THIS SOFTWARE IS PROVIDED BY THE NETBSD FOUNDATION, INC. AND CONTRIBUTORS
18 * ``AS IS'' AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED
19 * TO, THE IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR
20 * PURPOSE ARE DISCLAIMED.  IN NO EVENT SHALL THE FOUNDATION OR CONTRIBUTORS
21 * BE LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR
22 * CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF
23 * SUBSTITUTE GOODS OR SERVICES; LOSS OF USE, DATA, OR PROFITS; OR BUSINESS
24 * INTERRUPTION) HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN
25 * CONTRACT, STRICT LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE)
26 * ARISING IN ANY WAY OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE
27 * POSSIBILITY OF SUCH DAMAGE.
28 */
29
30#include <machine/asm.h>
31#include <sys/errno.h>
32
33RCSID("$NetBSD$")
34
35ENTRY(scalbnf, 0)
36	cvtfd	4(%ap), %r0
37	movl	8(%ap), %r2
38	jr	1f
39
40ALTENTRY(scalbnl)
41ENTRY(scalbn, 0)
42	movd	4(%ap), %r0
43	movl	12(%ap), %r2
44
45	/*
46	 * First see if the result will be within bounds
47	 */
481:	cmpl	%r2,$(256+55)
49	jgeq	overflow
50	cmpl	%r2,$-(256+55)
51	jleq	underflow
52
53	extzv	$7, $9, %r0, %r3	/* extract exponent */
54	jneq	5f			/* 0 * N is still 0 */
55	cmpw	%r3, $-256		/* is this a NAN? */
56	jgtr	2f			/*    nope */
57	ret
58
592:	movzbl	%r3, %r3		/* clear sign bit */
60	addl2	%r2, %r3		/* add to exponent */
61	jleq	underflow		/*   <= 0? */
62
633:	jbs	$8, %r3, overflow	/* did exponent overflow? */
64	insv	$7, $8, %r3, %r0	/* replace exponent */
65	ret
66
67underflow:
68	clrq	%r0			/* return 0.0 */
69	jr	copysign		/* or mabye -0.0 */
70overflow:
71	pushl	$ERANGE			/* overflow */
72	calls	$1, _C_LABEL(infnan)	/* possibly signal */
73copysign:
74	jbc	$15, 4(%ap), 5f		/* if clear, don't set sign-bit */
75	bisw2	$0x8000, %r0		/* preserve sign-bit */
765:
77	ret
78