1/*-
2 * Copyright (c) 2003, Steven G. Kargl
3 * All rights reserved.
4 *
5 * Redistribution and use in source and binary forms, with or without
6 * modification, are permitted provided that the following conditions
7 * are met:
8 * 1. Redistributions of source code must retain the above copyright
--- 11 unchanged lines hidden (view full) ---
20 * NOT LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES; LOSS OF USE,
21 * DATA, OR PROFITS; OR BUSINESS INTERRUPTION) HOWEVER CAUSED AND ON ANY
22 * THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT LIABILITY, OR TORT
23 * (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY OUT OF THE USE OF
24 * THIS SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF SUCH DAMAGE.
25 */
26
27#include <sys/cdefs.h>
28__FBSDID("$FreeBSD: head/lib/msun/src/s_round.c 133146 2004-08-05 01:44:55Z das $");
29
30#include <math.h>
31
32double
33round(double x)
34{
35 double t;
36
37 if (!isnormal(x))
38 return (x);
39
40 if (x >= 0.0) {
41 t = ceil(x);
42 if (t - x > 0.5)
43 t -= 1.0;
44 return (t);
45 } else {
46 t = ceil(-x);
47 if (t + x > 0.5)
48 t -= 1.0;
49 return (-t);
50 }
51}
2 * Copyright (c) 2003, Steven G. Kargl
3 * All rights reserved.
4 *
5 * Redistribution and use in source and binary forms, with or without
6 * modification, are permitted provided that the following conditions
7 * are met:
8 * 1. Redistributions of source code must retain the above copyright
--- 11 unchanged lines hidden (view full) ---
20 * NOT LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES; LOSS OF USE,
21 * DATA, OR PROFITS; OR BUSINESS INTERRUPTION) HOWEVER CAUSED AND ON ANY
22 * THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT LIABILITY, OR TORT
23 * (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY OUT OF THE USE OF
24 * THIS SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF SUCH DAMAGE.
25 */
26
27#include <sys/cdefs.h>
28__FBSDID("$FreeBSD: head/lib/msun/src/s_round.c 133146 2004-08-05 01:44:55Z das $");
29
30#include <math.h>
31
32double
33round(double x)
34{
35 double t;
36
37 if (!isnormal(x))
38 return (x);
39
40 if (x >= 0.0) {
41 t = ceil(x);
42 if (t - x > 0.5)
43 t -= 1.0;
44 return (t);
45 } else {
46 t = ceil(-x);
47 if (t + x > 0.5)
48 t -= 1.0;
49 return (-t);
50 }
51}