diff options
Diffstat (limited to 'REORG.TODO/sysdeps/m68k/m680x0')
187 files changed, 5417 insertions, 0 deletions
diff --git a/REORG.TODO/sysdeps/m68k/m680x0/Implies b/REORG.TODO/sysdeps/m68k/m680x0/Implies new file mode 100644 index 0000000000..abf356d9e8 --- /dev/null +++ b/REORG.TODO/sysdeps/m68k/m680x0/Implies @@ -0,0 +1 @@ +ieee754/ldbl-96 diff --git a/REORG.TODO/sysdeps/m68k/m680x0/Makefile b/REORG.TODO/sysdeps/m68k/m680x0/Makefile new file mode 100644 index 0000000000..582fa6f900 --- /dev/null +++ b/REORG.TODO/sysdeps/m68k/m680x0/Makefile @@ -0,0 +1,2 @@ +# The 68k `long double' is a distinct type we support. +long-double-fcts = yes diff --git a/REORG.TODO/sysdeps/m68k/m680x0/add_n.S b/REORG.TODO/sysdeps/m68k/m680x0/add_n.S new file mode 100644 index 0000000000..522ffcacf4 --- /dev/null +++ b/REORG.TODO/sysdeps/m68k/m680x0/add_n.S @@ -0,0 +1,82 @@ +/* mc68020 __mpn_add_n -- Add two limb vectors of the same length > 0 and store + sum in a third limb vector. + +Copyright (C) 1992-2017 Free Software Foundation, Inc. + +This file is part of the GNU MP Library. + +The GNU MP Library 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. + +The GNU MP 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 Lesser General Public +License for more details. + +You should have received a copy of the GNU Lesser General Public License +along with the GNU MP Library. If not, see <http://www.gnu.org/licenses/>. */ + +/* + INPUT PARAMETERS + res_ptr (sp + 4) + s1_ptr (sp + 8) + s2_ptr (sp + 16) + size (sp + 12) +*/ + +#include "sysdep.h" +#include "asm-syntax.h" + + TEXT +ENTRY(__mpn_add_n) +/* Save used registers on the stack. */ + movel R(d2),MEM_PREDEC(sp) + cfi_adjust_cfa_offset (4) + movel R(a2),MEM_PREDEC(sp) + cfi_adjust_cfa_offset (4) + cfi_rel_offset (R(d2), 4) + cfi_rel_offset (R(a2), 0) + +/* Copy the arguments to registers. Better use movem? */ + movel MEM_DISP(sp,12),R(a2) + movel MEM_DISP(sp,16),R(a0) + movel MEM_DISP(sp,20),R(a1) + movel MEM_DISP(sp,24),R(d2) + + eorw #1,R(d2) + lsrl #1,R(d2) + bcc L(L1) + subql #1,R(d2) /* clears cy as side effect */ + +L(Loop:) + movel MEM_POSTINC(a0),R(d0) + movel MEM_POSTINC(a1),R(d1) + addxl R(d1),R(d0) + movel R(d0),MEM_POSTINC(a2) +L(L1:) movel MEM_POSTINC(a0),R(d0) + movel MEM_POSTINC(a1),R(d1) + addxl R(d1),R(d0) + movel R(d0),MEM_POSTINC(a2) + + dbf R(d2),L(Loop) /* loop until 16 lsb of %4 == -1 */ + subxl R(d0),R(d0) /* d0 <= -cy; save cy as 0 or -1 in d0 */ + subl #0x10000,R(d2) + bcs L(L2) + addl R(d0),R(d0) /* restore cy */ + bra L(Loop) + +L(L2:) + negl R(d0) + +/* Restore used registers from stack frame. */ + movel MEM_POSTINC(sp),R(a2) + cfi_adjust_cfa_offset (-4) + cfi_restore (R(a2)) + movel MEM_POSTINC(sp),R(d2) + cfi_adjust_cfa_offset (-4) + cfi_restore (R(d2)) + + rts +END(__mpn_add_n) diff --git a/REORG.TODO/sysdeps/m68k/m680x0/bits/flt-eval-method.h b/REORG.TODO/sysdeps/m68k/m680x0/bits/flt-eval-method.h new file mode 100644 index 0000000000..3d28adcc96 --- /dev/null +++ b/REORG.TODO/sysdeps/m68k/m680x0/bits/flt-eval-method.h @@ -0,0 +1,25 @@ +/* Define __GLIBC_FLT_EVAL_METHOD. M68K version. + Copyright (C) 2016-2017 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 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. + + 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 + Lesser General Public License for more details. + + You should have received a copy of the GNU Lesser General Public + License along with the GNU C Library; if not, see + <http://www.gnu.org/licenses/>. */ + +#ifndef _MATH_H +# error "Never use <bits/flt-eval-method.h> directly; include <math.h> instead." +#endif + +/* The m68k FPUs evaluate all values in the 96-bit floating-point + format which is also available for the user as 'long double'. */ +#define __GLIBC_FLT_EVAL_METHOD 2 diff --git a/REORG.TODO/sysdeps/m68k/m680x0/bits/fp-logb.h b/REORG.TODO/sysdeps/m68k/m680x0/bits/fp-logb.h new file mode 100644 index 0000000000..ebde56d2e2 --- /dev/null +++ b/REORG.TODO/sysdeps/m68k/m680x0/bits/fp-logb.h @@ -0,0 +1,24 @@ +/* Define __FP_LOGB0_IS_MIN and __FP_LOGBNAN_IS_MIN. M68K version. + Copyright (C) 2016-2017 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 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. + + 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 + Lesser General Public License for more details. + + You should have received a copy of the GNU Lesser General Public + License along with the GNU C Library; if not, see + <http://www.gnu.org/licenses/>. */ + +#ifndef _MATH_H +# error "Never use <bits/fp-logb.h> directly; include <math.h> instead." +#endif + +#define __FP_LOGB0_IS_MIN 1 +#define __FP_LOGBNAN_IS_MIN 0 diff --git a/REORG.TODO/sysdeps/m68k/m680x0/bits/huge_vall.h b/REORG.TODO/sysdeps/m68k/m680x0/bits/huge_vall.h new file mode 100644 index 0000000000..14a4483582 --- /dev/null +++ b/REORG.TODO/sysdeps/m68k/m680x0/bits/huge_vall.h @@ -0,0 +1,41 @@ +/* `HUGE_VALL' constant for m68k (where it is infinity). + Used by <stdlib.h> and <math.h> functions for overflow. + Copyright (C) 1992-2017 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 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. + + 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 + Lesser General Public License for more details. + + You should have received a copy of the GNU Lesser General Public + License along with the GNU C Library. If not, see + <http://www.gnu.org/licenses/>. */ + +#ifndef _MATH_H +# error "Never use <bits/huge_val.h> directly; include <math.h> instead." +#endif + +#if __GNUC_PREREQ(3,3) +# define HUGE_VALL (__builtin_huge_vall ()) +#elif __GNUC_PREREQ(2,96) +# define HUGE_VALL (__extension__ 0x1.0p32767L) +#elif defined __GNUC__ + +# define HUGE_VALL \ + (__extension__ \ + ((union { unsigned long __l[3]; long double __ld; }) \ + { __l: { 0x7fff0000UL, 0x80000000UL, 0UL } }).__ld) + +#else /* not GCC */ + +static union { unsigned char __c[12]; long double __ld; } __huge_vall = + { { 0x7f, 0xff, 0, 0, 0x80, 0, 0, 0, 0, 0, 0, 0 } }; +# define HUGE_VALL (__huge_vall.__ld) + +#endif /* GCC 2.95. */ diff --git a/REORG.TODO/sysdeps/m68k/m680x0/fpu/bits/mathinline.h b/REORG.TODO/sysdeps/m68k/m680x0/fpu/bits/mathinline.h new file mode 100644 index 0000000000..8e6bdc4e27 --- /dev/null +++ b/REORG.TODO/sysdeps/m68k/m680x0/fpu/bits/mathinline.h @@ -0,0 +1,444 @@ +/* Definitions of inline math functions implemented by the m68881/2. + Copyright (C) 1991-2017 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 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. + + 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 + Lesser General Public License for more details. + + You should have received a copy of the GNU Lesser General Public + License along with the GNU C Library. If not, see + <http://www.gnu.org/licenses/>. */ + +#ifndef _MATH_H +# error "Never use <bits/mathinline.h> directly; include <math.h> instead." +#endif + +#ifndef __extern_inline +# define __MATH_INLINE __inline +#else +# define __MATH_INLINE __extern_inline +#endif + +#ifdef __GNUC__ + +#ifdef __USE_ISOC99 +/* GCC 3.1 and up have builtins that actually can be used. */ +# if !__GNUC_PREREQ (3,1) +/* ISO C99 defines some macros to perform unordered comparisons. The + m68k FPU supports this with special opcodes and we should use them. + These must not be inline functions since we have to be able to handle + all floating-point types. */ +# undef isgreater +# undef isgreaterequal +# undef isless +# undef islessequal +# undef islessgreater +# undef isunordered +# define isgreater(x, y) \ + __extension__ \ + ({ char __result; \ + __asm__ ("fcmp%.x %2,%1; fsogt %0" \ + : "=dm" (__result) : "f" (x), "f" (y)); \ + __result != 0; }) + +# define isgreaterequal(x, y) \ + __extension__ \ + ({ char __result; \ + __asm__ ("fcmp%.x %2,%1; fsoge %0" \ + : "=dm" (__result) : "f" (x), "f" (y)); \ + __result != 0; }) + +# define isless(x, y) \ + __extension__ \ + ({ char __result; \ + __asm__ ("fcmp%.x %2,%1; fsolt %0" \ + : "=dm" (__result) : "f" (x), "f" (y)); \ + __result != 0; }) + +# define islessequal(x, y) \ + __extension__ \ + ({ char __result; \ + __asm__ ("fcmp%.x %2,%1; fsole %0" \ + : "=dm" (__result) : "f" (x), "f" (y)); \ + __result != 0; }) + +# define islessgreater(x, y) \ + __extension__ \ + ({ char __result; \ + __asm__ ("fcmp%.x %2,%1; fsogl %0" \ + : "=dm" (__result) : "f" (x), "f" (y)); \ + __result != 0; }) + +# define isunordered(x, y) \ + __extension__ \ + ({ char __result; \ + __asm__ ("fcmp%.x %2,%1; fsun %0" \ + : "=dm" (__result) : "f" (x), "f" (y)); \ + __result != 0; }) +# endif /* GCC 3.1 */ + +/* Test for negative number. Used in the signbit() macro. */ +__MATH_INLINE int +__NTH (__signbitf (float __x)) +{ + __extension__ union { float __f; int __i; } __u = { __f: __x }; + return __u.__i < 0; +} +__MATH_INLINE int +__NTH (__signbit (double __x)) +{ + __extension__ union { double __d; int __i[2]; } __u = { __d: __x }; + return __u.__i[0] < 0; +} +__MATH_INLINE int +__NTH (__signbitl (long double __x)) +{ + __extension__ union { long double __d; int __i[3]; } __u = { __d: __x }; + return __u.__i[0] < 0; +} +#endif + + +#if (!defined __NO_MATH_INLINES && defined __OPTIMIZE__) \ + || defined __LIBC_INTERNAL_MATH_INLINES + +#ifdef __LIBC_INTERNAL_MATH_INLINES +/* This is used when defining the functions themselves. Define them with + __ names, and with `static inline' instead of `extern inline' so the + bodies will always be used, never an external function call. + Note: GCC 6 objects to __attribute__ ((__leaf__)) on static functions. */ +# define __m81_u(x) __CONCAT(__,x) +# define __m81_inline static __inline +# define __m81_nth(fn) __NTHNL (fn) +#else +# define __m81_u(x) x +# define __m81_inline __MATH_INLINE +# define __m81_nth(fn) __NTH (fn) +# define __M81_MATH_INLINES 1 +#endif + +/* Define a const math function. */ +#define __m81_defun(rettype, func, args) \ + __m81_inline rettype __attribute__((__const__)) \ + __m81_nth (__m81_u(func) args) + +/* Define the three variants of a math function that has a direct + implementation in the m68k fpu. FUNC is the name for C (which will be + suffixed with f and l for the float and long double version, resp). OP + is the name of the fpu operation (without leading f). */ + +#ifdef __USE_ISOC99 +# define __inline_mathop(func, op) \ + __inline_mathop1(double, func, op) \ + __inline_mathop1(float, __CONCAT(func,f), op) \ + __inline_mathop1(long double, __CONCAT(func,l), op) +#else +# define __inline_mathop(func, op) \ + __inline_mathop1(double, func, op) +#endif + +#define __inline_mathop1(float_type,func, op) \ + __m81_defun (float_type, func, (float_type __mathop_x)) \ + { \ + float_type __result; \ + __asm __volatile__ ("f" __STRING(op) "%.x %1, %0" \ + : "=f" (__result) : "f" (__mathop_x)); \ + return __result; \ + } + +__inline_mathop(__atan, atan) +__inline_mathop(__cos, cos) +__inline_mathop(__sin, sin) +__inline_mathop(__tan, tan) +__inline_mathop(__tanh, tanh) +__inline_mathop(__fabs, abs) + +#if defined __USE_XOPEN_EXTENDED || defined __USE_ISOC99 +__inline_mathop(__rint, int) +__inline_mathop(__expm1, etoxm1) +__inline_mathop(__log1p, lognp1) +#endif + +#ifdef __USE_MISC +__inline_mathop(__significand, getman) +#endif + +#ifdef __USE_ISOC99 +__inline_mathop(__trunc, intrz) +#endif + +#if !defined __NO_MATH_INLINES && defined __OPTIMIZE__ + +__inline_mathop(atan, atan) +__inline_mathop(tanh, tanh) + +# if defined __USE_XOPEN_EXTENDED || defined __USE_ISOC99 +__inline_mathop(rint, int) +__inline_mathop(log1p, lognp1) +# endif + +# ifdef __USE_MISC +__inline_mathop(significand, getman) +# endif + +# ifdef __USE_ISOC99 +__inline_mathop(trunc, intrz) +# endif + +#endif /* !__NO_MATH_INLINES && __OPTIMIZE__ */ + +/* This macro contains the definition for the rest of the inline + functions, using FLOAT_TYPE as the domain type and S as the suffix + for the function names. */ + +#define __inline_functions(float_type, s) \ +__m81_defun (float_type, __CONCAT(__floor,s), (float_type __x)) \ +{ \ + float_type __result; \ + unsigned long int __ctrl_reg; \ + __asm __volatile__ ("fmove%.l %!, %0" : "=dm" (__ctrl_reg)); \ + /* Set rounding towards negative infinity. */ \ + __asm __volatile__ ("fmove%.l %0, %!" : /* No outputs. */ \ + : "dmi" ((__ctrl_reg & ~0x10) | 0x20)); \ + /* Convert X to an integer, using -Inf rounding. */ \ + __asm __volatile__ ("fint%.x %1, %0" : "=f" (__result) : "f" (__x)); \ + /* Restore the previous rounding mode. */ \ + __asm __volatile__ ("fmove%.l %0, %!" : /* No outputs. */ \ + : "dmi" (__ctrl_reg)); \ + return __result; \ +} \ + \ +__m81_defun (float_type, __CONCAT(__ceil,s), (float_type __x)) \ +{ \ + float_type __result; \ + unsigned long int __ctrl_reg; \ + __asm __volatile__ ("fmove%.l %!, %0" : "=dm" (__ctrl_reg)); \ + /* Set rounding towards positive infinity. */ \ + __asm __volatile__ ("fmove%.l %0, %!" : /* No outputs. */ \ + : "dmi" (__ctrl_reg | 0x30)); \ + /* Convert X to an integer, using +Inf rounding. */ \ + __asm __volatile__ ("fint%.x %1, %0" : "=f" (__result) : "f" (__x)); \ + /* Restore the previous rounding mode. */ \ + __asm __volatile__ ("fmove%.l %0, %!" : /* No outputs. */ \ + : "dmi" (__ctrl_reg)); \ + return __result; \ +} + +__inline_functions(double,) +#ifdef __USE_ISOC99 +__inline_functions(float,f) +__inline_functions(long double,l) +#endif +#undef __inline_functions + +#ifdef __USE_MISC + +# define __inline_functions(float_type, s) \ +__m81_defun (int, __CONCAT(__isinf,s), (float_type __value)) \ +{ \ + /* There is no branch-condition for infinity, \ + so we must extract and examine the condition codes manually. */ \ + unsigned long int __fpsr; \ + __asm ("ftst%.x %1\n" \ + "fmove%.l %/fpsr, %0" : "=dm" (__fpsr) : "f" (__value)); \ + return (__fpsr & (2 << 24)) ? (__fpsr & (8 << 24) ? -1 : 1) : 0; \ +} \ + \ +__m81_defun (int, __CONCAT(__finite,s), (float_type __value)) \ +{ \ + /* There is no branch-condition for infinity, so we must extract and \ + examine the condition codes manually. */ \ + unsigned long int __fpsr; \ + __asm ("ftst%.x %1\n" \ + "fmove%.l %/fpsr, %0" : "=dm" (__fpsr) : "f" (__value)); \ + return (__fpsr & (3 << 24)) == 0; \ +} \ + \ +__m81_defun (float_type, __CONCAT(__scalbn,s), \ + (float_type __x, int __n)) \ +{ \ + float_type __result; \ + __asm __volatile__ ("fscale%.l %1, %0" : "=f" (__result) \ + : "dmi" (__n), "0" (__x)); \ + return __result; \ +} + +__inline_functions(double,) +__inline_functions(float,f) +__inline_functions(long double,l) +# undef __inline_functions + +#endif /* Use misc. */ + +#if defined __USE_MISC || defined __USE_XOPEN + +# define __inline_functions(float_type, s) \ +__m81_defun (int, __CONCAT(__isnan,s), (float_type __value)) \ +{ \ + char __result; \ + __asm ("ftst%.x %1\n" \ + "fsun %0" : "=dm" (__result) : "f" (__value)); \ + return __result; \ +} + +__inline_functions(double,) +# ifdef __USE_MISC +__inline_functions(float,f) +__inline_functions(long double,l) +# endif +# undef __inline_functions + +#endif + +#ifdef __USE_ISOC99 + +# define __inline_functions(float_type, s) \ +__m81_defun (float_type, __CONCAT(__scalbln,s), \ + (float_type __x, long int __n)) \ +{ \ + return __CONCAT(__scalbn,s) (__x, __n); \ +} \ + \ +__m81_defun (float_type, __CONCAT(__nearbyint,s), (float_type __x)) \ +{ \ + float_type __result; \ + unsigned long int __ctrl_reg; \ + __asm __volatile__ ("fmove%.l %!, %0" : "=dm" (__ctrl_reg)); \ + /* Temporarily disable the inexact exception. */ \ + __asm __volatile__ ("fmove%.l %0, %!" : /* No outputs. */ \ + : "dmi" (__ctrl_reg & ~0x200)); \ + __asm __volatile__ ("fint%.x %1, %0" : "=f" (__result) : "f" (__x)); \ + __asm __volatile__ ("fmove%.l %0, %!" : /* No outputs. */ \ + : "dmi" (__ctrl_reg)); \ + return __result; \ +} \ + \ +__m81_defun (long int, __CONCAT(__lrint,s), (float_type __x)) \ +{ \ + long int __result; \ + __asm __volatile__ ("fmove%.l %1, %0" : "=dm" (__result) : "f" (__x)); \ + return __result; \ +} + +__inline_functions (double,) +__inline_functions (float,f) +__inline_functions (long double,l) +# undef __inline_functions + +#endif /* Use ISO C9x */ + +#ifdef __USE_GNU + +# define __inline_functions(float_type, s) \ +__m81_inline void \ +__m81_nth (__m81_u(__CONCAT(__sincos,s)) \ + (float_type __x, float_type *__sinx, float_type *__cosx)) \ +{ \ + __asm __volatile__ ("fsincos%.x %2,%1:%0" \ + : "=f" (*__sinx), "=f" (*__cosx) : "f" (__x)); \ +} + +__inline_functions (double,) +__inline_functions (float,f) +__inline_functions (long double,l) +# undef __inline_functions + +#endif + +#if !defined __NO_MATH_INLINES && defined __OPTIMIZE__ + +/* Define inline versions of the user visible functions. */ + +# define __inline_forward_c(rettype, name, args1, args2) \ +__MATH_INLINE rettype __attribute__((__const__)) \ +__NTH (name args1) \ +{ \ + return __CONCAT(__,name) args2; \ +} + +# define __inline_forward(rettype, name, args1, args2) \ +__MATH_INLINE rettype __NTH (name args1) \ +{ \ + return __CONCAT(__,name) args2; \ +} + +__inline_forward_c(double,floor, (double __x), (__x)) +__inline_forward_c(double,ceil, (double __x), (__x)) +# ifdef __USE_MISC +# ifndef __USE_ISOC99 /* Conflict with macro of same name. */ +__inline_forward_c(int,isinf, (double __value), (__value)) +# endif +__inline_forward_c(int,finite, (double __value), (__value)) +__inline_forward_c(double,scalbn, (double __x, int __n), (__x, __n)) +# endif +# if defined __USE_MISC || defined __USE_XOPEN +# ifndef __USE_ISOC99 /* Conflict with macro of same name. */ +__inline_forward_c(int,isnan, (double __value), (__value)) +# endif +# endif +# ifdef __USE_ISOC99 +__inline_forward_c(double,scalbln, (double __x, long int __n), (__x, __n)) +__inline_forward_c(double,nearbyint, (double __value), (__value)) +__inline_forward_c(long int,lrint, (double __value), (__value)) +# endif +# ifdef __USE_GNU +__inline_forward(void,sincos, (double __x, double *__sinx, double *__cosx), + (__x, __sinx, __cosx)) +# endif + +# ifdef __USE_ISOC99 + +__inline_forward_c(float,floorf, (float __x), (__x)) +__inline_forward_c(float,ceilf, (float __x), (__x)) +# ifdef __USE_MISC +__inline_forward_c(int,isinff, (float __value), (__value)) +__inline_forward_c(int,finitef, (float __value), (__value)) +__inline_forward_c(float,scalbnf, (float __x, int __n), (__x, __n)) +__inline_forward_c(int,isnanf, (float __value), (__value)) +# endif +# ifdef __USE_ISOC99 +__inline_forward_c(float,scalblnf, (float __x, long int __n), (__x, __n)) +__inline_forward_c(float,nearbyintf, (float __value), (__value)) +__inline_forward_c(long int,lrintf, (float __value), (__value)) +# endif +# ifdef __USE_GNU +__inline_forward(void,sincosf, (float __x, float *__sinx, float *__cosx), + (__x, __sinx, __cosx)) +# endif + +__inline_forward_c(long double,floorl, (long double __x), (__x)) +__inline_forward_c(long double,ceill, (long double __x), (__x)) +# ifdef __USE_MISC +__inline_forward_c(int,isinfl, (long double __value), (__value)) +__inline_forward_c(int,finitel, (long double __value), (__value)) +__inline_forward_c(long double,scalbnl, (long double __x, int __n), (__x, __n)) +__inline_forward_c(int,isnanl, (long double __value), (__value)) +# endif +# ifdef __USE_ISOC99 +__inline_forward_c(long double,scalblnl, (long double __x, long int __n), + (__x, __n)) +__inline_forward_c(long double,nearbyintl, (long double __value), (__value)) +__inline_forward_c(long int,lrintl, (long double __value), (__value)) +# endif +# ifdef __USE_GNU +__inline_forward(void,sincosl, + (long double __x, long double *__sinx, long double *__cosx), + (__x, __sinx, __cosx)) +# endif + +#endif /* Use misc or ISO C99 */ + +#undef __inline_forward +#undef __inline_forward_c + +#endif /* !__NO_MATH_INLINES && __OPTIMIZE__ */ + +#endif +#endif /* GCC. */ diff --git a/REORG.TODO/sysdeps/m68k/m680x0/fpu/branred.c b/REORG.TODO/sysdeps/m68k/m680x0/fpu/branred.c new file mode 100644 index 0000000000..1cc8931700 --- /dev/null +++ b/REORG.TODO/sysdeps/m68k/m680x0/fpu/branred.c @@ -0,0 +1 @@ +/* Not needed. */ diff --git a/REORG.TODO/sysdeps/m68k/m680x0/fpu/doasin.c b/REORG.TODO/sysdeps/m68k/m680x0/fpu/doasin.c new file mode 100644 index 0000000000..1cc8931700 --- /dev/null +++ b/REORG.TODO/sysdeps/m68k/m680x0/fpu/doasin.c @@ -0,0 +1 @@ +/* Not needed. */ diff --git a/REORG.TODO/sysdeps/m68k/m680x0/fpu/dosincos.c b/REORG.TODO/sysdeps/m68k/m680x0/fpu/dosincos.c new file mode 100644 index 0000000000..1cc8931700 --- /dev/null +++ b/REORG.TODO/sysdeps/m68k/m680x0/fpu/dosincos.c @@ -0,0 +1 @@ +/* Not needed. */ diff --git a/REORG.TODO/sysdeps/m68k/m680x0/fpu/e_acos.c b/REORG.TODO/sysdeps/m68k/m680x0/fpu/e_acos.c new file mode 100644 index 0000000000..26eaa67242 --- /dev/null +++ b/REORG.TODO/sysdeps/m68k/m680x0/fpu/e_acos.c @@ -0,0 +1,37 @@ +/* Copyright (C) 1996-2017 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 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. + + 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 + Lesser General Public License for more details. + + You should have received a copy of the GNU Lesser General Public + License along with the GNU C Library. If not, see + <http://www.gnu.org/licenses/>. */ + +#include <math.h> +#include <math_private.h> +#include "mathimpl.h" + +#ifndef FUNC +# define FUNC __ieee754_acos +# define FUNC_FINITE __acos_finite +#endif +#ifndef float_type +# define float_type double +#endif + +float_type +FUNC (float_type x) +{ + return __m81_u(FUNC)(x); +} +#ifdef FUNC_FINITE +strong_alias (FUNC, FUNC_FINITE) +#endif diff --git a/REORG.TODO/sysdeps/m68k/m680x0/fpu/e_acosf.c b/REORG.TODO/sysdeps/m68k/m680x0/fpu/e_acosf.c new file mode 100644 index 0000000000..300fff601e --- /dev/null +++ b/REORG.TODO/sysdeps/m68k/m680x0/fpu/e_acosf.c @@ -0,0 +1,6 @@ +#ifndef FUNC +# define FUNC __ieee754_acosf +# define FUNC_FINITE __acosf_finite +#endif +#define float_type float +#include <e_acos.c> diff --git a/REORG.TODO/sysdeps/m68k/m680x0/fpu/e_acosl.c b/REORG.TODO/sysdeps/m68k/m680x0/fpu/e_acosl.c new file mode 100644 index 0000000000..35767f4e23 --- /dev/null +++ b/REORG.TODO/sysdeps/m68k/m680x0/fpu/e_acosl.c @@ -0,0 +1,6 @@ +#ifndef FUNC +# define FUNC __ieee754_acosl +# define FUNC_FINITE __acosl_finite +#endif +#define float_type long double +#include <e_acos.c> diff --git a/REORG.TODO/sysdeps/m68k/m680x0/fpu/e_asin.c b/REORG.TODO/sysdeps/m68k/m680x0/fpu/e_asin.c new file mode 100644 index 0000000000..271a1f23d2 --- /dev/null +++ b/REORG.TODO/sysdeps/m68k/m680x0/fpu/e_asin.c @@ -0,0 +1,3 @@ +#define FUNC __ieee754_asin +#define FUNC_FINITE __asin_finite +#include <e_acos.c> diff --git a/REORG.TODO/sysdeps/m68k/m680x0/fpu/e_asinf.c b/REORG.TODO/sysdeps/m68k/m680x0/fpu/e_asinf.c new file mode 100644 index 0000000000..ece691c6fb --- /dev/null +++ b/REORG.TODO/sysdeps/m68k/m680x0/fpu/e_asinf.c @@ -0,0 +1,3 @@ +#define FUNC __ieee754_asinf +#define FUNC_FINITE __asinf_finite +#include <e_acosf.c> diff --git a/REORG.TODO/sysdeps/m68k/m680x0/fpu/e_asinl.c b/REORG.TODO/sysdeps/m68k/m680x0/fpu/e_asinl.c new file mode 100644 index 0000000000..0c136f30eb --- /dev/null +++ b/REORG.TODO/sysdeps/m68k/m680x0/fpu/e_asinl.c @@ -0,0 +1,3 @@ +#define FUNC __ieee754_asinl +#define FUNC_FINITE __asinl_finite +#include <e_acosl.c> diff --git a/REORG.TODO/sysdeps/m68k/m680x0/fpu/e_atan2.c b/REORG.TODO/sysdeps/m68k/m680x0/fpu/e_atan2.c new file mode 100644 index 0000000000..dd7392f780 --- /dev/null +++ b/REORG.TODO/sysdeps/m68k/m680x0/fpu/e_atan2.c @@ -0,0 +1,104 @@ +/* Copyright (C) 1997-2017 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 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. + + 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 + Lesser General Public License for more details. + + You should have received a copy of the GNU Lesser General Public + License along with the GNU C Library. If not, see + <http://www.gnu.org/licenses/>. */ + +#include <math.h> +#include <math_private.h> +#include "mathimpl.h" + +#ifndef SUFF +#define SUFF +#endif +#ifndef float_type +#define float_type double +#endif + +#define CONCATX(a,b) __CONCAT(a,b) +#define s(name) CONCATX(name,SUFF) +#define m81(func) __m81_u(s(func)) + +float_type +s(__ieee754_atan2) (float_type y, float_type x) +{ + float_type pi, pi_2, z; + unsigned long y_cond, x_cond; + + __asm ("fmovecr%.x %#0, %0" : "=f" (pi)); + __asm ("fscale%.w %#-1, %0" : "=f" (pi_2) : "0" (pi)); + y_cond = __m81_test (y); + x_cond = __m81_test (x); + + if ((x_cond | y_cond) & __M81_COND_NAN) + z = x + y; + else if (y_cond & __M81_COND_ZERO) + { + if (x_cond & __M81_COND_NEG) + z = y_cond & __M81_COND_NEG ? -pi : pi; + else + z = y; + } + else if (x_cond & __M81_COND_INF) + { + if (y_cond & __M81_COND_INF) + { + float_type pi_4; + __asm ("fscale%.w %#-2, %0" : "=f" (pi_4) : "0" (pi)); + z = x_cond & __M81_COND_NEG ? 3 * pi_4 : pi_4; + } + else + z = x_cond & __M81_COND_NEG ? pi : 0; + if (y_cond & __M81_COND_NEG) + z = -z; + } + else if (y_cond & __M81_COND_INF) + z = y_cond & __M81_COND_NEG ? -pi_2 : pi_2; + else if (x_cond & __M81_COND_NEG) + { + if (y_cond & __M81_COND_NEG) + { + if (-x > -y) + z = -pi + m81(__atan) (y / x); + else + z = -pi_2 - m81(__atan) (x / y); + } + else + { + if (-x > y) + z = pi + m81(__atan) (y / x); + else + z = pi_2 - m81(__atan) (x / y); + } + } + else + { + if (y_cond & __M81_COND_NEG) + { + if (x > -y) + z = m81(__atan) (y / x); + else + z = -pi_2 - m81(__atan) (x / y); + } + else + { + if (x > y) + z = m81(__atan) (y / x); + else + z = pi_2 - m81(__atan) (x / y); + } + } + return z; +} +strong_alias (s(__ieee754_atan2), CONCATX (s (__atan2), _finite)) diff --git a/REORG.TODO/sysdeps/m68k/m680x0/fpu/e_atan2f.c b/REORG.TODO/sysdeps/m68k/m680x0/fpu/e_atan2f.c new file mode 100644 index 0000000000..a0c750a759 --- /dev/null +++ b/REORG.TODO/sysdeps/m68k/m680x0/fpu/e_atan2f.c @@ -0,0 +1,3 @@ +#define SUFF f +#define float_type float +#include <e_atan2.c> diff --git a/REORG.TODO/sysdeps/m68k/m680x0/fpu/e_atan2l.c b/REORG.TODO/sysdeps/m68k/m680x0/fpu/e_atan2l.c new file mode 100644 index 0000000000..426ca94baa --- /dev/null +++ b/REORG.TODO/sysdeps/m68k/m680x0/fpu/e_atan2l.c @@ -0,0 +1,3 @@ +#define SUFF l +#define float_type long double +#include <e_atan2.c> diff --git a/REORG.TODO/sysdeps/m68k/m680x0/fpu/e_atanh.c b/REORG.TODO/sysdeps/m68k/m680x0/fpu/e_atanh.c new file mode 100644 index 0000000000..f5f7e7b520 --- /dev/null +++ b/REORG.TODO/sysdeps/m68k/m680x0/fpu/e_atanh.c @@ -0,0 +1,3 @@ +#define FUNC __ieee754_atanh +#define FUNC_FINITE __atanh_finite +#include <e_acos.c> diff --git a/REORG.TODO/sysdeps/m68k/m680x0/fpu/e_atanhf.c b/REORG.TODO/sysdeps/m68k/m680x0/fpu/e_atanhf.c new file mode 100644 index 0000000000..cb464cbb7d --- /dev/null +++ b/REORG.TODO/sysdeps/m68k/m680x0/fpu/e_atanhf.c @@ -0,0 +1,3 @@ +#define FUNC __ieee754_atanhf +#define FUNC_FINITE __atanhf_finite +#include <e_acosf.c> diff --git a/REORG.TODO/sysdeps/m68k/m680x0/fpu/e_atanhl.c b/REORG.TODO/sysdeps/m68k/m680x0/fpu/e_atanhl.c new file mode 100644 index 0000000000..8f276f5fad --- /dev/null +++ b/REORG.TODO/sysdeps/m68k/m680x0/fpu/e_atanhl.c @@ -0,0 +1,3 @@ +#define FUNC __ieee754_atanhl +#define FUNC_FINITE __atanhl_finite +#include <e_acosl.c> diff --git a/REORG.TODO/sysdeps/m68k/m680x0/fpu/e_cosh.c b/REORG.TODO/sysdeps/m68k/m680x0/fpu/e_cosh.c new file mode 100644 index 0000000000..aaba11be3a --- /dev/null +++ b/REORG.TODO/sysdeps/m68k/m680x0/fpu/e_cosh.c @@ -0,0 +1,3 @@ +#define FUNC __ieee754_cosh +#define FUNC_FINITE __cosh_finite +#include <e_acos.c> diff --git a/REORG.TODO/sysdeps/m68k/m680x0/fpu/e_coshf.c b/REORG.TODO/sysdeps/m68k/m680x0/fpu/e_coshf.c new file mode 100644 index 0000000000..0884741c5c --- /dev/null +++ b/REORG.TODO/sysdeps/m68k/m680x0/fpu/e_coshf.c @@ -0,0 +1,3 @@ +#define FUNC __ieee754_coshf +#define FUNC_FINITE __coshf_finite +#include <e_acosf.c> diff --git a/REORG.TODO/sysdeps/m68k/m680x0/fpu/e_coshl.c b/REORG.TODO/sysdeps/m68k/m680x0/fpu/e_coshl.c new file mode 100644 index 0000000000..7fc9b7b91c --- /dev/null +++ b/REORG.TODO/sysdeps/m68k/m680x0/fpu/e_coshl.c @@ -0,0 +1,3 @@ +#define FUNC __ieee754_coshl +#define FUNC_FINITE __coshl_finite +#include <e_acosl.c> diff --git a/REORG.TODO/sysdeps/m68k/m680x0/fpu/e_exp.c b/REORG.TODO/sysdeps/m68k/m680x0/fpu/e_exp.c new file mode 100644 index 0000000000..9f228034a1 --- /dev/null +++ b/REORG.TODO/sysdeps/m68k/m680x0/fpu/e_exp.c @@ -0,0 +1,3 @@ +#define FUNC __ieee754_exp +#define FUNC_FINITE __exp_finite +#include <e_acos.c> diff --git a/REORG.TODO/sysdeps/m68k/m680x0/fpu/e_exp10.c b/REORG.TODO/sysdeps/m68k/m680x0/fpu/e_exp10.c new file mode 100644 index 0000000000..e66a80599c --- /dev/null +++ b/REORG.TODO/sysdeps/m68k/m680x0/fpu/e_exp10.c @@ -0,0 +1,3 @@ +#define FUNC __ieee754_exp10 +#define FUNC_FINITE __exp10_finite +#include <e_acos.c> diff --git a/REORG.TODO/sysdeps/m68k/m680x0/fpu/e_exp10f.c b/REORG.TODO/sysdeps/m68k/m680x0/fpu/e_exp10f.c new file mode 100644 index 0000000000..7a3bdd4186 --- /dev/null +++ b/REORG.TODO/sysdeps/m68k/m680x0/fpu/e_exp10f.c @@ -0,0 +1,3 @@ +#define FUNC __ieee754_exp10f +#define FUNC_FINITE __exp10f_finite +#include <e_acosf.c> diff --git a/REORG.TODO/sysdeps/m68k/m680x0/fpu/e_exp10l.c b/REORG.TODO/sysdeps/m68k/m680x0/fpu/e_exp10l.c new file mode 100644 index 0000000000..e2a6255a76 --- /dev/null +++ b/REORG.TODO/sysdeps/m68k/m680x0/fpu/e_exp10l.c @@ -0,0 +1,3 @@ +#define FUNC __ieee754_exp10l +#define FUNC_FINITE __exp10l_finite +#include <e_acosl.c> diff --git a/REORG.TODO/sysdeps/m68k/m680x0/fpu/e_exp2.c b/REORG.TODO/sysdeps/m68k/m680x0/fpu/e_exp2.c new file mode 100644 index 0000000000..35e29ae653 --- /dev/null +++ b/REORG.TODO/sysdeps/m68k/m680x0/fpu/e_exp2.c @@ -0,0 +1,3 @@ +#define FUNC __ieee754_exp2 +#define FUNC_FINITE __exp2_finite +#include <e_acos.c> diff --git a/REORG.TODO/sysdeps/m68k/m680x0/fpu/e_exp2f.c b/REORG.TODO/sysdeps/m68k/m680x0/fpu/e_exp2f.c new file mode 100644 index 0000000000..e57ddb5365 --- /dev/null +++ b/REORG.TODO/sysdeps/m68k/m680x0/fpu/e_exp2f.c @@ -0,0 +1,3 @@ +#define FUNC __ieee754_exp2f +#define FUNC_FINITE __exp2f_finite +#include <e_acosf.c> diff --git a/REORG.TODO/sysdeps/m68k/m680x0/fpu/e_exp2l.c b/REORG.TODO/sysdeps/m68k/m680x0/fpu/e_exp2l.c new file mode 100644 index 0000000000..c5854ae427 --- /dev/null +++ b/REORG.TODO/sysdeps/m68k/m680x0/fpu/e_exp2l.c @@ -0,0 +1,3 @@ +#define FUNC __ieee754_exp2l +#define FUNC_FINITE __exp2l_finite +#include <e_acosl.c> diff --git a/REORG.TODO/sysdeps/m68k/m680x0/fpu/e_expf.c b/REORG.TODO/sysdeps/m68k/m680x0/fpu/e_expf.c new file mode 100644 index 0000000000..af1e019504 --- /dev/null +++ b/REORG.TODO/sysdeps/m68k/m680x0/fpu/e_expf.c @@ -0,0 +1,3 @@ +#define FUNC __ieee754_expf +#define FUNC_FINITE __expf_finite +#include <e_acosf.c> diff --git a/REORG.TODO/sysdeps/m68k/m680x0/fpu/e_expl.c b/REORG.TODO/sysdeps/m68k/m680x0/fpu/e_expl.c new file mode 100644 index 0000000000..543eb0edf5 --- /dev/null +++ b/REORG.TODO/sysdeps/m68k/m680x0/fpu/e_expl.c @@ -0,0 +1,3 @@ +#define FUNC __ieee754_expl +#define FUNC_FINITE __expl_finite +#include <e_acosl.c> diff --git a/REORG.TODO/sysdeps/m68k/m680x0/fpu/e_fmod.c b/REORG.TODO/sysdeps/m68k/m680x0/fpu/e_fmod.c new file mode 100644 index 0000000000..cd8e37a55a --- /dev/null +++ b/REORG.TODO/sysdeps/m68k/m680x0/fpu/e_fmod.c @@ -0,0 +1,37 @@ +/* Copyright (C) 1996-2017 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 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. + + 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 + Lesser General Public License for more details. + + You should have received a copy of the GNU Lesser General Public + License along with the GNU C Library. If not, see + <http://www.gnu.org/licenses/>. */ + +#include <math.h> +#include <math_private.h> +#include "mathimpl.h" + +#ifndef FUNC +# define FUNC __ieee754_fmod +# define FUNC_FINITE __fmod_finite +#endif +#ifndef float_type +# define float_type double +#endif + +float_type +FUNC (float_type x, float_type y) +{ + return __m81_u(FUNC)(x, y); +} +#ifdef FUNC_FINITE +strong_alias (FUNC, FUNC_FINITE) +#endif diff --git a/REORG.TODO/sysdeps/m68k/m680x0/fpu/e_fmodf.c b/REORG.TODO/sysdeps/m68k/m680x0/fpu/e_fmodf.c new file mode 100644 index 0000000000..38995c3768 --- /dev/null +++ b/REORG.TODO/sysdeps/m68k/m680x0/fpu/e_fmodf.c @@ -0,0 +1,6 @@ +#ifndef FUNC +# define FUNC __ieee754_fmodf +# define FUNC_FINITE __fmodf_finite +#endif +#define float_type float +#include <e_fmod.c> diff --git a/REORG.TODO/sysdeps/m68k/m680x0/fpu/e_fmodl.c b/REORG.TODO/sysdeps/m68k/m680x0/fpu/e_fmodl.c new file mode 100644 index 0000000000..8d793b2de7 --- /dev/null +++ b/REORG.TODO/sysdeps/m68k/m680x0/fpu/e_fmodl.c @@ -0,0 +1,6 @@ +#ifndef FUNC +# define FUNC __ieee754_fmodl +# define FUNC_FINITE __fmodl_finite +#endif +#define float_type long double +#include <e_fmod.c> diff --git a/REORG.TODO/sysdeps/m68k/m680x0/fpu/e_ilogb.c b/REORG.TODO/sysdeps/m68k/m680x0/fpu/e_ilogb.c new file mode 100644 index 0000000000..0de5830372 --- /dev/null +++ b/REORG.TODO/sysdeps/m68k/m680x0/fpu/e_ilogb.c @@ -0,0 +1,47 @@ +/* Copyright (C) 1996-2017 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 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. + + 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 + Lesser General Public License for more details. + + You should have received a copy of the GNU Lesser General Public + License along with the GNU C Library. If not, see + <http://www.gnu.org/licenses/>. */ + +#include <math.h> +#include "mathimpl.h" + +#ifndef SUFF +#define SUFF +#endif +#ifndef float_type +#define float_type double +#endif + +#define CONCATX(a,b) __CONCAT(a,b) +#define s(name) CONCATX(name,SUFF) +#define m81(func) __m81_u(s(func)) + +int +s(__ieee754_ilogb) (float_type x) +{ + float_type result; + unsigned long x_cond; + + x_cond = __m81_test (x); + /* We must return consistent values for zero and NaN. */ + if (x_cond & __M81_COND_ZERO) + return FP_ILOGB0; + if (x_cond & (__M81_COND_NAN | __M81_COND_INF)) + return FP_ILOGBNAN; + + __asm ("fgetexp%.x %1, %0" : "=f" (result) : "f" (x)); + return (int) result; +} diff --git a/REORG.TODO/sysdeps/m68k/m680x0/fpu/e_ilogbf.c b/REORG.TODO/sysdeps/m68k/m680x0/fpu/e_ilogbf.c new file mode 100644 index 0000000000..34a8cd28e6 --- /dev/null +++ b/REORG.TODO/sysdeps/m68k/m680x0/fpu/e_ilogbf.c @@ -0,0 +1,3 @@ +#define SUFF f +#define float_type float +#include <e_ilogb.c> diff --git a/REORG.TODO/sysdeps/m68k/m680x0/fpu/e_ilogbl.c b/REORG.TODO/sysdeps/m68k/m680x0/fpu/e_ilogbl.c new file mode 100644 index 0000000000..b0e13af6b1 --- /dev/null +++ b/REORG.TODO/sysdeps/m68k/m680x0/fpu/e_ilogbl.c @@ -0,0 +1,3 @@ +#define SUFF l +#define float_type long double +#include <e_ilogb.c> diff --git a/REORG.TODO/sysdeps/m68k/m680x0/fpu/e_log.c b/REORG.TODO/sysdeps/m68k/m680x0/fpu/e_log.c new file mode 100644 index 0000000000..628c62b31a --- /dev/null +++ b/REORG.TODO/sysdeps/m68k/m680x0/fpu/e_log.c @@ -0,0 +1,3 @@ +#define FUNC __ieee754_log +#define FUNC_FINITE __log_finite +#include <e_acos.c> diff --git a/REORG.TODO/sysdeps/m68k/m680x0/fpu/e_log10.c b/REORG.TODO/sysdeps/m68k/m680x0/fpu/e_log10.c new file mode 100644 index 0000000000..78e0693e9d --- /dev/null +++ b/REORG.TODO/sysdeps/m68k/m680x0/fpu/e_log10.c @@ -0,0 +1,3 @@ +#define FUNC __ieee754_log10 +#define FUNC_FINITE __log10_finite +#include <e_acos.c> diff --git a/REORG.TODO/sysdeps/m68k/m680x0/fpu/e_log10f.c b/REORG.TODO/sysdeps/m68k/m680x0/fpu/e_log10f.c new file mode 100644 index 0000000000..452a75ec15 --- /dev/null +++ b/REORG.TODO/sysdeps/m68k/m680x0/fpu/e_log10f.c @@ -0,0 +1,3 @@ +#define FUNC __ieee754_log10f +#define FUNC_FINITE __log10f_finite +#include <e_acosf.c> diff --git a/REORG.TODO/sysdeps/m68k/m680x0/fpu/e_log10l.c b/REORG.TODO/sysdeps/m68k/m680x0/fpu/e_log10l.c new file mode 100644 index 0000000000..7c5dcb8e18 --- /dev/null +++ b/REORG.TODO/sysdeps/m68k/m680x0/fpu/e_log10l.c @@ -0,0 +1,3 @@ +#define FUNC __ieee754_log10l +#define FUNC_FINITE __log10l_finite +#include <e_acosl.c> diff --git a/REORG.TODO/sysdeps/m68k/m680x0/fpu/e_log2.c b/REORG.TODO/sysdeps/m68k/m680x0/fpu/e_log2.c new file mode 100644 index 0000000000..a00ddbdc96 --- /dev/null +++ b/REORG.TODO/sysdeps/m68k/m680x0/fpu/e_log2.c @@ -0,0 +1,3 @@ +#define FUNC __ieee754_log2 +#define FUNC_FINITE __log2_finite +#include <e_acos.c> diff --git a/REORG.TODO/sysdeps/m68k/m680x0/fpu/e_log2f.c b/REORG.TODO/sysdeps/m68k/m680x0/fpu/e_log2f.c new file mode 100644 index 0000000000..670b69c8a9 --- /dev/null +++ b/REORG.TODO/sysdeps/m68k/m680x0/fpu/e_log2f.c @@ -0,0 +1,3 @@ +#define FUNC __ieee754_log2f +#define FUNC_FINITE __log2f_finite +#include <e_acosf.c> diff --git a/REORG.TODO/sysdeps/m68k/m680x0/fpu/e_log2l.c b/REORG.TODO/sysdeps/m68k/m680x0/fpu/e_log2l.c new file mode 100644 index 0000000000..4a97a00235 --- /dev/null +++ b/REORG.TODO/sysdeps/m68k/m680x0/fpu/e_log2l.c @@ -0,0 +1,3 @@ +#define FUNC __ieee754_log2l +#define FUNC_FINITE __log2l_finite +#include <e_acosl.c> diff --git a/REORG.TODO/sysdeps/m68k/m680x0/fpu/e_logf.c b/REORG.TODO/sysdeps/m68k/m680x0/fpu/e_logf.c new file mode 100644 index 0000000000..1989a95431 --- /dev/null +++ b/REORG.TODO/sysdeps/m68k/m680x0/fpu/e_logf.c @@ -0,0 +1,3 @@ +#define FUNC __ieee754_logf +#define FUNC_FINITE __logf_finite +#include <e_acosf.c> diff --git a/REORG.TODO/sysdeps/m68k/m680x0/fpu/e_logl.c b/REORG.TODO/sysdeps/m68k/m680x0/fpu/e_logl.c new file mode 100644 index 0000000000..9ab842c3fc --- /dev/null +++ b/REORG.TODO/sysdeps/m68k/m680x0/fpu/e_logl.c @@ -0,0 +1,3 @@ +#define FUNC __ieee754_logl +#define FUNC_FINITE __logl_finite +#include <e_acosl.c> diff --git a/REORG.TODO/sysdeps/m68k/m680x0/fpu/e_pow.c b/REORG.TODO/sysdeps/m68k/m680x0/fpu/e_pow.c new file mode 100644 index 0000000000..f2dad685e3 --- /dev/null +++ b/REORG.TODO/sysdeps/m68k/m680x0/fpu/e_pow.c @@ -0,0 +1,126 @@ +/* Copyright (C) 1997-2017 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 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. + + 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 + Lesser General Public License for more details. + + You should have received a copy of the GNU Lesser General Public + License along with the GNU C Library. If not, see + <http://www.gnu.org/licenses/>. */ + +#include <math.h> +#include <math_private.h> +#include "mathimpl.h" + +#ifndef SUFF +#define SUFF +#endif +#ifndef float_type +#define float_type double +#endif + +#define CONCATX(a,b) __CONCAT(a,b) +#define s(name) CONCATX(name,SUFF) +#define m81(func) __m81_u(s(func)) + +float_type +s(__ieee754_pow) (float_type x, float_type y) +{ + float_type z; + float_type ax; + unsigned long x_cond, y_cond; + + y_cond = __m81_test (y); + if (y_cond & __M81_COND_ZERO) + return 1.0; + if (y_cond & __M81_COND_NAN) + return x == 1.0 ? x : x + y; + + x_cond = __m81_test (x); + if (x_cond & __M81_COND_NAN) + return x + y; + + if (y_cond & __M81_COND_INF) + { + ax = s(fabs) (x); + if (ax == 1.0) + return ax; + if (ax > 1.0) + return y_cond & __M81_COND_NEG ? 0 : y; + else + return y_cond & __M81_COND_NEG ? -y : 0; + } + + if (s(fabs) (y) == 1.0) + return y_cond & __M81_COND_NEG ? 1 / x : x; + + if (y == 2) + return x * x; + if (y == 0.5 && !(x_cond & __M81_COND_NEG)) + return m81(__ieee754_sqrt) (x); + + if (x == 10.0) + { + __asm ("ftentox%.x %1, %0" : "=f" (z) : "f" (y)); + return z; + } + if (x == 2.0) + { + __asm ("ftwotox%.x %1, %0" : "=f" (z) : "f" (y)); + return z; + } + + ax = s(fabs) (x); + if (x_cond & (__M81_COND_INF | __M81_COND_ZERO) || ax == 1.0) + { + z = ax; + if (y_cond & __M81_COND_NEG) + z = 1 / z; + if (x_cond & __M81_COND_NEG) + { + if (y != m81(__rint) (y)) + { + if (x == -1) + z = (z - z) / (z - z); + } + else + goto maybe_negate; + } + return z; + } + + if (x_cond & __M81_COND_NEG) + { + if (y == m81(__rint) (y)) + { + z = m81(__ieee754_exp) (y * m81(__ieee754_log) (-x)); + maybe_negate: + /* We always use the long double format, since y is already in + this format and rounding won't change the result. */ + { + int32_t exponent; + u_int32_t i0, i1; + GET_LDOUBLE_WORDS (exponent, i0, i1, y); + exponent = (exponent & 0x7fff) - 0x3fff; + if (exponent <= 31 + ? i0 & (1 << (31 - exponent)) + : (exponent <= 63 + && i1 & (1 << (63 - exponent)))) + z = -z; + } + } + else + z = (y - y) / (y - y); + } + else + z = m81(__ieee754_exp) (y * m81(__ieee754_log) (x)); + return z; +} +strong_alias (s(__ieee754_pow), CONCATX (s(__pow), _finite)) diff --git a/REORG.TODO/sysdeps/m68k/m680x0/fpu/e_powf.c b/REORG.TODO/sysdeps/m68k/m680x0/fpu/e_powf.c new file mode 100644 index 0000000000..379014355a --- /dev/null +++ b/REORG.TODO/sysdeps/m68k/m680x0/fpu/e_powf.c @@ -0,0 +1,3 @@ +#define SUFF f +#define float_type float +#include <e_pow.c> diff --git a/REORG.TODO/sysdeps/m68k/m680x0/fpu/e_powl.c b/REORG.TODO/sysdeps/m68k/m680x0/fpu/e_powl.c new file mode 100644 index 0000000000..f71fa34a26 --- /dev/null +++ b/REORG.TODO/sysdeps/m68k/m680x0/fpu/e_powl.c @@ -0,0 +1,3 @@ +#define SUFF l +#define float_type long double +#include <e_pow.c> diff --git a/REORG.TODO/sysdeps/m68k/m680x0/fpu/e_rem_pio2.c b/REORG.TODO/sysdeps/m68k/m680x0/fpu/e_rem_pio2.c new file mode 100644 index 0000000000..1347b0468c --- /dev/null +++ b/REORG.TODO/sysdeps/m68k/m680x0/fpu/e_rem_pio2.c @@ -0,0 +1,3 @@ +/* Empty. This file is only meant to avoid compiling the file with the + same name in the libm-ieee754 directory. The code is not used since + there is an assembler version for all users of this file. */ diff --git a/REORG.TODO/sysdeps/m68k/m680x0/fpu/e_rem_pio2f.c b/REORG.TODO/sysdeps/m68k/m680x0/fpu/e_rem_pio2f.c new file mode 100644 index 0000000000..1347b0468c --- /dev/null +++ b/REORG.TODO/sysdeps/m68k/m680x0/fpu/e_rem_pio2f.c @@ -0,0 +1,3 @@ +/* Empty. This file is only meant to avoid compiling the file with the + same name in the libm-ieee754 directory. The code is not used since + there is an assembler version for all users of this file. */ diff --git a/REORG.TODO/sysdeps/m68k/m680x0/fpu/e_rem_pio2l.c b/REORG.TODO/sysdeps/m68k/m680x0/fpu/e_rem_pio2l.c new file mode 100644 index 0000000000..1347b0468c --- /dev/null +++ b/REORG.TODO/sysdeps/m68k/m680x0/fpu/e_rem_pio2l.c @@ -0,0 +1,3 @@ +/* Empty. This file is only meant to avoid compiling the file with the + same name in the libm-ieee754 directory. The code is not used since + there is an assembler version for all users of this file. */ diff --git a/REORG.TODO/sysdeps/m68k/m680x0/fpu/e_remainder.c b/REORG.TODO/sysdeps/m68k/m680x0/fpu/e_remainder.c new file mode 100644 index 0000000000..f7732af8c5 --- /dev/null +++ b/REORG.TODO/sysdeps/m68k/m680x0/fpu/e_remainder.c @@ -0,0 +1,3 @@ +#define FUNC __ieee754_remainder +#define FUNC_FINITE __remainder_finite +#include <e_fmod.c> diff --git a/REORG.TODO/sysdeps/m68k/m680x0/fpu/e_remainderf.c b/REORG.TODO/sysdeps/m68k/m680x0/fpu/e_remainderf.c new file mode 100644 index 0000000000..94b53e7a80 --- /dev/null +++ b/REORG.TODO/sysdeps/m68k/m680x0/fpu/e_remainderf.c @@ -0,0 +1,3 @@ +#define FUNC __ieee754_remainderf +#define FUNC_FINITE __remainderf_finite +#include <e_fmodf.c> diff --git a/REORG.TODO/sysdeps/m68k/m680x0/fpu/e_remainderl.c b/REORG.TODO/sysdeps/m68k/m680x0/fpu/e_remainderl.c new file mode 100644 index 0000000000..d5b59607ab --- /dev/null +++ b/REORG.TODO/sysdeps/m68k/m680x0/fpu/e_remainderl.c @@ -0,0 +1,3 @@ +#define FUNC __ieee754_remainderl +#define FUNC_FINITE __remainderl_finite +#include <e_fmodl.c> diff --git a/REORG.TODO/sysdeps/m68k/m680x0/fpu/e_scalb.c b/REORG.TODO/sysdeps/m68k/m680x0/fpu/e_scalb.c new file mode 100644 index 0000000000..80bcade955 --- /dev/null +++ b/REORG.TODO/sysdeps/m68k/m680x0/fpu/e_scalb.c @@ -0,0 +1,60 @@ +/* Copyright (C) 1997-2017 Free Software Foundation, Inc. + This file is part of the GNU C Library. + Contributed by Andreas Schwab <schwab@issan.informatik.uni-dortmund.de>. + + The GNU C Library 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. + + 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 + Lesser General Public License for more details. + + You should have received a copy of the GNU Lesser General Public + License along with the GNU C Library. If not, see + <http://www.gnu.org/licenses/>. */ + +#include <math.h> +#include <math_private.h> +#include "mathimpl.h" + +#ifndef SUFF +#define SUFF +#endif +#ifndef float_type +#define float_type double +#endif + +#define CONCATX(a,b) __CONCAT(a,b) +#define s(name) CONCATX(name,SUFF) +#define m81(func) __m81_u(s(func)) + +float_type +s(__ieee754_scalb) (float_type x, float_type fn) +{ + float_type retval; + unsigned long x_cond = __m81_test (x); + unsigned long fn_cond = __m81_test (fn); + + if ((x_cond | fn_cond) & __M81_COND_NAN) + return x * fn; + + if (fn_cond & __M81_COND_INF) + { + if (!(fn_cond & __M81_COND_NEG)) + return x * fn; + else if (x_cond & __M81_COND_ZERO) + return x; + else + return x / -fn; + } + + if (m81(__rint) (fn) != fn) + return (x - x) / (x - x); + + __asm ("fscale%.x %1, %0" : "=f" (retval) : "f" (fn), "0" (x)); + return retval; +} +strong_alias (s(__ieee754_scalb), CONCATX (s(__scalb), _finite)) diff --git a/REORG.TODO/sysdeps/m68k/m680x0/fpu/e_scalbf.c b/REORG.TODO/sysdeps/m68k/m680x0/fpu/e_scalbf.c new file mode 100644 index 0000000000..7943571246 --- /dev/null +++ b/REORG.TODO/sysdeps/m68k/m680x0/fpu/e_scalbf.c @@ -0,0 +1,3 @@ +#define SUFF f +#define float_type float +#include <e_scalb.c> diff --git a/REORG.TODO/sysdeps/m68k/m680x0/fpu/e_scalbl.c b/REORG.TODO/sysdeps/m68k/m680x0/fpu/e_scalbl.c new file mode 100644 index 0000000000..35fb2dc0ed --- /dev/null +++ b/REORG.TODO/sysdeps/m68k/m680x0/fpu/e_scalbl.c @@ -0,0 +1,3 @@ +#define SUFF l +#define float_type long double +#include <e_scalb.c> diff --git a/REORG.TODO/sysdeps/m68k/m680x0/fpu/e_sinh.c b/REORG.TODO/sysdeps/m68k/m680x0/fpu/e_sinh.c new file mode 100644 index 0000000000..1e1b1c1180 --- /dev/null +++ b/REORG.TODO/sysdeps/m68k/m680x0/fpu/e_sinh.c @@ -0,0 +1,3 @@ +#define FUNC __ieee754_sinh +#define FUNC_FINITE __sinh_finite +#include <e_acos.c> diff --git a/REORG.TODO/sysdeps/m68k/m680x0/fpu/e_sinhf.c b/REORG.TODO/sysdeps/m68k/m680x0/fpu/e_sinhf.c new file mode 100644 index 0000000000..13c79f9d96 --- /dev/null +++ b/REORG.TODO/sysdeps/m68k/m680x0/fpu/e_sinhf.c @@ -0,0 +1,3 @@ +#define FUNC __ieee754_sinhf +#define FUNC_FINITE __sinhf_finite +#include <e_acosf.c> diff --git a/REORG.TODO/sysdeps/m68k/m680x0/fpu/e_sinhl.c b/REORG.TODO/sysdeps/m68k/m680x0/fpu/e_sinhl.c new file mode 100644 index 0000000000..47f0f24401 --- /dev/null +++ b/REORG.TODO/sysdeps/m68k/m680x0/fpu/e_sinhl.c @@ -0,0 +1,3 @@ +#define FUNC __ieee754_sinhl +#define FUNC_FINITE __sinhl_finite +#include <e_acosl.c> diff --git a/REORG.TODO/sysdeps/m68k/m680x0/fpu/e_sqrt.c b/REORG.TODO/sysdeps/m68k/m680x0/fpu/e_sqrt.c new file mode 100644 index 0000000000..9a250cbc25 --- /dev/null +++ b/REORG.TODO/sysdeps/m68k/m680x0/fpu/e_sqrt.c @@ -0,0 +1,3 @@ +#define FUNC __ieee754_sqrt +#define FUNC_FINITE __sqrt_finite +#include <e_acos.c> diff --git a/REORG.TODO/sysdeps/m68k/m680x0/fpu/e_sqrtf.c b/REORG.TODO/sysdeps/m68k/m680x0/fpu/e_sqrtf.c new file mode 100644 index 0000000000..372c3eb2f4 --- /dev/null +++ b/REORG.TODO/sysdeps/m68k/m680x0/fpu/e_sqrtf.c @@ -0,0 +1,3 @@ +#define FUNC __ieee754_sqrtf +#define FUNC_FINITE __sqrtf_finite +#include <e_acosf.c> diff --git a/REORG.TODO/sysdeps/m68k/m680x0/fpu/e_sqrtl.c b/REORG.TODO/sysdeps/m68k/m680x0/fpu/e_sqrtl.c new file mode 100644 index 0000000000..df5f44e6a7 --- /dev/null +++ b/REORG.TODO/sysdeps/m68k/m680x0/fpu/e_sqrtl.c @@ -0,0 +1,3 @@ +#define FUNC __ieee754_sqrtl +#define FUNC_FINITE __sqrtl_finite +#include <e_acosl.c> diff --git a/REORG.TODO/sysdeps/m68k/m680x0/fpu/fraiseexcpt.c b/REORG.TODO/sysdeps/m68k/m680x0/fpu/fraiseexcpt.c new file mode 100644 index 0000000000..08ba0a7ef5 --- /dev/null +++ b/REORG.TODO/sysdeps/m68k/m680x0/fpu/fraiseexcpt.c @@ -0,0 +1,83 @@ +/* Raise given exceptions. + Copyright (C) 1997-2017 Free Software Foundation, Inc. + This file is part of the GNU C Library. + Contributed by Andreas Schwab <schwab@issan.informatik.uni-dortmund.de> + + The GNU C Library 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. + + 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 + Lesser General Public License for more details. + + You should have received a copy of the GNU Lesser General Public + License along with the GNU C Library. If not, see + <http://www.gnu.org/licenses/>. */ + +#include <fenv.h> +#include <float.h> +#include <math.h> + +int +__feraiseexcept (int excepts) +{ + /* Raise exceptions represented by EXCEPTS. But we must raise only one + signal at a time. It is important that if the overflow/underflow + exception and the divide by zero exception are given at the same + time, the overflow/underflow exception follows the divide by zero + exception. */ + + /* First: invalid exception. */ + if (excepts & FE_INVALID) + { + /* One example of an invalid operation is 0 * Infinity. */ + double d = HUGE_VAL; + __asm__ __volatile__ ("fmul%.s %#0r0,%0; fnop" : "=f" (d) : "0" (d)); + } + + /* Next: division by zero. */ + if (excepts & FE_DIVBYZERO) + { + double d = 1.0; + __asm__ __volatile__ ("fdiv%.s %#0r0,%0; fnop" : "=f" (d) : "0" (d)); + } + + /* Next: overflow. */ + if (excepts & FE_OVERFLOW) + { + long double d = LDBL_MAX; + + __asm__ __volatile__ ("fmul%.x %0,%0; fnop" : "=f" (d) : "0" (d)); + } + + /* Next: underflow. */ + if (excepts & FE_UNDERFLOW) + { + long double d = -LDBL_MAX; + + __asm__ __volatile__ ("fetox%.x %0; fnop" : "=f" (d) : "0" (d)); + } + + /* Last: inexact. */ + if (excepts & FE_INEXACT) + { + long double d = 1.0; + __asm__ __volatile__ ("fdiv%.s %#0r3,%0; fnop" : "=f" (d) : "0" (d)); + } + + /* Success. */ + return 0; +} + +#include <shlib-compat.h> +#if SHLIB_COMPAT (libm, GLIBC_2_1, GLIBC_2_2) +strong_alias (__feraiseexcept, __old_feraiseexcept) +compat_symbol (libm, __old_feraiseexcept, feraiseexcept, GLIBC_2_1); +#endif + +libm_hidden_def (__feraiseexcept) +libm_hidden_ver (__feraiseexcept, feraiseexcept) +versioned_symbol (libm, __feraiseexcept, feraiseexcept, GLIBC_2_2); diff --git a/REORG.TODO/sysdeps/m68k/m680x0/fpu/halfulp.c b/REORG.TODO/sysdeps/m68k/m680x0/fpu/halfulp.c new file mode 100644 index 0000000000..1cc8931700 --- /dev/null +++ b/REORG.TODO/sysdeps/m68k/m680x0/fpu/halfulp.c @@ -0,0 +1 @@ +/* Not needed. */ diff --git a/REORG.TODO/sysdeps/m68k/m680x0/fpu/k_cosl.c b/REORG.TODO/sysdeps/m68k/m680x0/fpu/k_cosl.c new file mode 100644 index 0000000000..1cc8931700 --- /dev/null +++ b/REORG.TODO/sysdeps/m68k/m680x0/fpu/k_cosl.c @@ -0,0 +1 @@ +/* Not needed. */ diff --git a/REORG.TODO/sysdeps/m68k/m680x0/fpu/k_rem_pio2.c b/REORG.TODO/sysdeps/m68k/m680x0/fpu/k_rem_pio2.c new file mode 100644 index 0000000000..1347b0468c --- /dev/null +++ b/REORG.TODO/sysdeps/m68k/m680x0/fpu/k_rem_pio2.c @@ -0,0 +1,3 @@ +/* Empty. This file is only meant to avoid compiling the file with the + same name in the libm-ieee754 directory. The code is not used since + there is an assembler version for all users of this file. */ diff --git a/REORG.TODO/sysdeps/m68k/m680x0/fpu/k_rem_pio2f.c b/REORG.TODO/sysdeps/m68k/m680x0/fpu/k_rem_pio2f.c new file mode 100644 index 0000000000..1347b0468c --- /dev/null +++ b/REORG.TODO/sysdeps/m68k/m680x0/fpu/k_rem_pio2f.c @@ -0,0 +1,3 @@ +/* Empty. This file is only meant to avoid compiling the file with the + same name in the libm-ieee754 directory. The code is not used since + there is an assembler version for all users of this file. */ diff --git a/REORG.TODO/sysdeps/m68k/m680x0/fpu/k_sinl.c b/REORG.TODO/sysdeps/m68k/m680x0/fpu/k_sinl.c new file mode 100644 index 0000000000..1cc8931700 --- /dev/null +++ b/REORG.TODO/sysdeps/m68k/m680x0/fpu/k_sinl.c @@ -0,0 +1 @@ +/* Not needed. */ diff --git a/REORG.TODO/sysdeps/m68k/m680x0/fpu/k_tanf.c b/REORG.TODO/sysdeps/m68k/m680x0/fpu/k_tanf.c new file mode 100644 index 0000000000..1cc8931700 --- /dev/null +++ b/REORG.TODO/sysdeps/m68k/m680x0/fpu/k_tanf.c @@ -0,0 +1 @@ +/* Not needed. */ diff --git a/REORG.TODO/sysdeps/m68k/m680x0/fpu/k_tanl.c b/REORG.TODO/sysdeps/m68k/m680x0/fpu/k_tanl.c new file mode 100644 index 0000000000..1cc8931700 --- /dev/null +++ b/REORG.TODO/sysdeps/m68k/m680x0/fpu/k_tanl.c @@ -0,0 +1 @@ +/* Not needed. */ diff --git a/REORG.TODO/sysdeps/m68k/m680x0/fpu/libm-test-ulps b/REORG.TODO/sysdeps/m68k/m680x0/fpu/libm-test-ulps new file mode 100644 index 0000000000..2247a4af60 --- /dev/null +++ b/REORG.TODO/sysdeps/m68k/m680x0/fpu/libm-test-ulps @@ -0,0 +1,1602 @@ +# Begin of automatic generation + +# Maximal error of functions: +Function: "acosh": +double: 1 +float: 1 +idouble: 1 +ifloat: 1 +ildouble: 1 +ldouble: 1 + +Function: "acosh_downward": +double: 1 +float: 1 +idouble: 1 +ifloat: 1 +ildouble: 1 +ldouble: 1 + +Function: "acosh_towardzero": +double: 1 +float: 1 +idouble: 1 +ifloat: 1 +ildouble: 1 +ldouble: 1 + +Function: "acosh_upward": +double: 2 +float: 1 +idouble: 2 +ifloat: 1 +ildouble: 1 +ldouble: 1 + +Function: "asinh": +double: 1 +float: 1 +idouble: 1 +ifloat: 1 +ildouble: 1 +ldouble: 1 + +Function: "asinh_downward": +double: 2 +float: 2 +idouble: 2 +ifloat: 2 +ildouble: 2 +ldouble: 2 + +Function: "asinh_towardzero": +double: 2 +float: 2 +idouble: 2 +ifloat: 2 +ildouble: 2 +ldouble: 2 + +Function: "asinh_upward": +double: 2 +float: 2 +idouble: 2 +ifloat: 2 +ildouble: 3 +ldouble: 3 + +Function: "atan2": +float: 1 +ifloat: 1 +ildouble: 1 +ldouble: 1 + +Function: "atan2_downward": +double: 1 +float: 1 +idouble: 1 +ifloat: 1 +ildouble: 2 +ldouble: 2 + +Function: "atan2_towardzero": +double: 1 +float: 1 +idouble: 1 +ifloat: 1 +ildouble: 2 +ldouble: 2 + +Function: "atan2_upward": +double: 1 +float: 1 +idouble: 1 +ifloat: 1 +ildouble: 2 +ldouble: 2 + +Function: "atan_downward": +double: 1 +float: 1 +idouble: 1 +ifloat: 1 +ildouble: 1 +ldouble: 1 + +Function: "atan_towardzero": +float: 1 +ifloat: 1 +ildouble: 1 +ldouble: 1 + +Function: "atan_upward": +double: 1 +float: 1 +idouble: 1 +ifloat: 1 +ildouble: 1 +ldouble: 1 + +Function: "cabs": +double: 1 +idouble: 1 +ildouble: 1 +ldouble: 1 + +Function: "cabs_downward": +double: 1 +idouble: 1 +ildouble: 1 +ldouble: 1 + +Function: "cabs_towardzero": +double: 1 +idouble: 1 +ildouble: 1 +ldouble: 1 + +Function: "cabs_upward": +double: 1 +idouble: 1 +ildouble: 1 +ldouble: 1 + +Function: Real part of "cacos": +double: 1 +float: 2 +idouble: 1 +ifloat: 2 +ildouble: 1 +ldouble: 1 + +Function: Imaginary part of "cacos": +double: 1 +float: 1 +idouble: 1 +ifloat: 1 +ildouble: 2 +ldouble: 2 + +Function: Real part of "cacos_downward": +double: 2 +float: 2 +idouble: 2 +ifloat: 2 +ildouble: 2 +ldouble: 2 + +Function: Imaginary part of "cacos_downward": +double: 5 +float: 3 +idouble: 5 +ifloat: 3 +ildouble: 4 +ldouble: 4 + +Function: Real part of "cacos_towardzero": +double: 2 +float: 2 +idouble: 2 +ifloat: 2 +ildouble: 2 +ldouble: 2 + +Function: Imaginary part of "cacos_towardzero": +double: 5 +float: 3 +idouble: 5 +ifloat: 3 +ildouble: 4 +ldouble: 4 + +Function: Real part of "cacos_upward": +double: 2 +float: 2 +idouble: 2 +ifloat: 2 +ildouble: 3 +ldouble: 3 + +Function: Imaginary part of "cacos_upward": +double: 4 +float: 5 +idouble: 4 +ifloat: 5 +ildouble: 4 +ldouble: 4 + +Function: Real part of "cacosh": +double: 1 +float: 1 +idouble: 1 +ifloat: 1 +ildouble: 2 +ldouble: 2 + +Function: Imaginary part of "cacosh": +double: 1 +float: 2 +idouble: 1 +ifloat: 2 +ildouble: 1 +ldouble: 1 + +Function: Real part of "cacosh_downward": +double: 5 +float: 3 +idouble: 5 +ifloat: 3 +ildouble: 4 +ldouble: 4 + +Function: Imaginary part of "cacosh_downward": +double: 2 +float: 2 +idouble: 2 +ifloat: 2 +ildouble: 2 +ldouble: 2 + +Function: Real part of "cacosh_towardzero": +double: 5 +float: 3 +idouble: 5 +ifloat: 3 +ildouble: 4 +ldouble: 4 + +Function: Imaginary part of "cacosh_towardzero": +double: 2 +float: 2 +idouble: 2 +ifloat: 2 +ildouble: 2 +ldouble: 2 + +Function: Real part of "cacosh_upward": +double: 4 +float: 5 +idouble: 4 +ifloat: 5 +ildouble: 4 +ldouble: 4 + +Function: Imaginary part of "cacosh_upward": +double: 2 +float: 2 +idouble: 2 +ifloat: 2 +ildouble: 3 +ldouble: 3 + +Function: "carg": +float: 1 +ifloat: 1 +ildouble: 1 +ldouble: 1 + +Function: "carg_downward": +double: 1 +float: 1 +idouble: 1 +ifloat: 1 +ildouble: 1 +ldouble: 1 + +Function: "carg_towardzero": +float: 1 +ifloat: 1 +ildouble: 1 +ldouble: 1 + +Function: "carg_upward": +double: 1 +float: 1 +idouble: 1 +ifloat: 1 +ildouble: 1 +ldouble: 1 + +Function: Real part of "casin": +double: 1 +float: 1 +idouble: 1 +ifloat: 1 +ildouble: 1 +ldouble: 1 + +Function: Imaginary part of "casin": +double: 1 +float: 1 +idouble: 1 +ifloat: 1 +ildouble: 2 +ldouble: 2 + +Function: Real part of "casin_downward": +double: 2 +float: 2 +idouble: 2 +ifloat: 2 +ildouble: 2 +ldouble: 2 + +Function: Imaginary part of "casin_downward": +double: 5 +float: 3 +idouble: 5 +ifloat: 3 +ildouble: 4 +ldouble: 4 + +Function: Real part of "casin_towardzero": +double: 1 +float: 2 +idouble: 1 +ifloat: 2 +ildouble: 2 +ldouble: 2 + +Function: Imaginary part of "casin_towardzero": +double: 5 +float: 3 +idouble: 5 +ifloat: 3 +ildouble: 4 +ldouble: 4 + +Function: Real part of "casin_upward": +double: 2 +float: 2 +idouble: 2 +ifloat: 2 +ildouble: 2 +ldouble: 2 + +Function: Imaginary part of "casin_upward": +double: 4 +float: 5 +idouble: 4 +ifloat: 5 +ildouble: 4 +ldouble: 4 + +Function: Real part of "casinh": +double: 1 +float: 1 +idouble: 1 +ifloat: 1 +ildouble: 2 +ldouble: 2 + +Function: Imaginary part of "casinh": +double: 1 +float: 1 +idouble: 1 +ifloat: 1 +ildouble: 1 +ldouble: 1 + +Function: Real part of "casinh_downward": +double: 5 +float: 3 +idouble: 5 +ifloat: 3 +ildouble: 4 +ldouble: 4 + +Function: Imaginary part of "casinh_downward": +double: 2 +float: 2 +idouble: 2 +ifloat: 2 +ildouble: 2 +ldouble: 2 + +Function: Real part of "casinh_towardzero": +double: 5 +float: 3 +idouble: 5 +ifloat: 3 +ildouble: 4 +ldouble: 4 + +Function: Imaginary part of "casinh_towardzero": +double: 1 +float: 2 +idouble: 1 +ifloat: 2 +ildouble: 2 +ldouble: 2 + +Function: Real part of "casinh_upward": +double: 4 +float: 5 +idouble: 4 +ifloat: 5 +ildouble: 4 +ldouble: 4 + +Function: Imaginary part of "casinh_upward": +double: 2 +float: 2 +idouble: 2 +ifloat: 2 +ildouble: 2 +ldouble: 2 + +Function: Real part of "catan": +ildouble: 1 +ldouble: 1 + +Function: Imaginary part of "catan": +double: 1 +float: 1 +idouble: 1 +ifloat: 1 +ildouble: 1 +ldouble: 1 + +Function: Real part of "catan_downward": +double: 1 +float: 1 +idouble: 1 +ifloat: 1 +ildouble: 1 +ldouble: 1 + +Function: Imaginary part of "catan_downward": +double: 2 +float: 1 +idouble: 2 +ifloat: 1 +ildouble: 3 +ldouble: 3 + +Function: Real part of "catan_towardzero": +double: 1 +float: 1 +idouble: 1 +ifloat: 1 +ildouble: 1 +ldouble: 1 + +Function: Imaginary part of "catan_towardzero": +double: 2 +float: 1 +idouble: 2 +ifloat: 1 +ildouble: 3 +ldouble: 3 + +Function: Real part of "catan_upward": +double: 1 +float: 1 +idouble: 1 +ifloat: 1 +ildouble: 1 +ldouble: 1 + +Function: Imaginary part of "catan_upward": +double: 2 +float: 3 +idouble: 2 +ifloat: 3 +ildouble: 3 +ldouble: 3 + +Function: Real part of "catanh": +double: 1 +float: 1 +idouble: 1 +ifloat: 1 +ildouble: 1 +ldouble: 1 + +Function: Imaginary part of "catanh": +ildouble: 1 +ldouble: 1 + +Function: Real part of "catanh_downward": +double: 2 +float: 1 +idouble: 2 +ifloat: 1 +ildouble: 3 +ldouble: 3 + +Function: Imaginary part of "catanh_downward": +double: 1 +float: 1 +idouble: 1 +ifloat: 1 +ildouble: 1 +ldouble: 1 + +Function: Real part of "catanh_towardzero": +double: 2 +float: 1 +idouble: 2 +ifloat: 1 +ildouble: 3 +ldouble: 3 + +Function: Imaginary part of "catanh_towardzero": +double: 1 +float: 1 +idouble: 1 +ifloat: 1 +ildouble: 1 +ldouble: 1 + +Function: Real part of "catanh_upward": +double: 2 +float: 3 +idouble: 2 +ifloat: 3 +ildouble: 4 +ldouble: 4 + +Function: Imaginary part of "catanh_upward": +double: 1 +float: 1 +idouble: 1 +ifloat: 1 +ildouble: 1 +ldouble: 1 + +Function: "cbrt": +double: 1 +float: 1 +idouble: 1 +ifloat: 1 +ildouble: 1 +ldouble: 1 + +Function: "cbrt_downward": +double: 1 +float: 1 +idouble: 1 +ifloat: 1 +ildouble: 1 +ldouble: 1 + +Function: "cbrt_towardzero": +double: 1 +float: 1 +idouble: 1 +ifloat: 1 +ildouble: 1 +ldouble: 1 + +Function: "cbrt_upward": +double: 1 +float: 1 +idouble: 1 +ifloat: 1 +ildouble: 1 +ldouble: 1 + +Function: Real part of "ccos": +ildouble: 1 +ldouble: 1 + +Function: Imaginary part of "ccos": +ildouble: 1 +ldouble: 1 + +Function: Real part of "ccos_downward": +double: 1 +idouble: 1 +ildouble: 2 +ldouble: 2 + +Function: Imaginary part of "ccos_downward": +double: 1 +idouble: 1 +ildouble: 2 +ldouble: 2 + +Function: Real part of "ccos_towardzero": +double: 1 +idouble: 1 +ildouble: 2 +ldouble: 2 + +Function: Imaginary part of "ccos_towardzero": +double: 1 +idouble: 1 +ildouble: 2 +ldouble: 2 + +Function: Real part of "ccos_upward": +double: 1 +idouble: 1 +ildouble: 2 +ldouble: 2 + +Function: Imaginary part of "ccos_upward": +double: 1 +idouble: 1 +ildouble: 2 +ldouble: 2 + +Function: Imaginary part of "ccosh": +ildouble: 1 +ldouble: 1 + +Function: Real part of "ccosh_downward": +double: 1 +idouble: 1 +ildouble: 2 +ldouble: 2 + +Function: Imaginary part of "ccosh_downward": +double: 1 +idouble: 1 +ildouble: 2 +ldouble: 2 + +Function: Real part of "ccosh_towardzero": +double: 1 +idouble: 1 +ildouble: 2 +ldouble: 2 + +Function: Imaginary part of "ccosh_towardzero": +double: 1 +idouble: 1 +ildouble: 2 +ldouble: 2 + +Function: Real part of "ccosh_upward": +double: 1 +idouble: 1 +ildouble: 2 +ldouble: 2 + +Function: Imaginary part of "ccosh_upward": +double: 1 +idouble: 1 +ildouble: 2 +ldouble: 2 + +Function: Real part of "cexp": +ildouble: 1 +ldouble: 1 + +Function: Imaginary part of "cexp": +ildouble: 1 +ldouble: 1 + +Function: Real part of "cexp_downward": +double: 1 +float: 1 +idouble: 1 +ifloat: 1 +ildouble: 3 +ldouble: 3 + +Function: Imaginary part of "cexp_downward": +double: 1 +float: 1 +idouble: 1 +ifloat: 1 +ildouble: 3 +ldouble: 3 + +Function: Real part of "cexp_towardzero": +double: 1 +float: 1 +idouble: 1 +ifloat: 1 +ildouble: 3 +ldouble: 3 + +Function: Imaginary part of "cexp_towardzero": +double: 1 +float: 1 +idouble: 1 +ifloat: 1 +ildouble: 3 +ldouble: 3 + +Function: Real part of "cexp_upward": +double: 1 +idouble: 1 +ildouble: 2 +ldouble: 2 + +Function: Imaginary part of "cexp_upward": +double: 1 +idouble: 1 +ildouble: 2 +ldouble: 2 + +Function: Real part of "clog": +double: 3 +float: 2 +idouble: 3 +ifloat: 2 +ildouble: 3 +ldouble: 3 + +Function: Imaginary part of "clog": +double: 1 +float: 1 +idouble: 1 +ifloat: 1 +ildouble: 1 +ldouble: 1 + +Function: Real part of "clog10": +double: 2 +float: 2 +idouble: 2 +ifloat: 2 +ildouble: 3 +ldouble: 3 + +Function: Imaginary part of "clog10": +double: 1 +float: 1 +idouble: 1 +ifloat: 1 +ildouble: 2 +ldouble: 2 + +Function: Real part of "clog10_downward": +double: 6 +float: 4 +idouble: 6 +ifloat: 4 +ildouble: 7 +ldouble: 7 + +Function: Imaginary part of "clog10_downward": +double: 1 +float: 2 +idouble: 1 +ifloat: 2 +ildouble: 4 +ldouble: 4 + +Function: Real part of "clog10_towardzero": +double: 6 +float: 4 +idouble: 6 +ifloat: 4 +ildouble: 7 +ldouble: 7 + +Function: Imaginary part of "clog10_towardzero": +double: 1 +float: 2 +idouble: 1 +ifloat: 2 +ildouble: 4 +ldouble: 4 + +Function: Real part of "clog10_upward": +double: 7 +float: 5 +idouble: 7 +ifloat: 5 +ildouble: 6 +ldouble: 6 + +Function: Imaginary part of "clog10_upward": +double: 1 +float: 2 +idouble: 1 +ifloat: 2 +ildouble: 3 +ldouble: 3 + +Function: Real part of "clog_downward": +double: 7 +float: 5 +idouble: 7 +ifloat: 5 +ildouble: 8 +ldouble: 8 + +Function: Imaginary part of "clog_downward": +double: 1 +float: 1 +idouble: 1 +ifloat: 1 +ildouble: 2 +ldouble: 2 + +Function: Real part of "clog_towardzero": +double: 7 +float: 5 +idouble: 7 +ifloat: 5 +ildouble: 9 +ldouble: 9 + +Function: Imaginary part of "clog_towardzero": +double: 1 +float: 1 +idouble: 1 +ifloat: 1 +ildouble: 2 +ldouble: 2 + +Function: Real part of "clog_upward": +double: 8 +float: 5 +idouble: 8 +ifloat: 5 +ildouble: 7 +ldouble: 7 + +Function: Imaginary part of "clog_upward": +double: 1 +float: 1 +idouble: 1 +ifloat: 1 +ildouble: 2 +ldouble: 2 + +Function: "cos": +double: 1 +idouble: 1 + +Function: Real part of "cpow": +double: 1 +float: 3 +idouble: 1 +ifloat: 3 +ildouble: 3 +ldouble: 3 + +Function: Imaginary part of "cpow": +float: 5 +ifloat: 5 +ildouble: 1 +ldouble: 1 + +Function: Real part of "cpow_downward": +double: 5 +float: 8 +idouble: 5 +ifloat: 8 +ildouble: 7 +ldouble: 7 + +Function: Imaginary part of "cpow_downward": +double: 1 +float: 1 +idouble: 1 +ifloat: 1 +ildouble: 2 +ldouble: 2 + +Function: Real part of "cpow_towardzero": +double: 5 +float: 8 +idouble: 5 +ifloat: 8 +ildouble: 7 +ldouble: 7 + +Function: Imaginary part of "cpow_towardzero": +double: 1 +float: 1 +idouble: 1 +ifloat: 1 +ildouble: 2 +ldouble: 2 + +Function: Real part of "cpow_upward": +double: 4 +float: 1 +idouble: 4 +ifloat: 1 +ildouble: 2 +ldouble: 2 + +Function: Imaginary part of "cpow_upward": +double: 1 +float: 1 +idouble: 1 +ifloat: 1 +ildouble: 2 +ldouble: 2 + +Function: Real part of "csin": +ildouble: 1 +ldouble: 1 + +Function: Real part of "csin_downward": +double: 1 +float: 1 +idouble: 1 +ifloat: 1 +ildouble: 2 +ldouble: 2 + +Function: Imaginary part of "csin_downward": +double: 1 +float: 1 +idouble: 1 +ifloat: 1 +ildouble: 1 +ldouble: 1 + +Function: Real part of "csin_towardzero": +double: 1 +float: 1 +idouble: 1 +ifloat: 1 +ildouble: 2 +ldouble: 2 + +Function: Imaginary part of "csin_towardzero": +double: 1 +float: 1 +idouble: 1 +ifloat: 1 +ildouble: 1 +ldouble: 1 + +Function: Real part of "csin_upward": +double: 1 +float: 1 +idouble: 1 +ifloat: 1 +ildouble: 2 +ldouble: 2 + +Function: Imaginary part of "csin_upward": +double: 1 +float: 1 +idouble: 1 +ifloat: 1 +ildouble: 1 +ldouble: 1 + +Function: Real part of "csinh": +ildouble: 1 +ldouble: 1 + +Function: Real part of "csinh_downward": +double: 1 +float: 1 +idouble: 1 +ifloat: 1 +ildouble: 1 +ldouble: 1 + +Function: Imaginary part of "csinh_downward": +double: 1 +float: 1 +idouble: 1 +ifloat: 1 +ildouble: 2 +ldouble: 2 + +Function: Real part of "csinh_towardzero": +double: 1 +float: 1 +idouble: 1 +ifloat: 1 +ildouble: 1 +ldouble: 1 + +Function: Imaginary part of "csinh_towardzero": +double: 1 +float: 1 +idouble: 1 +ifloat: 1 +ildouble: 2 +ldouble: 2 + +Function: Real part of "csinh_upward": +double: 1 +float: 1 +idouble: 1 +ifloat: 1 +ildouble: 1 +ldouble: 1 + +Function: Imaginary part of "csinh_upward": +double: 1 +float: 1 +idouble: 1 +ifloat: 1 +ildouble: 2 +ldouble: 2 + +Function: Real part of "csqrt": +double: 1 +float: 1 +idouble: 1 +ifloat: 1 +ildouble: 2 +ldouble: 2 + +Function: Imaginary part of "csqrt": +double: 1 +float: 1 +idouble: 1 +ifloat: 1 +ildouble: 2 +ldouble: 2 + +Function: Real part of "csqrt_downward": +double: 1 +float: 2 +idouble: 1 +ifloat: 2 +ildouble: 4 +ldouble: 4 + +Function: Imaginary part of "csqrt_downward": +double: 1 +float: 2 +idouble: 1 +ifloat: 2 +ildouble: 3 +ldouble: 3 + +Function: Real part of "csqrt_towardzero": +double: 1 +float: 1 +idouble: 1 +ifloat: 1 +ildouble: 4 +ldouble: 4 + +Function: Imaginary part of "csqrt_towardzero": +double: 1 +float: 2 +idouble: 1 +ifloat: 2 +ildouble: 3 +ldouble: 3 + +Function: Real part of "csqrt_upward": +double: 1 +float: 2 +idouble: 1 +ifloat: 2 +ildouble: 5 +ldouble: 5 + +Function: Imaginary part of "csqrt_upward": +double: 2 +float: 2 +idouble: 2 +ifloat: 2 +ildouble: 3 +ldouble: 3 + +Function: Real part of "ctan": +double: 1 +float: 1 +idouble: 1 +ifloat: 1 +ildouble: 2 +ldouble: 2 + +Function: Imaginary part of "ctan": +double: 1 +float: 1 +idouble: 1 +ifloat: 1 +ildouble: 2 +ldouble: 2 + +Function: Real part of "ctan_downward": +double: 2 +float: 2 +idouble: 2 +ifloat: 2 +ildouble: 2 +ldouble: 2 + +Function: Imaginary part of "ctan_downward": +double: 3 +float: 2 +idouble: 3 +ifloat: 2 +ildouble: 2 +ldouble: 2 + +Function: Real part of "ctan_towardzero": +double: 1 +float: 2 +idouble: 1 +ifloat: 2 +ildouble: 2 +ldouble: 2 + +Function: Imaginary part of "ctan_towardzero": +double: 3 +float: 2 +idouble: 3 +ifloat: 2 +ildouble: 4 +ldouble: 4 + +Function: Real part of "ctan_upward": +double: 1 +float: 3 +idouble: 1 +ifloat: 3 +ildouble: 5 +ldouble: 5 + +Function: Imaginary part of "ctan_upward": +double: 1 +float: 1 +idouble: 1 +ifloat: 1 +ildouble: 7 +ldouble: 7 + +Function: Real part of "ctanh": +double: 1 +float: 1 +idouble: 1 +ifloat: 1 +ildouble: 2 +ldouble: 2 + +Function: Imaginary part of "ctanh": +double: 1 +float: 2 +idouble: 1 +ifloat: 2 +ildouble: 2 +ldouble: 2 + +Function: Real part of "ctanh_downward": +double: 3 +float: 2 +idouble: 3 +ifloat: 2 +ildouble: 2 +ldouble: 2 + +Function: Imaginary part of "ctanh_downward": +double: 2 +float: 2 +idouble: 2 +ifloat: 2 +ildouble: 2 +ldouble: 2 + +Function: Real part of "ctanh_towardzero": +double: 3 +float: 2 +idouble: 3 +ifloat: 2 +ildouble: 4 +ldouble: 4 + +Function: Imaginary part of "ctanh_towardzero": +double: 1 +float: 1 +idouble: 1 +ifloat: 1 +ildouble: 2 +ldouble: 2 + +Function: Real part of "ctanh_upward": +double: 1 +float: 2 +idouble: 1 +ifloat: 2 +ildouble: 7 +ldouble: 7 + +Function: Imaginary part of "ctanh_upward": +double: 1 +float: 3 +idouble: 1 +ifloat: 3 +ildouble: 5 +ldouble: 5 + +Function: "erf": +float: 1 +ifloat: 1 +ildouble: 1 +ldouble: 1 + +Function: "erf_downward": +float: 1 +ifloat: 1 + +Function: "erf_towardzero": +ildouble: 1 +ldouble: 1 + +Function: "erf_upward": +float: 1 +ifloat: 1 +ildouble: 1 +ldouble: 1 + +Function: "erfc": +float: 1 +ifloat: 1 +ildouble: 2 +ldouble: 2 + +Function: "erfc_downward": +double: 1 +float: 3 +idouble: 1 +ifloat: 3 +ildouble: 4 +ldouble: 4 + +Function: "erfc_towardzero": +double: 1 +float: 1 +idouble: 1 +ifloat: 1 +ildouble: 4 +ldouble: 4 + +Function: "erfc_upward": +double: 1 +float: 3 +idouble: 1 +ifloat: 3 +ildouble: 3 +ldouble: 3 + +Function: "exp2": +double: 1 +idouble: 1 + +Function: "gamma": +float: 1 +ifloat: 1 +ildouble: 2 +ldouble: 2 + +Function: "gamma_downward": +double: 1 +float: 1 +idouble: 1 +ifloat: 1 +ildouble: 3 +ldouble: 3 + +Function: "gamma_towardzero": +double: 1 +float: 1 +idouble: 1 +ifloat: 1 +ildouble: 3 +ldouble: 3 + +Function: "gamma_upward": +double: 1 +float: 1 +idouble: 1 +ifloat: 1 +ildouble: 2 +ldouble: 2 + +Function: "hypot": +double: 1 +idouble: 1 +ildouble: 1 +ldouble: 1 + +Function: "hypot_downward": +double: 1 +idouble: 1 +ildouble: 1 +ldouble: 1 + +Function: "hypot_towardzero": +double: 1 +idouble: 1 +ildouble: 1 +ldouble: 1 + +Function: "hypot_upward": +double: 1 +idouble: 1 +ildouble: 1 +ldouble: 1 + +Function: "j0": +double: 1 +float: 2 +idouble: 1 +ifloat: 2 +ildouble: 2 +ldouble: 2 + +Function: "j0_downward": +double: 4 +float: 3 +idouble: 4 +ifloat: 3 +ildouble: 4 +ldouble: 4 + +Function: "j0_towardzero": +double: 1 +float: 1 +idouble: 1 +ifloat: 1 +ildouble: 5 +ldouble: 5 + +Function: "j0_upward": +double: 2 +float: 3 +idouble: 2 +ifloat: 3 +ildouble: 5 +ldouble: 5 + +Function: "j1": +float: 2 +ifloat: 2 +ildouble: 1 +ldouble: 1 + +Function: "j1_downward": +double: 1 +float: 4 +idouble: 1 +ifloat: 4 +ildouble: 5 +ldouble: 5 + +Function: "j1_towardzero": +double: 1 +float: 2 +idouble: 1 +ifloat: 2 +ildouble: 4 +ldouble: 4 + +Function: "j1_upward": +double: 1 +float: 3 +idouble: 1 +ifloat: 3 +ildouble: 3 +ldouble: 3 + +Function: "jn": +double: 2 +float: 2 +idouble: 2 +ifloat: 2 +ildouble: 4 +ldouble: 4 + +Function: "jn_downward": +double: 2 +float: 4 +idouble: 2 +ifloat: 4 +ildouble: 5 +ldouble: 5 + +Function: "jn_towardzero": +double: 2 +float: 3 +idouble: 2 +ifloat: 3 +ildouble: 5 +ldouble: 5 + +Function: "jn_upward": +double: 2 +float: 3 +idouble: 2 +ifloat: 3 +ildouble: 5 +ldouble: 5 + +Function: "lgamma": +float: 1 +ifloat: 1 +ildouble: 2 +ldouble: 2 + +Function: "lgamma_downward": +double: 1 +float: 1 +idouble: 1 +ifloat: 1 +ildouble: 3 +ldouble: 3 + +Function: "lgamma_towardzero": +double: 1 +float: 1 +idouble: 1 +ifloat: 1 +ildouble: 3 +ldouble: 3 + +Function: "lgamma_upward": +double: 1 +float: 1 +idouble: 1 +ifloat: 1 +ildouble: 2 +ldouble: 2 + +Function: "pow": +double: 1 +float: 7 +idouble: 1 +ifloat: 7 +ildouble: 9 +ldouble: 9 + +Function: "pow_downward": +double: 1 +float: 7 +idouble: 1 +ifloat: 7 +ildouble: 9 +ldouble: 9 + +Function: "pow_towardzero": +double: 1 +float: 4 +idouble: 1 +ifloat: 4 +ildouble: 4 +ldouble: 4 + +Function: "pow_upward": +double: 1 +float: 2 +idouble: 1 +ifloat: 2 +ildouble: 8 +ldouble: 8 + +Function: "sin": +double: 1 +idouble: 1 + +Function: "tgamma": +double: 1 +float: 4 +idouble: 1 +ifloat: 4 +ildouble: 9 +ldouble: 9 + +Function: "tgamma_downward": +double: 2 +float: 3 +idouble: 2 +ifloat: 3 +ildouble: 9 +ldouble: 9 + +Function: "tgamma_towardzero": +double: 2 +float: 4 +idouble: 2 +ifloat: 4 +ildouble: 9 +ldouble: 9 + +Function: "tgamma_upward": +double: 1 +float: 4 +idouble: 1 +ifloat: 4 +ildouble: 9 +ldouble: 9 + +Function: "y0": +double: 1 +float: 1 +idouble: 1 +ifloat: 1 +ildouble: 1 +ldouble: 1 + +Function: "y0_downward": +double: 1 +float: 3 +idouble: 1 +ifloat: 3 +ildouble: 6 +ldouble: 6 + +Function: "y0_towardzero": +double: 1 +float: 2 +idouble: 1 +ifloat: 2 +ildouble: 5 +ldouble: 5 + +Function: "y0_upward": +double: 1 +float: 3 +idouble: 1 +ifloat: 3 +ildouble: 3 +ldouble: 3 + +Function: "y1": +double: 1 +float: 3 +idouble: 1 +ifloat: 3 +ildouble: 2 +ldouble: 2 + +Function: "y1_downward": +double: 4 +float: 3 +idouble: 4 +ifloat: 3 +ildouble: 6 +ldouble: 6 + +Function: "y1_towardzero": +double: 2 +float: 3 +idouble: 2 +ifloat: 3 +ildouble: 5 +ldouble: 5 + +Function: "y1_upward": +double: 2 +float: 2 +idouble: 2 +ifloat: 2 +ildouble: 7 +ldouble: 7 + +Function: "yn": +double: 2 +float: 3 +idouble: 2 +ifloat: 3 +ildouble: 4 +ldouble: 4 + +Function: "yn_downward": +double: 2 +float: 3 +idouble: 2 +ifloat: 3 +ildouble: 6 +ldouble: 6 + +Function: "yn_towardzero": +double: 3 +float: 3 +idouble: 3 +ifloat: 3 +ildouble: 5 +ldouble: 5 + +Function: "yn_upward": +double: 3 +float: 3 +idouble: 3 +ifloat: 3 +ildouble: 4 +ldouble: 4 + +# end of automatic generation diff --git a/REORG.TODO/sysdeps/m68k/m680x0/fpu/libm-test-ulps-name b/REORG.TODO/sysdeps/m68k/m680x0/fpu/libm-test-ulps-name new file mode 100644 index 0000000000..35c8a9e5f0 --- /dev/null +++ b/REORG.TODO/sysdeps/m68k/m680x0/fpu/libm-test-ulps-name @@ -0,0 +1 @@ +M68k diff --git a/REORG.TODO/sysdeps/m68k/m680x0/fpu/math_private.h b/REORG.TODO/sysdeps/m68k/m680x0/fpu/math_private.h new file mode 100644 index 0000000000..2b7909535e --- /dev/null +++ b/REORG.TODO/sysdeps/m68k/m680x0/fpu/math_private.h @@ -0,0 +1,20 @@ +#ifndef M68K_MATH_PRIVATE_H +#define M68K_MATH_PRIVATE_H 1 + +#define math_opt_barrier(x) \ +({ __typeof (x) __x; \ + __asm ("" : "=f" (__x) : "0" (x)); \ + __x; }) +#define math_force_eval(x) \ +do \ + { \ + __typeof (x) __x = (x); \ + if (sizeof (x) <= sizeof (double)) \ + __asm __volatile ("" : : "m" (__x)); \ + else \ + __asm __volatile ("" : : "f" (__x)); \ + } \ +while (0) + +#include_next <math_private.h> +#endif diff --git a/REORG.TODO/sysdeps/m68k/m680x0/fpu/mathimpl.h b/REORG.TODO/sysdeps/m68k/m680x0/fpu/mathimpl.h new file mode 100644 index 0000000000..653df96873 --- /dev/null +++ b/REORG.TODO/sysdeps/m68k/m680x0/fpu/mathimpl.h @@ -0,0 +1,93 @@ +/* Definitions of libc internal inline math functions implemented + by the m68881/2. + Copyright (C) 1991-2017 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 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. + + 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 + Lesser General Public License for more details. + + You should have received a copy of the GNU Lesser General Public + License along with the GNU C Library. If not, see + <http://www.gnu.org/licenses/>. */ + +/* This file contains the definitions of the inline math functions that + are only used internally inside libm, not visible to the user. */ + +__inline_mathop (__ieee754_acos, acos) +__inline_mathop (__ieee754_asin, asin) +__inline_mathop (__ieee754_cosh, cosh) +__inline_mathop (__ieee754_sinh, sinh) +__inline_mathop (__ieee754_exp, etox) +__inline_mathop (__ieee754_exp2, twotox) +__inline_mathop (__ieee754_exp10, tentox) +__inline_mathop (__ieee754_log10, log10) +__inline_mathop (__ieee754_log2, log2) +__inline_mathop (__ieee754_log, logn) +__inline_mathop (__ieee754_sqrt, sqrt) +__inline_mathop (__ieee754_atanh, atanh) + +__m81_defun (double, __ieee754_remainder, (double __x, double __y)) +{ + double __result; + __asm ("frem%.x %1, %0" : "=f" (__result) : "f" (__y), "0" (__x)); + return __result; +} + +__m81_defun (float, __ieee754_remainderf, (float __x, float __y)) +{ + float __result; + __asm ("frem%.x %1, %0" : "=f" (__result) : "f" (__y), "0" (__x)); + return __result; +} + +__m81_defun (long double, + __ieee754_remainderl, (long double __x, long double __y)) +{ + long double __result; + __asm ("frem%.x %1, %0" : "=f" (__result) : "f" (__y), "0" (__x)); + return __result; +} + +__m81_defun (double, __ieee754_fmod, (double __x, double __y)) +{ + double __result; + __asm ("fmod%.x %1, %0" : "=f" (__result) : "f" (__y), "0" (__x)); + return __result; +} + +__m81_defun (float, __ieee754_fmodf, (float __x, float __y)) +{ + float __result; + __asm ("fmod%.x %1, %0" : "=f" (__result) : "f" (__y), "0" (__x)); + return __result; +} + +__m81_defun (long double, + __ieee754_fmodl, (long double __x, long double __y)) +{ + long double __result; + __asm ("fmod%.x %1, %0" : "=f" (__result) : "f" (__y), "0" (__x)); + return __result; +} + +/* Get the m68881 condition codes, to quickly check multiple conditions. */ +static __inline__ unsigned long +__m81_test (long double __val) +{ + unsigned long __fpsr; + __asm ("ftst%.x %1; fmove%.l %/fpsr,%0" : "=dm" (__fpsr) : "f" (__val)); + return __fpsr; +} + +/* Bit values returned by __m81_test. */ +#define __M81_COND_NAN (1 << 24) +#define __M81_COND_INF (2 << 24) +#define __M81_COND_ZERO (4 << 24) +#define __M81_COND_NEG (8 << 24) diff --git a/REORG.TODO/sysdeps/m68k/m680x0/fpu/mpa.c b/REORG.TODO/sysdeps/m68k/m680x0/fpu/mpa.c new file mode 100644 index 0000000000..1cc8931700 --- /dev/null +++ b/REORG.TODO/sysdeps/m68k/m680x0/fpu/mpa.c @@ -0,0 +1 @@ +/* Not needed. */ diff --git a/REORG.TODO/sysdeps/m68k/m680x0/fpu/mpatan.c b/REORG.TODO/sysdeps/m68k/m680x0/fpu/mpatan.c new file mode 100644 index 0000000000..1cc8931700 --- /dev/null +++ b/REORG.TODO/sysdeps/m68k/m680x0/fpu/mpatan.c @@ -0,0 +1 @@ +/* Not needed. */ diff --git a/REORG.TODO/sysdeps/m68k/m680x0/fpu/mpatan2.c b/REORG.TODO/sysdeps/m68k/m680x0/fpu/mpatan2.c new file mode 100644 index 0000000000..1cc8931700 --- /dev/null +++ b/REORG.TODO/sysdeps/m68k/m680x0/fpu/mpatan2.c @@ -0,0 +1 @@ +/* Not needed. */ diff --git a/REORG.TODO/sysdeps/m68k/m680x0/fpu/mpexp.c b/REORG.TODO/sysdeps/m68k/m680x0/fpu/mpexp.c new file mode 100644 index 0000000000..1cc8931700 --- /dev/null +++ b/REORG.TODO/sysdeps/m68k/m680x0/fpu/mpexp.c @@ -0,0 +1 @@ +/* Not needed. */ diff --git a/REORG.TODO/sysdeps/m68k/m680x0/fpu/mplog.c b/REORG.TODO/sysdeps/m68k/m680x0/fpu/mplog.c new file mode 100644 index 0000000000..1cc8931700 --- /dev/null +++ b/REORG.TODO/sysdeps/m68k/m680x0/fpu/mplog.c @@ -0,0 +1 @@ +/* Not needed. */ diff --git a/REORG.TODO/sysdeps/m68k/m680x0/fpu/mpsqrt.c b/REORG.TODO/sysdeps/m68k/m680x0/fpu/mpsqrt.c new file mode 100644 index 0000000000..1cc8931700 --- /dev/null +++ b/REORG.TODO/sysdeps/m68k/m680x0/fpu/mpsqrt.c @@ -0,0 +1 @@ +/* Not needed. */ diff --git a/REORG.TODO/sysdeps/m68k/m680x0/fpu/mptan.c b/REORG.TODO/sysdeps/m68k/m680x0/fpu/mptan.c new file mode 100644 index 0000000000..1cc8931700 --- /dev/null +++ b/REORG.TODO/sysdeps/m68k/m680x0/fpu/mptan.c @@ -0,0 +1 @@ +/* Not needed. */ diff --git a/REORG.TODO/sysdeps/m68k/m680x0/fpu/s_atan.c b/REORG.TODO/sysdeps/m68k/m680x0/fpu/s_atan.c new file mode 100644 index 0000000000..55e26312fc --- /dev/null +++ b/REORG.TODO/sysdeps/m68k/m680x0/fpu/s_atan.c @@ -0,0 +1,34 @@ +/* Copyright (C) 1996-2017 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 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. + + 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 + Lesser General Public License for more details. + + You should have received a copy of the GNU Lesser General Public + License along with the GNU C Library. If not, see + <http://www.gnu.org/licenses/>. */ + +#include <math.h> + +#ifndef FUNC +#define FUNC atan +#endif +#ifndef float_type +#define float_type double +#endif + +#define __CONCATX(a,b) __CONCAT(a,b) + +float_type +__CONCATX(__,FUNC) (float_type x) +{ + return __m81_u(__CONCATX(__,FUNC))(x); +} +weak_alias (__CONCATX(__,FUNC), FUNC) diff --git a/REORG.TODO/sysdeps/m68k/m680x0/fpu/s_atanf.c b/REORG.TODO/sysdeps/m68k/m680x0/fpu/s_atanf.c new file mode 100644 index 0000000000..c98559a8ba --- /dev/null +++ b/REORG.TODO/sysdeps/m68k/m680x0/fpu/s_atanf.c @@ -0,0 +1,5 @@ +#ifndef FUNC +#define FUNC atanf +#endif +#define float_type float +#include <s_atan.c> diff --git a/REORG.TODO/sysdeps/m68k/m680x0/fpu/s_atanl.c b/REORG.TODO/sysdeps/m68k/m680x0/fpu/s_atanl.c new file mode 100644 index 0000000000..b7e608addd --- /dev/null +++ b/REORG.TODO/sysdeps/m68k/m680x0/fpu/s_atanl.c @@ -0,0 +1,5 @@ +#ifndef FUNC +#define FUNC atanl +#endif +#define float_type long double +#include <s_atan.c> diff --git a/REORG.TODO/sysdeps/m68k/m680x0/fpu/s_ccosh_template.c b/REORG.TODO/sysdeps/m68k/m680x0/fpu/s_ccosh_template.c new file mode 100644 index 0000000000..26d11fc740 --- /dev/null +++ b/REORG.TODO/sysdeps/m68k/m680x0/fpu/s_ccosh_template.c @@ -0,0 +1,68 @@ +/* Complex cosine hyperbole function. m68k fpu version + Copyright (C) 1997-2017 Free Software Foundation, Inc. + This file is part of the GNU C Library. + Contributed by Andreas Schwab <schwab@issan.informatik.uni-dortmund.de>. + + The GNU C Library 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. + + 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 + Lesser General Public License for more details. + + You should have received a copy of the GNU Lesser General Public + License along with the GNU C Library. If not, see + <http://www.gnu.org/licenses/>. */ + +#include <complex.h> +#include <math.h> +#include "mathimpl.h" + +#define s(name) M_SUF (name) +#define m81(func) __m81_u(s(func)) + +CFLOAT +s(__ccosh) (CFLOAT x) +{ + CFLOAT retval; + unsigned long ix_cond = __m81_test (__imag__ x); + + if ((ix_cond & (__M81_COND_INF|__M81_COND_NAN)) == 0) + { + /* Imaginary part is finite. */ + FLOAT sin_ix, cos_ix; + + __asm ("fsincos%.x %2,%1:%0" : "=f" (sin_ix), "=f" (cos_ix) + : "f" (__imag__ x)); + __real__ retval = cos_ix * m81(__ieee754_cosh) (__real__ x); + if (ix_cond & __M81_COND_ZERO) + __imag__ retval = (signbit (__real__ x) + ? -__imag__ x : __imag__ x); + else + __imag__ retval = sin_ix * m81(__ieee754_sinh) (__real__ x); + } + else + { + unsigned long rx_cond = __m81_test (__real__ x); + + if (rx_cond & __M81_COND_ZERO) + { + __real__ retval = __imag__ x - __imag__ x; + __imag__ retval = __real__ x; + } + else + { + if (rx_cond & __M81_COND_INF) + __real__ retval = s(fabs) (__real__ x); + else + __real__ retval = s(__nan) (""); + __imag__ retval = __imag__ x - __imag__ x; + } + } + + return retval; +} +weak_alias (s(__ccosh), s(ccosh)) diff --git a/REORG.TODO/sysdeps/m68k/m680x0/fpu/s_ceil.c b/REORG.TODO/sysdeps/m68k/m680x0/fpu/s_ceil.c new file mode 100644 index 0000000000..93d5ad72e5 --- /dev/null +++ b/REORG.TODO/sysdeps/m68k/m680x0/fpu/s_ceil.c @@ -0,0 +1,2 @@ +#define FUNC ceil +#include <s_atan.c> diff --git a/REORG.TODO/sysdeps/m68k/m680x0/fpu/s_ceilf.c b/REORG.TODO/sysdeps/m68k/m680x0/fpu/s_ceilf.c new file mode 100644 index 0000000000..b3ba6a5700 --- /dev/null +++ b/REORG.TODO/sysdeps/m68k/m680x0/fpu/s_ceilf.c @@ -0,0 +1,2 @@ +#define FUNC ceilf +#include <s_atanf.c> diff --git a/REORG.TODO/sysdeps/m68k/m680x0/fpu/s_ceill.c b/REORG.TODO/sysdeps/m68k/m680x0/fpu/s_ceill.c new file mode 100644 index 0000000000..2bf95b00c6 --- /dev/null +++ b/REORG.TODO/sysdeps/m68k/m680x0/fpu/s_ceill.c @@ -0,0 +1,2 @@ +#define FUNC ceill +#include <s_atanl.c> diff --git a/REORG.TODO/sysdeps/m68k/m680x0/fpu/s_cexp_template.c b/REORG.TODO/sysdeps/m68k/m680x0/fpu/s_cexp_template.c new file mode 100644 index 0000000000..17cf145685 --- /dev/null +++ b/REORG.TODO/sysdeps/m68k/m680x0/fpu/s_cexp_template.c @@ -0,0 +1,129 @@ +/* Complex exponential function. m68k fpu version + Copyright (C) 1997-2017 Free Software Foundation, Inc. + This file is part of the GNU C Library. + Contributed by Andreas Schwab <schwab@issan.informatik.uni-dortmund.de> + + The GNU C Library 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. + + 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 + Lesser General Public License for more details. + + You should have received a copy of the GNU Lesser General Public + License along with the GNU C Library. If not, see + <http://www.gnu.org/licenses/>. */ + +#include <float.h> +#include <complex.h> +#include <math.h> +#include "mathimpl.h" + +#define CONCATX(a,b) __CONCAT(a,b) +#define s(name) M_SUF (name) +#define m81(func) __m81_u(s(func)) + +CFLOAT +s(__cexp) (CFLOAT x) +{ + CFLOAT retval; + unsigned long ix_cond; + + ix_cond = __m81_test (__imag__ x); + + if ((ix_cond & (__M81_COND_NAN|__M81_COND_INF)) == 0) + { + /* Imaginary part is finite. */ + unsigned long rx_cond = __m81_test (__real__ x); + + if ((rx_cond & (__M81_COND_NAN|__M81_COND_INF)) == 0) + { + const int t = (int) ((LDBL_MAX_EXP - 1) * M_LN2l); + long double sin_ix, cos_ix, exp_val; + + __m81_u (__sincosl) (__imag__ x, &sin_ix, &cos_ix); + + if (__real__ x > t) + { + long double exp_t = __m81_u(__ieee754_expl) (t); + __real__ x -= t; + sin_ix *= exp_t; + cos_ix *= exp_t; + if (__real__ x > t) + { + __real__ x -= t; + sin_ix *= exp_t; + cos_ix *= exp_t; + } + } + + exp_val = __m81_u(__ieee754_expl) (__real__ x); + __real__ retval = exp_val * cos_ix; + if (ix_cond & __M81_COND_ZERO) + __imag__ retval = __imag__ x; + else + __imag__ retval = exp_val * sin_ix; + } + else + { + /* Compute the sign of the result. */ + long double remainder, pi_2; + int quadrant; + + if ((rx_cond & (__M81_COND_NAN|__M81_COND_NEG)) == __M81_COND_NEG) + __real__ retval = __imag__ retval = 0.0; + else + __real__ retval = __imag__ retval = __real__ x; + __asm ("fmovecr %#0,%0\n\tfscale%.w %#-1,%0" : "=f" (pi_2)); + __asm ("fmod%.x %2,%0\n\tfmove%.l %/fpsr,%1" + : "=f" (remainder), "=dm" (quadrant) + : "f" (pi_2), "0" (__imag__ x)); + quadrant = (quadrant >> 16) & 0x83; + if (quadrant & 0x80) + quadrant ^= 0x83; + switch (quadrant) + { + default: + break; + case 1: + __real__ retval = -__real__ retval; + break; + case 2: + __real__ retval = -__real__ retval; + case 3: + __imag__ retval = -__imag__ retval; + break; + } + if (ix_cond & __M81_COND_ZERO && (rx_cond & __M81_COND_NAN) == 0) + __imag__ retval = __imag__ x; + } + } + else + { + unsigned long rx_cond = __m81_test (__real__ x); + + if (rx_cond & __M81_COND_INF) + { + /* Real part is infinite. */ + if (rx_cond & __M81_COND_NEG) + { + __real__ retval = __imag__ retval = 0.0; + if (ix_cond & __M81_COND_NEG) + __imag__ retval = -__imag__ retval; + } + else + { + __real__ retval = __real__ x; + __imag__ retval = __imag__ x - __imag__ x; + } + } + else + __real__ retval = __imag__ retval = __imag__ x - __imag__ x; + } + + return retval; +} +weak_alias (s(__cexp), s(cexp)) diff --git a/REORG.TODO/sysdeps/m68k/m680x0/fpu/s_cos.c b/REORG.TODO/sysdeps/m68k/m680x0/fpu/s_cos.c new file mode 100644 index 0000000000..2f1adf77f5 --- /dev/null +++ b/REORG.TODO/sysdeps/m68k/m680x0/fpu/s_cos.c @@ -0,0 +1,2 @@ +#define FUNC cos +#include <s_sin.c> diff --git a/REORG.TODO/sysdeps/m68k/m680x0/fpu/s_cosf.c b/REORG.TODO/sysdeps/m68k/m680x0/fpu/s_cosf.c new file mode 100644 index 0000000000..b6d6accb4a --- /dev/null +++ b/REORG.TODO/sysdeps/m68k/m680x0/fpu/s_cosf.c @@ -0,0 +1,2 @@ +#define FUNC cosf +#include <s_sinf.c> diff --git a/REORG.TODO/sysdeps/m68k/m680x0/fpu/s_cosl.c b/REORG.TODO/sysdeps/m68k/m680x0/fpu/s_cosl.c new file mode 100644 index 0000000000..eb61521fdd --- /dev/null +++ b/REORG.TODO/sysdeps/m68k/m680x0/fpu/s_cosl.c @@ -0,0 +1,2 @@ +#define FUNC cosl +#include <s_sinl.c> diff --git a/REORG.TODO/sysdeps/m68k/m680x0/fpu/s_csin_template.c b/REORG.TODO/sysdeps/m68k/m680x0/fpu/s_csin_template.c new file mode 100644 index 0000000000..be72208ef7 --- /dev/null +++ b/REORG.TODO/sysdeps/m68k/m680x0/fpu/s_csin_template.c @@ -0,0 +1,59 @@ +/* Complex sine function. m68k fpu version + Copyright (C) 1997-2017 Free Software Foundation, Inc. + This file is part of the GNU C Library. + Contributed by Andreas Schwab <schwab@issan.informatik.uni-dortmund.de>. + + The GNU C Library 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. + + 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 + Lesser General Public License for more details. + + You should have received a copy of the GNU Lesser General Public + License along with the GNU C Library. If not, see + <http://www.gnu.org/licenses/>. */ + +#include <complex.h> +#include <math.h> +#include "mathimpl.h" + +#define s(name) M_SUF (name) +#define m81(func) __m81_u(s(func)) + +CFLOAT +s(__csin) (CFLOAT x) +{ + CFLOAT retval; + unsigned long rx_cond = __m81_test (__real__ x); + + if ((rx_cond & (__M81_COND_INF|__M81_COND_NAN)) == 0) + { + /* Real part is finite. */ + FLOAT sin_rx, cos_rx; + + __asm ("fsincos%.x %2,%1:%0" : "=f" (sin_rx), "=f" (cos_rx) + : "f" (__real__ x)); + if (rx_cond & __M81_COND_ZERO) + __real__ retval = __real__ x; + else + __real__ retval = sin_rx * m81(__ieee754_cosh) (__imag__ x); + __imag__ retval = cos_rx * m81(__ieee754_sinh) (__imag__ x); + } + else + { + unsigned long ix_cond = __m81_test (__imag__ x); + + __real__ retval = __real__ x - __real__ x; + if (ix_cond & (__M81_COND_ZERO|__M81_COND_INF|__M81_COND_NAN)) + __imag__ retval = __imag__ x; + else + __imag__ retval = __real__ retval; + } + + return retval; +} +weak_alias (s(__csin), s(csin)) diff --git a/REORG.TODO/sysdeps/m68k/m680x0/fpu/s_csinh_template.c b/REORG.TODO/sysdeps/m68k/m680x0/fpu/s_csinh_template.c new file mode 100644 index 0000000000..f5180459e6 --- /dev/null +++ b/REORG.TODO/sysdeps/m68k/m680x0/fpu/s_csinh_template.c @@ -0,0 +1,62 @@ +/* Complex sine hyperbole function. m68k fpu version + Copyright (C) 1997-2017 Free Software Foundation, Inc. + This file is part of the GNU C Library. + Contributed by Andreas Schwab <schwab@issan.informatik.uni-dortmund.de>. + + The GNU C Library 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. + + 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 + Lesser General Public License for more details. + + You should have received a copy of the GNU Lesser General Public + License along with the GNU C Library. If not, see + <http://www.gnu.org/licenses/>. */ + +#include <complex.h> +#include <math.h> +#include "mathimpl.h" + +#define CONCATX(a,b) __CONCAT(a,b) +#define s(name) M_SUF (name) +#define m81(func) __m81_u(s(func)) + +CFLOAT +s(__csinh) (CFLOAT x) +{ + CFLOAT retval; + unsigned long ix_cond; + + ix_cond = __m81_test (__imag__ x); + + if ((ix_cond & (__M81_COND_INF|__M81_COND_NAN)) == 0) + { + /* Imaginary part is finite. */ + FLOAT sin_ix, cos_ix; + + __asm ("fsincos%.x %2,%1:%0" : "=f" (sin_ix), "=f" (cos_ix) + : "f" (__imag__ x)); + __real__ retval = cos_ix * m81(__ieee754_sinh) (__real__ x); + if (ix_cond & __M81_COND_ZERO) + __imag__ retval = __imag__ x; + else + __imag__ retval = sin_ix * m81(__ieee754_cosh) (__real__ x); + } + else + { + unsigned long rx_cond = __m81_test (__real__ x); + + __imag__ retval = __imag__ x - __imag__ x; + if (rx_cond & (__M81_COND_ZERO|__M81_COND_INF|__M81_COND_NAN)) + __real__ retval = __real__ x; + else + __real__ retval = __imag__ retval; + } + + return retval; +} +weak_alias (s(__csinh), s(csinh)) diff --git a/REORG.TODO/sysdeps/m68k/m680x0/fpu/s_expm1.c b/REORG.TODO/sysdeps/m68k/m680x0/fpu/s_expm1.c new file mode 100644 index 0000000000..c969a7f029 --- /dev/null +++ b/REORG.TODO/sysdeps/m68k/m680x0/fpu/s_expm1.c @@ -0,0 +1,41 @@ +/* Copyright (C) 2012-2017 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 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. + + 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 + Lesser General Public License for more details. + + You should have received a copy of the GNU Lesser General Public + License along with the GNU C Library. If not, see + <http://www.gnu.org/licenses/>. */ + +#include <math.h> +#include <errno.h> +#include "mathimpl.h" + +#ifndef FUNC +# define FUNC expm1 +#endif +#ifndef float_type +# define float_type double +#endif +#ifndef o_threshold +# define o_threshold 7.09782712893383973096e+02 +#endif + +#define CONCATX(a,b) __CONCAT(a,b) + +float_type +CONCATX(__,FUNC) (float_type x) +{ + if ((__m81_test (x) & __M81_COND_INF) == 0 && isgreater (x, o_threshold)) + __set_errno (ERANGE); + return __m81_u(CONCATX(__, FUNC)) (x); +} +weak_alias (CONCATX(__, FUNC), FUNC) diff --git a/REORG.TODO/sysdeps/m68k/m680x0/fpu/s_expm1f.c b/REORG.TODO/sysdeps/m68k/m680x0/fpu/s_expm1f.c new file mode 100644 index 0000000000..2cd8ffcfeb --- /dev/null +++ b/REORG.TODO/sysdeps/m68k/m680x0/fpu/s_expm1f.c @@ -0,0 +1,4 @@ +#define FUNC expm1f +#define float_type float +#define o_threshold 8.8722831726e+01F +#include <s_expm1.c> diff --git a/REORG.TODO/sysdeps/m68k/m680x0/fpu/s_expm1l.c b/REORG.TODO/sysdeps/m68k/m680x0/fpu/s_expm1l.c new file mode 100644 index 0000000000..a2c80860ba --- /dev/null +++ b/REORG.TODO/sysdeps/m68k/m680x0/fpu/s_expm1l.c @@ -0,0 +1,6 @@ +#define FUNC expm1l +#define float_type long double +/* 400c0000b17217f7d1cf79ab */ +#define o_threshold 1.13565234062941439487914863093465101e+04L +#include <s_expm1.c> +libm_hidden_def (__expm1l) diff --git a/REORG.TODO/sysdeps/m68k/m680x0/fpu/s_fabs.c b/REORG.TODO/sysdeps/m68k/m680x0/fpu/s_fabs.c new file mode 100644 index 0000000000..1f0631e2ff --- /dev/null +++ b/REORG.TODO/sysdeps/m68k/m680x0/fpu/s_fabs.c @@ -0,0 +1,2 @@ +#define FUNC fabs +#include <s_atan.c> diff --git a/REORG.TODO/sysdeps/m68k/m680x0/fpu/s_fabsf.c b/REORG.TODO/sysdeps/m68k/m680x0/fpu/s_fabsf.c new file mode 100644 index 0000000000..8f9421998a --- /dev/null +++ b/REORG.TODO/sysdeps/m68k/m680x0/fpu/s_fabsf.c @@ -0,0 +1,2 @@ +#define FUNC fabsf +#include <s_atanf.c> diff --git a/REORG.TODO/sysdeps/m68k/m680x0/fpu/s_fabsl.c b/REORG.TODO/sysdeps/m68k/m680x0/fpu/s_fabsl.c new file mode 100644 index 0000000000..8ac14d5b8c --- /dev/null +++ b/REORG.TODO/sysdeps/m68k/m680x0/fpu/s_fabsl.c @@ -0,0 +1,2 @@ +#define FUNC fabsl +#include <s_atanl.c> diff --git a/REORG.TODO/sysdeps/m68k/m680x0/fpu/s_finite.c b/REORG.TODO/sysdeps/m68k/m680x0/fpu/s_finite.c new file mode 100644 index 0000000000..dafbd5901d --- /dev/null +++ b/REORG.TODO/sysdeps/m68k/m680x0/fpu/s_finite.c @@ -0,0 +1,2 @@ +#define FUNC finite +#include <s_isinf.c> diff --git a/REORG.TODO/sysdeps/m68k/m680x0/fpu/s_finitef.c b/REORG.TODO/sysdeps/m68k/m680x0/fpu/s_finitef.c new file mode 100644 index 0000000000..b81342e5c6 --- /dev/null +++ b/REORG.TODO/sysdeps/m68k/m680x0/fpu/s_finitef.c @@ -0,0 +1,2 @@ +#define FUNC finitef +#include <s_isinff.c> diff --git a/REORG.TODO/sysdeps/m68k/m680x0/fpu/s_finitel.c b/REORG.TODO/sysdeps/m68k/m680x0/fpu/s_finitel.c new file mode 100644 index 0000000000..bd346a220c --- /dev/null +++ b/REORG.TODO/sysdeps/m68k/m680x0/fpu/s_finitel.c @@ -0,0 +1,2 @@ +#define FUNC finitel +#include <s_isinfl.c> diff --git a/REORG.TODO/sysdeps/m68k/m680x0/fpu/s_floor.c b/REORG.TODO/sysdeps/m68k/m680x0/fpu/s_floor.c new file mode 100644 index 0000000000..e1219c602a --- /dev/null +++ b/REORG.TODO/sysdeps/m68k/m680x0/fpu/s_floor.c @@ -0,0 +1,2 @@ +#define FUNC floor +#include <s_atan.c> diff --git a/REORG.TODO/sysdeps/m68k/m680x0/fpu/s_floorf.c b/REORG.TODO/sysdeps/m68k/m680x0/fpu/s_floorf.c new file mode 100644 index 0000000000..f4f9b9a1d8 --- /dev/null +++ b/REORG.TODO/sysdeps/m68k/m680x0/fpu/s_floorf.c @@ -0,0 +1,2 @@ +#define FUNC floorf +#include <s_atanf.c> diff --git a/REORG.TODO/sysdeps/m68k/m680x0/fpu/s_floorl.c b/REORG.TODO/sysdeps/m68k/m680x0/fpu/s_floorl.c new file mode 100644 index 0000000000..2c1ffd7d2c --- /dev/null +++ b/REORG.TODO/sysdeps/m68k/m680x0/fpu/s_floorl.c @@ -0,0 +1,2 @@ +#define FUNC floorl +#include <s_atanl.c> diff --git a/REORG.TODO/sysdeps/m68k/m680x0/fpu/s_fpclassifyl.c b/REORG.TODO/sysdeps/m68k/m680x0/fpu/s_fpclassifyl.c new file mode 100644 index 0000000000..248a2d3ca4 --- /dev/null +++ b/REORG.TODO/sysdeps/m68k/m680x0/fpu/s_fpclassifyl.c @@ -0,0 +1,43 @@ +/* Return classification value corresponding to argument. m68k version. + Copyright (C) 1997-2017 Free Software Foundation, Inc. + This file is part of the GNU C Library. + Contributed by Ulrich Drepper <drepper@cygnus.com>, 1997. + Fixed for m68k by Andreas Schwab <schwab@suse.de>. + + The GNU C Library 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. + + 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 + Lesser General Public License for more details. + + You should have received a copy of the GNU Lesser General Public + License along with the GNU C Library. If not, see + <http://www.gnu.org/licenses/>. */ + +#include <math.h> + +#include <math_private.h> + + +int +__fpclassifyl (long double x) +{ + u_int32_t ex, hx, lx; + int retval = FP_NORMAL; + + GET_LDOUBLE_WORDS (ex, hx, lx, x); + ex &= 0x7fff; + if ((ex | hx | lx) == 0) + retval = FP_ZERO; + else if (ex == 0 && (hx & 0x80000000) == 0) + retval = FP_SUBNORMAL; + else if (ex == 0x7fff) + retval = ((hx & 0x7fffffff) | lx) != 0 ? FP_NAN : FP_INFINITE; + + return retval; +} +libm_hidden_def (__fpclassifyl) diff --git a/REORG.TODO/sysdeps/m68k/m680x0/fpu/s_frexp.c b/REORG.TODO/sysdeps/m68k/m680x0/fpu/s_frexp.c new file mode 100644 index 0000000000..009b1e3223 --- /dev/null +++ b/REORG.TODO/sysdeps/m68k/m680x0/fpu/s_frexp.c @@ -0,0 +1,53 @@ +/* Copyright (C) 1996-2017 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 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. + + 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 + Lesser General Public License for more details. + + You should have received a copy of the GNU Lesser General Public + License along with the GNU C Library. If not, see + <http://www.gnu.org/licenses/>. */ + +#include <math.h> + +#ifndef FUNC +#define FUNC frexp +#endif +#ifndef float_type +#define float_type double +#endif + +#define __CONCATX(a,b) __CONCAT(a,b) + +float_type +__CONCATX(__,FUNC) (float_type value, int *expptr) +{ + float_type mantissa, exponent; + int iexponent; + unsigned long fpsr; + + __asm ("ftst%.x %1\n" + "fmove%.l %/fpsr, %0" + : "=dm" (fpsr) : "f" (value)); + if (fpsr & (7 << 24)) + { + /* Not finite or zero. */ + *expptr = 0; + return value; + } + __asm ("fgetexp%.x %1, %0" : "=f" (exponent) : "f" (value)); + iexponent = (int) exponent + 1; + *expptr = iexponent; + __asm ("fscale%.l %2, %0" + : "=f" (mantissa) + : "0" (value), "dmi" (-iexponent)); + return mantissa; +} +weak_alias (__CONCATX(__,FUNC), FUNC) diff --git a/REORG.TODO/sysdeps/m68k/m680x0/fpu/s_frexpf.c b/REORG.TODO/sysdeps/m68k/m680x0/fpu/s_frexpf.c new file mode 100644 index 0000000000..893b6ad3cf --- /dev/null +++ b/REORG.TODO/sysdeps/m68k/m680x0/fpu/s_frexpf.c @@ -0,0 +1,3 @@ +#define FUNC frexpf +#define float_type float +#include <s_frexp.c> diff --git a/REORG.TODO/sysdeps/m68k/m680x0/fpu/s_frexpl.c b/REORG.TODO/sysdeps/m68k/m680x0/fpu/s_frexpl.c new file mode 100644 index 0000000000..7d1d3964e7 --- /dev/null +++ b/REORG.TODO/sysdeps/m68k/m680x0/fpu/s_frexpl.c @@ -0,0 +1,58 @@ +/* Copyright (C) 2003-2017 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 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. + + 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 + Lesser General Public License for more details. + + You should have received a copy of the GNU Lesser General Public + License along with the GNU C Library. If not, see + <http://www.gnu.org/licenses/>. */ + +#include <math.h> + +long double +__frexpl (long double value, int *expptr) +{ + long double mantissa, exponent; + int iexponent; + unsigned long fpsr; + + __asm ("ftst%.x %1\n" + "fmove%.l %/fpsr, %0" + : "=dm" (fpsr) : "f" (value)); + if (fpsr & (7 << 24)) + { + /* Not finite or zero. */ + *expptr = 0; + return value; + } + __asm ("fgetexp%.x %1, %0" : "=f" (exponent) : "f" (value)); + iexponent = (int) exponent + 1; + *expptr = iexponent; + /* Unnormalized numbers must be handled specially, otherwise fscale + results in overflow. */ + if (iexponent <= -16384) + { + value *= 0x1p16383L; + iexponent += 16383; + } + else if (iexponent >= 16384) + { + value *= 0x1p-16383L; + iexponent -= 16383; + } + + __asm ("fscale%.l %2, %0" + : "=f" (mantissa) + : "0" (value), "dmi" (-iexponent)); + return mantissa; +} + +weak_alias (__frexpl, frexpl) diff --git a/REORG.TODO/sysdeps/m68k/m680x0/fpu/s_isinf.c b/REORG.TODO/sysdeps/m68k/m680x0/fpu/s_isinf.c new file mode 100644 index 0000000000..8764bb7e78 --- /dev/null +++ b/REORG.TODO/sysdeps/m68k/m680x0/fpu/s_isinf.c @@ -0,0 +1,37 @@ +/* Copyright (C) 1996-2017 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 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. + + 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 + Lesser General Public License for more details. + + You should have received a copy of the GNU Lesser General Public + License along with the GNU C Library. If not, see + <http://www.gnu.org/licenses/>. */ + +#include <math.h> + +#ifndef FUNC +#define FUNC isinf +#endif +#ifndef float_type +#define float_type double +#endif + +#define __CONCATX(a,b) __CONCAT(a,b) + +int +__CONCATX(__,FUNC) (float_type x) +{ + return __m81_u(__CONCATX(__,FUNC))(x); +} + +#define hidden_defx(a) hidden_def(a) +hidden_defx(__CONCATX(__,FUNC)) +weak_alias (__CONCATX(__,FUNC), FUNC) diff --git a/REORG.TODO/sysdeps/m68k/m680x0/fpu/s_isinff.c b/REORG.TODO/sysdeps/m68k/m680x0/fpu/s_isinff.c new file mode 100644 index 0000000000..ebf4f2ae96 --- /dev/null +++ b/REORG.TODO/sysdeps/m68k/m680x0/fpu/s_isinff.c @@ -0,0 +1,5 @@ +#ifndef FUNC +#define FUNC isinff +#endif +#define float_type float +#include <s_isinf.c> diff --git a/REORG.TODO/sysdeps/m68k/m680x0/fpu/s_isinfl.c b/REORG.TODO/sysdeps/m68k/m680x0/fpu/s_isinfl.c new file mode 100644 index 0000000000..963725ad77 --- /dev/null +++ b/REORG.TODO/sysdeps/m68k/m680x0/fpu/s_isinfl.c @@ -0,0 +1,5 @@ +#ifndef FUNC +#define FUNC isinfl +#endif +#define float_type long double +#include <s_isinf.c> diff --git a/REORG.TODO/sysdeps/m68k/m680x0/fpu/s_isnan.c b/REORG.TODO/sysdeps/m68k/m680x0/fpu/s_isnan.c new file mode 100644 index 0000000000..151d6dc6b4 --- /dev/null +++ b/REORG.TODO/sysdeps/m68k/m680x0/fpu/s_isnan.c @@ -0,0 +1,2 @@ +#define FUNC isnan +#include <s_isinf.c> diff --git a/REORG.TODO/sysdeps/m68k/m680x0/fpu/s_isnanf.c b/REORG.TODO/sysdeps/m68k/m680x0/fpu/s_isnanf.c new file mode 100644 index 0000000000..667bca7235 --- /dev/null +++ b/REORG.TODO/sysdeps/m68k/m680x0/fpu/s_isnanf.c @@ -0,0 +1,2 @@ +#define FUNC isnanf +#include <s_isinff.c> diff --git a/REORG.TODO/sysdeps/m68k/m680x0/fpu/s_isnanl.c b/REORG.TODO/sysdeps/m68k/m680x0/fpu/s_isnanl.c new file mode 100644 index 0000000000..bbacb64f11 --- /dev/null +++ b/REORG.TODO/sysdeps/m68k/m680x0/fpu/s_isnanl.c @@ -0,0 +1,2 @@ +#define FUNC isnanl +#include <s_isinfl.c> diff --git a/REORG.TODO/sysdeps/m68k/m680x0/fpu/s_llrint.c b/REORG.TODO/sysdeps/m68k/m680x0/fpu/s_llrint.c new file mode 100644 index 0000000000..333f680947 --- /dev/null +++ b/REORG.TODO/sysdeps/m68k/m680x0/fpu/s_llrint.c @@ -0,0 +1,75 @@ +/* Round argument to nearest integral value according to current rounding + direction. + Copyright (C) 1997-2017 Free Software Foundation, Inc. + This file is part of the GNU C Library. + Contributed by Andreas Schwab <schwab@issan.informatik.uni-dortmund.de> + + The GNU C Library 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. + + 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 + Lesser General Public License for more details. + + You should have received a copy of the GNU Lesser General Public + License along with the GNU C Library. If not, see + <http://www.gnu.org/licenses/>. */ + +#include <math.h> +#include <math_private.h> +#include "mathimpl.h" + +long long int +__llrint (double x) +{ + int32_t e; + u_int32_t h, l, s; + long long int result; + + x = __m81_u(__rint) (x); + + /* We could use __fixxfdi from libgcc, but here we can take advantage of + the known floating point format. */ + EXTRACT_WORDS (h, l, x); + + e = ((h >> 20) & 0x7ff) - 0x3ff; + if (e < 0) + return 0; + s = h; + h &= 0xfffff; + h |= 0x100000; + + if (e < 63) + { + if (e > 52) + { + h <<= e - 52; + h |= l >> (84 - e); + l <<= e - 52; + result = ((long long int) h << 32) | l; + } + else if (e > 20) + { + l >>= 52 - e; + l |= h << (e - 20); + h >>= 52 - e; + result = ((long long int) h << 32) | l; + } + else + result = h >> (20 - e); + if (s & 0x80000000) + result = -result; + } + else + /* The number is too large or not finite. The standard leaves it + undefined what to return when the number is too large to fit in a + `long long int'. */ + result = -1LL; + + return result; +} + +weak_alias (__llrint, llrint) diff --git a/REORG.TODO/sysdeps/m68k/m680x0/fpu/s_llrintf.c b/REORG.TODO/sysdeps/m68k/m680x0/fpu/s_llrintf.c new file mode 100644 index 0000000000..b4367a1494 --- /dev/null +++ b/REORG.TODO/sysdeps/m68k/m680x0/fpu/s_llrintf.c @@ -0,0 +1,65 @@ +/* Round argument to nearest integral value according to current rounding + direction. + Copyright (C) 1997-2017 Free Software Foundation, Inc. + This file is part of the GNU C Library. + Contributed by Andreas Schwab <schwab@issan.informatik.uni-dortmund.de> + + The GNU C Library 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. + + 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 + Lesser General Public License for more details. + + You should have received a copy of the GNU Lesser General Public + License along with the GNU C Library. If not, see + <http://www.gnu.org/licenses/>. */ + +#include <math.h> +#include <math_private.h> +#include "mathimpl.h" + +long long int +__llrintf (float x) +{ + int32_t e; + u_int32_t i, s; + long long int result; + + x = __m81_u(__rintf) (x); + + GET_FLOAT_WORD (i, x); + + e = ((i >> 23) & 0xff) - 0x7f; + if (e < 0) + return 0; + s = i; + i &= 0x7fffff; + i |= 0x800000; + + if (e < 63) + { + if (e > 55) + result = (long long int) (i << (e - 55)) << 32; + else if (e > 31) + result = (((long long int) (i >> (55 - e)) << 32) | (i << (e - 23))); + else if (e > 23) + result = i << (e - 23); + else + result = i >> (23 - e); + if (s & 0x80000000) + result = -result; + } + else + /* The number is too large or not finite. The standard leaves it + undefined what to return when the number is too large to fit in a + `long long int'. */ + result = -1LL; + + return result; +} + +weak_alias (__llrintf, llrintf) diff --git a/REORG.TODO/sysdeps/m68k/m680x0/fpu/s_llrintl.c b/REORG.TODO/sysdeps/m68k/m680x0/fpu/s_llrintl.c new file mode 100644 index 0000000000..17a7d7d7b0 --- /dev/null +++ b/REORG.TODO/sysdeps/m68k/m680x0/fpu/s_llrintl.c @@ -0,0 +1,64 @@ +/* Round argument to nearest integral value according to current rounding + direction. + Copyright (C) 1997-2017 Free Software Foundation, Inc. + This file is part of the GNU C Library. + Contributed by Andreas Schwab <schwab@issan.informatik.uni-dortmund.de> + + The GNU C Library 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. + + 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 + Lesser General Public License for more details. + + You should have received a copy of the GNU Lesser General Public + License along with the GNU C Library. If not, see + <http://www.gnu.org/licenses/>. */ + +#include <math.h> +#include <math_private.h> +#include "mathimpl.h" + +long long int +__llrintl (long double x) +{ + int32_t e, s; + u_int32_t h, l; + long long int result; + + x = __m81_u(__rintl) (x); + + GET_LDOUBLE_WORDS (e, h, l, x); + + s = e; + e = (e & 0x7fff) - 0x3fff; + if (e < 0) + return 0; + + if (e < 63) + { + if (e > 31) + { + l >>= 63 - e; + l |= h << (e - 31); + h >>= 63 - e; + result = ((long long int) h << 32) | l; + } + else + result = h >> (31 - e); + if (s & 0x8000) + result = -result; + } + else + /* The number is too large or not finite. The standard leaves it + undefined what to return when the number is too large to fit in a + `long long int'. */ + result = -1LL; + + return result; +} + +weak_alias (__llrintl, llrintl) diff --git a/REORG.TODO/sysdeps/m68k/m680x0/fpu/s_log1p.c b/REORG.TODO/sysdeps/m68k/m680x0/fpu/s_log1p.c new file mode 100644 index 0000000000..082618df18 --- /dev/null +++ b/REORG.TODO/sysdeps/m68k/m680x0/fpu/s_log1p.c @@ -0,0 +1,4 @@ +#define FUNC log1p +#undef weak_alias +#define weak_alias(a,b) +#include <s_atan.c> diff --git a/REORG.TODO/sysdeps/m68k/m680x0/fpu/s_log1pf.c b/REORG.TODO/sysdeps/m68k/m680x0/fpu/s_log1pf.c new file mode 100644 index 0000000000..480c39519f --- /dev/null +++ b/REORG.TODO/sysdeps/m68k/m680x0/fpu/s_log1pf.c @@ -0,0 +1,4 @@ +#define FUNC log1pf +#undef weak_alias +#define weak_alias(a,b) +#include <s_atanf.c> diff --git a/REORG.TODO/sysdeps/m68k/m680x0/fpu/s_log1pl.c b/REORG.TODO/sysdeps/m68k/m680x0/fpu/s_log1pl.c new file mode 100644 index 0000000000..a4f34a4f86 --- /dev/null +++ b/REORG.TODO/sysdeps/m68k/m680x0/fpu/s_log1pl.c @@ -0,0 +1,4 @@ +#define FUNC log1pl +#undef weak_alias +#define weak_alias(a,b) +#include <s_atanl.c> diff --git a/REORG.TODO/sysdeps/m68k/m680x0/fpu/s_logbl.c b/REORG.TODO/sysdeps/m68k/m680x0/fpu/s_logbl.c new file mode 100644 index 0000000000..2586622bf7 --- /dev/null +++ b/REORG.TODO/sysdeps/m68k/m680x0/fpu/s_logbl.c @@ -0,0 +1,49 @@ +/* s_logbl.c -- long double version of s_logb.c. + * Conversion to long double by Ulrich Drepper, + * Cygnus Support, drepper@cygnus.com. + */ + +/* + * ==================================================== + * Copyright (C) 1993 by Sun Microsystems, Inc. All rights reserved. + * + * Developed at SunPro, a Sun Microsystems, Inc. business. + * Permission to use, copy, modify, and distribute this + * software is freely granted, provided that this notice + * is preserved. + * ==================================================== + */ + +/* + * long double logbl(x) + * IEEE 754 logb. Included to pass IEEE test suite. Not recommend. + * Use ilogb instead. + */ + +#include <math.h> +#include <math_private.h> + +long double +__logbl (long double x) +{ + int32_t es, lx, ix; + + GET_LDOUBLE_WORDS (es, ix, lx, x); + es &= 0x7fff; /* exponent */ + if ((es | ix | lx) == 0) + return -1.0 / fabsl (x); + if (es == 0x7fff) + return x * x; + if (es == 0) /* IEEE 754 logb */ + { + /* POSIX specifies that denormal number is treated as + though it were normalized. */ + if (ix == 0) + es = -(__builtin_clz (lx) + 32); + else + es = -__builtin_clz (ix); + } + return (long double) (es - 16383); +} + +weak_alias (__logbl, logbl) diff --git a/REORG.TODO/sysdeps/m68k/m680x0/fpu/s_lrint.c b/REORG.TODO/sysdeps/m68k/m680x0/fpu/s_lrint.c new file mode 100644 index 0000000000..d25a548dcd --- /dev/null +++ b/REORG.TODO/sysdeps/m68k/m680x0/fpu/s_lrint.c @@ -0,0 +1,37 @@ +/* Round argument to nearest integral value according to current rounding + direction. + Copyright (C) 1997-2017 Free Software Foundation, Inc. + This file is part of the GNU C Library. + Contributed by Andreas Schwab <schwab@issan.informatik.uni-dortmund.de> + + The GNU C Library 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. + + 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 + Lesser General Public License for more details. + + You should have received a copy of the GNU Lesser General Public + License along with the GNU C Library. If not, see + <http://www.gnu.org/licenses/>. */ + +#include <math.h> + +#ifndef suffix +#define suffix /*empty*/ +#endif +#ifndef float_type +#define float_type double +#endif + +#define CONCATX(a,b) __CONCAT(a,b) + +long int +CONCATX(__lrint,suffix) (float_type x) +{ + return __m81_u(CONCATX(__lrint,suffix)) (x); +} +weak_alias (CONCATX(__lrint,suffix), CONCATX(lrint,suffix)) diff --git a/REORG.TODO/sysdeps/m68k/m680x0/fpu/s_lrintf.c b/REORG.TODO/sysdeps/m68k/m680x0/fpu/s_lrintf.c new file mode 100644 index 0000000000..44924cb82f --- /dev/null +++ b/REORG.TODO/sysdeps/m68k/m680x0/fpu/s_lrintf.c @@ -0,0 +1,3 @@ +#define suffix f +#define float_type float +#include <s_lrint.c> diff --git a/REORG.TODO/sysdeps/m68k/m680x0/fpu/s_lrintl.c b/REORG.TODO/sysdeps/m68k/m680x0/fpu/s_lrintl.c new file mode 100644 index 0000000000..cd0bd23b8a --- /dev/null +++ b/REORG.TODO/sysdeps/m68k/m680x0/fpu/s_lrintl.c @@ -0,0 +1,3 @@ +#define suffix l +#define float_type long double +#include <s_lrint.c> diff --git a/REORG.TODO/sysdeps/m68k/m680x0/fpu/s_modf.c b/REORG.TODO/sysdeps/m68k/m680x0/fpu/s_modf.c new file mode 100644 index 0000000000..3c1a5a2749 --- /dev/null +++ b/REORG.TODO/sysdeps/m68k/m680x0/fpu/s_modf.c @@ -0,0 +1,53 @@ +/* Copyright (C) 1996-2017 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 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. + + 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 + Lesser General Public License for more details. + + You should have received a copy of the GNU Lesser General Public + License along with the GNU C Library. If not, see + <http://www.gnu.org/licenses/>. */ + +#include <math.h> +#include "mathimpl.h" + +#ifndef SUFF +#define SUFF +#endif +#ifndef float_type +#define float_type double +#endif + +#define CONCATX(a,b) __CONCAT(a,b) +#define s(name) CONCATX(name,SUFF) +#define m81(func) __m81_u(s(func)) + +float_type +s(__modf) (float_type x, float_type *iptr) +{ + float_type x_int, result; + unsigned long x_cond; + + __asm ("fintrz%.x %1, %0" : "=f" (x_int) : "f" (x)); + *iptr = x_int; + x_cond = __m81_test (x); + if (x_cond & __M81_COND_INF) + { + result = 0; + if (x_cond & __M81_COND_NEG) + result = -result; + } + else if (x_cond & __M81_COND_ZERO) + result = x; + else + result = x - x_int; + return result; +} +weak_alias (s(__modf), s(modf)) diff --git a/REORG.TODO/sysdeps/m68k/m680x0/fpu/s_modff.c b/REORG.TODO/sysdeps/m68k/m680x0/fpu/s_modff.c new file mode 100644 index 0000000000..0c44d7c0ab --- /dev/null +++ b/REORG.TODO/sysdeps/m68k/m680x0/fpu/s_modff.c @@ -0,0 +1,3 @@ +#define SUFF f +#define float_type float +#include <s_modf.c> diff --git a/REORG.TODO/sysdeps/m68k/m680x0/fpu/s_modfl.c b/REORG.TODO/sysdeps/m68k/m680x0/fpu/s_modfl.c new file mode 100644 index 0000000000..c7075b3ba9 --- /dev/null +++ b/REORG.TODO/sysdeps/m68k/m680x0/fpu/s_modfl.c @@ -0,0 +1,3 @@ +#define SUFF l +#define float_type long double +#include <s_modf.c> diff --git a/REORG.TODO/sysdeps/m68k/m680x0/fpu/s_nearbyint.c b/REORG.TODO/sysdeps/m68k/m680x0/fpu/s_nearbyint.c new file mode 100644 index 0000000000..b87f5e2160 --- /dev/null +++ b/REORG.TODO/sysdeps/m68k/m680x0/fpu/s_nearbyint.c @@ -0,0 +1,2 @@ +#define FUNC nearbyint +#include <s_atan.c> diff --git a/REORG.TODO/sysdeps/m68k/m680x0/fpu/s_nearbyintf.c b/REORG.TODO/sysdeps/m68k/m680x0/fpu/s_nearbyintf.c new file mode 100644 index 0000000000..70d08ab44c --- /dev/null +++ b/REORG.TODO/sysdeps/m68k/m680x0/fpu/s_nearbyintf.c @@ -0,0 +1,2 @@ +#define FUNC nearbyintf +#include <s_atanf.c> diff --git a/REORG.TODO/sysdeps/m68k/m680x0/fpu/s_nearbyintl.c b/REORG.TODO/sysdeps/m68k/m680x0/fpu/s_nearbyintl.c new file mode 100644 index 0000000000..230cd7784c --- /dev/null +++ b/REORG.TODO/sysdeps/m68k/m680x0/fpu/s_nearbyintl.c @@ -0,0 +1,2 @@ +#define FUNC nearbyintl +#include <s_atanl.c> diff --git a/REORG.TODO/sysdeps/m68k/m680x0/fpu/s_nextafterl.c b/REORG.TODO/sysdeps/m68k/m680x0/fpu/s_nextafterl.c new file mode 100644 index 0000000000..c46c0e76ec --- /dev/null +++ b/REORG.TODO/sysdeps/m68k/m680x0/fpu/s_nextafterl.c @@ -0,0 +1,108 @@ +/* s_nextafterl.c -- long double version of s_nextafter.c. + * Conversion to long double by Ulrich Drepper, + * Cygnus Support, drepper@cygnus.com. + * Fixed for m68k by Andreas Schwab <schwab@suse.de>. + */ + +/* + * ==================================================== + * Copyright (C) 1993 by Sun Microsystems, Inc. All rights reserved. + * + * Developed at SunPro, a Sun Microsystems, Inc. business. + * Permission to use, copy, modify, and distribute this + * software is freely granted, provided that this notice + * is preserved. + * ==================================================== + */ + +#if defined(LIBM_SCCS) && !defined(lint) +static char rcsid[] = "$NetBSD: $"; +#endif + +/* IEEE functions + * nextafterl(x,y) + * return the next machine floating-point number of x in the + * direction toward y. + * Special cases: + */ + +#include <errno.h> +#include <math.h> +#include <math_private.h> + +long double __nextafterl(long double x, long double y) +{ + int32_t ix,iy,esx,esy; + u_int32_t hx,hy,lx,ly; + + GET_LDOUBLE_WORDS(esx,hx,lx,x); + GET_LDOUBLE_WORDS(esy,hy,ly,y); + ix = esx&0x7fff; /* |x| */ + iy = esy&0x7fff; /* |y| */ + + if(((ix==0x7fff)&&((hx&0x7fffffff)|lx)!=0) || /* x is nan */ + ((iy==0x7fff)&&((hy&0x7fffffff)|ly)!=0)) /* y is nan */ + return x+y; + if(x==y) return y; /* x=y, return y */ + if((ix|hx|lx)==0) { /* x == 0 */ + SET_LDOUBLE_WORDS(x,esy&0x8000,0,1);/* return +-minsubnormal */ + y = x*x; + math_force_eval (y); /* raise underflow flag */ + return x; + } + if(esx>=0) { /* x > 0 */ + if(esx>esy||((esx==esy) && (hx>hy||((hx==hy)&&(lx>ly))))) { + /* x > y, x -= ulp */ + if(lx==0) { + if (ix != 0 && hx == 0x80000000) hx = 0; + if (hx==0) esx -= 1; + hx -= 1; + } + lx -= 1; + } else { /* x < y, x += ulp */ + lx += 1; + if(lx==0) { + hx += 1; + if (hx==0) { + hx = 0x80000000; + esx += 1; + } + } + } + } else { /* x < 0 */ + if(esy>=0||esx>esy||((esx==esy) && (hx>hy||((hx==hy)&&(lx>ly))))){ + /* x < y, x -= ulp */ + if(lx==0) { + if (ix != 0 && hx == 0x80000000) hx = 0; + if (hx==0) esx -= 1; + hx -= 1; + } + lx -= 1; + } else { /* x > y, x += ulp */ + lx += 1; + if(lx==0) { + hx += 1; + if (hx==0) { + hx = 0x80000000; + esx += 1; + } + } + } + } + esy = esx&0x7fff; + if(esy==0x7fff) { + long double u = x + x; /* overflow */ + math_force_eval (u); + __set_errno (ERANGE); + } + if(esy==0 && (hx & 0x80000000) == 0) { /* underflow */ + y = x*x; + math_force_eval (y); /* raise underflow flag */ + __set_errno (ERANGE); + } + SET_LDOUBLE_WORDS(x,esx,hx,lx); + return x; +} +weak_alias (__nextafterl, nextafterl) +strong_alias (__nextafterl, __nexttowardl) +weak_alias (__nextafterl, nexttowardl) diff --git a/REORG.TODO/sysdeps/m68k/m680x0/fpu/s_remquo.c b/REORG.TODO/sysdeps/m68k/m680x0/fpu/s_remquo.c new file mode 100644 index 0000000000..554751ca87 --- /dev/null +++ b/REORG.TODO/sysdeps/m68k/m680x0/fpu/s_remquo.c @@ -0,0 +1,46 @@ +/* Compute remainder and a congruent to the quotient. m68k fpu version + Copyright (C) 1997-2017 Free Software Foundation, Inc. + This file is part of the GNU C Library. + Contributed by Andreas Schwab <schwab@issan.informatik.uni-dortmund.de> + + The GNU C Library 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. + + 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 + Lesser General Public License for more details. + + You should have received a copy of the GNU Lesser General Public + License along with the GNU C Library. If not, see + <http://www.gnu.org/licenses/>. */ + +#include <math.h> + +#ifndef SUFF +#define SUFF +#endif +#ifndef float_type +#define float_type double +#endif + +#define CONCATX(a,b) __CONCAT(a,b) +#define s(name) CONCATX(name,SUFF) + +float_type +s(__remquo) (float_type x, float_type y, int *quo) +{ + float_type result; + int cquo, fpsr; + + __asm ("frem%.x %2,%0\n\tfmove%.l %/fpsr,%1" + : "=f" (result), "=dm" (fpsr) : "f" (y), "0" (x)); + cquo = (fpsr >> 16) & 0x7f; + if (fpsr & (1 << 23)) + cquo = -cquo; + *quo = cquo; + return result; +} +weak_alias (s(__remquo), s(remquo)) diff --git a/REORG.TODO/sysdeps/m68k/m680x0/fpu/s_remquof.c b/REORG.TODO/sysdeps/m68k/m680x0/fpu/s_remquof.c new file mode 100644 index 0000000000..8a292fc26c --- /dev/null +++ b/REORG.TODO/sysdeps/m68k/m680x0/fpu/s_remquof.c @@ -0,0 +1,3 @@ +#define SUFF f +#define float_type float +#include <s_remquo.c> diff --git a/REORG.TODO/sysdeps/m68k/m680x0/fpu/s_remquol.c b/REORG.TODO/sysdeps/m68k/m680x0/fpu/s_remquol.c new file mode 100644 index 0000000000..d236cfd1f9 --- /dev/null +++ b/REORG.TODO/sysdeps/m68k/m680x0/fpu/s_remquol.c @@ -0,0 +1,3 @@ +#define SUFF l +#define float_type long double +#include <s_remquo.c> diff --git a/REORG.TODO/sysdeps/m68k/m680x0/fpu/s_rint.c b/REORG.TODO/sysdeps/m68k/m680x0/fpu/s_rint.c new file mode 100644 index 0000000000..f0f18c7346 --- /dev/null +++ b/REORG.TODO/sysdeps/m68k/m680x0/fpu/s_rint.c @@ -0,0 +1,2 @@ +#define FUNC rint +#include <s_atan.c> diff --git a/REORG.TODO/sysdeps/m68k/m680x0/fpu/s_rintf.c b/REORG.TODO/sysdeps/m68k/m680x0/fpu/s_rintf.c new file mode 100644 index 0000000000..4e00cab0fb --- /dev/null +++ b/REORG.TODO/sysdeps/m68k/m680x0/fpu/s_rintf.c @@ -0,0 +1,2 @@ +#define FUNC rintf +#include <s_atanf.c> diff --git a/REORG.TODO/sysdeps/m68k/m680x0/fpu/s_rintl.c b/REORG.TODO/sysdeps/m68k/m680x0/fpu/s_rintl.c new file mode 100644 index 0000000000..305667b3a1 --- /dev/null +++ b/REORG.TODO/sysdeps/m68k/m680x0/fpu/s_rintl.c @@ -0,0 +1,2 @@ +#define FUNC rintl +#include <s_atanl.c> diff --git a/REORG.TODO/sysdeps/m68k/m680x0/fpu/s_scalbln.c b/REORG.TODO/sysdeps/m68k/m680x0/fpu/s_scalbln.c new file mode 100644 index 0000000000..1009713fbc --- /dev/null +++ b/REORG.TODO/sysdeps/m68k/m680x0/fpu/s_scalbln.c @@ -0,0 +1,2 @@ +/* Nothing to do. This function is the same as scalbn. So we define an + alias. */ diff --git a/REORG.TODO/sysdeps/m68k/m680x0/fpu/s_scalblnf.c b/REORG.TODO/sysdeps/m68k/m680x0/fpu/s_scalblnf.c new file mode 100644 index 0000000000..5e558c3540 --- /dev/null +++ b/REORG.TODO/sysdeps/m68k/m680x0/fpu/s_scalblnf.c @@ -0,0 +1,2 @@ +/* Nothing to do. This function is the same as scalbnf. So we define an + alias. */ diff --git a/REORG.TODO/sysdeps/m68k/m680x0/fpu/s_scalblnl.c b/REORG.TODO/sysdeps/m68k/m680x0/fpu/s_scalblnl.c new file mode 100644 index 0000000000..cda2ec11c8 --- /dev/null +++ b/REORG.TODO/sysdeps/m68k/m680x0/fpu/s_scalblnl.c @@ -0,0 +1,2 @@ +/* Nothing to do. This function is the same as scalbnl. So we define an + alias. */ diff --git a/REORG.TODO/sysdeps/m68k/m680x0/fpu/s_scalbn.c b/REORG.TODO/sysdeps/m68k/m680x0/fpu/s_scalbn.c new file mode 100644 index 0000000000..0a1ffdb734 --- /dev/null +++ b/REORG.TODO/sysdeps/m68k/m680x0/fpu/s_scalbn.c @@ -0,0 +1,52 @@ +/* Copyright (C) 1996-2017 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 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. + + 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 + Lesser General Public License for more details. + + You should have received a copy of the GNU Lesser General Public + License along with the GNU C Library. If not, see + <http://www.gnu.org/licenses/>. */ + +#define scalbln __no_scalbln_decl +#define scalblnf __no_scalblnf_decl +#define scalblnl __no_scalblnl_decl +#define __scalbln __no__scalbln_decl +#define __scalblnf __no__scalblnf_decl +#define __scalblnl __no__scalblnl_decl +#include <math.h> +#undef scalbln +#undef scalblnf +#undef scalblnl +#undef __scalbln +#undef __scalblnf +#undef __scalblnl + +#ifndef suffix +#define suffix /*empty*/ +#endif +#ifndef float_type +#define float_type double +#endif + +#define __CONCATX(a,b) __CONCAT(a,b) + +float_type +__CONCATX(__scalbn,suffix) (float_type x, int exp) +{ + return __m81_u(__CONCATX(__scalbn,suffix))(x, exp); +} +strong_alias (__CONCATX(__scalbn,suffix), __CONCATX(__scalbln,suffix)) + +#include <shlib-compat.h> +#if SHLIB_COMPAT (libc, GLIBC_2_1, GLIBC_2_20) +compat_symbol (libc, __CONCATX(__scalbn,suffix), __CONCATX(scalbln,suffix), + GLIBC_2_1); +#endif diff --git a/REORG.TODO/sysdeps/m68k/m680x0/fpu/s_scalbnf.c b/REORG.TODO/sysdeps/m68k/m680x0/fpu/s_scalbnf.c new file mode 100644 index 0000000000..547971836a --- /dev/null +++ b/REORG.TODO/sysdeps/m68k/m680x0/fpu/s_scalbnf.c @@ -0,0 +1,3 @@ +#define suffix f +#define float_type float +#include <s_scalbn.c> diff --git a/REORG.TODO/sysdeps/m68k/m680x0/fpu/s_scalbnl.c b/REORG.TODO/sysdeps/m68k/m680x0/fpu/s_scalbnl.c new file mode 100644 index 0000000000..874bafb04b --- /dev/null +++ b/REORG.TODO/sysdeps/m68k/m680x0/fpu/s_scalbnl.c @@ -0,0 +1,3 @@ +#define suffix l +#define float_type long double +#include <s_scalbn.c> diff --git a/REORG.TODO/sysdeps/m68k/m680x0/fpu/s_significand.c b/REORG.TODO/sysdeps/m68k/m680x0/fpu/s_significand.c new file mode 100644 index 0000000000..34d4ea3d14 --- /dev/null +++ b/REORG.TODO/sysdeps/m68k/m680x0/fpu/s_significand.c @@ -0,0 +1,2 @@ +#define FUNC significand +#include <s_atan.c> diff --git a/REORG.TODO/sysdeps/m68k/m680x0/fpu/s_significandf.c b/REORG.TODO/sysdeps/m68k/m680x0/fpu/s_significandf.c new file mode 100644 index 0000000000..4e769ca317 --- /dev/null +++ b/REORG.TODO/sysdeps/m68k/m680x0/fpu/s_significandf.c @@ -0,0 +1,2 @@ +#define FUNC significandf +#include <s_atanf.c> diff --git a/REORG.TODO/sysdeps/m68k/m680x0/fpu/s_significandl.c b/REORG.TODO/sysdeps/m68k/m680x0/fpu/s_significandl.c new file mode 100644 index 0000000000..8c6fc7e610 --- /dev/null +++ b/REORG.TODO/sysdeps/m68k/m680x0/fpu/s_significandl.c @@ -0,0 +1,2 @@ +#define FUNC significandl +#include <s_atanl.c> diff --git a/REORG.TODO/sysdeps/m68k/m680x0/fpu/s_sin.c b/REORG.TODO/sysdeps/m68k/m680x0/fpu/s_sin.c new file mode 100644 index 0000000000..953fb445cf --- /dev/null +++ b/REORG.TODO/sysdeps/m68k/m680x0/fpu/s_sin.c @@ -0,0 +1,38 @@ +/* Copyright (C) 2012-2017 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 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. + + 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 + Lesser General Public License for more details. + + You should have received a copy of the GNU Lesser General Public + License along with the GNU C Library. If not, see + <http://www.gnu.org/licenses/>. */ + +#include <math.h> +#include <errno.h> +#include "mathimpl.h" + +#ifndef FUNC +# define FUNC sin +#endif +#ifndef float_type +# define float_type double +#endif + +#define CONCATX(a,b) __CONCAT(a,b) + +float_type +CONCATX(__,FUNC) (float_type x) +{ + if (__m81_test (x) & __M81_COND_INF) + __set_errno (EDOM); + return __m81_u(CONCATX(__, FUNC)) (x); +} +weak_alias (CONCATX(__, FUNC), FUNC) diff --git a/REORG.TODO/sysdeps/m68k/m680x0/fpu/s_sincos.c b/REORG.TODO/sysdeps/m68k/m680x0/fpu/s_sincos.c new file mode 100644 index 0000000000..ca93d93333 --- /dev/null +++ b/REORG.TODO/sysdeps/m68k/m680x0/fpu/s_sincos.c @@ -0,0 +1,34 @@ +/* Copyright (C) 1997-2017 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 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. + + 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 + Lesser General Public License for more details. + + You should have received a copy of the GNU Lesser General Public + License along with the GNU C Library. If not, see + <http://www.gnu.org/licenses/>. */ + +#include <math.h> + +#ifndef FUNC +#define FUNC sincos +#endif +#ifndef float_type +#define float_type double +#endif + +#define CONCATX(a,b) __CONCAT(a,b) + +void +CONCATX(__,FUNC) (float_type x, float_type *sinx, float_type *cosx) +{ + __m81_u(CONCATX(__,FUNC))(x, sinx, cosx); +} +weak_alias (CONCATX(__,FUNC), FUNC) diff --git a/REORG.TODO/sysdeps/m68k/m680x0/fpu/s_sincosf.c b/REORG.TODO/sysdeps/m68k/m680x0/fpu/s_sincosf.c new file mode 100644 index 0000000000..7ee2ec6600 --- /dev/null +++ b/REORG.TODO/sysdeps/m68k/m680x0/fpu/s_sincosf.c @@ -0,0 +1,3 @@ +#define FUNC sincosf +#define float_type float +#include <s_sincos.c> diff --git a/REORG.TODO/sysdeps/m68k/m680x0/fpu/s_sincosl.c b/REORG.TODO/sysdeps/m68k/m680x0/fpu/s_sincosl.c new file mode 100644 index 0000000000..f998cc0977 --- /dev/null +++ b/REORG.TODO/sysdeps/m68k/m680x0/fpu/s_sincosl.c @@ -0,0 +1,3 @@ +#define FUNC sincosl +#define float_type long double +#include <s_sincos.c> diff --git a/REORG.TODO/sysdeps/m68k/m680x0/fpu/s_sinf.c b/REORG.TODO/sysdeps/m68k/m680x0/fpu/s_sinf.c new file mode 100644 index 0000000000..304f4f76f6 --- /dev/null +++ b/REORG.TODO/sysdeps/m68k/m680x0/fpu/s_sinf.c @@ -0,0 +1,5 @@ +#ifndef FUNC +# define FUNC sinf +#endif +#define float_type float +#include <s_sin.c> diff --git a/REORG.TODO/sysdeps/m68k/m680x0/fpu/s_sinl.c b/REORG.TODO/sysdeps/m68k/m680x0/fpu/s_sinl.c new file mode 100644 index 0000000000..31dac0bee7 --- /dev/null +++ b/REORG.TODO/sysdeps/m68k/m680x0/fpu/s_sinl.c @@ -0,0 +1,5 @@ +#ifndef FUNC +# define FUNC sinl +#endif +#define float_type long double +#include <s_sin.c> diff --git a/REORG.TODO/sysdeps/m68k/m680x0/fpu/s_tan.c b/REORG.TODO/sysdeps/m68k/m680x0/fpu/s_tan.c new file mode 100644 index 0000000000..dfbbbb2270 --- /dev/null +++ b/REORG.TODO/sysdeps/m68k/m680x0/fpu/s_tan.c @@ -0,0 +1,2 @@ +#define FUNC tan +#include <s_sin.c> diff --git a/REORG.TODO/sysdeps/m68k/m680x0/fpu/s_tanf.c b/REORG.TODO/sysdeps/m68k/m680x0/fpu/s_tanf.c new file mode 100644 index 0000000000..9ae66ddbf3 --- /dev/null +++ b/REORG.TODO/sysdeps/m68k/m680x0/fpu/s_tanf.c @@ -0,0 +1,2 @@ +#define FUNC tanf +#include <s_sinf.c> diff --git a/REORG.TODO/sysdeps/m68k/m680x0/fpu/s_tanh.c b/REORG.TODO/sysdeps/m68k/m680x0/fpu/s_tanh.c new file mode 100644 index 0000000000..ac2e7dbb79 --- /dev/null +++ b/REORG.TODO/sysdeps/m68k/m680x0/fpu/s_tanh.c @@ -0,0 +1,2 @@ +#define FUNC tanh +#include <s_atan.c> diff --git a/REORG.TODO/sysdeps/m68k/m680x0/fpu/s_tanhf.c b/REORG.TODO/sysdeps/m68k/m680x0/fpu/s_tanhf.c new file mode 100644 index 0000000000..1addaae4ff --- /dev/null +++ b/REORG.TODO/sysdeps/m68k/m680x0/fpu/s_tanhf.c @@ -0,0 +1,2 @@ +#define FUNC tanhf +#include <s_atanf.c> diff --git a/REORG.TODO/sysdeps/m68k/m680x0/fpu/s_tanhl.c b/REORG.TODO/sysdeps/m68k/m680x0/fpu/s_tanhl.c new file mode 100644 index 0000000000..6e997911f2 --- /dev/null +++ b/REORG.TODO/sysdeps/m68k/m680x0/fpu/s_tanhl.c @@ -0,0 +1,2 @@ +#define FUNC tanhl +#include <s_atanl.c> diff --git a/REORG.TODO/sysdeps/m68k/m680x0/fpu/s_tanl.c b/REORG.TODO/sysdeps/m68k/m680x0/fpu/s_tanl.c new file mode 100644 index 0000000000..27daf8f29b --- /dev/null +++ b/REORG.TODO/sysdeps/m68k/m680x0/fpu/s_tanl.c @@ -0,0 +1,2 @@ +#define FUNC tanl +#include <s_sinl.c> diff --git a/REORG.TODO/sysdeps/m68k/m680x0/fpu/s_trunc.c b/REORG.TODO/sysdeps/m68k/m680x0/fpu/s_trunc.c new file mode 100644 index 0000000000..96f29a776c --- /dev/null +++ b/REORG.TODO/sysdeps/m68k/m680x0/fpu/s_trunc.c @@ -0,0 +1,2 @@ +#define FUNC trunc +#include <s_atan.c> diff --git a/REORG.TODO/sysdeps/m68k/m680x0/fpu/s_truncf.c b/REORG.TODO/sysdeps/m68k/m680x0/fpu/s_truncf.c new file mode 100644 index 0000000000..44dca748ca --- /dev/null +++ b/REORG.TODO/sysdeps/m68k/m680x0/fpu/s_truncf.c @@ -0,0 +1,2 @@ +#define FUNC truncf +#include <s_atanf.c> diff --git a/REORG.TODO/sysdeps/m68k/m680x0/fpu/s_truncl.c b/REORG.TODO/sysdeps/m68k/m680x0/fpu/s_truncl.c new file mode 100644 index 0000000000..8d35777aed --- /dev/null +++ b/REORG.TODO/sysdeps/m68k/m680x0/fpu/s_truncl.c @@ -0,0 +1,2 @@ +#define FUNC truncl +#include <s_atanl.c> diff --git a/REORG.TODO/sysdeps/m68k/m680x0/fpu/sincos32.c b/REORG.TODO/sysdeps/m68k/m680x0/fpu/sincos32.c new file mode 100644 index 0000000000..1cc8931700 --- /dev/null +++ b/REORG.TODO/sysdeps/m68k/m680x0/fpu/sincos32.c @@ -0,0 +1 @@ +/* Not needed. */ diff --git a/REORG.TODO/sysdeps/m68k/m680x0/fpu/sincostab.c b/REORG.TODO/sysdeps/m68k/m680x0/fpu/sincostab.c new file mode 100644 index 0000000000..1cc8931700 --- /dev/null +++ b/REORG.TODO/sysdeps/m68k/m680x0/fpu/sincostab.c @@ -0,0 +1 @@ +/* Not needed. */ diff --git a/REORG.TODO/sysdeps/m68k/m680x0/fpu/slowexp.c b/REORG.TODO/sysdeps/m68k/m680x0/fpu/slowexp.c new file mode 100644 index 0000000000..1cc8931700 --- /dev/null +++ b/REORG.TODO/sysdeps/m68k/m680x0/fpu/slowexp.c @@ -0,0 +1 @@ +/* Not needed. */ diff --git a/REORG.TODO/sysdeps/m68k/m680x0/fpu/slowpow.c b/REORG.TODO/sysdeps/m68k/m680x0/fpu/slowpow.c new file mode 100644 index 0000000000..1cc8931700 --- /dev/null +++ b/REORG.TODO/sysdeps/m68k/m680x0/fpu/slowpow.c @@ -0,0 +1 @@ +/* Not needed. */ diff --git a/REORG.TODO/sysdeps/m68k/m680x0/fpu/t_exp.c b/REORG.TODO/sysdeps/m68k/m680x0/fpu/t_exp.c new file mode 100644 index 0000000000..fd37963b05 --- /dev/null +++ b/REORG.TODO/sysdeps/m68k/m680x0/fpu/t_exp.c @@ -0,0 +1 @@ +/* Empty. Not needed. */ diff --git a/REORG.TODO/sysdeps/m68k/m680x0/lshift.S b/REORG.TODO/sysdeps/m68k/m680x0/lshift.S new file mode 100644 index 0000000000..16b1b2f210 --- /dev/null +++ b/REORG.TODO/sysdeps/m68k/m680x0/lshift.S @@ -0,0 +1,168 @@ +/* mc68020 __mpn_lshift -- Shift left a low-level natural-number integer. + +Copyright (C) 1996-2017 Free Software Foundation, Inc. + +This file is part of the GNU MP Library. + +The GNU MP Library 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. + +The GNU MP 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 Lesser General Public +License for more details. + +You should have received a copy of the GNU Lesser General Public License +along with the GNU MP Library. If not, see <http://www.gnu.org/licenses/>. */ + +/* + INPUT PARAMETERS + res_ptr (sp + 4) + s_ptr (sp + 8) + s_size (sp + 16) + cnt (sp + 12) +*/ + +#include "sysdep.h" +#include "asm-syntax.h" + +#define res_ptr a1 +#define s_ptr a0 +#define s_size d6 +#define cnt d4 + + TEXT +ENTRY(__mpn_lshift) + +/* Save used registers on the stack. */ + moveml R(d2)-R(d6)/R(a2),MEM_PREDEC(sp) + cfi_adjust_cfa_offset (6*4) + cfi_rel_offset (R(d2), 0) + cfi_rel_offset (R(d3), 4) + cfi_rel_offset (R(d4), 8) + cfi_rel_offset (R(d5), 12) + cfi_rel_offset (R(d6), 16) + cfi_rel_offset (R(a2), 20) + +/* Copy the arguments to registers. */ + movel MEM_DISP(sp,28),R(res_ptr) + movel MEM_DISP(sp,32),R(s_ptr) + movel MEM_DISP(sp,36),R(s_size) + movel MEM_DISP(sp,40),R(cnt) + + moveql #1,R(d5) + cmpl R(d5),R(cnt) + bne L(Lnormal) + cmpl R(s_ptr),R(res_ptr) + bls L(Lspecial) /* jump if s_ptr >= res_ptr */ +#if (defined (__mc68020__) || defined (__NeXT__) || defined(mc68020)) + lea MEM_INDX1(s_ptr,s_size,l,4),R(a2) +#else /* not mc68020 */ + movel R(s_size),R(d0) + asll #2,R(d0) + lea MEM_INDX(s_ptr,d0,l),R(a2) +#endif + cmpl R(res_ptr),R(a2) + bls L(Lspecial) /* jump if res_ptr >= s_ptr + s_size */ + +L(Lnormal:) + moveql #32,R(d5) + subl R(cnt),R(d5) + +#if (defined (__mc68020__) || defined (__NeXT__) || defined(mc68020)) + lea MEM_INDX1(s_ptr,s_size,l,4),R(s_ptr) + lea MEM_INDX1(res_ptr,s_size,l,4),R(res_ptr) +#else /* not mc68000 */ + movel R(s_size),R(d0) + asll #2,R(d0) + addl R(s_size),R(s_ptr) + addl R(s_size),R(res_ptr) +#endif + movel MEM_PREDEC(s_ptr),R(d2) + movel R(d2),R(d0) + lsrl R(d5),R(d0) /* compute carry limb */ + + lsll R(cnt),R(d2) + movel R(d2),R(d1) + subql #1,R(s_size) + beq L(Lend) + lsrl #1,R(s_size) + bcs L(L1) + subql #1,R(s_size) + +L(Loop:) + movel MEM_PREDEC(s_ptr),R(d2) + movel R(d2),R(d3) + lsrl R(d5),R(d3) + orl R(d3),R(d1) + movel R(d1),MEM_PREDEC(res_ptr) + lsll R(cnt),R(d2) +L(L1:) + movel MEM_PREDEC(s_ptr),R(d1) + movel R(d1),R(d3) + lsrl R(d5),R(d3) + orl R(d3),R(d2) + movel R(d2),MEM_PREDEC(res_ptr) + lsll R(cnt),R(d1) + + dbf R(s_size),L(Loop) + subl #0x10000,R(s_size) + bcc L(Loop) + +L(Lend:) + movel R(d1),MEM_PREDEC(res_ptr) /* store least significant limb */ + +/* Restore used registers from stack frame. */ + moveml MEM_POSTINC(sp),R(d2)-R(d6)/R(a2) + cfi_remember_state + cfi_adjust_cfa_offset (-6*4) + cfi_restore (R(d2)) + cfi_restore (R(d3)) + cfi_restore (R(d4)) + cfi_restore (R(d5)) + cfi_restore (R(d6)) + cfi_restore (R(a2)) + rts + +/* We loop from least significant end of the arrays, which is only + permissible if the source and destination don't overlap, since the + function is documented to work for overlapping source and destination. */ + + cfi_restore_state +L(Lspecial:) + clrl R(d0) /* initialize carry */ + eorw #1,R(s_size) + lsrl #1,R(s_size) + bcc L(LL1) + subql #1,R(s_size) + +L(LLoop:) + movel MEM_POSTINC(s_ptr),R(d2) + addxl R(d2),R(d2) + movel R(d2),MEM_POSTINC(res_ptr) +L(LL1:) + movel MEM_POSTINC(s_ptr),R(d2) + addxl R(d2),R(d2) + movel R(d2),MEM_POSTINC(res_ptr) + + dbf R(s_size),L(LLoop) + addxl R(d0),R(d0) /* save cy in lsb */ + subl #0x10000,R(s_size) + bcs L(LLend) + lsrl #1,R(d0) /* restore cy */ + bra L(LLoop) + +L(LLend:) +/* Restore used registers from stack frame. */ + moveml MEM_POSTINC(sp),R(d2)-R(d6)/R(a2) + cfi_adjust_cfa_offset (-6*4) + cfi_restore (R(d2)) + cfi_restore (R(d3)) + cfi_restore (R(d4)) + cfi_restore (R(d5)) + cfi_restore (R(d6)) + cfi_restore (R(a2)) + rts +END(__mpn_lshift) diff --git a/REORG.TODO/sysdeps/m68k/m680x0/m68020/Makefile b/REORG.TODO/sysdeps/m68k/m680x0/m68020/Makefile new file mode 100644 index 0000000000..b17635467d --- /dev/null +++ b/REORG.TODO/sysdeps/m68k/m680x0/m68020/Makefile @@ -0,0 +1,3 @@ +ifeq ($(subdir),db2) +CPPFLAGS += -DHAVE_SPINLOCKS=1 -DHAVE_ASSEM_MC68020_GCC=1 +endif diff --git a/REORG.TODO/sysdeps/m68k/m680x0/m68020/addmul_1.S b/REORG.TODO/sysdeps/m68k/m680x0/m68020/addmul_1.S new file mode 100644 index 0000000000..53856c0846 --- /dev/null +++ b/REORG.TODO/sysdeps/m68k/m680x0/m68020/addmul_1.S @@ -0,0 +1,88 @@ +/* mc68020 __mpn_addmul_1 -- Multiply a limb vector with a limb and add + the result to a second limb vector. + +Copyright (C) 1992-2017 Free Software Foundation, Inc. + +This file is part of the GNU MP Library. + +The GNU MP Library 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. + +The GNU MP 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 Lesser General Public +License for more details. + +You should have received a copy of the GNU Lesser General Public License +along with the GNU MP Library. If not, see <http://www.gnu.org/licenses/>. */ + +/* + INPUT PARAMETERS + res_ptr (sp + 4) + s1_ptr (sp + 8) + s1_size (sp + 12) + s2_limb (sp + 16) +*/ + +#include "sysdep.h" +#include "asm-syntax.h" + + TEXT +ENTRY(__mpn_addmul_1) + +#define res_ptr a0 +#define s1_ptr a1 +#define s1_size d2 +#define s2_limb d4 + +/* Save used registers on the stack. */ + moveml R(d2)-R(d5),MEM_PREDEC(sp) + cfi_adjust_cfa_offset (4*4) + cfi_rel_offset (R(d2), 0) + cfi_rel_offset (R(d3), 4) + cfi_rel_offset (R(d4), 8) + cfi_rel_offset (R(d5), 12) + +/* Copy the arguments to registers. Better use movem? */ + movel MEM_DISP(sp,20),R(res_ptr) + movel MEM_DISP(sp,24),R(s1_ptr) + movel MEM_DISP(sp,28),R(s1_size) + movel MEM_DISP(sp,32),R(s2_limb) + + eorw #1,R(s1_size) + clrl R(d1) + clrl R(d5) + lsrl #1,R(s1_size) + bcc L(L1) + subql #1,R(s1_size) + subl R(d0),R(d0) /* (d0,cy) <= (0,0) */ + +L(Loop:) + movel MEM_POSTINC(s1_ptr),R(d3) + mulul R(s2_limb),R(d1):R(d3) + addxl R(d0),R(d3) + addxl R(d5),R(d1) + addl R(d3),MEM_POSTINC(res_ptr) +L(L1:) movel MEM_POSTINC(s1_ptr),R(d3) + mulul R(s2_limb),R(d0):R(d3) + addxl R(d1),R(d3) + addxl R(d5),R(d0) + addl R(d3),MEM_POSTINC(res_ptr) + + dbf R(s1_size),L(Loop) + addxl R(d5),R(d0) + subl #0x10000,R(s1_size) + bcc L(Loop) + +/* Restore used registers from stack frame. */ + moveml MEM_POSTINC(sp),R(d2)-R(d5) + cfi_adjust_cfa_offset (-4*4) + cfi_restore (R(d2)) + cfi_restore (R(d3)) + cfi_restore (R(d4)) + cfi_restore (R(d5)) + + rts +END(__mpn_addmul_1) diff --git a/REORG.TODO/sysdeps/m68k/m680x0/m68020/atomic-machine.h b/REORG.TODO/sysdeps/m68k/m680x0/m68020/atomic-machine.h new file mode 100644 index 0000000000..83238f51de --- /dev/null +++ b/REORG.TODO/sysdeps/m68k/m680x0/m68020/atomic-machine.h @@ -0,0 +1,261 @@ +/* Copyright (C) 2003-2017 Free Software Foundation, Inc. + This file is part of the GNU C Library. + Contributed by Andreas Schwab <schwab@suse.de>, 2003. + + The GNU C Library 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. + + 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 + Lesser General Public License for more details. + + You should have received a copy of the GNU Lesser General Public + License along with the GNU C Library. If not, see + <http://www.gnu.org/licenses/>. */ + +#include <stdint.h> + + +typedef int8_t atomic8_t; +typedef uint8_t uatomic8_t; +typedef int_fast8_t atomic_fast8_t; +typedef uint_fast8_t uatomic_fast8_t; + +typedef int16_t atomic16_t; +typedef uint16_t uatomic16_t; +typedef int_fast16_t atomic_fast16_t; +typedef uint_fast16_t uatomic_fast16_t; + +typedef int32_t atomic32_t; +typedef uint32_t uatomic32_t; +typedef int_fast32_t atomic_fast32_t; +typedef uint_fast32_t uatomic_fast32_t; + +typedef int64_t atomic64_t; +typedef uint64_t uatomic64_t; +typedef int_fast64_t atomic_fast64_t; +typedef uint_fast64_t uatomic_fast64_t; + +typedef intptr_t atomicptr_t; +typedef uintptr_t uatomicptr_t; +typedef intmax_t atomic_max_t; +typedef uintmax_t uatomic_max_t; + +#define __HAVE_64B_ATOMICS 1 +#define USE_ATOMIC_COMPILER_BUILTINS 0 + +/* XXX Is this actually correct? */ +#define ATOMIC_EXCHANGE_USES_CAS 1 + +#define __arch_compare_and_exchange_val_8_acq(mem, newval, oldval) \ + ({ __typeof (*(mem)) __ret; \ + __asm __volatile ("cas%.b %0,%2,%1" \ + : "=d" (__ret), "+m" (*(mem)) \ + : "d" (newval), "0" (oldval)); \ + __ret; }) + +#define __arch_compare_and_exchange_val_16_acq(mem, newval, oldval) \ + ({ __typeof (*(mem)) __ret; \ + __asm __volatile ("cas%.w %0,%2,%1" \ + : "=d" (__ret), "+m" (*(mem)) \ + : "d" (newval), "0" (oldval)); \ + __ret; }) + +#define __arch_compare_and_exchange_val_32_acq(mem, newval, oldval) \ + ({ __typeof (*(mem)) __ret; \ + __asm __volatile ("cas%.l %0,%2,%1" \ + : "=d" (__ret), "+m" (*(mem)) \ + : "d" (newval), "0" (oldval)); \ + __ret; }) + +# define __arch_compare_and_exchange_val_64_acq(mem, newval, oldval) \ + ({ __typeof (*(mem)) __ret; \ + __typeof (mem) __memp = (mem); \ + __asm __volatile ("cas2%.l %0:%R0,%1:%R1,(%2):(%3)" \ + : "=d" (__ret) \ + : "d" ((__typeof (*(mem))) (newval)), "r" (__memp), \ + "r" ((char *) __memp + 4), "0" (oldval) \ + : "memory"); \ + __ret; }) + +#define atomic_exchange_acq(mem, newvalue) \ + ({ __typeof (*(mem)) __result = *(mem); \ + if (sizeof (*(mem)) == 1) \ + __asm __volatile ("1: cas%.b %0,%2,%1;" \ + " jbne 1b" \ + : "=d" (__result), "+m" (*(mem)) \ + : "d" (newvalue), "0" (__result)); \ + else if (sizeof (*(mem)) == 2) \ + __asm __volatile ("1: cas%.w %0,%2,%1;" \ + " jbne 1b" \ + : "=d" (__result), "+m" (*(mem)) \ + : "d" (newvalue), "0" (__result)); \ + else if (sizeof (*(mem)) == 4) \ + __asm __volatile ("1: cas%.l %0,%2,%1;" \ + " jbne 1b" \ + : "=d" (__result), "+m" (*(mem)) \ + : "d" (newvalue), "0" (__result)); \ + else \ + { \ + __typeof (mem) __memp = (mem); \ + __asm __volatile ("1: cas2%.l %0:%R0,%1:%R1,(%2):(%3);" \ + " jbne 1b" \ + : "=d" (__result) \ + : "d" ((__typeof (*(mem))) (newvalue)), \ + "r" (__memp), "r" ((char *) __memp + 4), \ + "0" (__result) \ + : "memory"); \ + } \ + __result; }) + +#define atomic_exchange_and_add(mem, value) \ + ({ __typeof (*(mem)) __result = *(mem); \ + __typeof (*(mem)) __temp; \ + if (sizeof (*(mem)) == 1) \ + __asm __volatile ("1: move%.b %0,%2;" \ + " add%.b %3,%2;" \ + " cas%.b %0,%2,%1;" \ + " jbne 1b" \ + : "=d" (__result), "+m" (*(mem)), \ + "=&d" (__temp) \ + : "d" (value), "0" (__result)); \ + else if (sizeof (*(mem)) == 2) \ + __asm __volatile ("1: move%.w %0,%2;" \ + " add%.w %3,%2;" \ + " cas%.w %0,%2,%1;" \ + " jbne 1b" \ + : "=d" (__result), "+m" (*(mem)), \ + "=&d" (__temp) \ + : "d" (value), "0" (__result)); \ + else if (sizeof (*(mem)) == 4) \ + __asm __volatile ("1: move%.l %0,%2;" \ + " add%.l %3,%2;" \ + " cas%.l %0,%2,%1;" \ + " jbne 1b" \ + : "=d" (__result), "+m" (*(mem)), \ + "=&d" (__temp) \ + : "d" (value), "0" (__result)); \ + else \ + { \ + __typeof (mem) __memp = (mem); \ + __asm __volatile ("1: move%.l %0,%1;" \ + " move%.l %R0,%R1;" \ + " add%.l %R2,%R1;" \ + " addx%.l %2,%1;" \ + " cas2%.l %0:%R0,%1:%R1,(%3):(%4);" \ + " jbne 1b" \ + : "=d" (__result), "=&d" (__temp) \ + : "d" ((__typeof (*(mem))) (value)), "r" (__memp), \ + "r" ((char *) __memp + 4), "0" (__result) \ + : "memory"); \ + } \ + __result; }) + +#define atomic_add(mem, value) \ + (void) ({ if (sizeof (*(mem)) == 1) \ + __asm __volatile ("add%.b %1,%0" \ + : "+m" (*(mem)) \ + : "id" (value)); \ + else if (sizeof (*(mem)) == 2) \ + __asm __volatile ("add%.w %1,%0" \ + : "+m" (*(mem)) \ + : "id" (value)); \ + else if (sizeof (*(mem)) == 4) \ + __asm __volatile ("add%.l %1,%0" \ + : "+m" (*(mem)) \ + : "id" (value)); \ + else \ + { \ + __typeof (mem) __memp = (mem); \ + __typeof (*(mem)) __oldval = *__memp; \ + __typeof (*(mem)) __temp; \ + __asm __volatile ("1: move%.l %0,%1;" \ + " move%.l %R0,%R1;" \ + " add%.l %R2,%R1;" \ + " addx%.l %2,%1;" \ + " cas2%.l %0:%R0,%1:%R1,(%3):(%4);" \ + " jbne 1b" \ + : "=d" (__oldval), "=&d" (__temp) \ + : "d" ((__typeof (*(mem))) (value)), \ + "r" (__memp), "r" ((char *) __memp + 4), \ + "0" (__oldval) \ + : "memory"); \ + } \ + }) + +#define atomic_increment_and_test(mem) \ + ({ char __result; \ + if (sizeof (*(mem)) == 1) \ + __asm __volatile ("addq%.b %#1,%1; seq %0" \ + : "=dm" (__result), "+m" (*(mem))); \ + else if (sizeof (*(mem)) == 2) \ + __asm __volatile ("addq%.w %#1,%1; seq %0" \ + : "=dm" (__result), "+m" (*(mem))); \ + else if (sizeof (*(mem)) == 4) \ + __asm __volatile ("addq%.l %#1,%1; seq %0" \ + : "=dm" (__result), "+m" (*(mem))); \ + else \ + { \ + __typeof (mem) __memp = (mem); \ + __typeof (*(mem)) __oldval = *__memp; \ + __typeof (*(mem)) __temp; \ + __asm __volatile ("1: move%.l %1,%2;" \ + " move%.l %R1,%R2;" \ + " addq%.l %#1,%R2;" \ + " addx%.l %5,%2;" \ + " seq %0;" \ + " cas2%.l %1:%R1,%2:%R2,(%3):(%4);" \ + " jbne 1b" \ + : "=&dm" (__result), "=d" (__oldval), \ + "=&d" (__temp) \ + : "r" (__memp), "r" ((char *) __memp + 4), \ + "d" (0), "1" (__oldval) \ + : "memory"); \ + } \ + __result; }) + +#define atomic_decrement_and_test(mem) \ + ({ char __result; \ + if (sizeof (*(mem)) == 1) \ + __asm __volatile ("subq%.b %#1,%1; seq %0" \ + : "=dm" (__result), "+m" (*(mem))); \ + else if (sizeof (*(mem)) == 2) \ + __asm __volatile ("subq%.w %#1,%1; seq %0" \ + : "=dm" (__result), "+m" (*(mem))); \ + else if (sizeof (*(mem)) == 4) \ + __asm __volatile ("subq%.l %#1,%1; seq %0" \ + : "=dm" (__result), "+m" (*(mem))); \ + else \ + { \ + __typeof (mem) __memp = (mem); \ + __typeof (*(mem)) __oldval = *__memp; \ + __typeof (*(mem)) __temp; \ + __asm __volatile ("1: move%.l %1,%2;" \ + " move%.l %R1,%R2;" \ + " subq%.l %#1,%R2;" \ + " subx%.l %5,%2;" \ + " seq %0;" \ + " cas2%.l %1:%R1,%2:%R2,(%3):(%4);" \ + " jbne 1b" \ + : "=&dm" (__result), "=d" (__oldval), \ + "=&d" (__temp) \ + : "r" (__memp), "r" ((char *) __memp + 4), \ + "d" (0), "1" (__oldval) \ + : "memory"); \ + } \ + __result; }) + +#define atomic_bit_set(mem, bit) \ + __asm __volatile ("bfset %0{%1,#1}" \ + : "+m" (*(mem)) \ + : "di" (sizeof (*(mem)) * 8 - (bit) - 1)) + +#define atomic_bit_test_set(mem, bit) \ + ({ char __result; \ + __asm __volatile ("bfset %1{%2,#1}; sne %0" \ + : "=dm" (__result), "+m" (*(mem)) \ + : "di" (sizeof (*(mem)) * 8 - (bit) - 1)); \ + __result; }) diff --git a/REORG.TODO/sysdeps/m68k/m680x0/m68020/bits/string.h b/REORG.TODO/sysdeps/m68k/m680x0/m68020/bits/string.h new file mode 100644 index 0000000000..120afc0d79 --- /dev/null +++ b/REORG.TODO/sysdeps/m68k/m680x0/m68020/bits/string.h @@ -0,0 +1,24 @@ +/* Optimized, inlined string functions. m680x0 version, x >= 2. + Copyright (C) 1997-2017 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 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. + + 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 + Lesser General Public License for more details. + + You should have received a copy of the GNU Lesser General Public + License along with the GNU C Library. If not, see + <http://www.gnu.org/licenses/>. */ + +#ifndef _STRING_H +# error "Never use <bits/string.h> directly; include <string.h> instead." +#endif + +/* Use the unaligned string inline ABI. */ +#define _STRING_INLINE_unaligned 1 diff --git a/REORG.TODO/sysdeps/m68k/m680x0/m68020/mul_1.S b/REORG.TODO/sysdeps/m68k/m680x0/m68020/mul_1.S new file mode 100644 index 0000000000..ee871d2e44 --- /dev/null +++ b/REORG.TODO/sysdeps/m68k/m680x0/m68020/mul_1.S @@ -0,0 +1,83 @@ +/* mc68020 __mpn_mul_1 -- Multiply a limb vector with a limb and store + the result in a second limb vector. + +Copyright (C) 1992-2017 Free Software Foundation, Inc. + +This file is part of the GNU MP Library. + +The GNU MP Library 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. + +The GNU MP 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 Lesser General Public +License for more details. + +You should have received a copy of the GNU Lesser General Public License +along with the GNU MP Library. If not, see <http://www.gnu.org/licenses/>. */ + +/* + INPUT PARAMETERS + res_ptr (sp + 4) + s1_ptr (sp + 8) + s1_size (sp + 12) + s2_limb (sp + 16) +*/ + +#include "sysdep.h" +#include "asm-syntax.h" + + TEXT +ENTRY(__mpn_mul_1) + +#define res_ptr a0 +#define s1_ptr a1 +#define s1_size d2 +#define s2_limb d4 + +/* Save used registers on the stack. */ + moveml R(d2)-R(d4),MEM_PREDEC(sp) + cfi_adjust_cfa_offset (3*4) + cfi_rel_offset (R(d2), 0) + cfi_rel_offset (R(d3), 4) + cfi_rel_offset (R(d4), 8) + +/* Copy the arguments to registers. Better use movem? */ + movel MEM_DISP(sp,16),R(res_ptr) + movel MEM_DISP(sp,20),R(s1_ptr) + movel MEM_DISP(sp,24),R(s1_size) + movel MEM_DISP(sp,28),R(s2_limb) + + eorw #1,R(s1_size) + clrl R(d1) + lsrl #1,R(s1_size) + bcc L(L1) + subql #1,R(s1_size) + subl R(d0),R(d0) /* (d0,cy) <= (0,0) */ + +L(Loop:) + movel MEM_POSTINC(s1_ptr),R(d3) + mulul R(s2_limb),R(d1):R(d3) + addxl R(d0),R(d3) + movel R(d3),MEM_POSTINC(res_ptr) +L(L1:) movel MEM_POSTINC(s1_ptr),R(d3) + mulul R(s2_limb),R(d0):R(d3) + addxl R(d1),R(d3) + movel R(d3),MEM_POSTINC(res_ptr) + + dbf R(s1_size),L(Loop) + clrl R(d3) + addxl R(d3),R(d0) + subl #0x10000,R(s1_size) + bcc L(Loop) + +/* Restore used registers from stack frame. */ + moveml MEM_POSTINC(sp),R(d2)-R(d4) + cfi_adjust_cfa_offset (-3*4) + cfi_restore (R(d2)) + cfi_restore (R(d3)) + cfi_restore (R(d4)) + rts +END(__mpn_mul_1) diff --git a/REORG.TODO/sysdeps/m68k/m680x0/m68020/string_private.h b/REORG.TODO/sysdeps/m68k/m680x0/m68020/string_private.h new file mode 100644 index 0000000000..a35ce053bb --- /dev/null +++ b/REORG.TODO/sysdeps/m68k/m680x0/m68020/string_private.h @@ -0,0 +1,21 @@ +/* Define _STRING_ARCH_unaligned. m680x0 version, x >= 2. + Copyright (C) 2016-2017 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 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. + + 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 + Lesser General Public License for more details. + + You should have received a copy of the GNU Lesser General Public + License along with the GNU C Library; if not, see + <http://www.gnu.org/licenses/>. */ + +/* Tell the generic inline macros that unaligned memory access is + possible. */ +#define _STRING_ARCH_unaligned 1 diff --git a/REORG.TODO/sysdeps/m68k/m680x0/m68020/submul_1.S b/REORG.TODO/sysdeps/m68k/m680x0/m68020/submul_1.S new file mode 100644 index 0000000000..268e6bf38d --- /dev/null +++ b/REORG.TODO/sysdeps/m68k/m680x0/m68020/submul_1.S @@ -0,0 +1,88 @@ +/* mc68020 __mpn_submul_1 -- Multiply a limb vector with a limb and subtract + the result from a second limb vector. + +Copyright (C) 1992-2017 Free Software Foundation, Inc. + +This file is part of the GNU MP Library. + +The GNU MP Library 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. + +The GNU MP 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 Lesser General Public +License for more details. + +You should have received a copy of the GNU Lesser General Public License +along with the GNU MP Library. If not, see <http://www.gnu.org/licenses/>. */ + +/* + INPUT PARAMETERS + res_ptr (sp + 4) + s1_ptr (sp + 8) + s1_size (sp + 12) + s2_limb (sp + 16) +*/ + +#include "sysdep.h" +#include "asm-syntax.h" + + TEXT +ENTRY(__mpn_submul_1) + +#define res_ptr a0 +#define s1_ptr a1 +#define s1_size d2 +#define s2_limb d4 + +/* Save used registers on the stack. */ + moveml R(d2)-R(d5),MEM_PREDEC(sp) + cfi_adjust_cfa_offset (4*4) + cfi_rel_offset (R(d2), 0) + cfi_rel_offset (R(d3), 4) + cfi_rel_offset (R(d4), 8) + cfi_rel_offset (R(d5), 12) + +/* Copy the arguments to registers. Better use movem? */ + movel MEM_DISP(sp,20),R(res_ptr) + movel MEM_DISP(sp,24),R(s1_ptr) + movel MEM_DISP(sp,28),R(s1_size) + movel MEM_DISP(sp,32),R(s2_limb) + + eorw #1,R(s1_size) + clrl R(d1) + clrl R(d5) + lsrl #1,R(s1_size) + bcc L(L1) + subql #1,R(s1_size) + subl R(d0),R(d0) /* (d0,cy) <= (0,0) */ + +L(Loop:) + movel MEM_POSTINC(s1_ptr),R(d3) + mulul R(s2_limb),R(d1):R(d3) + addxl R(d0),R(d3) + addxl R(d5),R(d1) + subl R(d3),MEM_POSTINC(res_ptr) +L(L1:) movel MEM_POSTINC(s1_ptr),R(d3) + mulul R(s2_limb),R(d0):R(d3) + addxl R(d1),R(d3) + addxl R(d5),R(d0) + subl R(d3),MEM_POSTINC(res_ptr) + + dbf R(s1_size),L(Loop) + addxl R(d5),R(d0) + subl #0x10000,R(s1_size) + bcc L(Loop) + +/* Restore used registers from stack frame. */ + moveml MEM_POSTINC(sp),R(d2)-R(d5) + cfi_adjust_cfa_offset (-4*4) + cfi_restore (R(d2)) + cfi_restore (R(d3)) + cfi_restore (R(d4)) + cfi_restore (R(d5)) + + rts +END(__mpn_submul_1) diff --git a/REORG.TODO/sysdeps/m68k/m680x0/m68020/wordcopy.S b/REORG.TODO/sysdeps/m68k/m680x0/m68020/wordcopy.S new file mode 100644 index 0000000000..4fb1a4518f --- /dev/null +++ b/REORG.TODO/sysdeps/m68k/m680x0/m68020/wordcopy.S @@ -0,0 +1 @@ +/* Empty, not needed. */ diff --git a/REORG.TODO/sysdeps/m68k/m680x0/printf_fphex.c b/REORG.TODO/sysdeps/m68k/m680x0/printf_fphex.c new file mode 100644 index 0000000000..d021a090ab --- /dev/null +++ b/REORG.TODO/sysdeps/m68k/m680x0/printf_fphex.c @@ -0,0 +1,2 @@ +#define LONG_DOUBLE_DENORM_BIAS IEEE854_LONG_DOUBLE_BIAS +#include <sysdeps/ieee754/ldbl-96/printf_fphex.c> diff --git a/REORG.TODO/sysdeps/m68k/m680x0/rshift.S b/REORG.TODO/sysdeps/m68k/m680x0/rshift.S new file mode 100644 index 0000000000..911179a273 --- /dev/null +++ b/REORG.TODO/sysdeps/m68k/m680x0/rshift.S @@ -0,0 +1,167 @@ +/* mc68020 __mpn_rshift -- Shift right a low-level natural-number integer. + +Copyright (C) 1996-2017 Free Software Foundation, Inc. + +This file is part of the GNU MP Library. + +The GNU MP Library 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. + +The GNU MP 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 Lesser General Public +License for more details. + +You should have received a copy of the GNU Lesser General Public License +along with the GNU MP Library. If not, see <http://www.gnu.org/licenses/>. */ + +/* + INPUT PARAMETERS + res_ptr (sp + 4) + s_ptr (sp + 8) + s_size (sp + 16) + cnt (sp + 12) +*/ + +#include "sysdep.h" +#include "asm-syntax.h" + +#define res_ptr a1 +#define s_ptr a0 +#define s_size d6 +#define cnt d4 + + TEXT +ENTRY(__mpn_rshift) +/* Save used registers on the stack. */ + moveml R(d2)-R(d6)/R(a2),MEM_PREDEC(sp) + cfi_adjust_cfa_offset (6*4) + cfi_rel_offset (R(d2), 0) + cfi_rel_offset (R(d3), 4) + cfi_rel_offset (R(d4), 8) + cfi_rel_offset (R(d5), 12) + cfi_rel_offset (R(d6), 16) + cfi_rel_offset (R(a2), 20) + +/* Copy the arguments to registers. */ + movel MEM_DISP(sp,28),R(res_ptr) + movel MEM_DISP(sp,32),R(s_ptr) + movel MEM_DISP(sp,36),R(s_size) + movel MEM_DISP(sp,40),R(cnt) + + moveql #1,R(d5) + cmpl R(d5),R(cnt) + bne L(Lnormal) + cmpl R(res_ptr),R(s_ptr) + bls L(Lspecial) /* jump if res_ptr >= s_ptr */ +#if (defined (__mc68020__) || defined (__NeXT__) || defined(mc68020)) + lea MEM_INDX1(res_ptr,s_size,l,4),R(a2) +#else /* not mc68020 */ + movel R(s_size),R(d0) + asll #2,R(d0) + lea MEM_INDX(res_ptr,d0,l),R(a2) +#endif + cmpl R(s_ptr),R(a2) + bls L(Lspecial) /* jump if s_ptr >= res_ptr + s_size */ + +L(Lnormal:) + moveql #32,R(d5) + subl R(cnt),R(d5) + movel MEM_POSTINC(s_ptr),R(d2) + movel R(d2),R(d0) + lsll R(d5),R(d0) /* compute carry limb */ + + lsrl R(cnt),R(d2) + movel R(d2),R(d1) + subql #1,R(s_size) + beq L(Lend) + lsrl #1,R(s_size) + bcs L(L1) + subql #1,R(s_size) + +L(Loop:) + movel MEM_POSTINC(s_ptr),R(d2) + movel R(d2),R(d3) + lsll R(d5),R(d3) + orl R(d3),R(d1) + movel R(d1),MEM_POSTINC(res_ptr) + lsrl R(cnt),R(d2) +L(L1:) + movel MEM_POSTINC(s_ptr),R(d1) + movel R(d1),R(d3) + lsll R(d5),R(d3) + orl R(d3),R(d2) + movel R(d2),MEM_POSTINC(res_ptr) + lsrl R(cnt),R(d1) + + dbf R(s_size),L(Loop) + subl #0x10000,R(s_size) + bcc L(Loop) + +L(Lend:) + movel R(d1),MEM(res_ptr) /* store most significant limb */ + +/* Restore used registers from stack frame. */ + moveml MEM_POSTINC(sp),R(d2)-R(d6)/R(a2) + cfi_remember_state + cfi_adjust_cfa_offset (-6*4) + cfi_restore (R(d2)) + cfi_restore (R(d3)) + cfi_restore (R(d4)) + cfi_restore (R(d5)) + cfi_restore (R(d6)) + cfi_restore (R(a2)) + rts + +/* We loop from most significant end of the arrays, which is only + permissible if the source and destination don't overlap, since the + function is documented to work for overlapping source and destination. */ + + cfi_restore_state +L(Lspecial:) +#if (defined (__mc68020__) || defined (__NeXT__) || defined(mc68020)) + lea MEM_INDX1(s_ptr,s_size,l,4),R(s_ptr) + lea MEM_INDX1(res_ptr,s_size,l,4),R(res_ptr) +#else /* not mc68000 */ + movel R(s_size),R(d0) + asll #2,R(d0) + addl R(s_size),R(s_ptr) + addl R(s_size),R(res_ptr) +#endif + + clrl R(d0) /* initialize carry */ + eorw #1,R(s_size) + lsrl #1,R(s_size) + bcc L(LL1) + subql #1,R(s_size) + +L(LLoop:) + movel MEM_PREDEC(s_ptr),R(d2) + roxrl #1,R(d2) + movel R(d2),MEM_PREDEC(res_ptr) +L(LL1:) + movel MEM_PREDEC(s_ptr),R(d2) + roxrl #1,R(d2) + movel R(d2),MEM_PREDEC(res_ptr) + + dbf R(s_size),L(LLoop) + roxrl #1,R(d0) /* save cy in msb */ + subl #0x10000,R(s_size) + bcs L(LLend) + addl R(d0),R(d0) /* restore cy */ + bra L(LLoop) + +L(LLend:) +/* Restore used registers from stack frame. */ + moveml MEM_POSTINC(sp),R(d2)-R(d6)/R(a2) + cfi_adjust_cfa_offset (-6*4) + cfi_restore (R(d2)) + cfi_restore (R(d3)) + cfi_restore (R(d4)) + cfi_restore (R(d5)) + cfi_restore (R(d6)) + cfi_restore (R(a2)) + rts +END(__mpn_rshift) diff --git a/REORG.TODO/sysdeps/m68k/m680x0/strtold_l.c b/REORG.TODO/sysdeps/m68k/m680x0/strtold_l.c new file mode 100644 index 0000000000..481d9924c2 --- /dev/null +++ b/REORG.TODO/sysdeps/m68k/m680x0/strtold_l.c @@ -0,0 +1,2 @@ +#define DENORM_EXP (MIN_EXP - 1) +#include <sysdeps/ieee754/ldbl-96/strtold_l.c> diff --git a/REORG.TODO/sysdeps/m68k/m680x0/sub_n.S b/REORG.TODO/sysdeps/m68k/m680x0/sub_n.S new file mode 100644 index 0000000000..9836b88be7 --- /dev/null +++ b/REORG.TODO/sysdeps/m68k/m680x0/sub_n.S @@ -0,0 +1,82 @@ +/* mc68020 __mpn_sub_n -- Subtract two limb vectors of the same length > 0 and + store difference in a third limb vector. + +Copyright (C) 1992-2017 Free Software Foundation, Inc. + +This file is part of the GNU MP Library. + +The GNU MP Library 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. + +The GNU MP 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 Lesser General Public +License for more details. + +You should have received a copy of the GNU Lesser General Public License +along with the GNU MP Library. If not, see <http://www.gnu.org/licenses/>. */ + +/* + INPUT PARAMETERS + res_ptr (sp + 4) + s1_ptr (sp + 8) + s2_ptr (sp + 16) + size (sp + 12) +*/ + +#include "sysdep.h" +#include "asm-syntax.h" + + TEXT +ENTRY(__mpn_sub_n) +/* Save used registers on the stack. */ + movel R(d2),MEM_PREDEC(sp) + cfi_adjust_cfa_offset (4) + movel R(a2),MEM_PREDEC(sp) + cfi_adjust_cfa_offset (4) + cfi_rel_offset (R(d2), 4) + cfi_rel_offset (R(a2), 0) + +/* Copy the arguments to registers. Better use movem? */ + movel MEM_DISP(sp,12),R(a2) + movel MEM_DISP(sp,16),R(a0) + movel MEM_DISP(sp,20),R(a1) + movel MEM_DISP(sp,24),R(d2) + + eorw #1,R(d2) + lsrl #1,R(d2) + bcc L(L1) + subql #1,R(d2) /* clears cy as side effect */ + +L(Loop:) + movel MEM_POSTINC(a0),R(d0) + movel MEM_POSTINC(a1),R(d1) + subxl R(d1),R(d0) + movel R(d0),MEM_POSTINC(a2) +L(L1:) movel MEM_POSTINC(a0),R(d0) + movel MEM_POSTINC(a1),R(d1) + subxl R(d1),R(d0) + movel R(d0),MEM_POSTINC(a2) + + dbf R(d2),L(Loop) /* loop until 16 lsb of %4 == -1 */ + subxl R(d0),R(d0) /* d0 <= -cy; save cy as 0 or -1 in d0 */ + subl #0x10000,R(d2) + bcs L(L2) + addl R(d0),R(d0) /* restore cy */ + bra L(Loop) + +L(L2:) + negl R(d0) + +/* Restore used registers from stack frame. */ + movel MEM_POSTINC(sp),R(a2) + cfi_adjust_cfa_offset (-4) + cfi_restore (R(a2)) + movel MEM_POSTINC(sp),R(d2) + cfi_adjust_cfa_offset (-4) + cfi_restore (R(d2)) + + rts +END(__mpn_sub_n) diff --git a/REORG.TODO/sysdeps/m68k/m680x0/sysdep.h b/REORG.TODO/sysdeps/m68k/m680x0/sysdep.h new file mode 100644 index 0000000000..b9235fcb42 --- /dev/null +++ b/REORG.TODO/sysdeps/m68k/m680x0/sysdep.h @@ -0,0 +1,39 @@ +/* Assembler macros for m680x0. + Copyright (C) 2010-2017 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 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. + + 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 + Lesser General Public License for more details. + + You should have received a copy of the GNU Lesser General Public + License along with the GNU C Library. If not, see + <http://www.gnu.org/licenses/>. */ + +#include <sysdeps/m68k/sysdep.h> + +#ifdef __ASSEMBLER__ + +/* Perform operation OP with PC-relative SRC as the first operand and + DST as the second. TMP is available as a temporary if needed. */ +# define PCREL_OP(OP, SRC, DST, TMP) \ + OP SRC(%pc), DST + +/* Load the address of the GOT into register R. */ +# define LOAD_GOT(R) \ + lea _GLOBAL_OFFSET_TABLE_@GOTPC (%pc), R + +#else + +/* As above, but PC is the spelling of the PC register. We need this + so that the macro can be used in both normal and extended asms. */ +#define PCREL_OP(OP, SRC, DST, TMP, PC) \ + OP " " SRC "(" PC "), " DST + +#endif /* __ASSEMBLER__ */ |