Home | History | Annotate | Line # | Download | only in spmath
setovfl.c revision 1.6
      1  1.6       dsl /*	$NetBSD: setovfl.c,v 1.6 2009/03/14 15:36:07 dsl Exp $	*/
      2  1.1  fredette 
      3  1.1  fredette /*	$OpenBSD: setovfl.c,v 1.4 2001/03/29 03:58:19 mickey Exp $	*/
      4  1.1  fredette 
      5  1.1  fredette /*
      6  1.1  fredette  * Copyright 1996 1995 by Open Software Foundation, Inc.
      7  1.1  fredette  *              All Rights Reserved
      8  1.1  fredette  *
      9  1.1  fredette  * Permission to use, copy, modify, and distribute this software and
     10  1.1  fredette  * its documentation for any purpose and without fee is hereby granted,
     11  1.1  fredette  * provided that the above copyright notice appears in all copies and
     12  1.1  fredette  * that both the copyright notice and this permission notice appear in
     13  1.1  fredette  * supporting documentation.
     14  1.1  fredette  *
     15  1.1  fredette  * OSF DISCLAIMS ALL WARRANTIES WITH REGARD TO THIS SOFTWARE
     16  1.1  fredette  * INCLUDING ALL IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS
     17  1.1  fredette  * FOR A PARTICULAR PURPOSE.
     18  1.1  fredette  *
     19  1.1  fredette  * IN NO EVENT SHALL OSF BE LIABLE FOR ANY SPECIAL, INDIRECT, OR
     20  1.1  fredette  * CONSEQUENTIAL DAMAGES OR ANY DAMAGES WHATSOEVER RESULTING FROM
     21  1.1  fredette  * LOSS OF USE, DATA OR PROFITS, WHETHER IN ACTION OF CONTRACT,
     22  1.1  fredette  * NEGLIGENCE, OR OTHER TORTIOUS ACTION, ARISING OUT OF OR IN CONNECTION
     23  1.1  fredette  * WITH THE USE OR PERFORMANCE OF THIS SOFTWARE.
     24  1.1  fredette  *
     25  1.1  fredette  */
     26  1.1  fredette /*
     27  1.1  fredette  * pmk1.1
     28  1.1  fredette  */
     29  1.1  fredette /*
     30  1.1  fredette  * (c) Copyright 1986 HEWLETT-PACKARD COMPANY
     31  1.1  fredette  *
     32  1.1  fredette  * To anyone who acknowledges that this file is provided "AS IS"
     33  1.1  fredette  * without any express or implied warranty:
     34  1.1  fredette  *     permission to use, copy, modify, and distribute this file
     35  1.1  fredette  * for any purpose is hereby granted without fee, provided that
     36  1.1  fredette  * the above copyright notice and this notice appears in all
     37  1.1  fredette  * copies, and that the name of Hewlett-Packard Company not be
     38  1.1  fredette  * used in advertising or publicity pertaining to distribution
     39  1.1  fredette  * of the software without specific, written prior permission.
     40  1.1  fredette  * Hewlett-Packard Company makes no representations about the
     41  1.1  fredette  * suitability of this software for any purpose.
     42  1.1  fredette  */
     43  1.1  fredette 
     44  1.2     lukem #include <sys/cdefs.h>
     45  1.6       dsl __KERNEL_RCSID(0, "$NetBSD: setovfl.c,v 1.6 2009/03/14 15:36:07 dsl Exp $");
     46  1.1  fredette 
     47  1.1  fredette #include "../spmath/float.h"
     48  1.1  fredette #include "../spmath/sgl_float.h"
     49  1.1  fredette #include "../spmath/dbl_float.h"
     50  1.1  fredette 
     51  1.3      matt sgl_floating_point
     52  1.6       dsl sgl_setoverflow(unsigned int sign)
     53  1.1  fredette {
     54  1.1  fredette 	sgl_floating_point result;
     55  1.1  fredette 
     56  1.1  fredette 	/* set result to infinity or largest number */
     57  1.3      matt #if 0
     58  1.1  fredette 	switch (Rounding_mode()) {
     59  1.1  fredette 		case ROUNDPLUS:
     60  1.1  fredette 			if (sign) {
     61  1.1  fredette 				Sgl_setlargestnegative(result);
     62  1.3      matt 			} else {
     63  1.1  fredette 				Sgl_setinfinitypositive(result);
     64  1.1  fredette 			}
     65  1.1  fredette 			break;
     66  1.1  fredette 		case ROUNDMINUS:
     67  1.1  fredette 			if (sign==0) {
     68  1.1  fredette 				Sgl_setlargestpositive(result);
     69  1.3      matt 			} else {
     70  1.1  fredette 				Sgl_setinfinitynegative(result);
     71  1.1  fredette 			}
     72  1.1  fredette 			break;
     73  1.1  fredette 		case ROUNDNEAREST:
     74  1.1  fredette 			Sgl_setinfinity(result,sign);
     75  1.1  fredette 			break;
     76  1.3      matt 		default:
     77  1.1  fredette 		case ROUNDZERO:
     78  1.1  fredette 			Sgl_setlargest(result,sign);
     79  1.3      matt 			break;
     80  1.1  fredette 	}
     81  1.3      matt #endif
     82  1.3      matt 	Sgl_setinfinity(result,sign);
     83  1.1  fredette 	return(result);
     84  1.1  fredette }
     85  1.1  fredette 
     86  1.1  fredette dbl_floating_point dbl_setoverflow(sign)
     87  1.1  fredette 
     88  1.1  fredette unsigned int sign;
     89  1.1  fredette {
     90  1.1  fredette 	dbl_floating_point result;
     91  1.1  fredette 
     92  1.1  fredette 	/* set result to infinity or largest number */
     93  1.1  fredette 	/* ignore for now
     94  1.1  fredette 	switch (Rounding_mode()) {
     95  1.1  fredette 		case ROUNDPLUS:
     96  1.1  fredette 			if (sign) {
     97  1.1  fredette 				Dbl_setlargestnegative(result);
     98  1.1  fredette 			}
     99  1.1  fredette 			else {
    100  1.1  fredette 				Dbl_setinfinitypositive(result);
    101  1.1  fredette 			}
    102  1.1  fredette 			break;
    103  1.1  fredette 		case ROUNDMINUS:
    104  1.1  fredette 			if (sign==0) {
    105  1.1  fredette 				Dbl_setlargestpositive(result);
    106  1.1  fredette 			}
    107  1.1  fredette 			else {
    108  1.1  fredette 				Dbl_setinfinitynegative(result);
    109  1.1  fredette 			}
    110  1.1  fredette 			break;
    111  1.1  fredette 		case ROUNDNEAREST:
    112  1.1  fredette 			Dbl_setinfinity(result,sign);
    113  1.1  fredette 			break;
    114  1.1  fredette 		case ROUNDZERO:
    115  1.1  fredette 			Dbl_setlargest(result,sign);
    116  1.1  fredette 	}
    117  1.1  fredette 	*/
    118  1.5     skrll 	Dbl_setinfinity(result.wd0,result.wd1,sign);
    119  1.1  fredette 	return(result);
    120  1.1  fredette }
    121