Home | History | Annotate | Line # | Download | only in fpe
      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