+/* Copyright 2001 by Stephen L. Moshier <moshier@na-net.ornl.gov>
+
+ This program is free software: you can redistribute it and/or modify
+ it under the terms of the GNU General Public License as published by
+ the Free Software Foundation; either version 3 of the License, or
+ (at your option) any later version.
+
+ This program is distributed in the hope that it will be useful,
+ but WITHOUT ANY WARRANTY; without even the implied warranty of
+ MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
+ GNU General Public License for more details.
+
+ You should have received a copy of the GNU General Public License
+ along with this program. If not, see <http://www.gnu.org/licenses/>. */
+
+#include <config.h>
+
+/* Specification. */
+#include <math.h>
+
+#if HAVE_SAME_LONG_DOUBLE_AS_DOUBLE
+
+long double
+logl (long double x)
+{
+ return log (x);
+}
+
+#elif HAVE_LOGL
+
+long double
+logl (long double x)
+# undef logl
+{
+ /* Work around the OSF/1 5.1 bug. */
+ if (x == 0.0L)
+ /* Return -Infinity. */
+ return -1.0L / 0.0L;
+ return logl (x);
+}
+
+#else
+
+/* Code based on glibc/sysdeps/ieee754/ldbl-128/e_logl.c. */
+
+/* logll.c