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
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
|
/**CHeaderFile*****************************************************************
FileName [epd.h]
PackageName [epd]
Synopsis [The University of Colorado extended double precision package.]
Description [arithmetic functions with extended double precision.]
SeeAlso []
Author [In-Ho Moon]
Copyright [This file was created at the University of Colorado at
Boulder. The University of Colorado at Boulder makes no warranty
about the suitability of this software for any purpose. It is
presented on an AS IS basis.]
Revision [$Id: epd.h,v 1.1.1.1 2003/02/24 22:23:57 wjiang Exp $]
******************************************************************************/
#ifndef _EPD
#define _EPD
ABC_NAMESPACE_HEADER_START
/*---------------------------------------------------------------------------*/
/* Constant declarations */
/*---------------------------------------------------------------------------*/
#define EPD_MAX_BIN 1023
#define EPD_MAX_DEC 308
#define EPD_EXP_INF 0x7ff
/*---------------------------------------------------------------------------*/
/* Structure declarations */
/*---------------------------------------------------------------------------*/
/**Struct**********************************************************************
Synopsis [IEEE double struct.]
Description [IEEE double struct.]
SeeAlso []
******************************************************************************/
#ifdef EPD_BIG_ENDIAN
struct IeeeDoubleStruct { /* BIG_ENDIAN */
unsigned int sign: 1;
unsigned int exponent: 11;
unsigned int mantissa0: 20;
unsigned int mantissa1: 32;
};
#else
struct IeeeDoubleStruct { /* LITTLE_ENDIAN */
unsigned int mantissa1: 32;
unsigned int mantissa0: 20;
unsigned int exponent: 11;
unsigned int sign: 1;
};
#endif
/**Struct**********************************************************************
Synopsis [IEEE double NaN struct.]
Description [IEEE double NaN struct.]
SeeAlso []
******************************************************************************/
#ifdef EPD_BIG_ENDIAN
struct IeeeNanStruct { /* BIG_ENDIAN */
unsigned int sign: 1;
unsigned int exponent: 11;
unsigned int quiet_bit: 1;
unsigned int mantissa0: 19;
unsigned int mantissa1: 32;
};
#else
struct IeeeNanStruct { /* LITTLE_ENDIAN */
unsigned int mantissa1: 32;
unsigned int mantissa0: 19;
unsigned int quiet_bit: 1;
unsigned int exponent: 11;
unsigned int sign: 1;
};
#endif
/**Struct**********************************************************************
Synopsis [Extended precision double to keep very large value.]
Description [Extended precision double to keep very large value.]
SeeAlso []
******************************************************************************/
struct EpDoubleStruct {
union {
double value;
struct IeeeDoubleStruct bits;
struct IeeeNanStruct nan;
} type;
int exponent;
};
/*---------------------------------------------------------------------------*/
/* Type declarations */
/*---------------------------------------------------------------------------*/
typedef struct EpDoubleStruct EpDouble;
typedef struct IeeeDoubleStruct IeeeDouble;
typedef struct IeeeNanStruct IeeeNan;
/*---------------------------------------------------------------------------*/
/* Function prototypes */
/*---------------------------------------------------------------------------*/
EpDouble *EpdAlloc();
int EpdCmp(const char *key1, const char *key2);
void EpdFree(EpDouble *epd);
void EpdGetString(EpDouble *epd, char *str);
void EpdConvert(double value, EpDouble *epd);
void EpdMultiply(EpDouble *epd1, double value);
void EpdMultiply2(EpDouble *epd1, EpDouble *epd2);
void EpdMultiply2Decimal(EpDouble *epd1, EpDouble *epd2);
void EpdMultiply3(EpDouble *epd1, EpDouble *epd2, EpDouble *epd3);
void EpdMultiply3Decimal(EpDouble *epd1, EpDouble *epd2, EpDouble *epd3);
void EpdDivide(EpDouble *epd1, double value);
void EpdDivide2(EpDouble *epd1, EpDouble *epd2);
void EpdDivide3(EpDouble *epd1, EpDouble *epd2, EpDouble *epd3);
void EpdAdd(EpDouble *epd1, double value);
void EpdAdd2(EpDouble *epd1, EpDouble *epd2);
void EpdAdd3(EpDouble *epd1, EpDouble *epd2, EpDouble *epd3);
void EpdSubtract(EpDouble *epd1, double value);
void EpdSubtract2(EpDouble *epd1, EpDouble *epd2);
void EpdSubtract3(EpDouble *epd1, EpDouble *epd2, EpDouble *epd3);
void EpdPow2(int n, EpDouble *epd);
void EpdPow2Decimal(int n, EpDouble *epd);
void EpdNormalize(EpDouble *epd);
void EpdNormalizeDecimal(EpDouble *epd);
void EpdGetValueAndDecimalExponent(EpDouble *epd, double *value, int *exponent);
int EpdGetExponent(double value);
int EpdGetExponentDecimal(double value);
void EpdMakeInf(EpDouble *epd, int sign);
void EpdMakeZero(EpDouble *epd, int sign);
void EpdMakeNan(EpDouble *epd);
void EpdCopy(EpDouble *from, EpDouble *to);
int EpdIsInf(EpDouble *epd);
int EpdIsZero(EpDouble *epd);
int EpdIsNan(EpDouble *epd);
int EpdIsNanOrInf(EpDouble *epd);
int IsInfDouble(double value);
int IsNanDouble(double value);
int IsNanOrInfDouble(double value);
ABC_NAMESPACE_HEADER_END
#endif /* _EPD */
|
