ldexp_881.c revision 1.1.6.2 1 1.1.6.2 yamt /* $NetBSD: ldexp_881.c,v 1.1.6.2 2014/05/22 11:36:47 yamt Exp $ */
2 1.1.6.2 yamt
3 1.1.6.2 yamt /*-
4 1.1.6.2 yamt * Copyright (c) 1998 The NetBSD Foundation, Inc.
5 1.1.6.2 yamt * All rights reserved.
6 1.1.6.2 yamt *
7 1.1.6.2 yamt * This code is derived from software contributed to The NetBSD Foundation
8 1.1.6.2 yamt * by Charles M. Hannum.
9 1.1.6.2 yamt *
10 1.1.6.2 yamt * Redistribution and use in source and binary forms, with or without
11 1.1.6.2 yamt * modification, are permitted provided that the following conditions
12 1.1.6.2 yamt * are met:
13 1.1.6.2 yamt * 1. Redistributions of source code must retain the above copyright
14 1.1.6.2 yamt * notice, this list of conditions and the following disclaimer.
15 1.1.6.2 yamt * 2. Redistributions in binary form must reproduce the above copyright
16 1.1.6.2 yamt * notice, this list of conditions and the following disclaimer in the
17 1.1.6.2 yamt * documentation and/or other materials provided with the distribution.
18 1.1.6.2 yamt *
19 1.1.6.2 yamt * THIS SOFTWARE IS PROVIDED BY THE NETBSD FOUNDATION, INC. AND CONTRIBUTORS
20 1.1.6.2 yamt * ``AS IS'' AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED
21 1.1.6.2 yamt * TO, THE IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR
22 1.1.6.2 yamt * PURPOSE ARE DISCLAIMED. IN NO EVENT SHALL THE FOUNDATION OR CONTRIBUTORS
23 1.1.6.2 yamt * BE LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR
24 1.1.6.2 yamt * CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF
25 1.1.6.2 yamt * SUBSTITUTE GOODS OR SERVICES; LOSS OF USE, DATA, OR PROFITS; OR BUSINESS
26 1.1.6.2 yamt * INTERRUPTION) HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN
27 1.1.6.2 yamt * CONTRACT, STRICT LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE)
28 1.1.6.2 yamt * ARISING IN ANY WAY OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE
29 1.1.6.2 yamt * POSSIBILITY OF SUCH DAMAGE.
30 1.1.6.2 yamt */
31 1.1.6.2 yamt
32 1.1.6.2 yamt #include <sys/cdefs.h>
33 1.1.6.2 yamt #if defined(LIBC_SCCS) && !defined(lint)
34 1.1.6.2 yamt __RCSID("$NetBSD: ldexp_881.c,v 1.1.6.2 2014/05/22 11:36:47 yamt Exp $");
35 1.1.6.2 yamt #endif /* LIBC_SCCS and not lint */
36 1.1.6.2 yamt
37 1.1.6.2 yamt #include <sys/types.h>
38 1.1.6.2 yamt #include <machine/ieee.h>
39 1.1.6.2 yamt #include <math.h>
40 1.1.6.2 yamt
41 1.1.6.2 yamt /*
42 1.1.6.2 yamt * ldexp(value, exp): return value * (2 ** exp).
43 1.1.6.2 yamt */
44 1.1.6.2 yamt double
45 1.1.6.2 yamt ldexp(value, exp2)
46 1.1.6.2 yamt double value;
47 1.1.6.2 yamt int exp2;
48 1.1.6.2 yamt {
49 1.1.6.2 yamt double temp;
50 1.1.6.2 yamt
51 1.1.6.2 yamt __asm ("fscalel %2,%1"
52 1.1.6.2 yamt : "=f" (temp)
53 1.1.6.2 yamt : "0" (value), "g" (exp2));
54 1.1.6.2 yamt return (temp);
55 1.1.6.2 yamt }
56