summaryrefslogtreecommitdiff
path: root/math
diff options
context:
space:
mode:
authorRoland McGrath <roland@gnu.org>1995-02-18 01:27:10 +0000
committerRoland McGrath <roland@gnu.org>1995-02-18 01:27:10 +0000
commit28f540f45bbacd939bfd07f213bcad2bf730b1bf (patch)
tree15f07c4c43d635959c6afee96bde71fb1b3614ee /math
downloadglibc-28f540f45bbacd939bfd07f213bcad2bf730b1bf.tar.gz
initial import
Diffstat (limited to 'math')
-rw-r--r--math/.cvsignore4
-rw-r--r--math/Makefile39
-rw-r--r--math/finite.c31
-rw-r--r--math/math.h317
-rw-r--r--math/test-math.c148
5 files changed, 539 insertions, 0 deletions
diff --git a/math/.cvsignore b/math/.cvsignore
new file mode 100644
index 0000000000..1f69fd919a
--- /dev/null
+++ b/math/.cvsignore
@@ -0,0 +1,4 @@
+*.gz *.Z *.tar *.tgz
+=*
+TODO COPYING* AUTHORS copyr-* copying.*
+glibc-*
diff --git a/math/Makefile b/math/Makefile
new file mode 100644
index 0000000000..c86d44a072
--- /dev/null
+++ b/math/Makefile
@@ -0,0 +1,39 @@
+# Copyright (C) 1991, 1992, 1994, 1995 Free Software Foundation, Inc.
+# This file is part of the GNU C Library.
+
+# The GNU C Library is free software; you can redistribute it and/or
+# modify it under the terms of the GNU Library General Public License as
+# published by the Free Software Foundation; either version 2 of the
+# License, or (at your option) any later version.
+
+# The GNU C Library 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
+# Library General Public License for more details.
+
+# You should have received a copy of the GNU Library General Public
+# License along with the GNU C Library; see the file COPYING.LIB. If
+# not, write to the Free Software Foundation, Inc., 675 Mass Ave,
+# Cambridge, MA 02139, USA.
+
+#
+# Makefile for math.
+#
+subdir := math
+
+headers := math.h __math.h huge_val.h nan.h
+
+
+routines := acos asin atan cos sin tan cosh sinh tanh exp fabs ldexp \
+ log log10 floor sqrt fmod frexp pow atan2 ceil modf \
+ isinf isnan finite infnan copysign drem logb \
+ rint hypot cabs cbrt expm1 log1p acosh asinh atanh \
+ isinfl isnanl
+tests := # test-math
+install-lib := libm.a
+
+
+include ../Rules
+
+$(objpfx)libm.a: $(dep-dummy-lib); $(make-dummy-lib)
+lib: $(objpfx)libm.a
diff --git a/math/finite.c b/math/finite.c
new file mode 100644
index 0000000000..48f8f8011a
--- /dev/null
+++ b/math/finite.c
@@ -0,0 +1,31 @@
+/* Copyright (C) 1991, 1992, 1995 Free Software Foundation, Inc.
+This file is part of the GNU C Library.
+
+The GNU C Library is free software; you can redistribute it and/or
+modify it under the terms of the GNU Library General Public License as
+published by the Free Software Foundation; either version 2 of the
+License, or (at your option) any later version.
+
+The GNU C Library 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
+Library General Public License for more details.
+
+You should have received a copy of the GNU Library General Public
+License along with the GNU C Library; see the file COPYING.LIB. If
+not, write to the Free Software Foundation, Inc., 675 Mass Ave,
+Cambridge, MA 02139, USA. */
+
+#include <ansidecl.h>
+#include <math.h>
+
+#undef __finite
+
+/* Return nonzero if VALUE is finite and not NaN. */
+int
+DEFUN(__finite, (value), double value)
+{
+ return !__isinf (value) && !__isnan (value);
+}
+
+weak_alias (__finite, finite)
diff --git a/math/math.h b/math/math.h
new file mode 100644
index 0000000000..51ba57c088
--- /dev/null
+++ b/math/math.h
@@ -0,0 +1,317 @@
+/* Copyright (C) 1991, 1992, 1993 Free Software Foundation, Inc.
+This file is part of the GNU C Library.
+
+The GNU C Library is free software; you can redistribute it and/or
+modify it under the terms of the GNU Library General Public License as
+published by the Free Software Foundation; either version 2 of the
+License, or (at your option) any later version.
+
+The GNU C Library 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
+Library General Public License for more details.
+
+You should have received a copy of the GNU Library General Public
+License along with the GNU C Library; see the file COPYING.LIB. If
+not, write to the, 1992 Free Software Foundation, Inc., 675 Mass Ave,
+Cambridge, MA 02139, USA. */
+
+/*
+ * ANSI Standard: 4.5 MATHEMATICS <math.h>
+ */
+
+#ifndef _MATH_H
+
+#define _MATH_H 1
+#include <features.h>
+
+__BEGIN_DECLS
+
+#define __need_Emath
+#include <errno.h>
+
+/* Get machine-dependent HUGE_VAL value (returned on overflow). */
+#include <huge_val.h>
+
+/* Get machine-dependent NAN value (returned for some domain errors). */
+#ifdef __USE_GNU
+#include <nan.h>
+#endif
+
+
+/* Trigonometric functions. */
+
+/* Arc cosine of X. */
+extern double acos __P ((double __x)) __attribute__ ((__const__));
+/* Arc sine of X. */
+extern double asin __P ((double __x)) __attribute__ ((__const__));
+/* Arc tangent of X. */
+extern double atan __P ((double __x)) __attribute__ ((__const__));
+/* Arc tangent of Y/X. */
+extern double atan2 __P ((double __y, double __x)) __attribute__ ((__const__));
+
+/* Cosine of X. */
+extern double cos __P ((double __x)) __attribute__ ((__const__));
+/* Sine of X. */
+extern double sin __P ((double __x)) __attribute__ ((__const__));
+/* Tangent of X. */
+extern double tan __P ((double __x)) __attribute__ ((__const__));
+
+
+/* Hyperbolic functions. */
+
+/* Hyperbolic cosine of X. */
+extern double cosh __P ((double __x)) __attribute__ ((__const__));
+/* Hyperbolic sine of X. */
+extern double sinh __P ((double __x)) __attribute__ ((__const__));
+/* Hyperbolic tangent of X. */
+extern double tanh __P ((double __x)) __attribute__ ((__const__));
+
+#ifdef __USE_MISC
+/* Hyperbolic arc cosine of X. */
+extern double acosh __P ((double __x)) __attribute__ ((__const__));
+/* Hyperbolic arc sine of X. */
+extern double asinh __P ((double __x)) __attribute__ ((__const__));
+/* Hyperbolic arc tangent of X. */
+extern double atanh __P ((double __x)) __attribute__ ((__const__));
+#endif
+
+/* Exponential and logarithmic functions. */
+
+/* Exponentional function of X. */
+extern double exp __P ((double __x)) __attribute__ ((__const__));
+
+/* Break VALUE into a normalized fraction and an integral power of 2. */
+extern double frexp __P ((double __value, int *__exp));
+
+/* X times (two to the EXP power). */
+extern double ldexp __P ((double __x, int __exp)) __attribute__ ((__const__));
+
+/* Natural logarithm of X. */
+extern double log __P ((double __x)) __attribute__ ((__const__));
+
+/* Base-ten logarithm of X. */
+extern double log10 __P ((double __x)) __attribute__ ((__const__));
+
+#ifdef __USE_MISC
+/* Return exp(X) - 1. */
+extern double __expm1 __P ((double __x)) __attribute__ ((__const__));
+extern double expm1 __P ((double __x)) __attribute__ ((__const__));
+
+/* Return log(1 + X). */
+extern double log1p __P ((double __x)) __attribute__ ((__const__));
+#endif
+
+/* Break VALUE into integral and fractional parts. */
+extern double modf __P ((double __value, double *__iptr));
+
+
+/* Power functions. */
+
+/* Return X to the Y power. */
+extern double pow __P ((double __x, double __y)) __attribute__ ((__const__));
+
+/* Return the square root of X. */
+extern double sqrt __P ((double __x)) __attribute__ ((__const__));
+
+#ifdef __USE_MISC
+/* Return the cube root of X. */
+extern double cbrt __P ((double __x)) __attribute__ ((__const__));
+#endif
+
+
+/* Nearest integer, absolute value, and remainder functions. */
+
+/* Smallest integral value not less than X. */
+extern double ceil __P ((double __x)) __attribute__ ((__const__));
+
+/* Absolute value of X. */
+extern double fabs __P ((double __x)) __attribute__ ((__const__));
+
+/* Largest integer not greater than X. */
+extern double floor __P ((double __x)) __attribute__ ((__const__));
+
+/* Floating-point modulo remainder of X/Y. */
+extern double fmod __P ((double __x, double __y)) __attribute__ ((__const__));
+
+
+/* Return 0 if VALUE is finite or NaN, +1 if it
+ is +Infinity, -1 if it is -Infinity. */
+extern int __isinf __P ((double __value)) __attribute__ ((__const__));
+
+/* Return nonzero if VALUE is not a number. */
+extern int __isnan __P ((double __value)) __attribute__ ((__const__));
+
+/* Return nonzero if VALUE is finite and not NaN. */
+extern int __finite __P ((double __value)) __attribute__ ((__const__));
+#ifdef __OPTIMIZE__
+#define __finite(value) (!__isinf(value))
+#endif
+
+/* Deal with an infinite or NaN result.
+ If ERROR is ERANGE, result is +Inf;
+ if ERROR is - ERANGE, result is -Inf;
+ otherwise result is NaN.
+ This will set `errno' to either ERANGE or EDOM,
+ and may return an infinity or NaN, or may do something else. */
+extern double __infnan __P ((int __error));
+
+/* Return X with its signed changed to Y's. */
+extern double __copysign __P ((double __x, double __y))
+ __attribute__ ((__const__));
+
+/* Return X times (2 to the Nth power). */
+extern double __scalb __P ((double __x, int __n))
+ __attribute__ ((__const__));
+
+#ifdef __OPTIMIZE__
+#define __scalb(x, n) ldexp ((x), (n))
+#endif
+
+/* Return the remainder of X/Y. */
+extern double __drem __P ((double __x, double __y))
+ __attribute__ ((__const__));
+
+/* Return the base 2 signed integral exponent of X. */
+extern double __logb __P ((double __x)) __attribute__ ((__const__));
+
+#ifdef __USE_MISC
+
+/* Return the integer nearest X in the direction of the
+ prevailing rounding mode. */
+extern double __rint __P ((double __x)) __attribute__ ((__const__));
+extern double rint __P ((double __x)) __attribute__ ((__const__));
+
+/* Return `sqrt(X*X + Y*Y)'. */
+extern double hypot __P ((double __x, double __y)) __attribute__ ((__const__));
+
+struct __cabs_complex
+{
+ double __x, __y;
+};
+
+/* Return `sqrt(X*X + Y*Y)'. */
+extern double cabs __P ((struct __cabs_complex)) __attribute__ ((__const__));
+
+extern int isinf __P ((double __value)) __attribute__ ((__const__));
+extern int isnan __P ((double __value)) __attribute__ ((__const__));
+extern int finite __P ((double __value)) __attribute__ ((__const__));
+extern double infnan __P ((int __error)) __attribute__ ((__const__));
+extern double copysign __P ((double __x, double __y))
+ __attribute__ ((__const__));
+extern double scalb __P ((double __x, int __n)) __attribute__ ((__const__));
+extern double drem __P ((double __x, double __y)) __attribute__ ((__const__));
+extern double logb __P ((double __x)) __attribute__ ((__const__));
+
+#ifdef __OPTIMIZE__
+#define isinf(value) __isinf(value)
+#define isnan(value) __isnan(value)
+#define infnan(error) __infnan(error)
+#define finite(value) __finite(value)
+#define copysign(x, y) __copysign((x), (y))
+#define scalb(x, n) __scalb((x), (n))
+#define drem(x, y) __drem((x), (y))
+#define logb(x) __logb(x)
+#endif /* Optimizing. */
+
+#endif /* Use misc. */
+
+
+#if 0
+/* The "Future Library Directions" section of the
+ ANSI Standard reserves these as `float' and
+ `long double' versions of the above functions. */
+
+extern float acosf __P ((float __x)) __attribute__ ((__const__));
+extern float asinf __P ((float __x)) __attribute__ ((__const__));
+extern float atanf __P ((float __x)) __attribute__ ((__const__));
+extern float atan2f __P ((float __y, float __x)) __attribute__ ((__const__));
+extern float cosf __P ((float __x)) __attribute__ ((__const__));
+extern float sinf __P ((float __x)) __attribute__ ((__const__));
+extern float tanf __P ((float __x)) __attribute__ ((__const__));
+extern float coshf __P ((float __x)) __attribute__ ((__const__));
+extern float sinhf __P ((float __x)) __attribute__ ((__const__));
+extern float tanhf __P ((float __x)) __attribute__ ((__const__));
+extern float expf __P ((float __x)) __attribute__ ((__const__));
+extern float frexpf __P ((float __value, int *__exp));
+extern float ldexpf __P ((float __x, int __exp)) __attribute__ ((__const__));
+extern float logf __P ((float __x)) __attribute__ ((__const__));
+extern float log10f __P ((float __x)) __attribute__ ((__const__));
+extern float modff __P ((float __value, float *__iptr));
+extern float powf __P ((float __x, float __y)) __attribute__ ((__const__));
+extern float sqrtf __P ((float __x)) __attribute__ ((__const__));
+extern float ceilf __P ((float __x)) __attribute__ ((__const__));
+extern float fabsf __P ((float __x)) __attribute__ ((__const__));
+extern float floorf __P ((float __x)) __attribute__ ((__const__));
+extern float fmodf __P ((float __x, float __y)) __attribute__ ((__const__));
+
+extern __long_double_t acosl __P ((__long_double_t __x))
+ __attribute__ ((__const__));
+extern __long_double_t asinl __P ((__long_double_t __x))
+ __attribute__ ((__const__));
+extern __long_double_t atanl __P ((__long_double_t __x))
+ __attribute__ ((__const__));
+extern __long_double_t atan2l __P ((__long_double_t __y, __long_double_t __x))
+ __attribute__ ((__const__));
+extern __long_double_t cosl __P ((__long_double_t __x))
+ __attribute__ ((__const__));
+extern __long_double_t sinl __P ((__long_double_t __x))
+ __attribute__ ((__const__));
+extern __long_double_t tanl __P ((__long_double_t __x))
+ __attribute__ ((__const__));
+extern __long_double_t coshl __P ((__long_double_t __x))
+ __attribute__ ((__const__));
+extern __long_double_t sinhl __P ((__long_double_t __x))
+ __attribute__ ((__const__));
+extern __long_double_t tanhl __P ((__long_double_t __x))
+ __attribute__ ((__const__));
+extern __long_double_t expl __P ((__long_double_t __x))
+ __attribute__ ((__const__));
+extern __long_double_t frexpl __P ((__long_double_t __value, int *__exp));
+extern __long_double_t ldexpl __P ((__long_double_t __x, int __exp))
+ __attribute__ ((__const__));
+extern __long_double_t logl __P ((__long_double_t __x))
+ __attribute__ ((__const__));
+extern __long_double_t log10l __P ((__long_double_t __x))
+ __attribute__ ((__const__));
+extern __long_double_t modfl __P ((__long_double_t __value,
+ __long_double_t *__ip));
+extern __long_double_t powl __P ((__long_double_t __x, __long_double_t __y))
+ __attribute__ ((__const__));
+extern __long_double_t sqrtl __P ((__long_double_t __x))
+ __attribute__ ((__const__));
+extern __long_double_t ceill __P ((__long_double_t __x))
+ __attribute__ ((__const__));
+extern __long_double_t fabsl __P ((__long_double_t __x))
+ __attribute__ ((__const__));
+extern __long_double_t floorl __P ((__long_double_t __x))
+ __attribute__ ((__const__));
+extern __long_double_t fmodl __P ((__long_double_t __x, __long_double_t __y))
+ __attribute__ ((__const__));
+#endif /* 0 */
+
+/* Get machine-dependent inline versions (if there are any). */
+#include <__math.h>
+
+__END_DECLS
+
+
+#ifdef __USE_BSD
+/* Some useful constants. */
+#define M_E 2.7182818284590452354 /* e */
+#define M_LOG2E 1.4426950408889634074 /* log 2e */
+#define M_LOG10E 0.43429448190325182765 /* log 10e */
+#define M_LN2 0.69314718055994530942 /* log e2 */
+#define M_LN10 2.30258509299404568402 /* log e10 */
+#define M_PI 3.14159265358979323846 /* pi */
+#define M_PI_2 1.57079632679489661923 /* pi/2 */
+#define M_PI_4 0.78539816339744830962 /* pi/4 */
+#define M_1_PI 0.31830988618379067154 /* 1/pi */
+#define M_2_PI 0.63661977236758134308 /* 2/pi */
+#define M_2_SQRTPI 1.12837916709551257390 /* 2/sqrt(pi) */
+#define M_SQRT2 1.41421356237309504880 /* sqrt(2) */
+#define M_SQRT1_2 0.70710678118654752440 /* 1/sqrt(2) */
+#endif
+
+
+#endif /* math.h */
diff --git a/math/test-math.c b/math/test-math.c
new file mode 100644
index 0000000000..a82daed4e4
--- /dev/null
+++ b/math/test-math.c
@@ -0,0 +1,148 @@
+#include <ansidecl.h>
+#include <stdio.h>
+#include <math.h>
+#include <stdlib.h>
+
+void print_trig_stuff __P ((void));
+
+int
+DEFUN_VOID(main)
+{
+ CONST char str[] = "123.456";
+ double x,h,li,lr,a,lrr;
+
+ x = atof (str);
+
+ printf ("%g %g\n", x, pow (10.0, 3.0));
+
+ x = sinh(2.0);
+
+ printf("sinh(2.0) = %g\n", x);
+
+ x = sinh(3.0);
+
+ printf("sinh(3.0) = %g\n", x);
+
+ h = hypot(2.0,3.0);
+
+ printf("h=%g\n", h);
+
+ a = atan2(3.0, 2.0);
+
+ printf("atan2(3,2) = %g\n", a);
+
+ lr = pow(h,4.0);
+
+ printf("pow(%g,4.0) = %g\n", h, lr);
+
+ lrr = lr;
+
+ li = 4.0 * a;
+
+ lr = lr / exp(a*5.0);
+
+ printf("%g / exp(%g * 5) = %g\n", lrr, a, lr);
+
+ lrr = li;
+
+ li += 5.0 * log(h);
+
+ printf("%g + 5*log(%g) = %g\n", lrr, h, li);
+
+ printf("cos(%g) = %g, sin(%g) = %g\n", li, cos(li), li, sin(li));
+
+ x = drem(10.3435,6.2831852);
+
+ printf("drem(10.3435,6.2831852) = %g\n", x);
+
+ x = drem(-10.3435,6.2831852);
+
+ printf("drem(-10.3435,6.2831852) = %g\n", x);
+
+ x = drem(-10.3435,-6.2831852);
+
+ printf("drem(-10.3435,-6.2831852) = %g\n", x);
+
+ x = drem(10.3435,-6.2831852);
+
+ printf("drem(10.3435,-6.2831852) = %g\n", x);
+
+
+ printf("x%8.6gx\n", .5);
+ printf("x%-8.6gx\n", .5);
+ printf("x%6.6gx\n", .5);
+
+ {
+ double x = atof ("-1e-17-");
+ printf ("%g %c= %g %s!\n",
+ x,
+ x == -1e-17 ? '=' : '!',
+ -1e-17,
+ x == -1e-17 ? "Worked" : "Failed");
+ }
+
+ print_trig_stuff ();
+
+ return 0;
+}
+
+
+#define PI 3.14159265358979323846264338327
+
+const double RAD[5] = { 0, PI/2, PI, (3*PI)/2, 2*PI };
+const int DEG[5] = { 0, 90, 180, 360 };
+
+#define PRINT_IT_1_ARG(_func, _arg, _value) \
+ (_value) = (_func)((_arg)); \
+ if (errno) { \
+ errno = 0; \
+ printf("%s = ERROR %s\n", #_func, strerror(errno)); \
+ } else \
+ printf("%s(%g) = %g\n", #_func, _arg, (_value)); \
+
+#define PRINT_IT_2_ARG(_func, _arg1, _arg2, _value) \
+ (_value) = (_func)((_arg1),(_arg2)); \
+ if (errno) { \
+ errno = 0; \
+ printf("%s = ERROR %s\n", #_func, strerror(errno)); \
+ } else \
+ printf("%s(%g, %g) = %g\n", #_func, _arg1, _arg2, (_value)); \
+
+void
+DEFUN_VOID (print_trig_stuff)
+{
+ double value, arg1, arg2;
+ int i;
+
+ puts ("\n\nMath Test");
+
+ errno = 0; /* automatically reset on error condition */
+ for (i=0; i<4; i++)
+ {
+ PRINT_IT_1_ARG (sin, RAD[i], value);
+ PRINT_IT_1_ARG (cos, RAD[i], value);
+ PRINT_IT_1_ARG (tan, RAD[i], value);
+ PRINT_IT_1_ARG (asin, RAD[i], value);
+ PRINT_IT_1_ARG (acos, RAD[i], value);
+ PRINT_IT_1_ARG (atan, RAD[i], value);
+ PRINT_IT_2_ARG (atan2, RAD[i], -RAD[i % 4], value);
+ }
+
+ arg1 = 16;
+ arg2 = 3;
+ PRINT_IT_1_ARG (exp, arg1, value);
+ PRINT_IT_1_ARG (log, arg1, value);
+ PRINT_IT_1_ARG (log10, arg1, value);
+ PRINT_IT_2_ARG (pow, arg1, arg2, value);
+ PRINT_IT_1_ARG (sqrt, arg1, value);
+ PRINT_IT_1_ARG (cbrt, arg1, value);
+ PRINT_IT_2_ARG (hypot, arg1, arg2, value);
+ PRINT_IT_1_ARG (expm1, arg1, value);
+ PRINT_IT_1_ARG (log1p, arg1, value);
+ PRINT_IT_1_ARG (sinh, arg1, value);
+ PRINT_IT_1_ARG (cosh, arg1, value);
+ PRINT_IT_1_ARG (tanh, arg1, value);
+ PRINT_IT_1_ARG (asinh, arg1, value);
+ PRINT_IT_1_ARG (acosh, arg1, value);
+ PRINT_IT_1_ARG (atanh, arg1, value);
+}