1*c7933ddeSpatrick /* $OpenBSD: fabs.c,v 1.1 2017/01/11 18:09:24 patrick Exp $ */
2*c7933ddeSpatrick /*
3*c7933ddeSpatrick * Copyright (c) 2008 Martynas Venckus <martynas@openbsd.org>
4*c7933ddeSpatrick *
5*c7933ddeSpatrick * Permission to use, copy, modify, and distribute this software for any
6*c7933ddeSpatrick * purpose with or without fee is hereby granted, provided that the above
7*c7933ddeSpatrick * copyright notice and this permission notice appear in all copies.
8*c7933ddeSpatrick *
9*c7933ddeSpatrick * THE SOFTWARE IS PROVIDED "AS IS" AND THE AUTHOR DISCLAIMS ALL WARRANTIES
10*c7933ddeSpatrick * WITH REGARD TO THIS SOFTWARE INCLUDING ALL IMPLIED WARRANTIES OF
11*c7933ddeSpatrick * MERCHANTABILITY AND FITNESS. IN NO EVENT SHALL THE AUTHOR BE LIABLE FOR
12*c7933ddeSpatrick * ANY SPECIAL, DIRECT, INDIRECT, OR CONSEQUENTIAL DAMAGES OR ANY DAMAGES
13*c7933ddeSpatrick * WHATSOEVER RESULTING FROM LOSS OF USE, DATA OR PROFITS, WHETHER IN AN
14*c7933ddeSpatrick * ACTION OF CONTRACT, NEGLIGENCE OR OTHER TORTIOUS ACTION, ARISING OUT OF
15*c7933ddeSpatrick * OR IN CONNECTION WITH THE USE OR PERFORMANCE OF THIS SOFTWARE.
16*c7933ddeSpatrick */
17*c7933ddeSpatrick
18*c7933ddeSpatrick #include <sys/types.h>
19*c7933ddeSpatrick #include <machine/ieee.h>
20*c7933ddeSpatrick #include <math.h>
21*c7933ddeSpatrick
22*c7933ddeSpatrick /*
23*c7933ddeSpatrick * fabs(d) returns the absolute value of d.
24*c7933ddeSpatrick */
25*c7933ddeSpatrick double
fabs(double d)26*c7933ddeSpatrick fabs(double d)
27*c7933ddeSpatrick {
28*c7933ddeSpatrick struct ieee_double *p = (struct ieee_double *)&d;
29*c7933ddeSpatrick
30*c7933ddeSpatrick p->dbl_sign = 0;
31*c7933ddeSpatrick
32*c7933ddeSpatrick return(d);
33*c7933ddeSpatrick }
34*c7933ddeSpatrick
35*c7933ddeSpatrick __strong_alias(fabsl, fabs);
36