1/* $OpenBSD: round.c,v 1.3 2023/01/27 16:39:58 miod Exp $ */ 2 3/* Written by Michael Shalayeff, 2003, Public domain. */ 4 5#include <assert.h> 6#include <stdio.h> 7#include <stdlib.h> 8#include <string.h> 9#include <signal.h> 10#include <unistd.h> 11#include <math.h> 12#include <ieeefp.h> 13 14static void 15sigfpe(int sig, siginfo_t *si, void *v) 16{ 17 char buf[132]; 18 19 if (si) { 20 snprintf(buf, sizeof(buf), "sigfpe: addr=%p, code=%d\n", 21 si->si_addr, si->si_code); 22 write(1, buf, strlen(buf)); 23 } 24 _exit(1); 25} 26 27int 28main(int argc, char *argv[]) 29{ 30 struct sigaction sa; 31 32 memset(&sa, 0, sizeof(sa)); 33 sa.sa_sigaction = sigfpe; 34 sa.sa_flags = SA_SIGINFO; 35 sigaction(SIGFPE, &sa, NULL); 36 37 assert(round(8.6) == 9.); 38 assert(roundf(8.6F) == 9.); 39 assert(roundl(8.6L) == 9.); 40 assert(lround(8.6) == 9L); 41 assert(lroundf(8.6F) == 9L); 42 assert(lroundl(8.6L) == 9L); 43 assert(llround(8.6) == 9LL); 44 assert(llroundf(8.6F) == 9LL); 45 assert(llroundl(8.6L) == 9LL); 46 47 assert(lround(0.0) == 0L); 48 assert(lroundf(0.0F) == 0L); 49 assert(lroundl(0.0L) == 0L); 50 assert(lround(-0.0) == 0L); 51 assert(lroundf(-0.0F) == 0L); 52 assert(lroundl(-0.0L) == 0L); 53 54 /* Test proper rounding direction of halfway values. */ 55 assert(round(2.5) == 3.); 56 assert(roundf(2.5F) == 3.); 57 assert(roundl(2.5L) == 3.); 58 assert(round(-2.5) == -3.); 59 assert(roundf(-2.5F) == -3.); 60 assert(roundl(-2.5L) == -3.); 61 62 assert(llround(4503599627370496.0) == 4503599627370496LL); 63 assert(llroundf(4503599627370496.0F) == 4503599627370496LL); 64 assert(llroundl(4503599627370496.0L) == 4503599627370496LL); 65 assert(llround(-4503599627370496.0) == -4503599627370496LL); 66 assert(llroundf(-4503599627370496.0F) == -4503599627370496LL); 67 assert(llroundl(-4503599627370496.0L) == -4503599627370496LL); 68 69 assert(llround(0x7ffffffffffffc00.0p0) == 0x7ffffffffffffc00LL); 70 assert(llroundf(0x7fffff8000000000.0p0F) == 0x7fffff8000000000LL); 71 assert(llroundl(0x7fffff8000000000.0p0L) == 0x7fffff8000000000LL); 72 assert(llround(-0x8000000000000000.0p0) == -0x8000000000000000LL); 73 assert(llroundf(-0x8000000000000000.0p0F) == -0x8000000000000000LL); 74 assert(llroundl(-0x8000000000000000.0p0L) == -0x8000000000000000LL); 75 76 exit(0); 77} 78