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
14 static void
sigfpe(int sig,siginfo_t * si,void * v)15 sigfpe(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
27 int
main(int argc,char * argv[])28 main(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