1 /* 2 * Written by Michael Shalayeff. Public Domain 3 */ 4 5 #if defined(LIBM_SCCS) && !defined(lint) 6 static char rcsid[] = "$OpenBSD: s_ceilf.c,v 1.3 2003/01/16 19:17:33 mickey Exp $"; 7 #endif 8 9 #include <sys/types.h> 10 #include <machine/ieeefp.h> 11 #include "math.h" 12 13 float 14 ceilf(float x) 15 { 16 u_int64_t ofpsr, fpsr; 17 18 __asm__ __volatile__("fstds %%fr0,0(%1)" : "=m" (ofpsr) : "r" (&ofpsr)); 19 fpsr = (ofpsr & ~((u_int64_t)FP_RM << (9 + 32))) | 20 ((u_int64_t)FP_RP << (9 + 32)); 21 __asm__ __volatile__("fldds 0(%0), %%fr0" :: "r" (&fpsr)); 22 23 __asm__ __volatile__("frnd,sgl %0,%0" : "+f" (x)); 24 25 __asm__ __volatile__("fldds 0(%0), %%fr0" :: "r" (&ofpsr)); 26 return (x); 27 } 28