Mercurial > hg > octave-kai > gnulib-hg

annotate lib/log10f.c @ 17342:c75939cb6254

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
merge with default branch
author Michael Goffioul <michael.goffioul@gmail.com>
date Thu, 21 Feb 2013 14:57:31 +0000
parents e542fd46ad6f
children
Ignore whitespace changes - Everywhere: Within whitespace: At end of lines:
rev   line source
15866
215da3b412e0 New module 'log10f'.
Bruno Haible <bruno@clisp.org>
parents:
diff changeset
1 /* Base 10 logarithmic function.
17249
e542fd46ad6f maint: update all copyright year number ranges
Eric Blake <eblake@redhat.com>
parents: 16747
diff changeset
2 Copyright (C) 2011-2013 Free Software Foundation, Inc.
15866
215da3b412e0 New module 'log10f'.
Bruno Haible <bruno@clisp.org>
parents:
diff changeset
3
215da3b412e0 New module 'log10f'.
Bruno Haible <bruno@clisp.org>
parents:
diff changeset
4 This program is free software: you can redistribute it and/or modify
215da3b412e0 New module 'log10f'.
Bruno Haible <bruno@clisp.org>
parents:
diff changeset
5 it under the terms of the GNU General Public License as published by
215da3b412e0 New module 'log10f'.
Bruno Haible <bruno@clisp.org>
parents:
diff changeset
6 the Free Software Foundation; either version 3 of the License, or
215da3b412e0 New module 'log10f'.
Bruno Haible <bruno@clisp.org>
parents:
diff changeset
7 (at your option) any later version.
215da3b412e0 New module 'log10f'.
Bruno Haible <bruno@clisp.org>
parents:
diff changeset
8
215da3b412e0 New module 'log10f'.
Bruno Haible <bruno@clisp.org>
parents:
diff changeset
9 This program is distributed in the hope that it will be useful,
215da3b412e0 New module 'log10f'.
Bruno Haible <bruno@clisp.org>
parents:
diff changeset
10 but WITHOUT ANY WARRANTY; without even the implied warranty of
215da3b412e0 New module 'log10f'.
Bruno Haible <bruno@clisp.org>
parents:
diff changeset
11 MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
215da3b412e0 New module 'log10f'.
Bruno Haible <bruno@clisp.org>
parents:
diff changeset
12 GNU General Public License for more details.
215da3b412e0 New module 'log10f'.
Bruno Haible <bruno@clisp.org>
parents:
diff changeset
13
215da3b412e0 New module 'log10f'.
Bruno Haible <bruno@clisp.org>
parents:
diff changeset
14 You should have received a copy of the GNU General Public License
215da3b412e0 New module 'log10f'.
Bruno Haible <bruno@clisp.org>
parents:
diff changeset
15 along with this program. If not, see <http://www.gnu.org/licenses/>. */
215da3b412e0 New module 'log10f'.
Bruno Haible <bruno@clisp.org>
parents:
diff changeset
16
215da3b412e0 New module 'log10f'.
Bruno Haible <bruno@clisp.org>
parents:
diff changeset
17 #include <config.h>
215da3b412e0 New module 'log10f'.
Bruno Haible <bruno@clisp.org>
parents:
diff changeset
18
215da3b412e0 New module 'log10f'.
Bruno Haible <bruno@clisp.org>
parents:
diff changeset
19 /* Specification. */
215da3b412e0 New module 'log10f'.
Bruno Haible <bruno@clisp.org>
parents:
diff changeset
20 #include <math.h>
215da3b412e0 New module 'log10f'.
Bruno Haible <bruno@clisp.org>
parents:
diff changeset
21
215da3b412e0 New module 'log10f'.
Bruno Haible <bruno@clisp.org>
parents:
diff changeset
22 float
215da3b412e0 New module 'log10f'.
Bruno Haible <bruno@clisp.org>
parents:
diff changeset
23 log10f (float x)
16737
6828b991e3b9 log10f: Work around OSF/1 5.1 bug.
Bruno Haible <bruno@clisp.org>
parents: 16201
diff changeset
24 #undef log10f
15866
215da3b412e0 New module 'log10f'.
Bruno Haible <bruno@clisp.org>
parents:
diff changeset
25 {
16737
6828b991e3b9 log10f: Work around OSF/1 5.1 bug.
Bruno Haible <bruno@clisp.org>
parents: 16201
diff changeset
26 #if HAVE_LOG10F
16747
676471a1cb33 log10f-ieee: Work around test failure on NetBSD 5.1.
Bruno Haible <bruno@clisp.org>
parents: 16737
diff changeset
27 if (x <= 0.0f)
676471a1cb33 log10f-ieee: Work around test failure on NetBSD 5.1.
Bruno Haible <bruno@clisp.org>
parents: 16737
diff changeset
28 {
676471a1cb33 log10f-ieee: Work around test failure on NetBSD 5.1.
Bruno Haible <bruno@clisp.org>
parents: 16737
diff changeset
29 /* Work around the OSF/1 5.1 bug. */
676471a1cb33 log10f-ieee: Work around test failure on NetBSD 5.1.
Bruno Haible <bruno@clisp.org>
parents: 16737
diff changeset
30 if (x == 0.0f)
676471a1cb33 log10f-ieee: Work around test failure on NetBSD 5.1.
Bruno Haible <bruno@clisp.org>
parents: 16737
diff changeset
31 /* Return -Infinity. */
676471a1cb33 log10f-ieee: Work around test failure on NetBSD 5.1.
Bruno Haible <bruno@clisp.org>
parents: 16737
diff changeset
32 return -1.0f / 0.0f;
676471a1cb33 log10f-ieee: Work around test failure on NetBSD 5.1.
Bruno Haible <bruno@clisp.org>
parents: 16737
diff changeset
33 /* Work around the NetBSD 5.1 bug. */
676471a1cb33 log10f-ieee: Work around test failure on NetBSD 5.1.
Bruno Haible <bruno@clisp.org>
parents: 16737
diff changeset
34 else /* x < 0.0 */
676471a1cb33 log10f-ieee: Work around test failure on NetBSD 5.1.
Bruno Haible <bruno@clisp.org>
parents: 16737
diff changeset
35 /* Return NaN. */
676471a1cb33 log10f-ieee: Work around test failure on NetBSD 5.1.
Bruno Haible <bruno@clisp.org>
parents: 16737
diff changeset
36 return 0.0f / 0.0f;
676471a1cb33 log10f-ieee: Work around test failure on NetBSD 5.1.
Bruno Haible <bruno@clisp.org>
parents: 16737
diff changeset
37 }
16737
6828b991e3b9 log10f: Work around OSF/1 5.1 bug.
Bruno Haible <bruno@clisp.org>
parents: 16201
diff changeset
38 return log10f (x);
6828b991e3b9 log10f: Work around OSF/1 5.1 bug.
Bruno Haible <bruno@clisp.org>
parents: 16201
diff changeset
39 #else
15866
215da3b412e0 New module 'log10f'.
Bruno Haible <bruno@clisp.org>
parents:
diff changeset
40 return (float) log10 ((double) x);
16737
6828b991e3b9 log10f: Work around OSF/1 5.1 bug.
Bruno Haible <bruno@clisp.org>
parents: 16201
diff changeset
41 #endif
15866
215da3b412e0 New module 'log10f'.
Bruno Haible <bruno@clisp.org>
parents:
diff changeset
42 }