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