1 |
|
/* |
2 |
|
* fpu/mathlib.h - Floating-point math support library |
3 |
|
* |
4 |
< |
* Basilisk II (C) 1997-2001 Christian Bauer |
4 |
> |
* Basilisk II (C) 1997-2008 Christian Bauer |
5 |
|
* |
6 |
|
* MC68881/68040 fpu emulation |
7 |
|
* |
103 |
|
|
104 |
|
/* This is the IEEE 754 single-precision format. */ |
105 |
|
struct { |
106 |
< |
#if UAE_BYTE_ORDER == UAE_BIG_ENDIAN |
106 |
> |
#ifdef WORDS_BIGENDIAN |
107 |
|
unsigned int negative:1; |
108 |
|
unsigned int exponent:8; |
109 |
|
unsigned int mantissa:23; |
110 |
< |
#endif /* Big endian. */ |
111 |
< |
#if UAE_BYTE_ORDER == UAE_LITTLE_ENDIAN |
110 |
> |
#else |
111 |
|
unsigned int mantissa:23; |
112 |
|
unsigned int exponent:8; |
113 |
|
unsigned int negative:1; |
114 |
< |
#endif /* Little endian. */ |
114 |
> |
#endif |
115 |
|
} ieee; |
116 |
|
|
117 |
|
/* This format makes it easier to see if a NaN is a signalling NaN. */ |
118 |
|
struct { |
119 |
< |
#if UAE_BYTE_ORDER == UAE_BIG_ENDIAN |
119 |
> |
#ifdef WORDS_BIGENDIAN |
120 |
|
unsigned int negative:1; |
121 |
|
unsigned int exponent:8; |
122 |
|
unsigned int quiet_nan:1; |
123 |
|
unsigned int mantissa:22; |
124 |
< |
#endif /* Big endian. */ |
126 |
< |
#if UAE_BYTE_ORDER == UAE_LITTLE_ENDIAN |
124 |
> |
#else |
125 |
|
unsigned int mantissa:22; |
126 |
|
unsigned int quiet_nan:1; |
127 |
|
unsigned int exponent:8; |
128 |
|
unsigned int negative:1; |
129 |
< |
#endif /* Little endian. */ |
129 |
> |
#endif |
130 |
|
} ieee_nan; |
131 |
|
}; |
132 |
|
|
136 |
|
|
137 |
|
/* This is the IEEE 754 double-precision format. */ |
138 |
|
struct { |
139 |
< |
#if UAE_BYTE_ORDER == UAE_BIG_ENDIAN |
139 |
> |
#ifdef WORDS_BIGENDIAN |
140 |
|
unsigned int negative:1; |
141 |
|
unsigned int exponent:11; |
142 |
|
/* Together these comprise the mantissa. */ |
143 |
|
unsigned int mantissa0:20; |
144 |
|
unsigned int mantissa1:32; |
145 |
< |
#endif /* Big endian. */ |
146 |
< |
#if UAE_BYTE_ORDER == UAE_LITTLE_ENDIAN |
149 |
< |
# if UAE_FLOAT_WORD_ORDER == UAE_BIG_ENDIAN |
145 |
> |
#else |
146 |
> |
# if HOST_FLOAT_WORDS_BIG_ENDIAN |
147 |
|
unsigned int mantissa0:20; |
148 |
|
unsigned int exponent:11; |
149 |
|
unsigned int negative:1; |
155 |
|
unsigned int exponent:11; |
156 |
|
unsigned int negative:1; |
157 |
|
# endif |
158 |
< |
#endif /* Little endian. */ |
158 |
> |
#endif |
159 |
|
} ieee; |
160 |
|
|
161 |
|
/* This format makes it easier to see if a NaN is a signalling NaN. */ |
162 |
|
struct { |
163 |
< |
#if UAE_BYTE_ORDER == UAE_BIG_ENDIAN |
163 |
> |
#ifdef WORDS_BIGENDIAN |
164 |
|
unsigned int negative:1; |
165 |
|
unsigned int exponent:11; |
166 |
|
unsigned int quiet_nan:1; |
168 |
|
unsigned int mantissa0:19; |
169 |
|
unsigned int mantissa1:32; |
170 |
|
#else |
171 |
< |
# if UAE_FLOAT_WORD_ORDER == UAE_BIG_ENDIAN |
171 |
> |
# if HOST_FLOAT_WORDS_BIG_ENDIAN |
172 |
|
unsigned int mantissa0:19; |
173 |
|
unsigned int quiet_nan:1; |
174 |
|
unsigned int exponent:11; |
187 |
|
|
188 |
|
/* This format is used to extract the sign_exponent and mantissa parts only */ |
189 |
|
struct { |
190 |
< |
#if UAE_FLOAT_WORD_ORDER == UAE_BIG_ENDIAN |
190 |
> |
#if HOST_FLOAT_WORDS_BIG_ENDIAN |
191 |
|
unsigned int msw:32; |
192 |
|
unsigned int lsw:32; |
193 |
|
#else |
204 |
|
|
205 |
|
/* This is the IEEE 854 double-extended-precision format. */ |
206 |
|
struct { |
207 |
< |
#if UAE_BYTE_ORDER == UAE_BIG_ENDIAN |
207 |
> |
#ifdef WORDS_BIGENDIAN |
208 |
|
unsigned int negative:1; |
209 |
|
unsigned int exponent:15; |
210 |
|
unsigned int empty:16; |
211 |
|
unsigned int mantissa0:32; |
212 |
|
unsigned int mantissa1:32; |
213 |
< |
#endif |
214 |
< |
#if UAE_BYTE_ORDER == UAE_LITTLE_ENDIAN |
218 |
< |
# if UAE_FLOAT_WORD_ORDER == UAE_BIG_ENDIAN |
213 |
> |
#else |
214 |
> |
# if HOST_FLOAT_WORDS_BIG_ENDIAN |
215 |
|
unsigned int exponent:15; |
216 |
|
unsigned int negative:1; |
217 |
|
unsigned int empty:16; |
229 |
|
|
230 |
|
/* This is for NaNs in the IEEE 854 double-extended-precision format. */ |
231 |
|
struct { |
232 |
< |
#if UAE_BYTE_ORDER == UAE_BIG_ENDIAN |
232 |
> |
#ifdef WORDS_BIGENDIAN |
233 |
|
unsigned int negative:1; |
234 |
|
unsigned int exponent:15; |
235 |
|
unsigned int empty:16; |
237 |
|
unsigned int quiet_nan:1; |
238 |
|
unsigned int mantissa0:30; |
239 |
|
unsigned int mantissa1:32; |
240 |
< |
#endif |
241 |
< |
#if UAE_BYTE_ORDER == UAE_LITTLE_ENDIAN |
246 |
< |
# if UAE_FLOAT_WORD_ORDER == UAE_BIG_ENDIAN |
240 |
> |
#else |
241 |
> |
# if HOST_FLOAT_WORDS_BIG_ENDIAN |
242 |
|
unsigned int exponent:15; |
243 |
|
unsigned int negative:1; |
244 |
|
unsigned int empty:16; |
260 |
|
|
261 |
|
/* This format is used to extract the sign_exponent and mantissa parts only */ |
262 |
|
struct { |
263 |
< |
#if UAE_FLOAT_WORD_ORDER == UAE_BIG_ENDIAN |
263 |
> |
#if HOST_FLOAT_WORDS_BIG_ENDIAN |
264 |
|
unsigned int sign_exponent:16; |
265 |
|
unsigned int empty:16; |
266 |
|
unsigned int msw:32; |
282 |
|
|
283 |
|
/* This is the IEEE 854 quad-precision format. */ |
284 |
|
struct { |
285 |
< |
#if UAE_BYTE_ORDER == UAE_BIG_ENDIAN |
285 |
> |
#ifdef WORDS_BIGENDIAN |
286 |
|
unsigned int negative:1; |
287 |
|
unsigned int exponent:15; |
288 |
|
unsigned int mantissa0:16; |
301 |
|
|
302 |
|
/* This is for NaNs in the IEEE 854 quad-precision format. */ |
303 |
|
struct { |
304 |
< |
#if UAE_BYTE_ORDER == UAE_BIG_ENDIAN |
304 |
> |
#ifdef WORDS_BIGENDIAN |
305 |
|
unsigned int negative:1; |
306 |
|
unsigned int exponent:15; |
307 |
|
unsigned int quiet_nan:1; |
321 |
|
} ieee_nan; |
322 |
|
|
323 |
|
/* This format is used to extract the sign_exponent and mantissa parts only */ |
324 |
< |
#if UAE_FLOAT_WORD_ORDER == UAE_BIG_ENDIAN |
324 |
> |
#if HOST_FLOAT_WORDS_BIG_ENDIAN |
325 |
|
struct { |
326 |
|
uae_u64 msw; |
327 |
|
uae_u64 lsw; |
395 |
|
return (int)(((uae_u32)(se)) >> 16); |
396 |
|
#endif |
397 |
|
#else |
398 |
< |
#ifndef USE_LONG_DOUBLE |
404 |
< |
fp_declare_init_shape(sxp, r, double); |
405 |
< |
return (sxp->ieee_nan.exponent == FP_DOUBLE_EXP_MAX) |
406 |
< |
#else |
398 |
> |
#if USE_LONG_DOUBLE || USE_QUAD_DOUBLE |
399 |
|
fp_declare_init_shape(sxp, r, extended); |
400 |
|
return (sxp->ieee_nan.exponent == FP_EXTENDED_EXP_MAX) |
401 |
+ |
#else |
402 |
+ |
fp_declare_init_shape(sxp, r, double); |
403 |
+ |
return (sxp->ieee_nan.exponent == FP_DOUBLE_EXP_MAX) |
404 |
|
#endif |
405 |
|
&& (sxp->ieee_nan.mantissa0 != 0) |
406 |
|
&& (sxp->ieee_nan.mantissa1 != 0) |
454 |
|
return ~(lx >> 31) & (1 - (se >> 14)); |
455 |
|
#endif |
456 |
|
#else |
457 |
< |
#ifndef USE_LONG_DOUBLE |
463 |
< |
fp_declare_init_shape(sxp, r, double); |
464 |
< |
return (sxp->ieee_nan.exponent == FP_DOUBLE_EXP_MAX) |
465 |
< |
#else |
457 |
> |
#if USE_LONG_DOUBLE || USE_QUAD_DOUBLE |
458 |
|
fp_declare_init_shape(sxp, r, extended); |
459 |
|
return (sxp->ieee_nan.exponent == FP_EXTENDED_EXP_MAX) |
460 |
+ |
#else |
461 |
+ |
fp_declare_init_shape(sxp, r, double); |
462 |
+ |
return (sxp->ieee_nan.exponent == FP_DOUBLE_EXP_MAX) |
463 |
|
#endif |
464 |
|
&& (sxp->ieee_nan.mantissa0 == 0) |
465 |
|
&& (sxp->ieee_nan.mantissa1 == 0) |
476 |
|
|
477 |
|
PRIVATE inline bool FFPU fp_do_isneg(fpu_register const & r) |
478 |
|
{ |
479 |
< |
#ifndef USE_LONG_DOUBLE |
485 |
< |
fp_declare_init_shape(sxp, r, double); |
486 |
< |
#else |
479 |
> |
#if USE_LONG_DOUBLE || USE_QUAD_DOUBLE |
480 |
|
fp_declare_init_shape(sxp, r, extended); |
481 |
+ |
#else |
482 |
+ |
fp_declare_init_shape(sxp, r, double); |
483 |
|
#endif |
484 |
|
return sxp->ieee.negative; |
485 |
|
} |
490 |
|
PRIVATE inline bool FFPU fp_do_iszero(fpu_register const & r) |
491 |
|
{ |
492 |
|
// TODO: BRANCHES_ARE_EXPENSIVE |
493 |
< |
#ifndef USE_LONG_DOUBLE |
499 |
< |
fp_declare_init_shape(sxp, r, double); |
500 |
< |
#else |
493 |
> |
#if USE_LONG_DOUBLE || USE_QUAD_DOUBLE |
494 |
|
fp_declare_init_shape(sxp, r, extended); |
495 |
+ |
#else |
496 |
+ |
fp_declare_init_shape(sxp, r, double); |
497 |
|
#endif |
498 |
|
return (sxp->ieee.exponent == 0) |
499 |
|
&& (sxp->ieee.mantissa0 == 0) |
526 |
|
PRIVATE inline void FFPU make_nan(fpu_register & r) |
527 |
|
{ |
528 |
|
// FIXME: is that correct ? |
529 |
< |
#ifndef USE_LONG_DOUBLE |
535 |
< |
fp_declare_init_shape(sxp, r, double); |
536 |
< |
sxp->ieee.exponent = FP_DOUBLE_EXP_MAX; |
537 |
< |
sxp->ieee.mantissa0 = 0xfffff; |
538 |
< |
#else |
529 |
> |
#if USE_LONG_DOUBLE || USE_QUAD_DOUBLE |
530 |
|
fp_declare_init_shape(sxp, r, extended); |
531 |
|
sxp->ieee.exponent = FP_EXTENDED_EXP_MAX; |
532 |
|
sxp->ieee.mantissa0 = 0xffffffff; |
533 |
+ |
#else |
534 |
+ |
fp_declare_init_shape(sxp, r, double); |
535 |
+ |
sxp->ieee.exponent = FP_DOUBLE_EXP_MAX; |
536 |
+ |
sxp->ieee.mantissa0 = 0xfffff; |
537 |
|
#endif |
538 |
|
sxp->ieee.mantissa1 = 0xffffffff; |
539 |
|
#ifdef USE_QUAD_DOUBLE |
547 |
|
#if 1 |
548 |
|
r = +0.0; |
549 |
|
#else |
550 |
< |
#ifndef USE_LONG_DOUBLE |
556 |
< |
fp_declare_init_shape(sxp, r, double); |
557 |
< |
#else |
550 |
> |
#if USE_LONG_DOUBLE || USE_QUAD_DOUBLE |
551 |
|
fp_declare_init_shape(sxp, r, extended); |
552 |
+ |
#else |
553 |
+ |
fp_declare_init_shape(sxp, r, double); |
554 |
|
#endif |
555 |
|
sxp->ieee.negative = 0; |
556 |
|
sxp->ieee.exponent = 0; |
568 |
|
#if 1 |
569 |
|
r = -0.0; |
570 |
|
#else |
571 |
< |
#ifndef USE_LONG_DOUBLE |
577 |
< |
fp_declare_init_shape(sxp, r, double); |
578 |
< |
#else |
571 |
> |
#if USE_LONG_DOUBLE || USE_QUAD_DOUBLE |
572 |
|
fp_declare_init_shape(sxp, r, extended); |
573 |
+ |
#else |
574 |
+ |
fp_declare_init_shape(sxp, r, double); |
575 |
|
#endif |
576 |
|
sxp->ieee.negative = 1; |
577 |
|
sxp->ieee.exponent = 0; |
586 |
|
|
587 |
|
PRIVATE inline void FFPU make_inf_positive(fpu_register & r) |
588 |
|
{ |
589 |
< |
#ifndef USE_LONG_DOUBLE |
595 |
< |
fp_declare_init_shape(sxp, r, double); |
596 |
< |
sxp->ieee_nan.exponent = FP_DOUBLE_EXP_MAX; |
597 |
< |
#else |
589 |
> |
#if USE_LONG_DOUBLE || USE_QUAD_DOUBLE |
590 |
|
fp_declare_init_shape(sxp, r, extended); |
591 |
|
sxp->ieee_nan.exponent = FP_EXTENDED_EXP_MAX; |
592 |
+ |
#else |
593 |
+ |
fp_declare_init_shape(sxp, r, double); |
594 |
+ |
sxp->ieee_nan.exponent = FP_DOUBLE_EXP_MAX; |
595 |
|
#endif |
596 |
|
sxp->ieee_nan.negative = 0; |
597 |
|
sxp->ieee_nan.mantissa0 = 0; |
604 |
|
|
605 |
|
PRIVATE inline void FFPU make_inf_negative(fpu_register & r) |
606 |
|
{ |
607 |
< |
#ifndef USE_LONG_DOUBLE |
613 |
< |
fp_declare_init_shape(sxp, r, double); |
614 |
< |
sxp->ieee_nan.exponent = FP_DOUBLE_EXP_MAX; |
615 |
< |
#else |
607 |
> |
#if USE_LONG_DOUBLE || USE_QUAD_DOUBLE |
608 |
|
fp_declare_init_shape(sxp, r, extended); |
609 |
|
sxp->ieee_nan.exponent = FP_EXTENDED_EXP_MAX; |
610 |
+ |
#else |
611 |
+ |
fp_declare_init_shape(sxp, r, double); |
612 |
+ |
sxp->ieee_nan.exponent = FP_DOUBLE_EXP_MAX; |
613 |
|
#endif |
614 |
|
sxp->ieee_nan.negative = 1; |
615 |
|
sxp->ieee_nan.mantissa0 = 0; |
622 |
|
|
623 |
|
PRIVATE inline fpu_register FFPU fast_fgetexp(fpu_register const & r) |
624 |
|
{ |
625 |
< |
#ifndef USE_LONG_DOUBLE |
631 |
< |
fp_declare_init_shape(sxp, r, double); |
632 |
< |
return (sxp->ieee.exponent - FP_DOUBLE_EXP_BIAS); |
633 |
< |
#else |
625 |
> |
#if USE_LONG_DOUBLE || USE_QUAD_DOUBLE |
626 |
|
fp_declare_init_shape(sxp, r, extended); |
627 |
|
return (sxp->ieee.exponent - FP_EXTENDED_EXP_BIAS); |
628 |
+ |
#else |
629 |
+ |
fp_declare_init_shape(sxp, r, double); |
630 |
+ |
return (sxp->ieee.exponent - FP_DOUBLE_EXP_BIAS); |
631 |
|
#endif |
632 |
|
} |
633 |
|
|
634 |
|
// Normalize to range 1..2 |
635 |
|
PRIVATE inline void FFPU fast_remove_exponent(fpu_register & r) |
636 |
|
{ |
637 |
< |
#ifndef USE_LONG_DOUBLE |
643 |
< |
fp_declare_init_shape(sxp, r, double); |
644 |
< |
sxp->ieee.exponent = FP_DOUBLE_EXP_BIAS; |
645 |
< |
#else |
637 |
> |
#if USE_LONG_DOUBLE || USE_QUAD_DOUBLE |
638 |
|
fp_declare_init_shape(sxp, r, extended); |
639 |
|
sxp->ieee.exponent = FP_EXTENDED_EXP_BIAS; |
640 |
+ |
#else |
641 |
+ |
fp_declare_init_shape(sxp, r, double); |
642 |
+ |
sxp->ieee.exponent = FP_DOUBLE_EXP_BIAS; |
643 |
|
#endif |
644 |
|
} |
645 |
|
|
647 |
|
// of the source and destination operands. |
648 |
|
PRIVATE inline uae_u32 FFPU get_quotient_sign(fpu_register const & ra, fpu_register const & rb) |
649 |
|
{ |
650 |
< |
#ifndef USE_LONG_DOUBLE |
656 |
< |
fp_declare_init_shape(sap, ra, double); |
657 |
< |
fp_declare_init_shape(sbp, rb, double); |
658 |
< |
#else |
650 |
> |
#if USE_LONG_DOUBLE || USE_QUAD_DOUBLE |
651 |
|
fp_declare_init_shape(sap, ra, extended); |
652 |
|
fp_declare_init_shape(sbp, rb, extended); |
653 |
+ |
#else |
654 |
+ |
fp_declare_init_shape(sap, ra, double); |
655 |
+ |
fp_declare_init_shape(sbp, rb, double); |
656 |
|
#endif |
657 |
|
return ((sap->ieee.negative ^ sbp->ieee.negative) ? FPSR_QUOTIENT_SIGN : 0); |
658 |
|
} |
661 |
|
/* --- Math functions --- */ |
662 |
|
/* -------------------------------------------------------------------------- */ |
663 |
|
|
664 |
< |
#if FPU_USE_ISO_C99 && USE_LONG_DOUBLE |
665 |
< |
# define fp_log logl |
666 |
< |
# define fp_log10 log10l |
667 |
< |
# define fp_exp expl |
668 |
< |
# define fp_pow powl |
669 |
< |
# define fp_fabs fabsl |
670 |
< |
# define fp_sqrt sqrtl |
671 |
< |
# define fp_sin sinl |
672 |
< |
# define fp_cos cosl |
673 |
< |
# define fp_tan tanl |
674 |
< |
# define fp_sinh sinhl |
675 |
< |
# define fp_cosh coshl |
676 |
< |
# define fp_tanh tanhl |
677 |
< |
# define fp_asin asinl |
678 |
< |
# define fp_acos acosl |
679 |
< |
# define fp_atan atanl |
680 |
< |
# define fp_asinh asinhl |
681 |
< |
# define fp_acosh acoshl |
682 |
< |
# define fp_atanh atanhl |
683 |
< |
# define fp_floor floorl |
684 |
< |
# define fp_ceil ceill |
685 |
< |
#else |
664 |
> |
#if FPU_USE_ISO_C99 && (USE_LONG_DOUBLE || USE_QUAD_DOUBLE) |
665 |
> |
# ifdef HAVE_LOGL |
666 |
> |
# define fp_log logl |
667 |
> |
# endif |
668 |
> |
# ifdef HAVE_LOG10L |
669 |
> |
# define fp_log10 log10l |
670 |
> |
# endif |
671 |
> |
# ifdef HAVE_EXPL |
672 |
> |
# define fp_exp expl |
673 |
> |
# endif |
674 |
> |
# ifdef HAVE_POWL |
675 |
> |
# define fp_pow powl |
676 |
> |
# endif |
677 |
> |
# ifdef HAVE_FABSL |
678 |
> |
# define fp_fabs fabsl |
679 |
> |
# endif |
680 |
> |
# ifdef HAVE_SQRTL |
681 |
> |
# define fp_sqrt sqrtl |
682 |
> |
# endif |
683 |
> |
# ifdef HAVE_SINL |
684 |
> |
# define fp_sin sinl |
685 |
> |
# endif |
686 |
> |
# ifdef HAVE_COSL |
687 |
> |
# define fp_cos cosl |
688 |
> |
# endif |
689 |
> |
# ifdef HAVE_TANL |
690 |
> |
# define fp_tan tanl |
691 |
> |
# endif |
692 |
> |
# ifdef HAVE_SINHL |
693 |
> |
# define fp_sinh sinhl |
694 |
> |
# endif |
695 |
> |
# ifdef HAVE_COSHL |
696 |
> |
# define fp_cosh coshl |
697 |
> |
# endif |
698 |
> |
# ifdef HAVE_TANHL |
699 |
> |
# define fp_tanh tanhl |
700 |
> |
# endif |
701 |
> |
# ifdef HAVE_ASINL |
702 |
> |
# define fp_asin asinl |
703 |
> |
# endif |
704 |
> |
# ifdef HAVE_ACOSL |
705 |
> |
# define fp_acos acosl |
706 |
> |
# endif |
707 |
> |
# ifdef HAVE_ATANL |
708 |
> |
# define fp_atan atanl |
709 |
> |
# endif |
710 |
> |
# ifdef HAVE_ASINHL |
711 |
> |
# define fp_asinh asinhl |
712 |
> |
# endif |
713 |
> |
# ifdef HAVE_ACOSHL |
714 |
> |
# define fp_acosh acoshl |
715 |
> |
# endif |
716 |
> |
# ifdef HAVE_ATANHL |
717 |
> |
# define fp_atanh atanhl |
718 |
> |
# endif |
719 |
> |
# ifdef HAVE_FLOORL |
720 |
> |
# define fp_floor floorl |
721 |
> |
# endif |
722 |
> |
# ifdef HAVE_CEILL |
723 |
> |
# define fp_ceil ceill |
724 |
> |
# endif |
725 |
> |
#endif |
726 |
> |
|
727 |
> |
#ifndef fp_log |
728 |
|
# define fp_log log |
729 |
+ |
#endif |
730 |
+ |
#ifndef fp_log10 |
731 |
|
# define fp_log10 log10 |
732 |
+ |
#endif |
733 |
+ |
#ifndef fp_exp |
734 |
|
# define fp_exp exp |
735 |
+ |
#endif |
736 |
+ |
#ifndef fp_pow |
737 |
|
# define fp_pow pow |
738 |
+ |
#endif |
739 |
+ |
#ifndef fp_fabs |
740 |
|
# define fp_fabs fabs |
741 |
+ |
#endif |
742 |
+ |
#ifndef fp_sqrt |
743 |
|
# define fp_sqrt sqrt |
744 |
+ |
#endif |
745 |
+ |
#ifndef fp_sin |
746 |
|
# define fp_sin sin |
747 |
+ |
#endif |
748 |
+ |
#ifndef fp_cos |
749 |
|
# define fp_cos cos |
750 |
+ |
#endif |
751 |
+ |
#ifndef fp_tan |
752 |
|
# define fp_tan tan |
753 |
+ |
#endif |
754 |
+ |
#ifndef fp_sinh |
755 |
|
# define fp_sinh sinh |
756 |
+ |
#endif |
757 |
+ |
#ifndef fp_cosh |
758 |
|
# define fp_cosh cosh |
759 |
+ |
#endif |
760 |
+ |
#ifndef fp_tanh |
761 |
|
# define fp_tanh tanh |
762 |
+ |
#endif |
763 |
+ |
#ifndef fp_asin |
764 |
|
# define fp_asin asin |
765 |
+ |
#endif |
766 |
+ |
#ifndef fp_acos |
767 |
|
# define fp_acos acos |
768 |
+ |
#endif |
769 |
+ |
#ifndef fp_atan |
770 |
|
# define fp_atan atan |
771 |
+ |
#endif |
772 |
+ |
#ifndef fp_asinh |
773 |
|
# define fp_asinh asinh |
774 |
+ |
#endif |
775 |
+ |
#ifndef fp_acosh |
776 |
|
# define fp_acosh acosh |
777 |
+ |
#endif |
778 |
+ |
#ifndef fp_atanh |
779 |
|
# define fp_atanh atanh |
780 |
+ |
#endif |
781 |
+ |
#ifndef fp_floor |
782 |
|
# define fp_floor floor |
783 |
+ |
#endif |
784 |
+ |
#ifndef fp_ceil |
785 |
|
# define fp_ceil ceil |
786 |
|
#endif |
787 |
|
|
788 |
< |
#if defined(FPU_IEEE) && defined(X86_ASSEMBLY) |
788 |
> |
#if defined(FPU_IEEE) && defined(USE_X87_ASSEMBLY) |
789 |
|
// Assembly optimized support functions. Taken from glibc 2.2.2 |
790 |
|
|
791 |
|
#undef fp_log |
923 |
|
|
924 |
|
PRIVATE inline fpu_extended fp_do_sgn1(fpu_extended x) |
925 |
|
{ |
926 |
+ |
#if USE_LONG_DOUBLE || USE_QUAD_DOUBLE |
927 |
|
fp_declare_init_shape(sxp, x, extended); |
928 |
|
sxp->ieee_nan.exponent = FP_EXTENDED_EXP_MAX; |
929 |
|
sxp->ieee_nan.one = 1; |
930 |
+ |
#else |
931 |
+ |
fp_declare_init_shape(sxp, x, double); |
932 |
+ |
sxp->ieee_nan.exponent = FP_DOUBLE_EXP_MAX; |
933 |
+ |
#endif |
934 |
|
sxp->ieee_nan.quiet_nan = 0; |
935 |
|
sxp->ieee_nan.mantissa0 = 0; |
936 |
|
sxp->ieee_nan.mantissa1 = 0; |
1110 |
|
|
1111 |
|
DEFINE_ROUND_FUNC(nearest, 0x000) |
1112 |
|
|
1113 |
< |
#endif /* X86_ASSEMBLY */ |
1113 |
> |
#endif /* USE_X87_ASSEMBLY */ |
1114 |
|
|
1115 |
|
#ifndef fp_round_to_minus_infinity |
1116 |
|
#define fp_round_to_minus_infinity(x) fp_floor(x) |