1 /* $OpenBSD: s_nearbyint.c,v 1.1 2011/04/28 18:05:39 martynas Exp $ */ 2 3 /* 4 * Copyright (c) 2011 Martynas Venckus <martynas@openbsd.org> 5 * 6 * Permission to use, copy, modify, and distribute this software for any 7 * purpose with or without fee is hereby granted, provided that the above 8 * copyright notice and this permission notice appear in all copies. 9 * 10 * THE SOFTWARE IS PROVIDED "AS IS" AND THE AUTHOR DISCLAIMS ALL WARRANTIES 11 * WITH REGARD TO THIS SOFTWARE INCLUDING ALL IMPLIED WARRANTIES OF 12 * MERCHANTABILITY AND FITNESS. IN NO EVENT SHALL THE AUTHOR BE LIABLE FOR 13 * ANY SPECIAL, DIRECT, INDIRECT, OR CONSEQUENTIAL DAMAGES OR ANY DAMAGES 14 * WHATSOEVER RESULTING FROM LOSS OF USE, DATA OR PROFITS, WHETHER IN AN 15 * ACTION OF CONTRACT, NEGLIGENCE OR OTHER TORTIOUS ACTION, ARISING OUT OF 16 * OR IN CONNECTION WITH THE USE OR PERFORMANCE OF THIS SOFTWARE. 17 */ 18 19 #include <fenv.h> 20 #include <math.h> 21 22 double 23 nearbyint(double x) 24 { 25 fenv_t fenv; 26 double res; 27 28 feholdexcept(&fenv); 29 res = rint(x); 30 fesetenv(&fenv); 31 32 return (res); 33 } 34 35 float 36 nearbyintf(float x) 37 { 38 fenv_t fenv; 39 float res; 40 41 feholdexcept(&fenv); 42 res = rintf(x); 43 fesetenv(&fenv); 44 45 return (res); 46 } 47 48 long double 49 nearbyintl(long double x) 50 { 51 fenv_t fenv; 52 long double res; 53 54 feholdexcept(&fenv); 55 res = rintl(x); 56 fesetenv(&fenv); 57 58 return (res); 59 } 60