line |
stmt |
bran |
cond |
sub |
pod |
time |
code |
1
|
|
|
|
|
|
|
## This file generated by InlineX::C2XS (version 0.22) using Inline::C (version 0.5501) |
2
|
|
|
|
|
|
|
package Math::Decimal128; |
3
|
|
|
|
|
|
|
|
4
|
30
|
|
|
30
|
|
18266
|
use 5.006; |
|
30
|
|
|
|
|
110
|
|
5
|
|
|
|
|
|
|
|
6
|
30
|
|
|
30
|
|
167
|
use warnings; |
|
30
|
|
|
|
|
64
|
|
|
30
|
|
|
|
|
719
|
|
7
|
30
|
|
|
30
|
|
151
|
use strict; |
|
30
|
|
|
|
|
66
|
|
|
30
|
|
|
|
|
1616
|
|
8
|
|
|
|
|
|
|
|
9
|
|
|
|
|
|
|
require Exporter; |
10
|
|
|
|
|
|
|
*import = \&Exporter::import; |
11
|
|
|
|
|
|
|
require DynaLoader; |
12
|
|
|
|
|
|
|
|
13
|
|
|
|
|
|
|
$Math::Decimal128::VERSION = '0.09'; |
14
|
|
|
|
|
|
|
|
15
|
30
|
|
|
30
|
|
14490
|
use subs qw(DEC128_MAX DEC128_MIN); |
|
30
|
|
|
|
|
633
|
|
|
30
|
|
|
|
|
139
|
|
16
|
|
|
|
|
|
|
|
17
|
|
|
|
|
|
|
use overload |
18
|
30
|
|
|
|
|
466
|
'+' => \&_overload_add, |
19
|
|
|
|
|
|
|
'*' => \&_overload_mul, |
20
|
|
|
|
|
|
|
'-' => \&_overload_sub, |
21
|
|
|
|
|
|
|
'/' => \&_overload_div, |
22
|
|
|
|
|
|
|
'+=' => \&_overload_add_eq, |
23
|
|
|
|
|
|
|
'*=' => \&_overload_mul_eq, |
24
|
|
|
|
|
|
|
'-=' => \&_overload_sub_eq, |
25
|
|
|
|
|
|
|
'/=' => \&_overload_div_eq, |
26
|
|
|
|
|
|
|
'""' => \&_overload_string, |
27
|
|
|
|
|
|
|
'==' => \&_overload_equiv, |
28
|
|
|
|
|
|
|
'!=' => \&_overload_not_equiv, |
29
|
|
|
|
|
|
|
'<' => \&_overload_lt, |
30
|
|
|
|
|
|
|
'>' => \&_overload_gt, |
31
|
|
|
|
|
|
|
'<=' => \&_overload_lte, |
32
|
|
|
|
|
|
|
'>=' => \&_overload_gte, |
33
|
|
|
|
|
|
|
'<=>' => \&_overload_spaceship, |
34
|
|
|
|
|
|
|
'=' => \&_overload_copy, |
35
|
|
|
|
|
|
|
'!' => \&_overload_not, |
36
|
|
|
|
|
|
|
'bool' => \&_overload_true, |
37
|
|
|
|
|
|
|
'abs' => \&_overload_abs, |
38
|
|
|
|
|
|
|
'++' => \&_overload_inc, |
39
|
|
|
|
|
|
|
'--' => \&_overload_dec, |
40
|
|
|
|
|
|
|
'int' => \&_overload_int, |
41
|
|
|
|
|
|
|
'neg' => \&_overload_neg, |
42
|
30
|
|
|
30
|
|
35568
|
; |
|
30
|
|
|
|
|
27814
|
|
43
|
|
|
|
|
|
|
|
44
|
|
|
|
|
|
|
DynaLoader::bootstrap Math::Decimal128 $Math::Decimal128::VERSION; |
45
|
|
|
|
|
|
|
|
46
|
|
|
|
|
|
|
@Math::Decimal128::EXPORT = (); |
47
|
|
|
|
|
|
|
@Math::Decimal128::EXPORT_OK = qw( |
48
|
|
|
|
|
|
|
NaND128 InfD128 ZeroD128 UnityD128 Exp10l NVtoD128 UVtoD128 IVtoD128 PVtoD128 STRtoD128 |
49
|
|
|
|
|
|
|
have_strtod128 D128toNV assignNaNl assignInfl D128toME DPDtoD128 |
50
|
|
|
|
|
|
|
D128toD128 D128toD128 is_NaND128 is_InfD128 is_ZeroD128 DEC128_MAX DEC128_MIN |
51
|
|
|
|
|
|
|
assignMEl d128_bytes MEtoD128 hex2binl decode_d128 decode_bidl decode_dpdl d128_fmt |
52
|
|
|
|
|
|
|
get_expl get_signl PVtoMEl MEtoPVl |
53
|
|
|
|
|
|
|
D128toFSTR D128toRSTR |
54
|
|
|
|
|
|
|
assignIVl assignUVl assignNVl assignD128 assignPVl assignDPDl |
55
|
|
|
|
|
|
|
); |
56
|
|
|
|
|
|
|
|
57
|
|
|
|
|
|
|
%Math::Decimal128::EXPORT_TAGS = (all => [qw( |
58
|
|
|
|
|
|
|
NaND128 InfD128 ZeroD128 UnityD128 Exp10l NVtoD128 UVtoD128 IVtoD128 PVtoD128 STRtoD128 |
59
|
|
|
|
|
|
|
have_strtod128 D128toNV assignNaNl assignInfl D128toME DPDtoD128 |
60
|
|
|
|
|
|
|
D128toD128 D128toD128 is_NaND128 is_InfD128 is_ZeroD128 DEC128_MAX DEC128_MIN |
61
|
|
|
|
|
|
|
assignMEl d128_bytes MEtoD128 hex2binl decode_d128 decode_bidl decode_dpdl d128_fmt |
62
|
|
|
|
|
|
|
get_expl get_signl PVtoMEl MEtoPVl |
63
|
|
|
|
|
|
|
D128toFSTR D128toRSTR |
64
|
|
|
|
|
|
|
assignIVl assignUVl assignNVl assignPVl assignDPDl assignD128 |
65
|
|
|
|
|
|
|
)]); |
66
|
|
|
|
|
|
|
|
67
|
|
|
|
|
|
|
####################################################################### |
68
|
|
|
|
|
|
|
####################################################################### |
69
|
|
|
|
|
|
|
|
70
|
|
|
|
|
|
|
$Math::Decimal128::nan_str = unpack("a*", pack( "B*", '011111' . ('0' x 122))); |
71
|
|
|
|
|
|
|
$Math::Decimal128::ninf_str = unpack("a*", pack( "B*", '11111' . ('0' x 123))); |
72
|
|
|
|
|
|
|
$Math::Decimal128::pinf_str = unpack("a*", pack( "B*", '01111' . ('0' x 123))); |
73
|
|
|
|
|
|
|
$Math::Decimal128::fmt = d128_fmt(); |
74
|
|
|
|
|
|
|
$Math::Decimal128::NNW = 0; # set to 1 to enable a non-numeric warning whenever |
75
|
|
|
|
|
|
|
# a string containing any non-numeric characters is |
76
|
|
|
|
|
|
|
# treated as a number. |
77
|
|
|
|
|
|
|
|
78
|
|
|
|
|
|
|
####################################################################### |
79
|
|
|
|
|
|
|
####################################################################### |
80
|
|
|
|
|
|
|
|
81
|
|
|
|
|
|
|
%Math::Decimal128::dpd_encode = $Math::Decimal128::fmt eq 'DPD' ? ( |
82
|
|
|
|
|
|
|
'0000000000' => '000', '0000000001' => '001', '0000000010' => '002', '0000000011' => '003', |
83
|
|
|
|
|
|
|
'0000000100' => '004', '0000000101' => '005', '0000000110' => '006', '0000000111' => '007', |
84
|
|
|
|
|
|
|
'0000001000' => '008', '0000001001' => '009', '0000010000' => '010', '0000010001' => '011', |
85
|
|
|
|
|
|
|
'0000010010' => '012', '0000010011' => '013', '0000010100' => '014', '0000010101' => '015', |
86
|
|
|
|
|
|
|
'0000010110' => '016', '0000010111' => '017', '0000011000' => '018', '0000011001' => '019', |
87
|
|
|
|
|
|
|
'0000100000' => '020', '0000100001' => '021', '0000100010' => '022', '0000100011' => '023', |
88
|
|
|
|
|
|
|
'0000100100' => '024', '0000100101' => '025', '0000100110' => '026', '0000100111' => '027', |
89
|
|
|
|
|
|
|
'0000101000' => '028', '0000101001' => '029', '0000110000' => '030', '0000110001' => '031', |
90
|
|
|
|
|
|
|
'0000110010' => '032', '0000110011' => '033', '0000110100' => '034', '0000110101' => '035', |
91
|
|
|
|
|
|
|
'0000110110' => '036', '0000110111' => '037', '0000111000' => '038', '0000111001' => '039', |
92
|
|
|
|
|
|
|
'0001000000' => '040', '0001000001' => '041', '0001000010' => '042', '0001000011' => '043', |
93
|
|
|
|
|
|
|
'0001000100' => '044', '0001000101' => '045', '0001000110' => '046', '0001000111' => '047', |
94
|
|
|
|
|
|
|
'0001001000' => '048', '0001001001' => '049', '0001010000' => '050', '0001010001' => '051', |
95
|
|
|
|
|
|
|
'0001010010' => '052', '0001010011' => '053', '0001010100' => '054', '0001010101' => '055', |
96
|
|
|
|
|
|
|
'0001010110' => '056', '0001010111' => '057', '0001011000' => '058', '0001011001' => '059', |
97
|
|
|
|
|
|
|
'0001100000' => '060', '0001100001' => '061', '0001100010' => '062', '0001100011' => '063', |
98
|
|
|
|
|
|
|
'0001100100' => '064', '0001100101' => '065', '0001100110' => '066', '0001100111' => '067', |
99
|
|
|
|
|
|
|
'0001101000' => '068', '0001101001' => '069', '0001110000' => '070', '0001110001' => '071', |
100
|
|
|
|
|
|
|
'0001110010' => '072', '0001110011' => '073', '0001110100' => '074', '0001110101' => '075', |
101
|
|
|
|
|
|
|
'0001110110' => '076', '0001110111' => '077', '0001111000' => '078', '0001111001' => '079', |
102
|
|
|
|
|
|
|
'0000001010' => '080', '0000001011' => '081', '0000101010' => '082', '0000101011' => '083', |
103
|
|
|
|
|
|
|
'0001001010' => '084', '0001001011' => '085', '0001101010' => '086', '0001101011' => '087', |
104
|
|
|
|
|
|
|
'0001001110' => '088', '0001001111' => '089', '0000011010' => '090', '0000011011' => '091', |
105
|
|
|
|
|
|
|
'0000111010' => '092', '0000111011' => '093', '0001011010' => '094', '0001011011' => '095', |
106
|
|
|
|
|
|
|
'0001111010' => '096', '0001111011' => '097', '0001011110' => '098', '0001011111' => '099', |
107
|
|
|
|
|
|
|
'0010000000' => '100', '0010000001' => '101', '0010000010' => '102', '0010000011' => '103', |
108
|
|
|
|
|
|
|
'0010000100' => '104', '0010000101' => '105', '0010000110' => '106', '0010000111' => '107', |
109
|
|
|
|
|
|
|
'0010001000' => '108', '0010001001' => '109', '0010010000' => '110', '0010010001' => '111', |
110
|
|
|
|
|
|
|
'0010010010' => '112', '0010010011' => '113', '0010010100' => '114', '0010010101' => '115', |
111
|
|
|
|
|
|
|
'0010010110' => '116', '0010010111' => '117', '0010011000' => '118', '0010011001' => '119', |
112
|
|
|
|
|
|
|
'0010100000' => '120', '0010100001' => '121', '0010100010' => '122', '0010100011' => '123', |
113
|
|
|
|
|
|
|
'0010100100' => '124', '0010100101' => '125', '0010100110' => '126', '0010100111' => '127', |
114
|
|
|
|
|
|
|
'0010101000' => '128', '0010101001' => '129', '0010110000' => '130', '0010110001' => '131', |
115
|
|
|
|
|
|
|
'0010110010' => '132', '0010110011' => '133', '0010110100' => '134', '0010110101' => '135', |
116
|
|
|
|
|
|
|
'0010110110' => '136', '0010110111' => '137', '0010111000' => '138', '0010111001' => '139', |
117
|
|
|
|
|
|
|
'0011000000' => '140', '0011000001' => '141', '0011000010' => '142', '0011000011' => '143', |
118
|
|
|
|
|
|
|
'0011000100' => '144', '0011000101' => '145', '0011000110' => '146', '0011000111' => '147', |
119
|
|
|
|
|
|
|
'0011001000' => '148', '0011001001' => '149', '0011010000' => '150', '0011010001' => '151', |
120
|
|
|
|
|
|
|
'0011010010' => '152', '0011010011' => '153', '0011010100' => '154', '0011010101' => '155', |
121
|
|
|
|
|
|
|
'0011010110' => '156', '0011010111' => '157', '0011011000' => '158', '0011011001' => '159', |
122
|
|
|
|
|
|
|
'0011100000' => '160', '0011100001' => '161', '0011100010' => '162', '0011100011' => '163', |
123
|
|
|
|
|
|
|
'0011100100' => '164', '0011100101' => '165', '0011100110' => '166', '0011100111' => '167', |
124
|
|
|
|
|
|
|
'0011101000' => '168', '0011101001' => '169', '0011110000' => '170', '0011110001' => '171', |
125
|
|
|
|
|
|
|
'0011110010' => '172', '0011110011' => '173', '0011110100' => '174', '0011110101' => '175', |
126
|
|
|
|
|
|
|
'0011110110' => '176', '0011110111' => '177', '0011111000' => '178', '0011111001' => '179', |
127
|
|
|
|
|
|
|
'0010001010' => '180', '0010001011' => '181', '0010101010' => '182', '0010101011' => '183', |
128
|
|
|
|
|
|
|
'0011001010' => '184', '0011001011' => '185', '0011101010' => '186', '0011101011' => '187', |
129
|
|
|
|
|
|
|
'0011001110' => '188', '0011001111' => '189', '0010011010' => '190', '0010011011' => '191', |
130
|
|
|
|
|
|
|
'0010111010' => '192', '0010111011' => '193', '0011011010' => '194', '0011011011' => '195', |
131
|
|
|
|
|
|
|
'0011111010' => '196', '0011111011' => '197', '0011011110' => '198', '0011011111' => '199', |
132
|
|
|
|
|
|
|
'0100000000' => '200', '0100000001' => '201', '0100000010' => '202', '0100000011' => '203', |
133
|
|
|
|
|
|
|
'0100000100' => '204', '0100000101' => '205', '0100000110' => '206', '0100000111' => '207', |
134
|
|
|
|
|
|
|
'0100001000' => '208', '0100001001' => '209', '0100010000' => '210', '0100010001' => '211', |
135
|
|
|
|
|
|
|
'0100010010' => '212', '0100010011' => '213', '0100010100' => '214', '0100010101' => '215', |
136
|
|
|
|
|
|
|
'0100010110' => '216', '0100010111' => '217', '0100011000' => '218', '0100011001' => '219', |
137
|
|
|
|
|
|
|
'0100100000' => '220', '0100100001' => '221', '0100100010' => '222', '0100100011' => '223', |
138
|
|
|
|
|
|
|
'0100100100' => '224', '0100100101' => '225', '0100100110' => '226', '0100100111' => '227', |
139
|
|
|
|
|
|
|
'0100101000' => '228', '0100101001' => '229', '0100110000' => '230', '0100110001' => '231', |
140
|
|
|
|
|
|
|
'0100110010' => '232', '0100110011' => '233', '0100110100' => '234', '0100110101' => '235', |
141
|
|
|
|
|
|
|
'0100110110' => '236', '0100110111' => '237', '0100111000' => '238', '0100111001' => '239', |
142
|
|
|
|
|
|
|
'0101000000' => '240', '0101000001' => '241', '0101000010' => '242', '0101000011' => '243', |
143
|
|
|
|
|
|
|
'0101000100' => '244', '0101000101' => '245', '0101000110' => '246', '0101000111' => '247', |
144
|
|
|
|
|
|
|
'0101001000' => '248', '0101001001' => '249', '0101010000' => '250', '0101010001' => '251', |
145
|
|
|
|
|
|
|
'0101010010' => '252', '0101010011' => '253', '0101010100' => '254', '0101010101' => '255', |
146
|
|
|
|
|
|
|
'0101010110' => '256', '0101010111' => '257', '0101011000' => '258', '0101011001' => '259', |
147
|
|
|
|
|
|
|
'0101100000' => '260', '0101100001' => '261', '0101100010' => '262', '0101100011' => '263', |
148
|
|
|
|
|
|
|
'0101100100' => '264', '0101100101' => '265', '0101100110' => '266', '0101100111' => '267', |
149
|
|
|
|
|
|
|
'0101101000' => '268', '0101101001' => '269', '0101110000' => '270', '0101110001' => '271', |
150
|
|
|
|
|
|
|
'0101110010' => '272', '0101110011' => '273', '0101110100' => '274', '0101110101' => '275', |
151
|
|
|
|
|
|
|
'0101110110' => '276', '0101110111' => '277', '0101111000' => '278', '0101111001' => '279', |
152
|
|
|
|
|
|
|
'0100001010' => '280', '0100001011' => '281', '0100101010' => '282', '0100101011' => '283', |
153
|
|
|
|
|
|
|
'0101001010' => '284', '0101001011' => '285', '0101101010' => '286', '0101101011' => '287', |
154
|
|
|
|
|
|
|
'0101001110' => '288', '0101001111' => '289', '0100011010' => '290', '0100011011' => '291', |
155
|
|
|
|
|
|
|
'0100111010' => '292', '0100111011' => '293', '0101011010' => '294', '0101011011' => '295', |
156
|
|
|
|
|
|
|
'0101111010' => '296', '0101111011' => '297', '0101011110' => '298', '0101011111' => '299', |
157
|
|
|
|
|
|
|
'0110000000' => '300', '0110000001' => '301', '0110000010' => '302', '0110000011' => '303', |
158
|
|
|
|
|
|
|
'0110000100' => '304', '0110000101' => '305', '0110000110' => '306', '0110000111' => '307', |
159
|
|
|
|
|
|
|
'0110001000' => '308', '0110001001' => '309', '0110010000' => '310', '0110010001' => '311', |
160
|
|
|
|
|
|
|
'0110010010' => '312', '0110010011' => '313', '0110010100' => '314', '0110010101' => '315', |
161
|
|
|
|
|
|
|
'0110010110' => '316', '0110010111' => '317', '0110011000' => '318', '0110011001' => '319', |
162
|
|
|
|
|
|
|
'0110100000' => '320', '0110100001' => '321', '0110100010' => '322', '0110100011' => '323', |
163
|
|
|
|
|
|
|
'0110100100' => '324', '0110100101' => '325', '0110100110' => '326', '0110100111' => '327', |
164
|
|
|
|
|
|
|
'0110101000' => '328', '0110101001' => '329', '0110110000' => '330', '0110110001' => '331', |
165
|
|
|
|
|
|
|
'0110110010' => '332', '0110110011' => '333', '0110110100' => '334', '0110110101' => '335', |
166
|
|
|
|
|
|
|
'0110110110' => '336', '0110110111' => '337', '0110111000' => '338', '0110111001' => '339', |
167
|
|
|
|
|
|
|
'0111000000' => '340', '0111000001' => '341', '0111000010' => '342', '0111000011' => '343', |
168
|
|
|
|
|
|
|
'0111000100' => '344', '0111000101' => '345', '0111000110' => '346', '0111000111' => '347', |
169
|
|
|
|
|
|
|
'0111001000' => '348', '0111001001' => '349', '0111010000' => '350', '0111010001' => '351', |
170
|
|
|
|
|
|
|
'0111010010' => '352', '0111010011' => '353', '0111010100' => '354', '0111010101' => '355', |
171
|
|
|
|
|
|
|
'0111010110' => '356', '0111010111' => '357', '0111011000' => '358', '0111011001' => '359', |
172
|
|
|
|
|
|
|
'0111100000' => '360', '0111100001' => '361', '0111100010' => '362', '0111100011' => '363', |
173
|
|
|
|
|
|
|
'0111100100' => '364', '0111100101' => '365', '0111100110' => '366', '0111100111' => '367', |
174
|
|
|
|
|
|
|
'0111101000' => '368', '0111101001' => '369', '0111110000' => '370', '0111110001' => '371', |
175
|
|
|
|
|
|
|
'0111110010' => '372', '0111110011' => '373', '0111110100' => '374', '0111110101' => '375', |
176
|
|
|
|
|
|
|
'0111110110' => '376', '0111110111' => '377', '0111111000' => '378', '0111111001' => '379', |
177
|
|
|
|
|
|
|
'0110001010' => '380', '0110001011' => '381', '0110101010' => '382', '0110101011' => '383', |
178
|
|
|
|
|
|
|
'0111001010' => '384', '0111001011' => '385', '0111101010' => '386', '0111101011' => '387', |
179
|
|
|
|
|
|
|
'0111001110' => '388', '0111001111' => '389', '0110011010' => '390', '0110011011' => '391', |
180
|
|
|
|
|
|
|
'0110111010' => '392', '0110111011' => '393', '0111011010' => '394', '0111011011' => '395', |
181
|
|
|
|
|
|
|
'0111111010' => '396', '0111111011' => '397', '0111011110' => '398', '0111011111' => '399', |
182
|
|
|
|
|
|
|
'1000000000' => '400', '1000000001' => '401', '1000000010' => '402', '1000000011' => '403', |
183
|
|
|
|
|
|
|
'1000000100' => '404', '1000000101' => '405', '1000000110' => '406', '1000000111' => '407', |
184
|
|
|
|
|
|
|
'1000001000' => '408', '1000001001' => '409', '1000010000' => '410', '1000010001' => '411', |
185
|
|
|
|
|
|
|
'1000010010' => '412', '1000010011' => '413', '1000010100' => '414', '1000010101' => '415', |
186
|
|
|
|
|
|
|
'1000010110' => '416', '1000010111' => '417', '1000011000' => '418', '1000011001' => '419', |
187
|
|
|
|
|
|
|
'1000100000' => '420', '1000100001' => '421', '1000100010' => '422', '1000100011' => '423', |
188
|
|
|
|
|
|
|
'1000100100' => '424', '1000100101' => '425', '1000100110' => '426', '1000100111' => '427', |
189
|
|
|
|
|
|
|
'1000101000' => '428', '1000101001' => '429', '1000110000' => '430', '1000110001' => '431', |
190
|
|
|
|
|
|
|
'1000110010' => '432', '1000110011' => '433', '1000110100' => '434', '1000110101' => '435', |
191
|
|
|
|
|
|
|
'1000110110' => '436', '1000110111' => '437', '1000111000' => '438', '1000111001' => '439', |
192
|
|
|
|
|
|
|
'1001000000' => '440', '1001000001' => '441', '1001000010' => '442', '1001000011' => '443', |
193
|
|
|
|
|
|
|
'1001000100' => '444', '1001000101' => '445', '1001000110' => '446', '1001000111' => '447', |
194
|
|
|
|
|
|
|
'1001001000' => '448', '1001001001' => '449', '1001010000' => '450', '1001010001' => '451', |
195
|
|
|
|
|
|
|
'1001010010' => '452', '1001010011' => '453', '1001010100' => '454', '1001010101' => '455', |
196
|
|
|
|
|
|
|
'1001010110' => '456', '1001010111' => '457', '1001011000' => '458', '1001011001' => '459', |
197
|
|
|
|
|
|
|
'1001100000' => '460', '1001100001' => '461', '1001100010' => '462', '1001100011' => '463', |
198
|
|
|
|
|
|
|
'1001100100' => '464', '1001100101' => '465', '1001100110' => '466', '1001100111' => '467', |
199
|
|
|
|
|
|
|
'1001101000' => '468', '1001101001' => '469', '1001110000' => '470', '1001110001' => '471', |
200
|
|
|
|
|
|
|
'1001110010' => '472', '1001110011' => '473', '1001110100' => '474', '1001110101' => '475', |
201
|
|
|
|
|
|
|
'1001110110' => '476', '1001110111' => '477', '1001111000' => '478', '1001111001' => '479', |
202
|
|
|
|
|
|
|
'1000001010' => '480', '1000001011' => '481', '1000101010' => '482', '1000101011' => '483', |
203
|
|
|
|
|
|
|
'1001001010' => '484', '1001001011' => '485', '1001101010' => '486', '1001101011' => '487', |
204
|
|
|
|
|
|
|
'1001001110' => '488', '1001001111' => '489', '1000011010' => '490', '1000011011' => '491', |
205
|
|
|
|
|
|
|
'1000111010' => '492', '1000111011' => '493', '1001011010' => '494', '1001011011' => '495', |
206
|
|
|
|
|
|
|
'1001111010' => '496', '1001111011' => '497', '1001011110' => '498', '1001011111' => '499', |
207
|
|
|
|
|
|
|
'1010000000' => '500', '1010000001' => '501', '1010000010' => '502', '1010000011' => '503', |
208
|
|
|
|
|
|
|
'1010000100' => '504', '1010000101' => '505', '1010000110' => '506', '1010000111' => '507', |
209
|
|
|
|
|
|
|
'1010001000' => '508', '1010001001' => '509', '1010010000' => '510', '1010010001' => '511', |
210
|
|
|
|
|
|
|
'1010010010' => '512', '1010010011' => '513', '1010010100' => '514', '1010010101' => '515', |
211
|
|
|
|
|
|
|
'1010010110' => '516', '1010010111' => '517', '1010011000' => '518', '1010011001' => '519', |
212
|
|
|
|
|
|
|
'1010100000' => '520', '1010100001' => '521', '1010100010' => '522', '1010100011' => '523', |
213
|
|
|
|
|
|
|
'1010100100' => '524', '1010100101' => '525', '1010100110' => '526', '1010100111' => '527', |
214
|
|
|
|
|
|
|
'1010101000' => '528', '1010101001' => '529', '1010110000' => '530', '1010110001' => '531', |
215
|
|
|
|
|
|
|
'1010110010' => '532', '1010110011' => '533', '1010110100' => '534', '1010110101' => '535', |
216
|
|
|
|
|
|
|
'1010110110' => '536', '1010110111' => '537', '1010111000' => '538', '1010111001' => '539', |
217
|
|
|
|
|
|
|
'1011000000' => '540', '1011000001' => '541', '1011000010' => '542', '1011000011' => '543', |
218
|
|
|
|
|
|
|
'1011000100' => '544', '1011000101' => '545', '1011000110' => '546', '1011000111' => '547', |
219
|
|
|
|
|
|
|
'1011001000' => '548', '1011001001' => '549', '1011010000' => '550', '1011010001' => '551', |
220
|
|
|
|
|
|
|
'1011010010' => '552', '1011010011' => '553', '1011010100' => '554', '1011010101' => '555', |
221
|
|
|
|
|
|
|
'1011010110' => '556', '1011010111' => '557', '1011011000' => '558', '1011011001' => '559', |
222
|
|
|
|
|
|
|
'1011100000' => '560', '1011100001' => '561', '1011100010' => '562', '1011100011' => '563', |
223
|
|
|
|
|
|
|
'1011100100' => '564', '1011100101' => '565', '1011100110' => '566', '1011100111' => '567', |
224
|
|
|
|
|
|
|
'1011101000' => '568', '1011101001' => '569', '1011110000' => '570', '1011110001' => '571', |
225
|
|
|
|
|
|
|
'1011110010' => '572', '1011110011' => '573', '1011110100' => '574', '1011110101' => '575', |
226
|
|
|
|
|
|
|
'1011110110' => '576', '1011110111' => '577', '1011111000' => '578', '1011111001' => '579', |
227
|
|
|
|
|
|
|
'1010001010' => '580', '1010001011' => '581', '1010101010' => '582', '1010101011' => '583', |
228
|
|
|
|
|
|
|
'1011001010' => '584', '1011001011' => '585', '1011101010' => '586', '1011101011' => '587', |
229
|
|
|
|
|
|
|
'1011001110' => '588', '1011001111' => '589', '1010011010' => '590', '1010011011' => '591', |
230
|
|
|
|
|
|
|
'1010111010' => '592', '1010111011' => '593', '1011011010' => '594', '1011011011' => '595', |
231
|
|
|
|
|
|
|
'1011111010' => '596', '1011111011' => '597', '1011011110' => '598', '1011011111' => '599', |
232
|
|
|
|
|
|
|
'1100000000' => '600', '1100000001' => '601', '1100000010' => '602', '1100000011' => '603', |
233
|
|
|
|
|
|
|
'1100000100' => '604', '1100000101' => '605', '1100000110' => '606', '1100000111' => '607', |
234
|
|
|
|
|
|
|
'1100001000' => '608', '1100001001' => '609', '1100010000' => '610', '1100010001' => '611', |
235
|
|
|
|
|
|
|
'1100010010' => '612', '1100010011' => '613', '1100010100' => '614', '1100010101' => '615', |
236
|
|
|
|
|
|
|
'1100010110' => '616', '1100010111' => '617', '1100011000' => '618', '1100011001' => '619', |
237
|
|
|
|
|
|
|
'1100100000' => '620', '1100100001' => '621', '1100100010' => '622', '1100100011' => '623', |
238
|
|
|
|
|
|
|
'1100100100' => '624', '1100100101' => '625', '1100100110' => '626', '1100100111' => '627', |
239
|
|
|
|
|
|
|
'1100101000' => '628', '1100101001' => '629', '1100110000' => '630', '1100110001' => '631', |
240
|
|
|
|
|
|
|
'1100110010' => '632', '1100110011' => '633', '1100110100' => '634', '1100110101' => '635', |
241
|
|
|
|
|
|
|
'1100110110' => '636', '1100110111' => '637', '1100111000' => '638', '1100111001' => '639', |
242
|
|
|
|
|
|
|
'1101000000' => '640', '1101000001' => '641', '1101000010' => '642', '1101000011' => '643', |
243
|
|
|
|
|
|
|
'1101000100' => '644', '1101000101' => '645', '1101000110' => '646', '1101000111' => '647', |
244
|
|
|
|
|
|
|
'1101001000' => '648', '1101001001' => '649', '1101010000' => '650', '1101010001' => '651', |
245
|
|
|
|
|
|
|
'1101010010' => '652', '1101010011' => '653', '1101010100' => '654', '1101010101' => '655', |
246
|
|
|
|
|
|
|
'1101010110' => '656', '1101010111' => '657', '1101011000' => '658', '1101011001' => '659', |
247
|
|
|
|
|
|
|
'1101100000' => '660', '1101100001' => '661', '1101100010' => '662', '1101100011' => '663', |
248
|
|
|
|
|
|
|
'1101100100' => '664', '1101100101' => '665', '1101100110' => '666', '1101100111' => '667', |
249
|
|
|
|
|
|
|
'1101101000' => '668', '1101101001' => '669', '1101110000' => '670', '1101110001' => '671', |
250
|
|
|
|
|
|
|
'1101110010' => '672', '1101110011' => '673', '1101110100' => '674', '1101110101' => '675', |
251
|
|
|
|
|
|
|
'1101110110' => '676', '1101110111' => '677', '1101111000' => '678', '1101111001' => '679', |
252
|
|
|
|
|
|
|
'1100001010' => '680', '1100001011' => '681', '1100101010' => '682', '1100101011' => '683', |
253
|
|
|
|
|
|
|
'1101001010' => '684', '1101001011' => '685', '1101101010' => '686', '1101101011' => '687', |
254
|
|
|
|
|
|
|
'1101001110' => '688', '1101001111' => '689', '1100011010' => '690', '1100011011' => '691', |
255
|
|
|
|
|
|
|
'1100111010' => '692', '1100111011' => '693', '1101011010' => '694', '1101011011' => '695', |
256
|
|
|
|
|
|
|
'1101111010' => '696', '1101111011' => '697', '1101011110' => '698', '1101011111' => '699', |
257
|
|
|
|
|
|
|
'1110000000' => '700', '1110000001' => '701', '1110000010' => '702', '1110000011' => '703', |
258
|
|
|
|
|
|
|
'1110000100' => '704', '1110000101' => '705', '1110000110' => '706', '1110000111' => '707', |
259
|
|
|
|
|
|
|
'1110001000' => '708', '1110001001' => '709', '1110010000' => '710', '1110010001' => '711', |
260
|
|
|
|
|
|
|
'1110010010' => '712', '1110010011' => '713', '1110010100' => '714', '1110010101' => '715', |
261
|
|
|
|
|
|
|
'1110010110' => '716', '1110010111' => '717', '1110011000' => '718', '1110011001' => '719', |
262
|
|
|
|
|
|
|
'1110100000' => '720', '1110100001' => '721', '1110100010' => '722', '1110100011' => '723', |
263
|
|
|
|
|
|
|
'1110100100' => '724', '1110100101' => '725', '1110100110' => '726', '1110100111' => '727', |
264
|
|
|
|
|
|
|
'1110101000' => '728', '1110101001' => '729', '1110110000' => '730', '1110110001' => '731', |
265
|
|
|
|
|
|
|
'1110110010' => '732', '1110110011' => '733', '1110110100' => '734', '1110110101' => '735', |
266
|
|
|
|
|
|
|
'1110110110' => '736', '1110110111' => '737', '1110111000' => '738', '1110111001' => '739', |
267
|
|
|
|
|
|
|
'1111000000' => '740', '1111000001' => '741', '1111000010' => '742', '1111000011' => '743', |
268
|
|
|
|
|
|
|
'1111000100' => '744', '1111000101' => '745', '1111000110' => '746', '1111000111' => '747', |
269
|
|
|
|
|
|
|
'1111001000' => '748', '1111001001' => '749', '1111010000' => '750', '1111010001' => '751', |
270
|
|
|
|
|
|
|
'1111010010' => '752', '1111010011' => '753', '1111010100' => '754', '1111010101' => '755', |
271
|
|
|
|
|
|
|
'1111010110' => '756', '1111010111' => '757', '1111011000' => '758', '1111011001' => '759', |
272
|
|
|
|
|
|
|
'1111100000' => '760', '1111100001' => '761', '1111100010' => '762', '1111100011' => '763', |
273
|
|
|
|
|
|
|
'1111100100' => '764', '1111100101' => '765', '1111100110' => '766', '1111100111' => '767', |
274
|
|
|
|
|
|
|
'1111101000' => '768', '1111101001' => '769', '1111110000' => '770', '1111110001' => '771', |
275
|
|
|
|
|
|
|
'1111110010' => '772', '1111110011' => '773', '1111110100' => '774', '1111110101' => '775', |
276
|
|
|
|
|
|
|
'1111110110' => '776', '1111110111' => '777', '1111111000' => '778', '1111111001' => '779', |
277
|
|
|
|
|
|
|
'1110001010' => '780', '1110001011' => '781', '1110101010' => '782', '1110101011' => '783', |
278
|
|
|
|
|
|
|
'1111001010' => '784', '1111001011' => '785', '1111101010' => '786', '1111101011' => '787', |
279
|
|
|
|
|
|
|
'1111001110' => '788', '1111001111' => '789', '1110011010' => '790', '1110011011' => '791', |
280
|
|
|
|
|
|
|
'1110111010' => '792', '1110111011' => '793', '1111011010' => '794', '1111011011' => '795', |
281
|
|
|
|
|
|
|
'1111111010' => '796', '1111111011' => '797', '1111011110' => '798', '1111011111' => '799', |
282
|
|
|
|
|
|
|
'0000001100' => '800', '0000001101' => '801', '0100001100' => '802', '0100001101' => '803', |
283
|
|
|
|
|
|
|
'1000001100' => '804', '1000001101' => '805', '1100001100' => '806', '1100001101' => '807', |
284
|
|
|
|
|
|
|
'0000101110' => '808', '0000101111' => '809', '0000011100' => '810', '0000011101' => '811', |
285
|
|
|
|
|
|
|
'0100011100' => '812', '0100011101' => '813', '1000011100' => '814', '1000011101' => '815', |
286
|
|
|
|
|
|
|
'1100011100' => '816', '1100011101' => '817', '0000111110' => '818', '0000111111' => '819', |
287
|
|
|
|
|
|
|
'0000101100' => '820', '0000101101' => '821', '0100101100' => '822', '0100101101' => '823', |
288
|
|
|
|
|
|
|
'1000101100' => '824', '1000101101' => '825', '1100101100' => '826', '1100101101' => '827', |
289
|
|
|
|
|
|
|
'0100101110' => '828', '0100101111' => '829', '0000111100' => '830', '0000111101' => '831', |
290
|
|
|
|
|
|
|
'0100111100' => '832', '0100111101' => '833', '1000111100' => '834', '1000111101' => '835', |
291
|
|
|
|
|
|
|
'1100111100' => '836', '1100111101' => '837', '0100111110' => '838', '0100111111' => '839', |
292
|
|
|
|
|
|
|
'0001001100' => '840', '0001001101' => '841', '0101001100' => '842', '0101001101' => '843', |
293
|
|
|
|
|
|
|
'1001001100' => '844', '1001001101' => '845', '1101001100' => '846', '1101001101' => '847', |
294
|
|
|
|
|
|
|
'1000101110' => '848', '1000101111' => '849', '0001011100' => '850', '0001011101' => '851', |
295
|
|
|
|
|
|
|
'0101011100' => '852', '0101011101' => '853', '1001011100' => '854', '1001011101' => '855', |
296
|
|
|
|
|
|
|
'1101011100' => '856', '1101011101' => '857', '1000111110' => '858', '1000111111' => '859', |
297
|
|
|
|
|
|
|
'0001101100' => '860', '0001101101' => '861', '0101101100' => '862', '0101101101' => '863', |
298
|
|
|
|
|
|
|
'1001101100' => '864', '1001101101' => '865', '1101101100' => '866', '1101101101' => '867', |
299
|
|
|
|
|
|
|
'1100101110' => '868', '1100101111' => '869', '0001111100' => '870', '0001111101' => '871', |
300
|
|
|
|
|
|
|
'0101111100' => '872', '0101111101' => '873', '1001111100' => '874', '1001111101' => '875', |
301
|
|
|
|
|
|
|
'1101111100' => '876', '1101111101' => '877', '1100111110' => '878', '1100111111' => '879', |
302
|
|
|
|
|
|
|
'0000001110' => '880', '0000001111' => '881', '0100001110' => '882', '0100001111' => '883', |
303
|
|
|
|
|
|
|
'1000001110' => '884', '1000001111' => '885', '1100001110' => '886', '1100001111' => '887', |
304
|
|
|
|
|
|
|
'0001101110' => '888', '0001101111' => '889', '0000011110' => '890', '0000011111' => '891', |
305
|
|
|
|
|
|
|
'0100011110' => '892', '0100011111' => '893', '1000011110' => '894', '1000011111' => '895', |
306
|
|
|
|
|
|
|
'1100011110' => '896', '1100011111' => '897', '0001111110' => '898', '0001111111' => '899', |
307
|
|
|
|
|
|
|
'0010001100' => '900', '0010001101' => '901', '0110001100' => '902', '0110001101' => '903', |
308
|
|
|
|
|
|
|
'1010001100' => '904', '1010001101' => '905', '1110001100' => '906', '1110001101' => '907', |
309
|
|
|
|
|
|
|
'0010101110' => '908', '0010101111' => '909', '0010011100' => '910', '0010011101' => '911', |
310
|
|
|
|
|
|
|
'0110011100' => '912', '0110011101' => '913', '1010011100' => '914', '1010011101' => '915', |
311
|
|
|
|
|
|
|
'1110011100' => '916', '1110011101' => '917', '0010111110' => '918', '0010111111' => '919', |
312
|
|
|
|
|
|
|
'0010101100' => '920', '0010101101' => '921', '0110101100' => '922', '0110101101' => '923', |
313
|
|
|
|
|
|
|
'1010101100' => '924', '1010101101' => '925', '1110101100' => '926', '1110101101' => '927', |
314
|
|
|
|
|
|
|
'0110101110' => '928', '0110101111' => '929', '0010111100' => '930', '0010111101' => '931', |
315
|
|
|
|
|
|
|
'0110111100' => '932', '0110111101' => '933', '1010111100' => '934', '1010111101' => '935', |
316
|
|
|
|
|
|
|
'1110111100' => '936', '1110111101' => '937', '0110111110' => '938', '0110111111' => '939', |
317
|
|
|
|
|
|
|
'0011001100' => '940', '0011001101' => '941', '0111001100' => '942', '0111001101' => '943', |
318
|
|
|
|
|
|
|
'1011001100' => '944', '1011001101' => '945', '1111001100' => '946', '1111001101' => '947', |
319
|
|
|
|
|
|
|
'1010101110' => '948', '1010101111' => '949', '0011011100' => '950', '0011011101' => '951', |
320
|
|
|
|
|
|
|
'0111011100' => '952', '0111011101' => '953', '1011011100' => '954', '1011011101' => '955', |
321
|
|
|
|
|
|
|
'1111011100' => '956', '1111011101' => '957', '1010111110' => '958', '1010111111' => '959', |
322
|
|
|
|
|
|
|
'0011101100' => '960', '0011101101' => '961', '0111101100' => '962', '0111101101' => '963', |
323
|
|
|
|
|
|
|
'1011101100' => '964', '1011101101' => '965', '1111101100' => '966', '1111101101' => '967', |
324
|
|
|
|
|
|
|
'1110101110' => '968', '1110101111' => '969', '0011111100' => '970', '0011111101' => '971', |
325
|
|
|
|
|
|
|
'0111111100' => '972', '0111111101' => '973', '1011111100' => '974', '1011111101' => '975', |
326
|
|
|
|
|
|
|
'1111111100' => '976', '1111111101' => '977', '1110111110' => '978', '1110111111' => '979', |
327
|
|
|
|
|
|
|
'0010001110' => '980', '0010001111' => '981', '0110001110' => '982', '0110001111' => '983', |
328
|
|
|
|
|
|
|
'1010001110' => '984', '1010001111' => '985', '1110001110' => '986', '1110001111' => '987', |
329
|
|
|
|
|
|
|
'0011101110' => '988', '0011101111' => '989', '0010011110' => '990', '0010011111' => '991', |
330
|
|
|
|
|
|
|
'0110011110' => '992', '0110011111' => '993', '1010011110' => '994', '1010011111' => '995', |
331
|
|
|
|
|
|
|
'1110011110' => '996', '1110011111' => '997', '0011111110' => '998', '0011111111' => '999', |
332
|
|
|
|
|
|
|
) : (); |
333
|
|
|
|
|
|
|
|
334
|
|
|
|
|
|
|
####################################################################### |
335
|
|
|
|
|
|
|
####################################################################### |
336
|
|
|
|
|
|
|
|
337
|
|
|
|
|
|
|
# %Math::Decimal128::dpd_decode is simply %Math::Decimal128::dpd_encode |
338
|
|
|
|
|
|
|
# with the keys and values interchanged. |
339
|
|
|
|
|
|
|
|
340
|
|
|
|
|
|
|
for my $key(keys(%Math::Decimal128::dpd_encode)) { |
341
|
|
|
|
|
|
|
$Math::Decimal128::dpd_decode{$Math::Decimal128::dpd_encode{$key}} = $key; |
342
|
|
|
|
|
|
|
} |
343
|
|
|
|
|
|
|
|
344
|
|
|
|
|
|
|
####################################################################### |
345
|
|
|
|
|
|
|
####################################################################### |
346
|
|
|
|
|
|
|
|
347
|
|
|
|
|
|
|
%Math::Decimal128::bid_decode = $Math::Decimal128::fmt eq 'BID' ? ( |
348
|
|
|
|
|
|
|
0 => MEtoD128('1' . ('0' x 33), 0), 1 => MEtoD128('1' . ('0' x 32), 0), |
349
|
|
|
|
|
|
|
2 => MEtoD128('1' . ('0' x 31), 0), 3 => MEtoD128('1' . ('0' x 30), 0), |
350
|
|
|
|
|
|
|
4 => MEtoD128('1' . ('0' x 29), 0), 5 => MEtoD128('1' . ('0' x 28), 0), |
351
|
|
|
|
|
|
|
6 => MEtoD128('1' . ('0' x 27), 0), 7 => MEtoD128('1' . ('0' x 26), 0), |
352
|
|
|
|
|
|
|
8 => MEtoD128('1' . ('0' x 25), 0), 9 => MEtoD128('1' . ('0' x 24), 0), |
353
|
|
|
|
|
|
|
10 => MEtoD128('1' . ('0' x 23), 0), 11 => MEtoD128('1' . ('0' x 22), 0), |
354
|
|
|
|
|
|
|
12 => MEtoD128('1' . ('0' x 21), 0), 13 => MEtoD128('1' . ('0' x 20), 0), |
355
|
|
|
|
|
|
|
14 => MEtoD128('1' . ('0' x 19), 0), 15 => MEtoD128('1' . ('0' x 18), 0), |
356
|
|
|
|
|
|
|
16 => MEtoD128('1' . ('0' x 17), 0), 17 => MEtoD128('1' . ('0' x 16), 0), |
357
|
|
|
|
|
|
|
18 => MEtoD128('1' . ('0' x 15), 0), 19 => MEtoD128('1' . ('0' x 14), 0), |
358
|
|
|
|
|
|
|
20 => MEtoD128('1' . ('0' x 13), 0), 21 => MEtoD128('1' . ('0' x 12), 0), |
359
|
|
|
|
|
|
|
22 => MEtoD128('1' . ('0' x 11), 0), 23 => MEtoD128('1' . ('0' x 10), 0), |
360
|
|
|
|
|
|
|
24 => MEtoD128('1' . ('0' x 9), 0), 25 => MEtoD128('1' . ('0' x 8), 0), |
361
|
|
|
|
|
|
|
26 => MEtoD128('1' . ('0' x 7), 0), 27 => MEtoD128('1' . ('0' x 6), 0), |
362
|
|
|
|
|
|
|
28 => MEtoD128('1' . ('0' x 5), 0), 29 => MEtoD128('1' . ('0' x 4), 0), |
363
|
|
|
|
|
|
|
30 => MEtoD128('1' . ('0' x 3), 0), 31 => MEtoD128('1' . ('0' x 2), 0), |
364
|
|
|
|
|
|
|
32 => MEtoD128('1' . ('0' x 1), 0), 33 => MEtoD128('1', 0) |
365
|
|
|
|
|
|
|
) : (); |
366
|
|
|
|
|
|
|
|
367
|
|
|
|
|
|
|
####################################################################### |
368
|
|
|
|
|
|
|
####################################################################### |
369
|
|
|
|
|
|
|
# Used only wrt BID encoding |
370
|
|
|
|
|
|
|
|
371
|
|
|
|
|
|
|
sub _decode_mant { |
372
|
189580
|
|
|
189580
|
|
299923
|
my $val = shift; |
373
|
189580
|
|
|
|
|
273479
|
my $ret = ''; |
374
|
189580
|
|
|
|
|
369468
|
for my $i(0 .. 33) { |
375
|
6445720
|
|
|
|
|
8595357
|
my $count = 0; |
376
|
6445720
|
100
|
|
|
|
14558900
|
if($val > 0) { |
377
|
4180101
|
|
|
|
|
10745581
|
while($val >= $Math::Decimal128::bid_decode{$i}) { |
378
|
7371369
|
|
|
|
|
14396331
|
$val -= $Math::Decimal128::bid_decode{$i}; |
379
|
7371369
|
|
|
|
|
17697047
|
$count++; |
380
|
|
|
|
|
|
|
} |
381
|
|
|
|
|
|
|
} |
382
|
6445720
|
|
|
|
|
10944016
|
$ret .= $count; |
383
|
|
|
|
|
|
|
} |
384
|
189580
|
|
|
|
|
611747
|
return $ret; |
385
|
|
|
|
|
|
|
} |
386
|
|
|
|
|
|
|
|
387
|
|
|
|
|
|
|
####################################################################### |
388
|
|
|
|
|
|
|
####################################################################### |
389
|
|
|
|
|
|
|
|
390
|
30
|
|
|
30
|
0
|
10889
|
sub dl_load_flags {0} # Prevent DynaLoader from complaining and croaking |
391
|
|
|
|
|
|
|
|
392
|
|
|
|
|
|
|
####################################################################### |
393
|
|
|
|
|
|
|
####################################################################### |
394
|
|
|
|
|
|
|
|
395
|
|
|
|
|
|
|
sub _overload_string { |
396
|
60
|
|
|
60
|
|
3721
|
my @ret = D128toME($_[0]); |
397
|
60
|
100
|
66
|
|
|
393
|
if(is_InfD128($_[0]) || !$_[0]) {return $ret[0]} |
|
12
|
|
|
|
|
45
|
|
398
|
48
|
|
|
|
|
177
|
return $ret[0] . 'e' . $ret[1]; |
399
|
|
|
|
|
|
|
} |
400
|
|
|
|
|
|
|
|
401
|
|
|
|
|
|
|
####################################################################### |
402
|
|
|
|
|
|
|
####################################################################### |
403
|
|
|
|
|
|
|
|
404
|
|
|
|
|
|
|
sub _overload_int { |
405
|
47732
|
100
|
100
|
47732
|
|
495492
|
if(is_NaND128($_[0]) || is_InfD128($_[0]) || is_ZeroD128($_[0])) {return $_[0]} |
|
333
|
|
100
|
|
|
990
|
|
406
|
47399
|
|
|
|
|
109130
|
my($man, $exp) = D128toME($_[0]); |
407
|
47399
|
100
|
|
|
|
129771
|
if($exp >= 0) {return $_[0]} |
|
23522
|
|
|
|
|
71134
|
|
408
|
23877
|
|
|
|
|
36745
|
my $man_length = length($man); |
409
|
23877
|
100
|
|
|
|
68556
|
$man_length-- if $man =~ /^\-/; |
410
|
23877
|
100
|
|
|
|
54557
|
if(-$exp >= $man_length) { # -1 <= $_[0] <= 1 |
411
|
23862
|
|
|
|
|
59841
|
my $z = ZeroD128(1); |
412
|
23862
|
100
|
|
|
|
73761
|
if($_[0] < $z) {return ZeroD128(-1)} # return -0 |
|
23861
|
|
|
|
|
110211
|
|
413
|
1
|
|
|
|
|
4
|
return $z; # return 0 |
414
|
|
|
|
|
|
|
} |
415
|
|
|
|
|
|
|
|
416
|
15
|
|
|
|
|
28
|
substr($man, $exp, -$exp, ''); |
417
|
15
|
|
|
|
|
33
|
return MEtoD128($man, 0); |
418
|
|
|
|
|
|
|
} |
419
|
|
|
|
|
|
|
|
420
|
|
|
|
|
|
|
####################################################################### |
421
|
|
|
|
|
|
|
####################################################################### |
422
|
|
|
|
|
|
|
|
423
|
|
|
|
|
|
|
sub new { |
424
|
|
|
|
|
|
|
|
425
|
|
|
|
|
|
|
# This function caters for 2 possibilities: |
426
|
|
|
|
|
|
|
# 1) that 'new' has been called OOP style - in which |
427
|
|
|
|
|
|
|
# case there will be a maximum of 2 args |
428
|
|
|
|
|
|
|
# 2) that 'new' has been called as a function - in |
429
|
|
|
|
|
|
|
# which case there will be a maximum of 1 arg. |
430
|
|
|
|
|
|
|
# If there are no args, then we just want to return a |
431
|
|
|
|
|
|
|
# Math::Decimal128 object that's a NaN. |
432
|
|
|
|
|
|
|
|
433
|
309406
|
100
|
|
309406
|
0
|
5016703
|
if(!@_) {return NaND128()} |
|
1
|
|
|
|
|
5
|
|
434
|
|
|
|
|
|
|
|
435
|
309405
|
100
|
|
|
|
660697
|
if(@_ > 3) {die "More than 3 arguments supplied to new()"} |
|
2
|
|
|
|
|
17
|
|
436
|
|
|
|
|
|
|
|
437
|
|
|
|
|
|
|
# If 'new' has been called OOP style, the first arg is the string |
438
|
|
|
|
|
|
|
# "Math::Decimal128" which we don't need - so let's remove it. However, |
439
|
|
|
|
|
|
|
# if the first arg is a Math::Decimal128 object (which is a possibility), |
440
|
|
|
|
|
|
|
# then we'll get a fatal error when we check it for equivalence to |
441
|
|
|
|
|
|
|
# the string "Math::Decimal128". So we first need to check that it's not |
442
|
|
|
|
|
|
|
# an object - which we'll do by using the ref() function: |
443
|
309403
|
100
|
100
|
|
|
1266314
|
if(!ref($_[0]) && $_[0] eq "Math::Decimal128") { |
444
|
309392
|
|
|
|
|
479049
|
shift; |
445
|
309392
|
100
|
|
|
|
600062
|
if(!@_) {return NaND128()} |
|
11
|
|
|
|
|
216
|
|
446
|
|
|
|
|
|
|
} |
447
|
|
|
|
|
|
|
|
448
|
|
|
|
|
|
|
# @_ can now contain max of 2 vals - the mantissa and exponent. |
449
|
|
|
|
|
|
|
# If @_ == 1 then it contains the value. |
450
|
309392
|
100
|
|
|
|
635745
|
if(@_ > 2) {die "Too many arguments supplied to new() - expected no more than 2"} |
|
1
|
|
|
|
|
5
|
|
451
|
|
|
|
|
|
|
|
452
|
309391
|
100
|
|
|
|
602920
|
if(@_ == 2) {return MEtoD128(shift, shift)} |
|
309258
|
|
|
|
|
604993
|
|
453
|
|
|
|
|
|
|
|
454
|
133
|
|
|
|
|
233
|
my $arg = shift; |
455
|
133
|
|
|
|
|
364
|
my $type = _itsa($arg); |
456
|
|
|
|
|
|
|
|
457
|
133
|
100
|
|
|
|
324
|
if($type == 1) { # UV |
458
|
13
|
|
|
|
|
102
|
return UVtoD128($arg); |
459
|
|
|
|
|
|
|
} |
460
|
|
|
|
|
|
|
|
461
|
120
|
100
|
|
|
|
278
|
if($type == 2) { # IV |
462
|
15
|
|
|
|
|
100
|
return IVtoD128($arg); |
463
|
|
|
|
|
|
|
} |
464
|
|
|
|
|
|
|
|
465
|
105
|
100
|
|
|
|
238
|
if($type == 3) { # NV |
466
|
3
|
|
|
|
|
19
|
die "new() cannot be used to assign an NV - use NVtoD128() instead"; |
467
|
|
|
|
|
|
|
} |
468
|
|
|
|
|
|
|
|
469
|
102
|
100
|
|
|
|
235
|
if($type == 4) { # PV |
470
|
74
|
50
|
|
|
|
227
|
return STRtoD128($arg) if have_strtod128(); |
471
|
74
|
|
|
|
|
744
|
return PVtoD128($arg); |
472
|
|
|
|
|
|
|
} |
473
|
|
|
|
|
|
|
|
474
|
28
|
100
|
|
|
|
64
|
if($type == 34) { # Math::Decimal128 object |
475
|
25
|
|
|
|
|
132
|
return D128toD128($arg); |
476
|
|
|
|
|
|
|
} |
477
|
|
|
|
|
|
|
|
478
|
3
|
|
|
|
|
13
|
die "Bad argument given to new"; |
479
|
|
|
|
|
|
|
} |
480
|
|
|
|
|
|
|
|
481
|
|
|
|
|
|
|
####################################################################### |
482
|
|
|
|
|
|
|
####################################################################### |
483
|
|
|
|
|
|
|
|
484
|
|
|
|
|
|
|
sub D128toME { |
485
|
179923
|
|
|
179923
|
0
|
601542
|
my($ret1, $ret2) = split /e/i, decode_d128($_[0]); |
486
|
179923
|
100
|
|
|
|
478739
|
$ret2 = 0 unless defined $ret2; |
487
|
179923
|
100
|
|
|
|
798614
|
$ret2 = 0 if is_ZeroD128($_[0]); |
488
|
179923
|
|
|
|
|
501371
|
return ($ret1, $ret2); |
489
|
|
|
|
|
|
|
} |
490
|
|
|
|
|
|
|
|
491
|
|
|
|
|
|
|
####################################################################### |
492
|
|
|
|
|
|
|
####################################################################### |
493
|
|
|
|
|
|
|
|
494
|
|
|
|
|
|
|
sub MEtoD128 { |
495
|
|
|
|
|
|
|
# Check that 2 args are supplied |
496
|
373640
|
50
|
|
373640
|
0
|
38332955
|
die "MEtoD128 takes 2 args" if @_ != 2; |
497
|
|
|
|
|
|
|
|
498
|
373640
|
|
|
|
|
704475
|
my ($arg1, $arg2) = (shift, shift); |
499
|
|
|
|
|
|
|
|
500
|
373640
|
100
|
|
|
|
1051849
|
die "Invalid 1st arg ($arg1) to MEtoD128" if $arg1 =~ /[^0-9\-\+]/; |
501
|
373639
|
50
|
|
|
|
808835
|
die "Invalid 2nd arg ($arg2) to MEtoD128" if $arg2 =~ /[^0-9\-\+]/; |
502
|
|
|
|
|
|
|
|
503
|
373639
|
|
|
|
|
611185
|
$arg1 =~ s/^\+//; |
504
|
373639
|
|
|
|
|
521935
|
$arg2 =~ s/^\+//; |
505
|
373639
|
100
|
|
|
|
1019612
|
my $sign = $arg1 =~ s/^\-// ? '-' : ''; |
506
|
|
|
|
|
|
|
|
507
|
373639
|
|
|
|
|
600954
|
my $len_1 = length($arg1); |
508
|
|
|
|
|
|
|
|
509
|
373639
|
100
|
100
|
|
|
1441277
|
if($len_1 > 34 || $arg2 < -6176) { |
510
|
11279
|
100
|
|
|
|
23535
|
die "$arg1 exceeds _Decimal128 precision. It needs to be shortened to no more than 34 decimal digits" |
511
|
|
|
|
|
|
|
if $len_1 > 34; |
512
|
11278
|
|
|
|
|
22384
|
($arg1, $arg2) = _round_as_needed($arg1, $arg2); |
513
|
|
|
|
|
|
|
} |
514
|
|
|
|
|
|
|
|
515
|
|
|
|
|
|
|
# split $arg1 into segments that don't exceed 16 digits. |
516
|
|
|
|
|
|
|
|
517
|
373638
|
|
|
|
|
594254
|
my($msd, $nsd, $lsd); |
518
|
|
|
|
|
|
|
|
519
|
|
|
|
|
|
|
{ |
520
|
30
|
|
|
30
|
|
97329
|
no warnings 'substr'; |
|
30
|
|
|
|
|
75
|
|
|
30
|
|
|
|
|
13135
|
|
|
373638
|
|
|
|
|
511718
|
|
521
|
373638
|
|
100
|
|
|
1181881
|
$msd = substr($arg1, -34, 10) || '0'; |
522
|
373638
|
|
100
|
|
|
1015855
|
$nsd = substr($arg1, -24, 12) || '0'; |
523
|
373638
|
|
100
|
|
|
839734
|
$lsd = substr($arg1, -12, 12) || '0'; |
524
|
|
|
|
|
|
|
} |
525
|
|
|
|
|
|
|
|
526
|
373638
|
|
|
|
|
32456206
|
return _MEtoD128($sign . $msd, $sign . $nsd, $sign . $lsd, $arg2); |
527
|
|
|
|
|
|
|
|
528
|
|
|
|
|
|
|
} |
529
|
|
|
|
|
|
|
|
530
|
|
|
|
|
|
|
####################################################################### |
531
|
|
|
|
|
|
|
####################################################################### |
532
|
|
|
|
|
|
|
|
533
|
|
|
|
|
|
|
# Values such as (d, -6178), (dd, -6179), (ddd, -6180), etc evaluate to zero. |
534
|
|
|
|
|
|
|
# But values such as (dddd, -6178), (ddd, -6179), (dddddddd, -6180), etc may be non-zero. |
535
|
|
|
|
|
|
|
# In such cases we'll remove the ignored (trailing) digits, rounding the leading |
536
|
|
|
|
|
|
|
# digits to nearest - tied to even for midway cases. |
537
|
|
|
|
|
|
|
|
538
|
|
|
|
|
|
|
sub _round_as_needed { |
539
|
11278
|
|
|
11278
|
|
21842
|
my($sign, $man, $exp) = ('', shift, shift); |
540
|
|
|
|
|
|
|
|
541
|
11278
|
50
|
|
|
|
24022
|
if($man =~ /^\-/) { |
542
|
0
|
|
|
|
|
0
|
$man =~ s/^\-//; |
543
|
0
|
|
|
|
|
0
|
$sign = '-'; |
544
|
|
|
|
|
|
|
} |
545
|
|
|
|
|
|
|
|
546
|
11278
|
|
|
|
|
17180
|
my $length = length $man; |
547
|
11278
|
|
|
|
|
17255
|
my $maxlen = -6176 - $exp; |
548
|
|
|
|
|
|
|
|
549
|
11278
|
100
|
|
|
|
23302
|
if($length >= $maxlen) { |
550
|
4286
|
|
|
|
|
7883
|
my $rounder = substr($man, $length - $maxlen); # The trailing (ignored) digits |
551
|
4286
|
100
|
|
|
|
9478
|
$man = $length > $maxlen ? substr($man, 0, $length - $maxlen) |
552
|
|
|
|
|
|
|
: '0'; |
553
|
4286
|
|
|
|
|
5892
|
my $roundup = 0; |
554
|
4286
|
100
|
|
|
|
11157
|
$roundup = 1 if substr($rounder, 0, 1) > 5; |
555
|
4286
|
100
|
100
|
|
|
11523
|
$roundup = 1 if ((substr($rounder, 0, 1) == 5) && |
|
|
|
66
|
|
|
|
|
556
|
|
|
|
|
|
|
((substr($rounder, 1) =~ /[1-9]/) || (substr($man, -1, 1) %2 == 1))); |
557
|
|
|
|
|
|
|
|
558
|
4286
|
100
|
|
|
|
9297
|
$man++ if $roundup; |
559
|
4286
|
|
|
|
|
6909
|
$exp += $maxlen; # Removal of trailing digits moved the implied |
560
|
|
|
|
|
|
|
# decimal point $maxlen places to the left |
561
|
|
|
|
|
|
|
} |
562
|
|
|
|
|
|
|
|
563
|
11278
|
|
|
|
|
30540
|
return ($sign . $man, $exp); |
564
|
|
|
|
|
|
|
} |
565
|
|
|
|
|
|
|
|
566
|
|
|
|
|
|
|
####################################################################### |
567
|
|
|
|
|
|
|
####################################################################### |
568
|
|
|
|
|
|
|
|
569
|
|
|
|
|
|
|
sub assignMEl { |
570
|
|
|
|
|
|
|
# Check that 3 args are supplied |
571
|
13
|
50
|
|
13
|
0
|
513
|
die "assignMEl takes 3 args" if @_ != 3; |
572
|
|
|
|
|
|
|
|
573
|
13
|
|
|
|
|
22
|
my $arg1 = shift; |
574
|
13
|
|
|
|
|
20
|
my $arg2 = shift; |
575
|
13
|
|
|
|
|
20
|
my $arg3 = shift; |
576
|
|
|
|
|
|
|
|
577
|
13
|
100
|
|
|
|
52
|
die "Invalid 1st arg ($arg1) to assignMEl" if _itsa($arg1) != 34; |
578
|
11
|
50
|
|
|
|
33
|
die "Invalid 2nd arg ($arg2) to assignMEl" if $arg2 =~ /[^0-9\-]/; |
579
|
11
|
50
|
|
|
|
31
|
die "Invalid 3rd arg ($arg3) to assignMEl" if $arg3 =~ /[^0-9\-]/; |
580
|
|
|
|
|
|
|
|
581
|
11
|
|
|
|
|
20
|
my $len_2 = length($arg2); |
582
|
11
|
100
|
|
|
|
31
|
my $sign = $arg2 =~ /^\-/ ? '-' : ''; |
583
|
11
|
100
|
|
|
|
27
|
if($sign) { |
584
|
4
|
|
|
|
|
8
|
$len_2--; |
585
|
4
|
|
|
|
|
11
|
$arg2 =~ s/^\-//; |
586
|
|
|
|
|
|
|
} |
587
|
|
|
|
|
|
|
|
588
|
11
|
50
|
33
|
|
|
50
|
if($len_2 > 34 || $arg3 < -6176) { |
589
|
0
|
0
|
|
|
|
0
|
die "$arg2 exceeds _Decimal128 precision. It needs to be shortened to no more than 34 decimal digits" |
590
|
|
|
|
|
|
|
if $len_2 > 34; |
591
|
0
|
|
|
|
|
0
|
($arg2, $arg3) = _round_as_needed($arg2, $arg3); |
592
|
|
|
|
|
|
|
} |
593
|
|
|
|
|
|
|
|
594
|
11
|
|
|
|
|
19
|
my($msd, $nsd, $lsd); |
595
|
|
|
|
|
|
|
|
596
|
|
|
|
|
|
|
{ |
597
|
30
|
|
|
30
|
|
226
|
no warnings 'substr'; |
|
30
|
|
|
|
|
71
|
|
|
30
|
|
|
|
|
85685
|
|
|
11
|
|
|
|
|
17
|
|
598
|
11
|
|
100
|
|
|
38
|
$msd = substr($arg2, -34, 10) || '0'; |
599
|
11
|
|
100
|
|
|
40
|
$nsd = substr($arg2, -24, 12) || '0'; |
600
|
11
|
|
100
|
|
|
37
|
$lsd = substr($arg2, -12, 12) || '0'; |
601
|
|
|
|
|
|
|
} |
602
|
|
|
|
|
|
|
|
603
|
11
|
|
|
|
|
1460
|
return _assignME($arg1, $sign . $msd, $sign . $nsd, $sign . $lsd, $arg3); |
604
|
|
|
|
|
|
|
|
605
|
|
|
|
|
|
|
} |
606
|
|
|
|
|
|
|
|
607
|
|
|
|
|
|
|
####################################################################### |
608
|
|
|
|
|
|
|
####################################################################### |
609
|
|
|
|
|
|
|
|
610
|
|
|
|
|
|
|
sub _sci2me { |
611
|
0
|
|
|
0
|
|
0
|
my @ret = split /e/i, $_[0]; |
612
|
0
|
|
|
|
|
0
|
chop $ret[0] while $ret[0] =~ /0\b/; |
613
|
0
|
|
|
|
|
0
|
my @adj = split /\./, $ret[0]; |
614
|
0
|
0
|
|
|
|
0
|
my $adj = defined $adj[1] ? length($adj[1]) |
615
|
|
|
|
|
|
|
: 0; |
616
|
0
|
|
|
|
|
0
|
$ret[0] =~ s/\.//; |
617
|
0
|
|
|
|
|
0
|
$ret[1] += $_[1] - $adj; |
618
|
|
|
|
|
|
|
|
619
|
0
|
|
|
|
|
0
|
return @ret; |
620
|
|
|
|
|
|
|
} |
621
|
|
|
|
|
|
|
|
622
|
|
|
|
|
|
|
####################################################################### |
623
|
|
|
|
|
|
|
####################################################################### |
624
|
|
|
|
|
|
|
|
625
|
8
|
|
|
8
|
|
137
|
sub DEC128_MAX {return _DEC128_MAX()} |
626
|
9
|
|
|
9
|
|
96
|
sub DEC128_MIN {return _DEC128_MIN()} |
627
|
|
|
|
|
|
|
|
628
|
|
|
|
|
|
|
####################################################################### |
629
|
|
|
|
|
|
|
####################################################################### |
630
|
|
|
|
|
|
|
|
631
|
|
|
|
|
|
|
sub d128_bytes { |
632
|
199416
|
|
|
199416
|
0
|
1714689
|
my @ret = _d128_bytes($_[0]); |
633
|
199416
|
|
|
|
|
876408
|
return join '', @ret; |
634
|
|
|
|
|
|
|
} |
635
|
|
|
|
|
|
|
|
636
|
|
|
|
|
|
|
####################################################################### |
637
|
|
|
|
|
|
|
####################################################################### |
638
|
|
|
|
|
|
|
|
639
|
|
|
|
|
|
|
sub hex2binl { |
640
|
194480
|
|
|
194480
|
0
|
773653
|
my $ret = unpack("B*", (pack "H*", $_[0])); |
641
|
194480
|
|
|
|
|
323503
|
my $len = length $ret; |
642
|
194480
|
50
|
|
|
|
419131
|
die "hex2binl() yielded $len bits" if $len != 128; |
643
|
194480
|
|
|
|
|
362661
|
return $ret; |
644
|
|
|
|
|
|
|
} |
645
|
|
|
|
|
|
|
|
646
|
|
|
|
|
|
|
####################################################################### |
647
|
|
|
|
|
|
|
####################################################################### |
648
|
|
|
|
|
|
|
|
649
|
|
|
|
|
|
|
sub d128_fmt { |
650
|
37
|
|
|
37
|
0
|
325
|
my $d128 = MEtoD128('1234567890123456789012345678901234', 0); |
651
|
|
|
|
|
|
|
# DPD: 2608134B9C1E28E56F3C127177823534 |
652
|
|
|
|
|
|
|
# BID: 30403CDE6FFF9732DE825CD07E96AFF2 |
653
|
37
|
50
|
|
|
|
164
|
return 'DPD' if d128_bytes($d128) =~ /534$/i; |
654
|
37
|
50
|
|
|
|
137
|
return 'BID' if d128_bytes($d128) =~ /FF2$/i; |
655
|
0
|
|
|
|
|
0
|
return 'Unknown'; |
656
|
|
|
|
|
|
|
} |
657
|
|
|
|
|
|
|
|
658
|
|
|
|
|
|
|
####################################################################### |
659
|
|
|
|
|
|
|
####################################################################### |
660
|
|
|
|
|
|
|
|
661
|
|
|
|
|
|
|
sub decode_dpdl { |
662
|
|
|
|
|
|
|
# Takes the Math::Decimal128 object as its arg. |
663
|
|
|
|
|
|
|
# Decodes Densely Packed Decimal formatting of the Decimal128 value. |
664
|
0
|
|
|
0
|
0
|
0
|
my $binstring = hex2binl(d128_bytes($_[0])); |
665
|
0
|
|
|
|
|
0
|
my @first = decode_dpd_1st($binstring); |
666
|
0
|
0
|
0
|
|
|
0
|
return ($first[0] . $first[1]) if ($first[1] =~ /inf/i || $first[1] =~ /nan/i); |
667
|
0
|
|
|
|
|
0
|
my $mantissa = $first[1] . decode_dpd_2nd($binstring); |
668
|
|
|
|
|
|
|
|
669
|
|
|
|
|
|
|
# Remove leading zeroes from the mantissa |
670
|
0
|
|
|
|
|
0
|
$mantissa =~ s/^0+//; |
671
|
0
|
0
|
|
|
|
0
|
if($mantissa eq '') {$mantissa = '0'} |
|
0
|
|
|
|
|
0
|
|
672
|
|
|
|
|
|
|
else { |
673
|
|
|
|
|
|
|
# Remove trailing zeroes |
674
|
0
|
|
|
|
|
0
|
while($mantissa =~ /0$/) { |
675
|
0
|
|
|
|
|
0
|
$mantissa =~ s/0$//; |
676
|
0
|
|
|
|
|
0
|
$first[2]++; |
677
|
|
|
|
|
|
|
} |
678
|
|
|
|
|
|
|
} |
679
|
|
|
|
|
|
|
|
680
|
0
|
|
|
|
|
0
|
my $ret = $first[0] . $mantissa . 'e' . $first[2]; |
681
|
|
|
|
|
|
|
|
682
|
|
|
|
|
|
|
} |
683
|
|
|
|
|
|
|
|
684
|
|
|
|
|
|
|
####################################################################### |
685
|
|
|
|
|
|
|
####################################################################### |
686
|
|
|
|
|
|
|
|
687
|
|
|
|
|
|
|
sub decode_dpd_1st{ |
688
|
|
|
|
|
|
|
# Takes the entire binary string as its arg. |
689
|
0
|
0
|
|
0
|
0
|
0
|
die "Argument to decode_dpd_1st is wrong size (", length($_[0]), ")" |
690
|
|
|
|
|
|
|
if length($_[0]) != 128; |
691
|
0
|
|
|
|
|
0
|
my $leading_bits = 18; |
692
|
0
|
|
|
|
|
0
|
my $trailing_bits = 110; |
693
|
0
|
|
|
|
|
0
|
my $msd; # significand's most siginificant digit |
694
|
|
|
|
|
|
|
my $exp; # exponent |
695
|
0
|
|
|
|
|
0
|
my $keep = substr($_[0], 0, $leading_bits); |
696
|
0
|
0
|
|
|
|
0
|
my $sign = substr($keep, 0, 1) ? '-' : ''; |
697
|
0
|
0
|
|
|
|
0
|
return ('','nan') if substr($keep, 1, 5) eq '11111'; |
698
|
0
|
0
|
|
|
|
0
|
if(substr($keep, 1 ,5) eq '11110') {return ($sign, 'inf')} |
|
0
|
|
|
|
|
0
|
|
699
|
0
|
|
|
|
|
0
|
my $pre = substr($keep, 1, 2); |
700
|
0
|
0
|
0
|
|
|
0
|
if($pre eq '00' || $pre eq '01' || $pre eq '10') { |
|
|
|
0
|
|
|
|
|
701
|
0
|
|
|
|
|
0
|
$msd = oct('0b0' . substr($keep, 3, 3)); |
702
|
0
|
|
|
|
|
0
|
$exp = oct('0b' . $pre . substr($keep, 6, 12)) - 6176; |
703
|
0
|
|
|
|
|
0
|
return ($sign, $msd, $exp); |
704
|
|
|
|
|
|
|
} |
705
|
0
|
|
|
|
|
0
|
$pre = substr($keep, 1, 4); |
706
|
0
|
0
|
0
|
|
|
0
|
if($pre eq '1100' || $pre eq '1101' || $pre eq '1110') { |
|
|
|
0
|
|
|
|
|
707
|
0
|
|
|
|
|
0
|
$exp = oct('0b' . substr($pre, 2, 2) . substr($keep, 6, 12)) - 6176; |
708
|
0
|
|
|
|
|
0
|
$msd = oct('0b' . '100' . substr($keep, 5, 1)); |
709
|
0
|
|
|
|
|
0
|
return ($sign, $msd, $exp); |
710
|
|
|
|
|
|
|
} |
711
|
0
|
|
|
|
|
0
|
die "decode_dpd_1st function failed to parse its argument ($_[0])"; |
712
|
|
|
|
|
|
|
} |
713
|
|
|
|
|
|
|
|
714
|
|
|
|
|
|
|
####################################################################### |
715
|
|
|
|
|
|
|
####################################################################### |
716
|
|
|
|
|
|
|
|
717
|
|
|
|
|
|
|
sub decode_dpd_2nd { |
718
|
|
|
|
|
|
|
# Takes the entire binary string as its arg. |
719
|
0
|
0
|
|
0
|
0
|
0
|
die "Argument to decode_dpd_2nd is wrong size (", length($_[0]), ")" |
720
|
|
|
|
|
|
|
if length($_[0]) != 128; |
721
|
0
|
|
|
|
|
0
|
my $leading_bits = 18; |
722
|
0
|
|
|
|
|
0
|
my $trailing_bits = 110; |
723
|
0
|
|
|
|
|
0
|
my $keep = substr($_[0], $leading_bits, $trailing_bits); |
724
|
0
|
|
|
|
|
0
|
my $ret = ''; |
725
|
0
|
|
|
|
|
0
|
for my $i(0, 10, 20, 30, 40, 50, 60, 70, 80, 90, 100) { |
726
|
0
|
|
|
|
|
0
|
my $key = substr($keep, $i, 10); |
727
|
0
|
|
|
|
|
0
|
$ret .= $Math::Decimal128::dpd_encode{$key}; |
728
|
|
|
|
|
|
|
} |
729
|
0
|
|
|
|
|
0
|
return $ret; |
730
|
|
|
|
|
|
|
} |
731
|
|
|
|
|
|
|
|
732
|
|
|
|
|
|
|
####################################################################### |
733
|
|
|
|
|
|
|
####################################################################### |
734
|
|
|
|
|
|
|
|
735
|
|
|
|
|
|
|
sub decode_bidl { |
736
|
|
|
|
|
|
|
# Takes a Math::Decimal128 object as its arg. |
737
|
|
|
|
|
|
|
# Decodes Binary Integer Decimal formatting of the _Decimal128 value. |
738
|
|
|
|
|
|
|
|
739
|
184756
|
|
|
184756
|
0
|
369162
|
my $keep = hex2binl(d128_bytes($_[0])); |
740
|
184756
|
50
|
|
|
|
440862
|
die "Base 2 representation is wrong size (", length($keep), ")" |
741
|
|
|
|
|
|
|
if length($keep) != 128; |
742
|
184756
|
|
|
|
|
259314
|
my $leading_bits = 17; |
743
|
184756
|
|
|
|
|
250045
|
my $trailing_bits = 111; |
744
|
184756
|
|
|
|
|
276456
|
my @mantissa; |
745
|
|
|
|
|
|
|
my $exp; # exponent |
746
|
184756
|
100
|
|
|
|
388214
|
my $sign = substr($keep, 0, 1) ? '-' : ''; |
747
|
184756
|
100
|
|
|
|
413163
|
return 'nan' if substr($keep, 1, 5) eq '11111'; |
748
|
184743
|
100
|
|
|
|
402270
|
if(substr($keep, 1 ,5) eq '11110') {return $sign . 'inf'} |
|
25
|
|
|
|
|
135
|
|
749
|
184718
|
|
|
|
|
283591
|
my $pre = substr($keep, 1, 2); |
750
|
184718
|
50
|
100
|
|
|
758849
|
if($pre eq '00' || $pre eq '01' || $pre eq '10') { |
|
|
|
66
|
|
|
|
|
751
|
184718
|
|
|
|
|
414326
|
$exp = oct('0b' . substr($keep, 1, 14)) - 6176; |
752
|
184718
|
|
|
|
|
4059812
|
@mantissa = reverse(split(//, '0' . substr($keep, 15, 113))); |
753
|
184718
|
|
|
|
|
3700474
|
my $m = _bid_mant(\@mantissa); # $m is a Math::Decimal128 object |
754
|
184718
|
|
|
|
|
427827
|
my $mantissa = _decode_mant($m); |
755
|
184718
|
|
|
|
|
602931
|
$mantissa =~ s/^0+//; |
756
|
184718
|
100
|
|
|
|
600128
|
if($mantissa !~ /[1-9]/) { $mantissa = '0'} |
|
925
|
|
|
|
|
1613
|
|
757
|
|
|
|
|
|
|
else { |
758
|
183793
|
|
|
|
|
524082
|
while($mantissa =~ /0$/) { |
759
|
2198318
|
|
|
|
|
4283603
|
$mantissa =~ s/0$//; |
760
|
2198318
|
|
|
|
|
5193366
|
$exp++; |
761
|
|
|
|
|
|
|
} |
762
|
|
|
|
|
|
|
} |
763
|
184718
|
|
|
|
|
1802439
|
return $sign . $mantissa . 'e' . $exp; |
764
|
|
|
|
|
|
|
} |
765
|
0
|
|
|
|
|
0
|
$pre = substr($keep, 1, 4); |
766
|
0
|
0
|
0
|
|
|
0
|
if($pre eq '1100' || $pre eq '1101' || $pre eq '1110') { |
|
|
|
0
|
|
|
|
|
767
|
0
|
|
|
|
|
0
|
$exp = oct('0b' . substr($keep, 3, 14)) - 6176; |
768
|
0
|
|
|
|
|
0
|
@mantissa = reverse(split(//,'100' . substr($keep, 17, 111))); |
769
|
0
|
|
|
|
|
0
|
my $m = _bid_mant(\@mantissa); # $m is a Math::Decimal128 object |
770
|
0
|
|
|
|
|
0
|
my $mantissa = _decode_mant($m); |
771
|
0
|
|
|
|
|
0
|
$mantissa =~ s/^0+//; |
772
|
0
|
0
|
|
|
|
0
|
if($mantissa !~ /[1-9]/) { $mantissa = '0'} |
|
0
|
|
|
|
|
0
|
|
773
|
|
|
|
|
|
|
else { |
774
|
0
|
|
|
|
|
0
|
while($mantissa =~ /0$/) { |
775
|
0
|
|
|
|
|
0
|
$mantissa =~ s/0$//; |
776
|
0
|
|
|
|
|
0
|
$exp++; |
777
|
|
|
|
|
|
|
} |
778
|
|
|
|
|
|
|
} |
779
|
0
|
|
|
|
|
0
|
return $sign . $mantissa . 'e' . $exp; |
780
|
|
|
|
|
|
|
} |
781
|
0
|
|
|
|
|
0
|
die "decode_bid function failed to parse its argument ($_[0])"; |
782
|
|
|
|
|
|
|
} |
783
|
|
|
|
|
|
|
|
784
|
|
|
|
|
|
|
####################################################################### |
785
|
|
|
|
|
|
|
####################################################################### |
786
|
|
|
|
|
|
|
# Now using XSub of same name |
787
|
|
|
|
|
|
|
# sub PVtoD128 { |
788
|
|
|
|
|
|
|
# |
789
|
|
|
|
|
|
|
# my($arg1, $arg2) = PVtoMEl($_[0]); |
790
|
|
|
|
|
|
|
# |
791
|
|
|
|
|
|
|
# if($arg1 =~ /inf|nan/i) { |
792
|
|
|
|
|
|
|
# $arg1 =~ /nan/i ? return NaND128() |
793
|
|
|
|
|
|
|
# : $arg1 =~ /^\-/ ? return InfD128(-1) |
794
|
|
|
|
|
|
|
# : return InfD128(1); |
795
|
|
|
|
|
|
|
# } |
796
|
|
|
|
|
|
|
# |
797
|
|
|
|
|
|
|
# return MEtoD128($arg1, $arg2); |
798
|
|
|
|
|
|
|
#} |
799
|
|
|
|
|
|
|
|
800
|
|
|
|
|
|
|
####################################################################### |
801
|
|
|
|
|
|
|
####################################################################### |
802
|
|
|
|
|
|
|
# Now using XSub of same name |
803
|
|
|
|
|
|
|
# sub assignPVl { |
804
|
|
|
|
|
|
|
# |
805
|
|
|
|
|
|
|
# my($arg1, $arg2) = PVtoMEl($_[1]); |
806
|
|
|
|
|
|
|
# if($arg1 =~ /inf|nan/i) { |
807
|
|
|
|
|
|
|
# $arg1 =~ /nan/i ? assignNaNl($_[0]) |
808
|
|
|
|
|
|
|
# : $arg1 =~ /^\-/ ? assignInfl($_[0], -1) |
809
|
|
|
|
|
|
|
# : assignInfl($_[0], 1); |
810
|
|
|
|
|
|
|
# } |
811
|
|
|
|
|
|
|
# else { |
812
|
|
|
|
|
|
|
# assignMEl($_[0], $arg1, $arg2); |
813
|
|
|
|
|
|
|
# } |
814
|
|
|
|
|
|
|
#} |
815
|
|
|
|
|
|
|
|
816
|
|
|
|
|
|
|
####################################################################### |
817
|
|
|
|
|
|
|
####################################################################### |
818
|
|
|
|
|
|
|
|
819
|
|
|
|
|
|
|
sub PVtoMEl { |
820
|
|
|
|
|
|
|
|
821
|
0
|
|
|
0
|
0
|
0
|
my($arg1, $arg2) = split /e/i, $_[0]; |
822
|
|
|
|
|
|
|
|
823
|
0
|
0
|
|
|
|
0
|
if($arg1 =~ /^(\-|\+)?inf|^(\-|\+)?nan/i) { |
824
|
0
|
|
|
|
|
0
|
return ($arg1, 0); |
825
|
|
|
|
|
|
|
} |
826
|
|
|
|
|
|
|
|
827
|
0
|
|
|
|
|
0
|
_sanitise_args($arg1, $arg2); |
828
|
0
|
|
|
|
|
0
|
return ($arg1, $arg2); |
829
|
|
|
|
|
|
|
} |
830
|
|
|
|
|
|
|
|
831
|
|
|
|
|
|
|
####################################################################### |
832
|
|
|
|
|
|
|
####################################################################### |
833
|
|
|
|
|
|
|
|
834
|
|
|
|
|
|
|
sub MEtoPVl { |
835
|
16
|
|
|
16
|
0
|
143
|
my $arg1 = shift; |
836
|
16
|
100
|
|
|
|
56
|
if($arg1 =~ /^(\-|\+)?inf|^(\-|\+)?nan/i) { |
837
|
8
|
|
|
|
|
21
|
$arg1 =~ s/\+//; |
838
|
8
|
|
|
|
|
18
|
return $arg1; |
839
|
|
|
|
|
|
|
} |
840
|
|
|
|
|
|
|
|
841
|
8
|
|
|
|
|
12
|
my $arg2 = shift; |
842
|
8
|
|
|
|
|
18
|
return $arg1 . 'e' . $arg2; |
843
|
|
|
|
|
|
|
} |
844
|
|
|
|
|
|
|
|
845
|
|
|
|
|
|
|
####################################################################### |
846
|
|
|
|
|
|
|
####################################################################### |
847
|
|
|
|
|
|
|
|
848
|
|
|
|
|
|
|
sub _sanitise_args { |
849
|
0
|
0
|
|
0
|
|
0
|
$_[1] = 0 unless defined $_[1]; |
850
|
0
|
|
|
|
|
0
|
$_[0] =~ s/\.0+$//; |
851
|
0
|
|
|
|
|
0
|
my @split = split /\./, $_[0]; |
852
|
0
|
0
|
|
|
|
0
|
$split[1] = '' unless defined $split[1]; |
853
|
0
|
|
|
|
|
0
|
$_[1] -= length($split[1]); |
854
|
0
|
|
|
|
|
0
|
$_[0] =~ s/\.//; |
855
|
0
|
|
|
|
|
0
|
$_[0] =~ s/^0+//; |
856
|
|
|
|
|
|
|
} |
857
|
|
|
|
|
|
|
|
858
|
|
|
|
|
|
|
####################################################################### |
859
|
|
|
|
|
|
|
####################################################################### |
860
|
|
|
|
|
|
|
|
861
|
|
|
|
|
|
|
sub get_expl { |
862
|
9724
|
|
|
9724
|
0
|
63483
|
my $keep = hex2binl(d128_bytes($_[0])); |
863
|
9724
|
|
|
|
|
19192
|
my $pre = substr($keep, 1, 2); |
864
|
9724
|
50
|
|
|
|
20615
|
if($Math::Decimal128::fmt eq 'DPD') { |
865
|
0
|
0
|
0
|
|
|
0
|
if($pre eq '00' || $pre eq '01' || $pre eq '10') { |
|
|
|
0
|
|
|
|
|
866
|
0
|
|
|
|
|
0
|
return oct('0b' . $pre . substr($keep, 6, 12)) - 6176; |
867
|
|
|
|
|
|
|
} |
868
|
|
|
|
|
|
|
else { |
869
|
0
|
|
|
|
|
0
|
return oct('0b' . substr($pre, 2, 2) . substr($keep, 6, 12)) - 6176; |
870
|
|
|
|
|
|
|
} |
871
|
|
|
|
|
|
|
} |
872
|
|
|
|
|
|
|
else { |
873
|
9724
|
50
|
100
|
|
|
40082
|
if($pre eq '00' || $pre eq '01' || $pre eq '10') { |
|
|
|
66
|
|
|
|
|
874
|
9724
|
|
|
|
|
29646
|
return oct('0b' . substr($keep, 1, 14)) - 6176; |
875
|
|
|
|
|
|
|
} |
876
|
|
|
|
|
|
|
else { |
877
|
0
|
|
|
|
|
0
|
return oct('0b' . substr($keep, 3, 14)) - 6176; |
878
|
|
|
|
|
|
|
} |
879
|
|
|
|
|
|
|
} |
880
|
|
|
|
|
|
|
} |
881
|
|
|
|
|
|
|
|
882
|
|
|
|
|
|
|
####################################################################### |
883
|
|
|
|
|
|
|
####################################################################### |
884
|
|
|
|
|
|
|
|
885
|
|
|
|
|
|
|
sub get_signl { |
886
|
4862
|
100
|
|
4862
|
0
|
20487
|
return '-' if hex(substr(d128_bytes($_[0]), 0, 1)) >= 8; |
887
|
2448
|
|
|
|
|
6392
|
return '+'; |
888
|
|
|
|
|
|
|
} |
889
|
|
|
|
|
|
|
|
890
|
|
|
|
|
|
|
####################################################################### |
891
|
|
|
|
|
|
|
####################################################################### |
892
|
|
|
|
|
|
|
|
893
|
|
|
|
|
|
|
sub DPDtoD128 { |
894
|
|
|
|
|
|
|
# Usable only where DPD format is in use. |
895
|
|
|
|
|
|
|
# Converts the 128-bit string returned by _MEtoBINSTR into |
896
|
|
|
|
|
|
|
# a Math::Decimal128 object set to the value encoded by the |
897
|
|
|
|
|
|
|
# the 128-bit string. This is all done without having to calculate |
898
|
|
|
|
|
|
|
# the actual value - and is typically ~25 times quicker than |
899
|
|
|
|
|
|
|
# MEtoD128. |
900
|
6
|
|
|
6
|
0
|
156
|
my($man, $exp) = (shift, shift); |
901
|
6
|
|
|
|
|
13
|
my $arg = _MEtoBINSTR($man, $exp); |
902
|
6
|
|
|
|
|
35
|
return _DPDtoD128(unpack("a*", pack( "B*", $arg))); |
903
|
|
|
|
|
|
|
} |
904
|
|
|
|
|
|
|
|
905
|
|
|
|
|
|
|
####################################################################### |
906
|
|
|
|
|
|
|
####################################################################### |
907
|
|
|
|
|
|
|
|
908
|
|
|
|
|
|
|
sub assignDPDl { |
909
|
6
|
|
|
6
|
0
|
68
|
_assignDPD($_[0], unpack("a*", pack("B*", _MEtoBINSTR($_[1], $_[2])))); |
910
|
|
|
|
|
|
|
} |
911
|
|
|
|
|
|
|
|
912
|
|
|
|
|
|
|
####################################################################### |
913
|
|
|
|
|
|
|
####################################################################### |
914
|
|
|
|
|
|
|
|
915
|
|
|
|
|
|
|
sub _MEtoBINSTR { |
916
|
|
|
|
|
|
|
# Converts (mantissa, exponent) strings to DPD encoded 128-bit string - without |
917
|
|
|
|
|
|
|
# the need to actually calculate the value. |
918
|
12
|
|
|
12
|
|
23
|
my $man = shift; |
919
|
|
|
|
|
|
|
|
920
|
12
|
50
|
|
|
|
54
|
if($man =~ /^(\-|\+)?inf|^(\-|\+)?nan/i) { |
921
|
12
|
100
|
|
|
|
75
|
$man =~ /\-inf/i ? return '11111' . ('0' x 123) |
|
|
100
|
|
|
|
|
|
922
|
|
|
|
|
|
|
: $man =~ /^(\-|\+)?nan/i ? return '011111' . ('0' x 122) |
923
|
|
|
|
|
|
|
: return '01111' . ('0' x 123); |
924
|
|
|
|
|
|
|
} |
925
|
|
|
|
|
|
|
|
926
|
0
|
|
|
|
|
0
|
my $exp = shift; |
927
|
|
|
|
|
|
|
|
928
|
|
|
|
|
|
|
# Determine the sign, and remove it. |
929
|
0
|
0
|
|
|
|
0
|
my $sign = $man =~ /^\-/ ? '1' : '0'; |
930
|
0
|
|
|
|
|
0
|
$man =~ s/[\+\-]//; |
931
|
0
|
0
|
|
|
|
0
|
die "_MEtoBINSTR has been passed (probably from DPDtoBINSTR) an illegal mantissa" |
932
|
|
|
|
|
|
|
if $man =~ /[^0-9]/; |
933
|
|
|
|
|
|
|
|
934
|
|
|
|
|
|
|
# Remove leading zeroes, and return zero (of appropriate sign) |
935
|
|
|
|
|
|
|
# if we're left with the empty string. |
936
|
0
|
|
|
|
|
0
|
$man =~ s/^0+//; |
937
|
0
|
0
|
|
|
|
0
|
return $sign . '01000011111111111' . ('0' x 110) unless $man; |
938
|
|
|
|
|
|
|
|
939
|
|
|
|
|
|
|
# Fill the mantissa with 34 digits - by zero padding the end. |
940
|
0
|
|
|
|
|
0
|
my $add_zeroes = 34 - length($man); |
941
|
0
|
|
|
|
|
0
|
$man .= '0' x $add_zeroes; |
942
|
0
|
|
|
|
|
0
|
$exp -= $add_zeroes; |
943
|
|
|
|
|
|
|
|
944
|
0
|
0
|
0
|
|
|
0
|
if(length($man) > 34 || $exp < -6176) { |
945
|
0
|
0
|
|
|
|
0
|
die "$man exceeds _Decimal128 precision. It needs to be shortened to no more than 34 decimal digits" |
946
|
|
|
|
|
|
|
if length($man) > 34; |
947
|
0
|
|
|
|
|
0
|
($man, $exp) = _round_as_needed($man, $exp); |
948
|
|
|
|
|
|
|
} |
949
|
|
|
|
|
|
|
|
950
|
|
|
|
|
|
|
# Return 0 if $exp is still less that -6176. |
951
|
0
|
0
|
|
|
|
0
|
return $sign . '01000011111111111' . ('0' x 110) if $exp < -6176; |
952
|
|
|
|
|
|
|
|
953
|
|
|
|
|
|
|
# Return -inf/inf if value is infinite |
954
|
0
|
0
|
|
|
|
0
|
if($exp > 6111) { |
955
|
0
|
0
|
|
|
|
0
|
return $sign . '1111' . ('0' x 123) if (length($man) + $exp) > 6145; |
956
|
|
|
|
|
|
|
} |
957
|
|
|
|
|
|
|
|
958
|
0
|
|
|
|
|
0
|
$man = '0' . $man while length($man) < 34; |
959
|
|
|
|
|
|
|
|
960
|
|
|
|
|
|
|
# The last 110 bits encode the last 33 digits. |
961
|
0
|
|
|
|
|
0
|
my $last_33_dig = substr($man, 1, 33); |
962
|
0
|
|
|
|
|
0
|
my $last_110_bits; |
963
|
0
|
|
|
|
|
0
|
for(my $i = 0; $i < 31; $i += 3) { |
964
|
0
|
|
|
|
|
0
|
$last_110_bits .= $Math::Decimal128::dpd_decode{substr($last_33_dig, $i, 3)} |
965
|
|
|
|
|
|
|
} |
966
|
|
|
|
|
|
|
|
967
|
0
|
|
|
|
|
0
|
my $len = length($last_110_bits); |
968
|
0
|
0
|
|
|
|
0
|
die "Wrong bitsize ($len != 110) in MEtoBINSTR()" if $len != 110; |
969
|
|
|
|
|
|
|
|
970
|
0
|
|
|
|
|
0
|
my $leading_digit = substr($man, 0, 1); # ie the msd (most siginificant digit). |
971
|
0
|
|
|
|
|
0
|
my $exp_base_2 = sprintf "%014b", $exp + 6176; |
972
|
|
|
|
|
|
|
|
973
|
|
|
|
|
|
|
# The encoding of the exponent and msd depends upon the value of the msd. |
974
|
|
|
|
|
|
|
# If it's 0..7, it's done one way; if it's 8 or 9 it's done th'other way. |
975
|
0
|
0
|
|
|
|
0
|
if($leading_digit < 8) { |
976
|
0
|
|
|
|
|
0
|
my $leading_digit_bits = sprintf "%03b", $leading_digit; |
977
|
0
|
|
|
|
|
0
|
substr($exp_base_2, 2, 0, $leading_digit_bits); |
978
|
|
|
|
|
|
|
} |
979
|
|
|
|
|
|
|
else { |
980
|
0
|
0
|
|
|
|
0
|
my $leading_digit_bit = $leading_digit == 8 ? '0' : '1'; |
981
|
0
|
|
|
|
|
0
|
$exp_base_2 = '11' . substr($exp_base_2, 0, 2) . $leading_digit_bit . substr($exp_base_2, 2, 12); |
982
|
|
|
|
|
|
|
} |
983
|
|
|
|
|
|
|
|
984
|
0
|
|
|
|
|
0
|
$len = length($exp_base_2); |
985
|
0
|
0
|
|
|
|
0
|
die "Exponent component length is wrong ($len != 17) in MEtoBINSTR()" if $len != 17; |
986
|
|
|
|
|
|
|
|
987
|
0
|
|
|
|
|
0
|
return $sign . $exp_base_2 . $last_110_bits; |
988
|
|
|
|
|
|
|
} |
989
|
|
|
|
|
|
|
|
990
|
|
|
|
|
|
|
####################################################################### |
991
|
|
|
|
|
|
|
####################################################################### |
992
|
|
|
|
|
|
|
|
993
|
|
|
|
|
|
|
sub D128toFSTR { |
994
|
|
|
|
|
|
|
# Converts the argument (M::D128 object) to a string in floating point |
995
|
|
|
|
|
|
|
# format - as distinct from scientific notation. |
996
|
1741
|
|
|
1741
|
0
|
7191
|
my($m, $e) = D128toME($_[0]); |
997
|
1741
|
100
|
|
|
|
5906
|
return 'nan' if is_NaND128($_[0]); |
998
|
1730
|
100
|
|
|
|
6516
|
if(is_InfD128($_[0])) { |
999
|
22
|
100
|
|
|
|
113
|
return 'inf' if is_InfD128($_[0])> 0; |
1000
|
11
|
|
|
|
|
47
|
return '-inf'; |
1001
|
|
|
|
|
|
|
} |
1002
|
1708
|
100
|
|
|
|
6413
|
return $m . '0' x $e if $e >= 0; |
1003
|
1062
|
|
|
|
|
2049
|
my($len, $sign) = (length $m, ''); |
1004
|
1062
|
|
|
|
|
2509
|
$m =~ s/^\-//; |
1005
|
1062
|
100
|
|
|
|
2577
|
if($len != length $m) { |
1006
|
534
|
|
|
|
|
815
|
$len--; |
1007
|
534
|
|
|
|
|
771
|
$sign = '-'; |
1008
|
|
|
|
|
|
|
} |
1009
|
1062
|
100
|
|
|
|
2434
|
if($len + $e > 0) { |
1010
|
273
|
|
|
|
|
580
|
substr($m, $e, 0, '.'); |
1011
|
273
|
|
|
|
|
758
|
return $sign . $m; |
1012
|
|
|
|
|
|
|
} |
1013
|
789
|
100
|
|
|
|
1822
|
if($len + $e < 0) { |
1014
|
752
|
|
|
|
|
2652
|
return $sign . '0.' . '0' x -($len + $e) . $m; |
1015
|
|
|
|
|
|
|
} |
1016
|
37
|
|
|
|
|
113
|
return $sign . '0.' . $m; |
1017
|
|
|
|
|
|
|
} |
1018
|
|
|
|
|
|
|
|
1019
|
|
|
|
|
|
|
####################################################################### |
1020
|
|
|
|
|
|
|
####################################################################### |
1021
|
|
|
|
|
|
|
|
1022
|
|
|
|
|
|
|
sub D128toRSTR { |
1023
|
|
|
|
|
|
|
# As for D128toFSTR, but rounds the string to the no. of |
1024
|
|
|
|
|
|
|
# decimal places specified by the second arg. |
1025
|
|
|
|
|
|
|
|
1026
|
123
|
100
|
|
123
|
0
|
935
|
die "2nd arg to D128toRSTR() must be greater than zero" |
1027
|
|
|
|
|
|
|
unless $_[1] >= 0; |
1028
|
|
|
|
|
|
|
|
1029
|
122
|
100
|
|
|
|
294
|
my $dp = $_[1] ? '.' : ''; |
1030
|
|
|
|
|
|
|
|
1031
|
122
|
|
|
|
|
290
|
my $str = D128toFSTR($_[0]); |
1032
|
|
|
|
|
|
|
|
1033
|
122
|
100
|
|
|
|
478
|
return $str if $str =~ /n/i; # inf/nan |
1034
|
|
|
|
|
|
|
|
1035
|
92
|
100
|
|
|
|
312
|
return $str . "$dp" . '0' x $_[1] unless $str =~ /\./; |
1036
|
|
|
|
|
|
|
|
1037
|
84
|
|
|
|
|
265
|
my($leading, $trailing) = split /\./, $str; |
1038
|
84
|
|
|
|
|
167
|
my $len_trail = length $trailing ; |
1039
|
|
|
|
|
|
|
|
1040
|
84
|
100
|
|
|
|
298
|
return $str if ($_[1] == $len_trail); |
1041
|
|
|
|
|
|
|
|
1042
|
66
|
100
|
|
|
|
175
|
if(length($trailing) <= $_[1]) { |
1043
|
14
|
|
|
|
|
37
|
$trailing .= '0' x ($_[1] - length($trailing)); |
1044
|
14
|
|
|
|
|
54
|
return $leading . "$dp" . $trailing; |
1045
|
|
|
|
|
|
|
} |
1046
|
|
|
|
|
|
|
|
1047
|
|
|
|
|
|
|
# $len_trail > specified number of decimal places ($_[1]). |
1048
|
|
|
|
|
|
|
# We need to round (to nearest, ties to even) from here on. |
1049
|
|
|
|
|
|
|
|
1050
|
52
|
100
|
100
|
|
|
365
|
return $leading . "$dp" . substr($trailing, 0, $_[1]) |
|
|
|
66
|
|
|
|
|
1051
|
|
|
|
|
|
|
if (substr($trailing, $_[1], 1) <= 4) || |
1052
|
|
|
|
|
|
|
(substr($trailing, $_[1]) =~ /^5(0+)?$/ && substr($trailing, $_[1] - 1, 1) % 2 == 0); |
1053
|
|
|
|
|
|
|
|
1054
|
34
|
|
|
|
|
72
|
my $to_inc = substr($trailing, 0, $_[1]); |
1055
|
|
|
|
|
|
|
|
1056
|
34
|
|
|
|
|
83
|
my $carry = _increment($to_inc); # $carry will either be mt string or '1'. If '1', then we |
1057
|
|
|
|
|
|
|
# also need to increment $leading. |
1058
|
|
|
|
|
|
|
|
1059
|
34
|
100
|
|
|
|
111
|
return $leading . "$dp" . $to_inc |
1060
|
|
|
|
|
|
|
if $carry eq ''; |
1061
|
|
|
|
|
|
|
|
1062
|
24
|
|
|
|
|
48
|
my($sign, $len_lead) = ('', length($leading)); |
1063
|
|
|
|
|
|
|
|
1064
|
24
|
|
|
|
|
62
|
$leading =~ s/^\-//; |
1065
|
|
|
|
|
|
|
|
1066
|
24
|
100
|
|
|
|
68
|
if($len_lead != length($leading)) { |
1067
|
12
|
|
|
|
|
23
|
$sign = '-'; |
1068
|
12
|
|
|
|
|
19
|
$len_lead--; |
1069
|
|
|
|
|
|
|
} |
1070
|
|
|
|
|
|
|
|
1071
|
24
|
|
|
|
|
49
|
$carry = _increment($leading); |
1072
|
|
|
|
|
|
|
|
1073
|
24
|
|
|
|
|
81
|
return $sign . $carry . $leading . "$dp" . $to_inc; |
1074
|
|
|
|
|
|
|
} |
1075
|
|
|
|
|
|
|
|
1076
|
|
|
|
|
|
|
####################################################################### |
1077
|
|
|
|
|
|
|
####################################################################### |
1078
|
|
|
|
|
|
|
|
1079
|
|
|
|
|
|
|
sub _increment { |
1080
|
58
|
|
|
58
|
|
97
|
my $carry = 1; |
1081
|
58
|
|
|
|
|
100
|
my $len = length($_[0]) * -1; |
1082
|
|
|
|
|
|
|
|
1083
|
58
|
|
|
|
|
143
|
for(my $offset = -1; $offset >= $len; $offset--) { |
1084
|
52
|
|
|
|
|
140
|
substr($_[0], $offset, 1) = (substr($_[0], $offset, 1) + 1) % 10; |
1085
|
52
|
100
|
|
|
|
161
|
if(substr($_[0], $offset, 1) ne '0') { |
1086
|
30
|
|
|
|
|
57
|
$carry = ''; |
1087
|
30
|
|
|
|
|
54
|
last; |
1088
|
|
|
|
|
|
|
} |
1089
|
|
|
|
|
|
|
} |
1090
|
|
|
|
|
|
|
|
1091
|
58
|
|
|
|
|
120
|
return $carry; |
1092
|
|
|
|
|
|
|
} |
1093
|
|
|
|
|
|
|
|
1094
|
|
|
|
|
|
|
####################################################################### |
1095
|
|
|
|
|
|
|
####################################################################### |
1096
|
|
|
|
|
|
|
|
1097
|
|
|
|
|
|
|
*decode_d128 = $Math::Decimal128::fmt eq 'DPD' ? \&decode_dpdl : \&decode_bidl; |
1098
|
|
|
|
|
|
|
|
1099
|
|
|
|
|
|
|
####################################################################### |
1100
|
|
|
|
|
|
|
####################################################################### |
1101
|
|
|
|
|
|
|
|
1102
|
|
|
|
|
|
|
1; |
1103
|
|
|
|
|
|
|
|
1104
|
|
|
|
|
|
|
__END__ |