Module Name: src
Committed By: joerg
Date: Wed Aug 3 14:13:08 UTC 2011
Modified Files:
src/distrib/sets/lists/comp: mi
src/lib/libm: Makefile
src/lib/libm/man: ieee_test.3
src/lib/libm/src: s_logb.c
Added Files:
src/lib/libm/arch/i387: s_logbl.S
src/lib/libm/src: s_logbl.c
Log Message:
Add logbl(3).
To generate a diff of this commit:
cvs rdiff -u -r1.1652 -r1.1653 src/distrib/sets/lists/comp/mi
cvs rdiff -u -r1.115 -r1.116 src/lib/libm/Makefile
cvs rdiff -u -r0 -r1.1 src/lib/libm/arch/i387/s_logbl.S
cvs rdiff -u -r1.10 -r1.11 src/lib/libm/man/ieee_test.3
cvs rdiff -u -r1.11 -r1.12 src/lib/libm/src/s_logb.c
cvs rdiff -u -r0 -r1.1 src/lib/libm/src/s_logbl.c
Please note that diffs are not public domain; they are subject to the
copyright notices on the relevant files.
Modified files:
Index: src/distrib/sets/lists/comp/mi
diff -u src/distrib/sets/lists/comp/mi:1.1652 src/distrib/sets/lists/comp/mi:1.1653
--- src/distrib/sets/lists/comp/mi:1.1652 Wed Aug 3 04:11:15 2011
+++ src/distrib/sets/lists/comp/mi Wed Aug 3 14:13:07 2011
@@ -1,4 +1,4 @@
-# $NetBSD: mi,v 1.1652 2011/08/03 04:11:15 manu Exp $
+# $NetBSD: mi,v 1.1653 2011/08/03 14:13:07 joerg Exp $
#
# Note: don't delete entries from here - mark them as "obsolete" instead.
#
@@ -7682,6 +7682,7 @@
./usr/share/man/cat3/log2f.0 comp-c-catman .cat
./usr/share/man/cat3/logb.0 comp-c-catman .cat
./usr/share/man/cat3/logbf.0 comp-c-catman .cat
+./usr/share/man/cat3/logbl.0 comp-c-catman .cat
./usr/share/man/cat3/logf.0 comp-c-catman .cat
./usr/share/man/cat3/login.0 comp-c-catman .cat
./usr/share/man/cat3/login_cap.0 comp-c-catman .cat
@@ -13758,6 +13759,7 @@
./usr/share/man/html3/log2f.html comp-c-htmlman html
./usr/share/man/html3/logb.html comp-c-htmlman html
./usr/share/man/html3/logbf.html comp-c-htmlman html
+./usr/share/man/html3/logbl.html comp-c-htmlman html
./usr/share/man/html3/logf.html comp-c-htmlman html
./usr/share/man/html3/login.html comp-c-htmlman html
./usr/share/man/html3/login_cap.html comp-c-htmlman html
@@ -19817,6 +19819,7 @@
./usr/share/man/man3/log2f.3 comp-c-man .man
./usr/share/man/man3/logb.3 comp-c-man .man
./usr/share/man/man3/logbf.3 comp-c-man .man
+./usr/share/man/man3/logbl.3 comp-c-man .man
./usr/share/man/man3/logf.3 comp-c-man .man
./usr/share/man/man3/login.3 comp-c-man .man
./usr/share/man/man3/login_cap.3 comp-c-man .man
Index: src/lib/libm/Makefile
diff -u src/lib/libm/Makefile:1.115 src/lib/libm/Makefile:1.116
--- src/lib/libm/Makefile:1.115 Thu Jul 28 22:32:28 2011
+++ src/lib/libm/Makefile Wed Aug 3 14:13:07 2011
@@ -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 @@
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 @@
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+=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
Index: src/lib/libm/man/ieee_test.3
diff -u src/lib/libm/man/ieee_test.3:1.10 src/lib/libm/man/ieee_test.3:1.11
--- src/lib/libm/man/ieee_test.3:1.10 Thu Aug 7 16:44:48 2003
+++ src/lib/libm/man/ieee_test.3 Wed Aug 3 14:13:07 2011
@@ -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
Index: src/lib/libm/src/s_logb.c
diff -u src/lib/libm/src/s_logb.c:1.11 src/lib/libm/src/s_logb.c:1.12
--- src/lib/libm/src/s_logb.c:1.11 Sun May 26 22:01:57 2002
+++ src/lib/libm/src/s_logb.c Wed Aug 3 14:13:07 2011
@@ -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 @@
#include "math.h"
#include "math_private.h"
+#ifndef __HAVE_LONG_DOUBLE
+__strong_alias(logbl,logb)
+#endif
+
double
logb(double x)
{
Added files:
Index: src/lib/libm/arch/i387/s_logbl.S
diff -u /dev/null src/lib/libm/arch/i387/s_logbl.S:1.1
--- /dev/null Wed Aug 3 14:13:08 2011
+++ src/lib/libm/arch/i387/s_logbl.S Wed Aug 3 14:13:07 2011
@@ -0,0 +1,16 @@
+/*
+ * Written by J.T. Conklin <j...@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
Index: src/lib/libm/src/s_logbl.c
diff -u /dev/null src/lib/libm/src/s_logbl.c:1.1
--- /dev/null Wed Aug 3 14:13:08 2011
+++ src/lib/libm/src/s_logbl.c Wed Aug 3 14:13:07 2011
@@ -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
