diff options
Diffstat (limited to 'math')
| -rw-r--r-- | math/.cvsignore | 4 | ||||
| -rw-r--r-- | math/Makefile | 39 | ||||
| -rw-r--r-- | math/finite.c | 31 | ||||
| -rw-r--r-- | math/math.h | 317 | ||||
| -rw-r--r-- | math/test-math.c | 148 |
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); +} |