Home | History | Annotate | Line # | Download | only in gen
fpsetround.c revision 1.11.2.1
      1  1.11.2.1     yamt /* $NetBSD: fpsetround.c,v 1.11.2.1 2012/04/17 00:05:12 yamt Exp $ */
      2       1.1      cgd 
      3       1.1      cgd /*
      4       1.1      cgd  * Copyright (c) 1995 Christopher G. Demetriou
      5       1.1      cgd  * All rights reserved.
      6       1.3      cgd  *
      7       1.1      cgd  * Redistribution and use in source and binary forms, with or without
      8       1.1      cgd  * modification, are permitted provided that the following conditions
      9       1.1      cgd  * are met:
     10       1.1      cgd  * 1. Redistributions of source code must retain the above copyright
     11       1.1      cgd  *    notice, this list of conditions and the following disclaimer.
     12       1.1      cgd  * 2. Redistributions in binary form must reproduce the above copyright
     13       1.1      cgd  *    notice, this list of conditions and the following disclaimer in the
     14       1.1      cgd  *    documentation and/or other materials provided with the distribution.
     15       1.1      cgd  * 3. All advertising materials mentioning features or use of this software
     16       1.1      cgd  *    must display the following acknowledgement:
     17       1.3      cgd  *          This product includes software developed for the
     18       1.7     salo  *          NetBSD Project.  See http://www.NetBSD.org/ for
     19       1.3      cgd  *          information about NetBSD.
     20       1.1      cgd  * 4. The name of the author may not be used to endorse or promote products
     21       1.3      cgd  *    derived from this software without specific prior written permission.
     22       1.3      cgd  *
     23       1.1      cgd  * THIS SOFTWARE IS PROVIDED BY THE AUTHOR ``AS IS'' AND ANY EXPRESS OR
     24       1.1      cgd  * IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE IMPLIED WARRANTIES
     25       1.1      cgd  * OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE ARE DISCLAIMED.
     26       1.1      cgd  * IN NO EVENT SHALL THE AUTHOR BE LIABLE FOR ANY DIRECT, INDIRECT,
     27       1.1      cgd  * INCIDENTAL, SPECIAL, EXEMPLARY, OR CONSEQUENTIAL DAMAGES (INCLUDING, BUT
     28       1.1      cgd  * NOT LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES; LOSS OF USE,
     29       1.1      cgd  * DATA, OR PROFITS; OR BUSINESS INTERRUPTION) HOWEVER CAUSED AND ON ANY
     30       1.1      cgd  * THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT LIABILITY, OR TORT
     31       1.1      cgd  * (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY OUT OF THE USE OF
     32       1.1      cgd  * THIS SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF SUCH DAMAGE.
     33       1.3      cgd  *
     34       1.3      cgd  * <<Id: LICENSE,v 1.2 2000/06/14 15:57:33 cgd Exp>>
     35       1.1      cgd  */
     36       1.1      cgd 
     37       1.6  thorpej #include <sys/cdefs.h>
     38       1.9    lukem #if defined(LIBC_SCCS) && !defined(lint)
     39  1.11.2.1     yamt __RCSID("$NetBSD: fpsetround.c,v 1.11.2.1 2012/04/17 00:05:12 yamt Exp $");
     40       1.9    lukem #endif /* LIBC_SCCS and not lint */
     41       1.6  thorpej 
     42       1.6  thorpej #include "namespace.h"
     43       1.6  thorpej 
     44       1.1      cgd #include <sys/types.h>
     45       1.1      cgd #include <ieeefp.h>
     46       1.6  thorpej 
     47       1.6  thorpej #ifdef __weak_alias
     48       1.6  thorpej __weak_alias(fpsetround,_fpsetround)
     49       1.6  thorpej #endif
     50       1.1      cgd 
     51       1.1      cgd fp_rnd
     52      1.11     matt fpsetround(fp_rnd rnd_dir)
     53       1.1      cgd {
     54      1.11     matt 	union {
     55      1.11     matt 		double d;
     56      1.11     matt 		uint64_t u64;
     57      1.11     matt 	} fpcrval;
     58      1.11     matt 	fp_rnd old;
     59       1.1      cgd 
     60      1.11     matt 	__asm("excb; mf_fpcr %0; excb" : "=f" (fpcrval.d));
     61  1.11.2.1     yamt 	old = (fp_rnd)(fpcrval.u64 >> 58) & 0x3;
     62       1.1      cgd 
     63      1.11     matt 	rnd_dir ^= old;
     64       1.1      cgd 
     65      1.11     matt 	fpcrval.u64 ^= (long)rnd_dir << 58;
     66       1.1      cgd 
     67      1.11     matt 	__asm("excb; mt_fpcr %0; excb" : : "f" (fpcrval.d));
     68      1.11     matt 
     69      1.11     matt 	return old;
     70       1.1      cgd }
     71