diff options
| author | joerg <joerg@NetBSD.org> | 2011-08-03 14:13:07 +0000 |
|---|---|---|
| committer | joerg <joerg@NetBSD.org> | 2011-08-03 14:13:07 +0000 |
| commit | f73a9c72a8ea0a66b7bfa17fd400c20a45b6da2c (patch) | |
| tree | 108a805eafcf1919e8378fddf09a46b381c9065b /lib/libm | |
| parent | 7c2bed2aa85a8ae35827ce5b121e3df125194305 (diff) | |
Add logbl(3).
Diffstat (limited to 'lib/libm')
| -rw-r--r-- | lib/libm/Makefile | 12 | ||||
| -rw-r--r-- | lib/libm/arch/i387/s_logbl.S | 16 | ||||
| -rw-r--r-- | lib/libm/man/ieee_test.3 | 7 | ||||
| -rw-r--r-- | lib/libm/src/s_logb.c | 6 | ||||
| -rw-r--r-- | lib/libm/src/s_logbl.c | 76 |
5 files changed, 109 insertions, 8 deletions
diff --git a/lib/libm/Makefile b/lib/libm/Makefile index adad1b7e036..e1a31fcf84e 100644 --- a/lib/libm/Makefile +++ b/lib/libm/Makefile @@ -1,4 +1,4 @@ -# $NetBSD: Makefile,v 1.115 2011/07/28 22:32:28 joerg Exp $ +# $NetBSD: Makefile,v 1.116 2011/08/03 14:13:07 joerg Exp $ # # @(#)Makefile 5.1beta 93/09/24 # @@ -68,8 +68,9 @@ ARCH_SRCS = e_acos.S e_asin.S e_atan2.S e_exp.S e_expf.S e_fmod.S e_log.S \ e_logf.S e_log10.S e_log10f.S e_log2.S e_log2f.S e_remainder.S \ e_remainderf.S e_scalb.S e_scalbf.S e_sqrt.S e_sqrtf.S s_atan.S \ s_atanf.S s_ceil.S s_ceilf.S s_copysign.S s_copysignf.S s_cos.S \ - s_cosf.S s_finite.S s_finitef.S s_floor.S s_floorf.S s_ilogb.S \ - s_ilogbf.S s_ilogbl.S s_log1p.S s_log1pf.S s_logb.S s_logbf.S \ + s_cosf.S s_finite.S s_finitef.S s_floor.S s_floorf.S \ + s_ilogb.S s_ilogbf.S s_ilogbl.S s_log1p.S s_log1pf.S \ + s_logb.S s_logbf.S s_logbl.S \ s_rint.S s_rintf.S s_scalbn.S s_scalbnf.S s_significand.S \ s_significandf.S s_sin.S s_sinf.S s_tan.S s_tanf.S lrint.S # do not pick up the i387 asm version, it is incorrect @@ -139,7 +140,8 @@ COMMON_SRCS+= e_acos.c e_acosf.c e_acosh.c e_acoshf.c e_asin.c e_asinf.c \ s_finite.c s_finitef.c \ s_floor.c s_floorf.c s_frexpf.c s_ilogb.c s_ilogbf.c s_ilogbl.c \ s_isinff.c s_isnanf.c s_ldexpf.c s_lib_version.c s_log1p.c \ - s_log1pf.c s_logb.c s_logbf.c s_matherr.c s_modff.c s_nextafter.c \ + s_log1pf.c s_logb.c s_logbf.c s_logbl.c \ + s_matherr.c s_modff.c s_nextafter.c \ s_nextafterf.c s_remquo.c s_remquof.c s_rint.c s_rintf.c s_round.c s_roundf.c s_scalbn.c \ s_scalbnf.c s_scalbnl.c s_signgam.c s_significand.c s_significandf.c s_sin.c \ s_sinf.c s_tan.c s_tanf.c s_tanh.c s_tanhf.c s_trunc.c s_truncf.c \ @@ -254,7 +256,7 @@ MLINKS+=fmod.3 fmodf.3 MLINKS+=hypot.3 hypotf.3 MLINKS+=ieee.3 finite.3 ieee.3 finitef.3 \ ieee.3 scalbn.3 ieee.3 scalbnf.3 ieee.3 scalbnl.3 -MLINKS+=ieee_test.3 logb.3 ieee_test.3 logbf.3 +MLINKS+=ieee_test.3 logb.3 ieee_test.3 logbf.3 ieee_test.3 logbl.3 MLINKS+=ieee_test.3 scalb.3 ieee_test.3 scalbf.3 MLINKS+=ieee_test.3 significand.3 ieee_test.3 significandf.3 MLINKS+=ilogb.3 ilogbf.3 ilogb.3 ilogbl.3 diff --git a/lib/libm/arch/i387/s_logbl.S b/lib/libm/arch/i387/s_logbl.S new file mode 100644 index 00000000000..2048a967e12 --- /dev/null +++ b/lib/libm/arch/i387/s_logbl.S @@ -0,0 +1,16 @@ +/* + * Written by J.T. Conklin <jtc@NetBSD.org>. + * Public domain. + */ + +#include <machine/asm.h> + +#include "abi.h" + +RCSID("$NetBSD: s_logbl.S,v 1.1 2011/08/03 14:13:07 joerg Exp $") + +ENTRY(logbl) + fldt ARG_LONG_DOUBLE_ONE + fxtract + fstp %st + ret diff --git a/lib/libm/man/ieee_test.3 b/lib/libm/man/ieee_test.3 index c82e4cf6acb..4194a29313c 100644 --- a/lib/libm/man/ieee_test.3 +++ b/lib/libm/man/ieee_test.3 @@ -26,14 +26,15 @@ .\" SUCH DAMAGE. .\" .\" from: @(#)ieee.3 6.4 (Berkeley) 5/6/91 -.\" $NetBSD: ieee_test.3,v 1.10 2003/08/07 16:44:48 agc Exp $ +.\" $NetBSD: ieee_test.3,v 1.11 2011/08/03 14:13:07 joerg Exp $ .\" -.Dd March 10, 1994 +.Dd August 3, 2011 .Dt IEEE_TEST 3 .Os .Sh NAME .Nm logb , .Nm logbf , +.Nm logbl , .Nm scalb , .Nm scalbf , .Nm significand , @@ -47,6 +48,8 @@ .Fn logb "double x" .Ft float .Fn logbf "float x" +.Ft long double +.Fn logbl "long double x" .Ft double .Fn scalb "double x" "double n" .Ft float diff --git a/lib/libm/src/s_logb.c b/lib/libm/src/s_logb.c index 9a77e49e6dc..6d89e5428a2 100644 --- a/lib/libm/src/s_logb.c +++ b/lib/libm/src/s_logb.c @@ -12,7 +12,7 @@ #include <sys/cdefs.h> #if defined(LIBM_SCCS) && !defined(lint) -__RCSID("$NetBSD: s_logb.c,v 1.11 2002/05/26 22:01:57 wiz Exp $"); +__RCSID("$NetBSD: s_logb.c,v 1.12 2011/08/03 14:13:07 joerg Exp $"); #endif /* @@ -24,6 +24,10 @@ __RCSID("$NetBSD: s_logb.c,v 1.11 2002/05/26 22:01:57 wiz Exp $"); #include "math.h" #include "math_private.h" +#ifndef __HAVE_LONG_DOUBLE +__strong_alias(logbl,logb) +#endif + double logb(double x) { diff --git a/lib/libm/src/s_logbl.c b/lib/libm/src/s_logbl.c new file mode 100644 index 00000000000..42899905156 --- /dev/null +++ b/lib/libm/src/s_logbl.c @@ -0,0 +1,76 @@ +/* $NetBSD: s_logbl.c,v 1.1 2011/08/03 14:13:07 joerg Exp $ */ + +/*- + * Copyright (c) 2011 The NetBSD Foundation, Inc. + * All rights reserved. + * + * This code is derived from software contributed to The NetBSD Foundation + * by Joerg Sonnenberger. + * + * Redistribution and use in source and binary forms, with or without + * modification, are permitted provided that the following conditions + * are met: + * 1. Redistributions of source code must retain the above copyright + * notice, this list of conditions and the following disclaimer. + * 2. Redistributions in binary form must reproduce the above copyright + * notice, this list of conditions and the following disclaimer in the + * documentation and/or other materials provided with the distribution. + * + * THIS SOFTWARE IS PROVIDED BY THE NETBSD FOUNDATION, INC. AND CONTRIBUTORS + * ``AS IS'' AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED + * TO, THE IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR + * PURPOSE ARE DISCLAIMED. IN NO EVENT SHALL THE FOUNDATION OR CONTRIBUTORS + * BE LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR + * CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF + * SUBSTITUTE GOODS OR SERVICES; LOSS OF USE, DATA, OR PROFITS; OR BUSINESS + * INTERRUPTION) HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN + * CONTRACT, STRICT LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE) + * ARISING IN ANY WAY OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE + * POSSIBILITY OF SUCH DAMAGE. + */ + +#include <sys/cdefs.h> +__RCSID("$NetBSD: s_logbl.c,v 1.1 2011/08/03 14:13:07 joerg Exp $"); + +#include "namespace.h" + +#include <float.h> +#include <math.h> +#include <machine/ieee.h> + +#ifdef __HAVE_LONG_DOUBLE + +#if LDBL_MANT_DIG == 64 +#define FROM_UNDERFLOW 0x1p65L +#elif LDBL_MANT_DIG == 113 +#define FROM_UNDERFLOW 0x1p114L +#else +#error Unsupported long double format +#endif + +long double +logbl(long double x) +{ + union ieee_ext_u u; + + if (x == 0.0L) + return -1.0L / fabsl(x); /* -HUGE_VALL + exception */ + + u.extu_ld = x; + + if (u.extu_ext.ext_exp == EXT_EXP_INFNAN) + return fabsl(x); /* NaN or +Inf */ + + if (u.extu_ext.ext_exp == 0) { + /* + * Scale denormalized numbers slightly, + * so that they are normal. + */ + u.extu_ld *= FROM_UNDERFLOW; + return u.extu_ext.ext_exp - EXT_EXP_BIAS - LDBL_MANT_DIG - 1; + } + return u.extu_ext.ext_exp - EXT_EXP_BIAS; + +} + +#endif |
