Mercurial > gnulib

annotate lib/logf.c @ 40186:8964917f9574

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
autoupdate
author Karl Berry <karl@freefriends.org>
date Mon, 18 Feb 2019 08:02:49 -0800
parents b06060465f09
children
Ignore whitespace changes - Everywhere: Within whitespace: At end of lines:
rev   line source
15863
c2fa2febdb1f New module 'logf'.
Bruno Haible <bruno@clisp.org>
parents:
diff changeset
1 /* Natural logarithmic function.
40057
b06060465f09 maint: Run 'make update-copyright'
Paul Eggert <eggert@cs.ucla.edu>
parents: 19484
diff changeset
2 Copyright (C) 2011-2019 Free Software Foundation, Inc.
15863
c2fa2febdb1f New module 'logf'.
Bruno Haible <bruno@clisp.org>
parents:
diff changeset
3
c2fa2febdb1f New module 'logf'.
Bruno Haible <bruno@clisp.org>
parents:
diff changeset
4 This program is free software: you can redistribute it and/or modify
c2fa2febdb1f New module 'logf'.
Bruno Haible <bruno@clisp.org>
parents:
diff changeset
5 it under the terms of the GNU General Public License as published by
c2fa2febdb1f New module 'logf'.
Bruno Haible <bruno@clisp.org>
parents:
diff changeset
6 the Free Software Foundation; either version 3 of the License, or
c2fa2febdb1f New module 'logf'.
Bruno Haible <bruno@clisp.org>
parents:
diff changeset
7 (at your option) any later version.
c2fa2febdb1f New module 'logf'.
Bruno Haible <bruno@clisp.org>
parents:
diff changeset
8
c2fa2febdb1f New module 'logf'.
Bruno Haible <bruno@clisp.org>
parents:
diff changeset
9 This program is distributed in the hope that it will be useful,
c2fa2febdb1f New module 'logf'.
Bruno Haible <bruno@clisp.org>
parents:
diff changeset
10 but WITHOUT ANY WARRANTY; without even the implied warranty of
c2fa2febdb1f New module 'logf'.
Bruno Haible <bruno@clisp.org>
parents:
diff changeset
11 MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
c2fa2febdb1f New module 'logf'.
Bruno Haible <bruno@clisp.org>
parents:
diff changeset
12 GNU General Public License for more details.
c2fa2febdb1f New module 'logf'.
Bruno Haible <bruno@clisp.org>
parents:
diff changeset
13
c2fa2febdb1f New module 'logf'.
Bruno Haible <bruno@clisp.org>
parents:
diff changeset
14 You should have received a copy of the GNU General Public License
19190
9759915b2aca all: prefer https: URLs
Paul Eggert <eggert@cs.ucla.edu>
parents: 18626
diff changeset
15 along with this program. If not, see <https://www.gnu.org/licenses/>. */
15863
c2fa2febdb1f New module 'logf'.
Bruno Haible <bruno@clisp.org>
parents:
diff changeset
16
c2fa2febdb1f New module 'logf'.
Bruno Haible <bruno@clisp.org>
parents:
diff changeset
17 #include <config.h>
c2fa2febdb1f New module 'logf'.
Bruno Haible <bruno@clisp.org>
parents:
diff changeset
18
c2fa2febdb1f New module 'logf'.
Bruno Haible <bruno@clisp.org>
parents:
diff changeset
19 /* Specification. */
c2fa2febdb1f New module 'logf'.
Bruno Haible <bruno@clisp.org>
parents:
diff changeset
20 #include <math.h>
c2fa2febdb1f New module 'logf'.
Bruno Haible <bruno@clisp.org>
parents:
diff changeset
21
c2fa2febdb1f New module 'logf'.
Bruno Haible <bruno@clisp.org>
parents:
diff changeset
22 float
c2fa2febdb1f New module 'logf'.
Bruno Haible <bruno@clisp.org>
parents:
diff changeset
23 logf (float x)
16632
2450f24b2cbc logf: Work around OSF/1 5.1 bug.
Bruno Haible <bruno@clisp.org>
parents: 16201
diff changeset
24 #undef logf
15863
c2fa2febdb1f New module 'logf'.
Bruno Haible <bruno@clisp.org>
parents:
diff changeset
25 {
16632
2450f24b2cbc logf: Work around OSF/1 5.1 bug.
Bruno Haible <bruno@clisp.org>
parents: 16201
diff changeset
26 #if HAVE_LOGF
16645
d7a6433b4967 logf-ieee: Work around test failure on NetBSD 5.1.
Bruno Haible <bruno@clisp.org>
parents: 16632
diff changeset
27 if (x <= 0.0f)
d7a6433b4967 logf-ieee: Work around test failure on NetBSD 5.1.
Bruno Haible <bruno@clisp.org>
parents: 16632
diff changeset
28 {
d7a6433b4967 logf-ieee: Work around test failure on NetBSD 5.1.
Bruno Haible <bruno@clisp.org>
parents: 16632
diff changeset
29 /* Work around the OSF/1 5.1 bug. */
d7a6433b4967 logf-ieee: Work around test failure on NetBSD 5.1.
Bruno Haible <bruno@clisp.org>
parents: 16632
diff changeset
30 if (x == 0.0f)
d7a6433b4967 logf-ieee: Work around test failure on NetBSD 5.1.
Bruno Haible <bruno@clisp.org>
parents: 16632
diff changeset
31 /* Return -Infinity. */
d7a6433b4967 logf-ieee: Work around test failure on NetBSD 5.1.
Bruno Haible <bruno@clisp.org>
parents: 16632
diff changeset
32 return -1.0f / 0.0f;
d7a6433b4967 logf-ieee: Work around test failure on NetBSD 5.1.
Bruno Haible <bruno@clisp.org>
parents: 16632
diff changeset
33 /* Work around the NetBSD 5.1 bug. */
d7a6433b4967 logf-ieee: Work around test failure on NetBSD 5.1.
Bruno Haible <bruno@clisp.org>
parents: 16632
diff changeset
34 else /* x < 0.0 */
d7a6433b4967 logf-ieee: Work around test failure on NetBSD 5.1.
Bruno Haible <bruno@clisp.org>
parents: 16632
diff changeset
35 /* Return NaN. */
d7a6433b4967 logf-ieee: Work around test failure on NetBSD 5.1.
Bruno Haible <bruno@clisp.org>
parents: 16632
diff changeset
36 return 0.0f / 0.0f;
d7a6433b4967 logf-ieee: Work around test failure on NetBSD 5.1.
Bruno Haible <bruno@clisp.org>
parents: 16632
diff changeset
37 }
16632
2450f24b2cbc logf: Work around OSF/1 5.1 bug.
Bruno Haible <bruno@clisp.org>
parents: 16201
diff changeset
38 return logf (x);
2450f24b2cbc logf: Work around OSF/1 5.1 bug.
Bruno Haible <bruno@clisp.org>
parents: 16201
diff changeset
39 #else
15863
c2fa2febdb1f New module 'logf'.
Bruno Haible <bruno@clisp.org>
parents:
diff changeset
40 return (float) log ((double) x);
16632
2450f24b2cbc logf: Work around OSF/1 5.1 bug.
Bruno Haible <bruno@clisp.org>
parents: 16201
diff changeset
41 #endif
15863
c2fa2febdb1f New module 'logf'.
Bruno Haible <bruno@clisp.org>
parents:
diff changeset
42 }