fpsetround.c revision 1.1 
11.1Schs/*	$NetBSD: fpsetround.c,v 1.1 2004/07/24 19:09:29 chs Exp $	*/
21.1Schs
31.1Schs/*	$OpenBSD: fpsetround.c,v 1.3 2002/10/21 18:41:05 mickey Exp $	*/
41.1Schs
51.1Schs/*
61.1Schs * Written by Miodrag Vallat.  Public domain
71.1Schs */
81.1Schs
91.1Schs#include <sys/types.h>
101.1Schs#include <ieeefp.h>
111.1Schs
121.1Schsfp_rnd
131.1Schsfpsetround(fp_rnd rnd_dir)
141.1Schs{
151.1Schs	uint64_t fpsr;
161.1Schs	fp_rnd old;
171.1Schs
181.1Schs	__asm__ __volatile__("fstd %%fr0,0(%1)" : "=m" (fpsr) : "r" (&fpsr));
191.1Schs	old = (fpsr >> 41) & 0x03;
201.1Schs	fpsr = (fpsr & 0xfffff9ff00000000LL) |
211.1Schs	    ((uint64_t)(rnd_dir & 0x03) << 41);
221.1Schs	__asm__ __volatile__("fldd 0(%0),%%fr0" : : "r" (&fpsr));
231.1Schs	return (old);
241.1Schs}
25