about summary refs log tree commit diff
path: root/REORG.TODO/sysdeps/ieee754/dbl-64/MathLib.h
diff options
context:
space:
mode:
Diffstat (limited to 'REORG.TODO/sysdeps/ieee754/dbl-64/MathLib.h')
-rw-r--r--REORG.TODO/sysdeps/ieee754/dbl-64/MathLib.h100
1 files changed, 100 insertions, 0 deletions
diff --git a/REORG.TODO/sysdeps/ieee754/dbl-64/MathLib.h b/REORG.TODO/sysdeps/ieee754/dbl-64/MathLib.h
new file mode 100644
index 0000000000..ffda70a413
--- /dev/null
+++ b/REORG.TODO/sysdeps/ieee754/dbl-64/MathLib.h
@@ -0,0 +1,100 @@
+/*
+ * IBM Accurate Mathematical Library
+ * Written by International Business Machines Corp.
+ * Copyright (C) 2001-2017 Free Software Foundation, Inc.
+ *
+ * This program is free software; you can redistribute it and/or modify
+ * it under the terms of the GNU Lesser General Public License as published by
+ * the Free Software Foundation; either version 2.1 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 Lesser General Public License for more details.
+ *
+ * You should have received a copy of the GNU Lesser General Public License
+ * along with this program; if not, see <http://www.gnu.org/licenses/>.
+ */
+
+/********************************************************************/
+/* Ultimate math functions. Each function computes the exact        */
+/* theoretical value of its argument rounded to nearest or even.    */
+/*                                                                  */
+/* Assumption: Machine arithmetic operations are performed in       */
+/* round nearest mode of IEEE 754 standard.                         */
+/********************************************************************/
+
+#ifndef UMATH_LIB
+#define UMATH_LIB
+/********************************************************************/
+/* Function changes the precision mode to IEEE 754 double precision */
+/* and the rounding mode to nearest or even.                        */
+/* It returns the original status of these modes.                   */
+/* See further explanations of usage in DPChange.h                  */
+/********************************************************************/
+unsigned short Init_Lib (void);
+
+/********************************************************************/
+/* Function that changes the precision and rounding modes to the    */
+/* specified by the argument received. See further explanations in  */
+/* DPChange.h                                                       */
+/********************************************************************/
+void Exit_Lib (unsigned short);
+
+
+/* The  asin() function calculates the arc sine of its argument.    */
+/* The  function returns the arc sine in radians                    */
+/* (between -PI/2 and PI/2).                                        */
+/* If the argument is greater than 1 or less than -1 it returns     */
+/* a NaN.                                                           */
+double uasin (double);
+
+
+/* The  acos() function calculates the arc cosine of its argument.  */
+/* The  function returns the arc cosine in radians                  */
+/* (between -PI/2 and PI/2).                                        */
+/* If the argument is greater than 1 or less than -1 it returns     */
+/* a NaN.                                                           */
+double uacos (double);
+
+/* The  atan() function calculates the arctanget of its argument.   */
+/* The  function returns the arc tangent in radians                 */
+/* (between -PI/2 and PI/2).                                        */
+double uatan (double);
+
+
+/* The uatan2() function calculates the arc tangent of the two arguments x   */
+/* and y (x is the right argument and y is the left one).The signs of both   */
+/* arguments are used to determine the quadrant of the result.               */
+/* The function returns the result in radians, which is between -PI and PI   */
+double uatan2 (double, double);
+
+/* Compute log(x). The base of log is e (natural logarithm)         */
+double ulog (double);
+
+/* Compute e raised to the power of argument x.                     */
+double uexp (double);
+
+/* Compute sin(x). The argument x is assumed to be given in radians.*/
+double usin (double);
+
+/* Compute cos(x). The argument x is assumed to be given in radians.*/
+double ucos (double);
+
+/* Compute tan(x). The argument x is assumed to be given in radians.*/
+double utan (double);
+
+/* Compute the square root of non-negative argument x.              */
+/* If x is negative the returned value is NaN.                      */
+double usqrt (double);
+
+/* Compute x raised to the power of y, where x is the left argument */
+/* and y is the right argument. The function returns a NaN if x<0.  */
+/* If x equals zero it returns -inf                                 */
+double upow (double, double);
+
+/* Computing x mod y, where x is the left argument and y is the     */
+/* right one.                                                       */
+double uremainder (double, double);
+#endif