1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
|
/* Copyright (C) 1993 Aladdin Enterprises. All rights reserved.
This file is part of GNU Ghostscript.
GNU Ghostscript is distributed in the hope that it will be useful, but
WITHOUT ANY WARRANTY. No author or distributor accepts responsibility to
anyone for the consequences of using it or for whether it serves any
particular purpose or works at all, unless he says so in writing. Refer
to the GNU Ghostscript General Public License for full details.
*/
/* gxfarith.h */
/* Floating point arithmetic macros for Ghostscript library */
#include "gconfigv.h" /* for USE_FPU */
#include "gxarith.h"
/*
* These macros replace the ones in gxarith.h on machines
* that are likely to have very slow floating point.
*
* None of these macros would be necessary if compilers had a clue
* about generating good floating point comparisons on machines with
* slow (or no) floating point hardware.
*/
# if USE_FPU <= 0 && arch_floats_are_IEEE && (arch_sizeof_float == arch_sizeof_int || arch_sizeof_float == arch_sizeof_long)
# if arch_sizeof_float == arch_sizeof_int
# define _f_int_t int /* not typedef, see _f_uint_t */
# else /* arch_sizeof_float == arch_sizeof_long */
# define _f_int_t long
# endif
# define _f_uint_t unsigned _f_int_t
# define _f_as_int(f) *(_f_int_t *)(&(f))
# define _f_as_uint(f) *(_f_uint_t *)(&(f))
# if arch_sizeof_double == arch_sizeof_int
# define _d_int_t int
# else
# if arch_sizeof_double == arch_sizeof_long
# define _d_int_t long
# endif
# endif
# define _d_uint_t unsigned _d_int_t
# define _d_as_int(f) *(_d_int_t *)(&(d))
# define _d_as_uint(f) *(_d_uint_t *)(&(d))
# define _ftest(v,f,n)\
(sizeof(v)==sizeof(float)?(f):(n))
# ifdef _d_int_t
# define _fdtest(v,f,d,n)\
(sizeof(v)==sizeof(float)?(f):sizeof(v)==sizeof(double)?(d):(n))
# else
# define _fdtest(v,f,d,n)\
_ftest(v,f,n)
# endif
# undef is_fzero
# define is_fzero(f) /* must handle both +0 and -0 */\
_fdtest(f, (_f_as_int(f) << 1) == 0, (_d_as_int(f) << 1) == 0,\
(f) == 0.0)
# undef is_fzero2
# define is_fzero2(f1,f2)\
(sizeof(f1) == sizeof(float) && sizeof(f2) == sizeof(float) ?\
((_f_as_int(f1) | _f_as_int(f2)) << 1) == 0 :\
(f1) == 0.0 && (f2) == 0.0)
# undef is_fneg
# if arch_is_big_endian
# define _is_fnegb(f) (*(byte *)&(f) >= 0x80)
# else
# define _is_fnegb(f) (((byte *)&(f))[sizeof(f) - 1] >= 0x80)
# endif
# if arch_sizeof_float == arch_sizeof_int
# define is_fneg(f)\
(sizeof(f) == sizeof(float) ? _f_as_int(f) < 0 :\
_is_fnegb(f))
# else
# define is_fneg(f) _is_fnegb(f)
# endif
# define IEEE_expt 0x7f800000 /* IEEE exponent mask */
# define IEEE_f1 0x3f800000 /* IEEE 1.0 */
# undef is_fge1
# if arch_sizeof_float == arch_sizeof_int
# define is_fge1(f)\
(sizeof(f) == sizeof(float) ?\
(_f_as_int(f)) >= IEEE_f1 :\
(f) >= 1.0)
# else /* arch_sizeof_float == arch_sizeof_long */
# define is_fge1(f)\
(sizeof(f) == sizeof(float) ?\
(int)(_f_as_int(f) >> 16) >= (IEEE_f1 >> 16) :\
(f) >= 1.0)
# endif
# undef f_fits_in_ubits
# undef f_fits_in_bits
# define _f_bits(n) (4.0 * (1L << ((n) - 2)))
# define f_fits_in_ubits(f, n)\
_ftest(f, _f_as_uint(f) < (_f_uint_t)IEEE_f1 + ((_f_uint_t)(n) << 23),\
(f) >= 0 && (f) < _f_bits(n))
# define f_fits_in_bits(f, n)\
_ftest(f, (_f_as_uint(f) & IEEE_expt) < IEEE_f1 + ((_f_uint_t)((n)-1) << 23),\
(f) >= -_f_bits((n)-1) && (f) < _f_bits((n)-1))
# endif /* USE_FPU <= 0 & ... */
|
