1*84d9c625SLionel Sambuc /* $NetBSD: fabs_ieee754.c,v 1.3 2013/11/29 22:15:26 joerg Exp $ */
22fe8fb19SBen Gras
32fe8fb19SBen Gras /*-
42fe8fb19SBen Gras * Copyright (c) 2003 The NetBSD Foundation, Inc.
52fe8fb19SBen Gras * All rights reserved.
62fe8fb19SBen Gras *
72fe8fb19SBen Gras * This code is derived from software contributed to The NetBSD Foundation
82fe8fb19SBen Gras * by Klaus Klein.
92fe8fb19SBen Gras *
102fe8fb19SBen Gras * Redistribution and use in source and binary forms, with or without
112fe8fb19SBen Gras * modification, are permitted provided that the following conditions
122fe8fb19SBen Gras * are met:
132fe8fb19SBen Gras * 1. Redistributions of source code must retain the above copyright
142fe8fb19SBen Gras * notice, this list of conditions and the following disclaimer.
152fe8fb19SBen Gras * 2. Redistributions in binary form must reproduce the above copyright
162fe8fb19SBen Gras * notice, this list of conditions and the following disclaimer in the
172fe8fb19SBen Gras * documentation and/or other materials provided with the distribution.
182fe8fb19SBen Gras *
192fe8fb19SBen Gras * THIS SOFTWARE IS PROVIDED BY THE NETBSD FOUNDATION, INC. AND CONTRIBUTORS
202fe8fb19SBen Gras * ``AS IS'' AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED
212fe8fb19SBen Gras * TO, THE IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR
222fe8fb19SBen Gras * PURPOSE ARE DISCLAIMED. IN NO EVENT SHALL THE FOUNDATION OR CONTRIBUTORS
232fe8fb19SBen Gras * BE LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR
242fe8fb19SBen Gras * CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF
252fe8fb19SBen Gras * SUBSTITUTE GOODS OR SERVICES; LOSS OF USE, DATA, OR PROFITS; OR BUSINESS
262fe8fb19SBen Gras * INTERRUPTION) HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN
272fe8fb19SBen Gras * CONTRACT, STRICT LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE)
282fe8fb19SBen Gras * ARISING IN ANY WAY OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE
292fe8fb19SBen Gras * POSSIBILITY OF SUCH DAMAGE.
302fe8fb19SBen Gras */
312fe8fb19SBen Gras
322fe8fb19SBen Gras #include <sys/cdefs.h>
332fe8fb19SBen Gras #if defined(LIBC_SCCS) && !defined(lint)
34*84d9c625SLionel Sambuc __RCSID("$NetBSD: fabs_ieee754.c,v 1.3 2013/11/29 22:15:26 joerg Exp $");
352fe8fb19SBen Gras #endif
362fe8fb19SBen Gras
372fe8fb19SBen Gras #include <machine/ieee.h>
382fe8fb19SBen Gras #include <math.h>
392fe8fb19SBen Gras
402fe8fb19SBen Gras /*
412fe8fb19SBen Gras * 7.12.3.6 fabs - return the absolute value of the argument
422fe8fb19SBen Gras * IEEE 754 double-precision version
432fe8fb19SBen Gras */
44*84d9c625SLionel Sambuc #ifndef __HAVE_LONG_DOUBLE
__strong_alias(fabsl,fabs)45*84d9c625SLionel Sambuc __strong_alias(fabsl, fabs)
46*84d9c625SLionel Sambuc #endif
47*84d9c625SLionel Sambuc
482fe8fb19SBen Gras double
492fe8fb19SBen Gras fabs(double x)
502fe8fb19SBen Gras {
512fe8fb19SBen Gras union ieee_double_u u;
522fe8fb19SBen Gras
532fe8fb19SBen Gras u.dblu_d = x;
542fe8fb19SBen Gras
552fe8fb19SBen Gras u.dblu_dbl.dbl_sign = 0;
562fe8fb19SBen Gras
572fe8fb19SBen Gras return u.dblu_d;
582fe8fb19SBen Gras }
59