fpu_getexp.c revision 1.5 1 1.5 isaki /* $NetBSD: fpu_getexp.c,v 1.5 2011/07/18 07:44:30 isaki Exp $ */
2 1.1 briggs
3 1.1 briggs /*
4 1.1 briggs * Copyright (c) 1995 Ken Nakata
5 1.1 briggs * All rights reserved.
6 1.1 briggs *
7 1.1 briggs * Redistribution and use in source and binary forms, with or without
8 1.1 briggs * modification, are permitted provided that the following conditions
9 1.1 briggs * are met:
10 1.1 briggs * 1. Redistributions of source code must retain the above copyright
11 1.1 briggs * notice, this list of conditions and the following disclaimer.
12 1.1 briggs * 2. Redistributions in binary form must reproduce the above copyright
13 1.1 briggs * notice, this list of conditions and the following disclaimer in the
14 1.1 briggs * documentation and/or other materials provided with the distribution.
15 1.1 briggs * 3. Neither the name of the author nor the names of its contributors
16 1.1 briggs * may be used to endorse or promote products derived from this software
17 1.1 briggs * without specific prior written permission.
18 1.1 briggs *
19 1.1 briggs * THIS SOFTWARE IS PROVIDED BY THE AUTHOR AND CONTRIBUTORS ``AS IS'' AND
20 1.1 briggs * ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE
21 1.1 briggs * IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE
22 1.1 briggs * ARE DISCLAIMED. IN NO EVENT SHALL THE AUTHOR OR CONTRIBUTORS BE LIABLE
23 1.1 briggs * FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR CONSEQUENTIAL
24 1.1 briggs * DAMAGES (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS
25 1.1 briggs * OR SERVICES; LOSS OF USE, DATA, OR PROFITS; OR BUSINESS INTERRUPTION)
26 1.1 briggs * HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT
27 1.1 briggs * LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY
28 1.1 briggs * OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF
29 1.1 briggs * SUCH DAMAGE.
30 1.1 briggs *
31 1.1 briggs * @(#)fpu_getexp.c 10/8/95
32 1.1 briggs */
33 1.2 lukem
34 1.2 lukem #include <sys/cdefs.h>
35 1.5 isaki __KERNEL_RCSID(0, "$NetBSD: fpu_getexp.c,v 1.5 2011/07/18 07:44:30 isaki Exp $");
36 1.1 briggs
37 1.1 briggs #include <sys/types.h>
38 1.1 briggs
39 1.1 briggs #include "fpu_emulate.h"
40 1.1 briggs
41 1.1 briggs struct fpn *
42 1.4 dsl fpu_getexp(struct fpemu *fe)
43 1.1 briggs {
44 1.5 isaki struct fpn *fp = &fe->fe_f2;
45 1.1 briggs
46 1.5 isaki fe->fe_fpsr &= ~FPSR_EXCP; /* clear all exceptions */
47 1.1 briggs
48 1.5 isaki if (fp->fp_class == FPC_INF) {
49 1.5 isaki fp = fpu_newnan(fe);
50 1.5 isaki fe->fe_fpsr |= FPSR_OPERR;
51 1.5 isaki } else if (fp->fp_class == FPC_NUM) { /* a number */
52 1.5 isaki fpu_explode(fe, &fe->fe_f3, FTYPE_LNG, &fp->fp_exp);
53 1.5 isaki fp = &fe->fe_f3;
54 1.5 isaki } else if (fp->fp_class == FPC_SNAN) { /* signaling NaN */
55 1.5 isaki fe->fe_fpsr |= FPSR_SNAN;
56 1.5 isaki } /* else if fp == zero or fp == quiet NaN, return itself */
57 1.5 isaki return fp;
58 1.1 briggs }
59 1.1 briggs
60 1.1 briggs struct fpn *
61 1.4 dsl fpu_getman(struct fpemu *fe)
62 1.1 briggs {
63 1.5 isaki struct fpn *fp = &fe->fe_f2;
64 1.1 briggs
65 1.5 isaki fe->fe_fpsr &= ~FPSR_EXCP; /* clear all exceptions */
66 1.1 briggs
67 1.5 isaki if (fp->fp_class == FPC_INF) {
68 1.5 isaki fp = fpu_newnan(fe);
69 1.5 isaki fe->fe_fpsr |= FPSR_OPERR;
70 1.5 isaki } else if (fp->fp_class == FPC_NUM) { /* a number */
71 1.5 isaki fp->fp_exp = 0;
72 1.5 isaki } else if (fp->fp_class == FPC_SNAN) { /* signaling NaN */
73 1.5 isaki fe->fe_fpsr |= FPSR_SNAN;
74 1.5 isaki } /* else if fp == zero or fp == quiet NaN, return itself */
75 1.5 isaki return fp;
76 1.1 briggs }
77