├── AUTHORS
├── CONTRIBUTING.md
├── CONTRIBUTORS
├── LICENSE
├── Makefile
├── README
└── dumpfp.cc
/AUTHORS:
--------------------------------------------------------------------------------
1 | # Copyright Holders.
2 |
3 | Google, Inc.
4 |
--------------------------------------------------------------------------------
/CONTRIBUTING.md:
--------------------------------------------------------------------------------
1 | ## Signing the CLA
2 |
3 | Please sign the [Google Contributor License Agreement
4 | (CLA)](https://cla.developers.google.com/)
5 | before sending pull requests. For any code changes to be
6 | accepted, the CLA must be signed. It's a quick process, I
7 | promise!
8 |
--------------------------------------------------------------------------------
/CONTRIBUTORS:
--------------------------------------------------------------------------------
1 | # People who have contributed code but are not copyright holders.
2 |
3 | Joshua Haberman
4 |
--------------------------------------------------------------------------------
/LICENSE:
--------------------------------------------------------------------------------
1 |
2 | Apache License
3 | Version 2.0, January 2004
4 | http://www.apache.org/licenses/
5 |
6 | TERMS AND CONDITIONS FOR USE, REPRODUCTION, AND DISTRIBUTION
7 |
8 | 1. Definitions.
9 |
10 | "License" shall mean the terms and conditions for use, reproduction,
11 | and distribution as defined by Sections 1 through 9 of this document.
12 |
13 | "Licensor" shall mean the copyright owner or entity authorized by
14 | the copyright owner that is granting the License.
15 |
16 | "Legal Entity" shall mean the union of the acting entity and all
17 | other entities that control, are controlled by, or are under common
18 | control with that entity. For the purposes of this definition,
19 | "control" means (i) the power, direct or indirect, to cause the
20 | direction or management of such entity, whether by contract or
21 | otherwise, or (ii) ownership of fifty percent (50%) or more of the
22 | outstanding shares, or (iii) beneficial ownership of such entity.
23 |
24 | "You" (or "Your") shall mean an individual or Legal Entity
25 | exercising permissions granted by this License.
26 |
27 | "Source" form shall mean the preferred form for making modifications,
28 | including but not limited to software source code, documentation
29 | source, and configuration files.
30 |
31 | "Object" form shall mean any form resulting from mechanical
32 | transformation or translation of a Source form, including but
33 | not limited to compiled object code, generated documentation,
34 | and conversions to other media types.
35 |
36 | "Work" shall mean the work of authorship, whether in Source or
37 | Object form, made available under the License, as indicated by a
38 | copyright notice that is included in or attached to the work
39 | (an example is provided in the Appendix below).
40 |
41 | "Derivative Works" shall mean any work, whether in Source or Object
42 | form, that is based on (or derived from) the Work and for which the
43 | editorial revisions, annotations, elaborations, or other modifications
44 | represent, as a whole, an original work of authorship. For the purposes
45 | of this License, Derivative Works shall not include works that remain
46 | separable from, or merely link (or bind by name) to the interfaces of,
47 | the Work and Derivative Works thereof.
48 |
49 | "Contribution" shall mean any work of authorship, including
50 | the original version of the Work and any modifications or additions
51 | to that Work or Derivative Works thereof, that is intentionally
52 | submitted to Licensor for inclusion in the Work by the copyright owner
53 | or by an individual or Legal Entity authorized to submit on behalf of
54 | the copyright owner. For the purposes of this definition, "submitted"
55 | means any form of electronic, verbal, or written communication sent
56 | to the Licensor or its representatives, including but not limited to
57 | communication on electronic mailing lists, source code control systems,
58 | and issue tracking systems that are managed by, or on behalf of, the
59 | Licensor for the purpose of discussing and improving the Work, but
60 | excluding communication that is conspicuously marked or otherwise
61 | designated in writing by the copyright owner as "Not a Contribution."
62 |
63 | "Contributor" shall mean Licensor and any individual or Legal Entity
64 | on behalf of whom a Contribution has been received by Licensor and
65 | subsequently incorporated within the Work.
66 |
67 | 2. Grant of Copyright License. Subject to the terms and conditions of
68 | this License, each Contributor hereby grants to You a perpetual,
69 | worldwide, non-exclusive, no-charge, royalty-free, irrevocable
70 | copyright license to reproduce, prepare Derivative Works of,
71 | publicly display, publicly perform, sublicense, and distribute the
72 | Work and such Derivative Works in Source or Object form.
73 |
74 | 3. Grant of Patent License. Subject to the terms and conditions of
75 | this License, each Contributor hereby grants to You a perpetual,
76 | worldwide, non-exclusive, no-charge, royalty-free, irrevocable
77 | (except as stated in this section) patent license to make, have made,
78 | use, offer to sell, sell, import, and otherwise transfer the Work,
79 | where such license applies only to those patent claims licensable
80 | by such Contributor that are necessarily infringed by their
81 | Contribution(s) alone or by combination of their Contribution(s)
82 | with the Work to which such Contribution(s) was submitted. If You
83 | institute patent litigation against any entity (including a
84 | cross-claim or counterclaim in a lawsuit) alleging that the Work
85 | or a Contribution incorporated within the Work constitutes direct
86 | or contributory patent infringement, then any patent licenses
87 | granted to You under this License for that Work shall terminate
88 | as of the date such litigation is filed.
89 |
90 | 4. Redistribution. You may reproduce and distribute copies of the
91 | Work or Derivative Works thereof in any medium, with or without
92 | modifications, and in Source or Object form, provided that You
93 | meet the following conditions:
94 |
95 | (a) You must give any other recipients of the Work or
96 | Derivative Works a copy of this License; and
97 |
98 | (b) You must cause any modified files to carry prominent notices
99 | stating that You changed the files; and
100 |
101 | (c) You must retain, in the Source form of any Derivative Works
102 | that You distribute, all copyright, patent, trademark, and
103 | attribution notices from the Source form of the Work,
104 | excluding those notices that do not pertain to any part of
105 | the Derivative Works; and
106 |
107 | (d) If the Work includes a "NOTICE" text file as part of its
108 | distribution, then any Derivative Works that You distribute must
109 | include a readable copy of the attribution notices contained
110 | within such NOTICE file, excluding those notices that do not
111 | pertain to any part of the Derivative Works, in at least one
112 | of the following places: within a NOTICE text file distributed
113 | as part of the Derivative Works; within the Source form or
114 | documentation, if provided along with the Derivative Works; or,
115 | within a display generated by the Derivative Works, if and
116 | wherever such third-party notices normally appear. The contents
117 | of the NOTICE file are for informational purposes only and
118 | do not modify the License. You may add Your own attribution
119 | notices within Derivative Works that You distribute, alongside
120 | or as an addendum to the NOTICE text from the Work, provided
121 | that such additional attribution notices cannot be construed
122 | as modifying the License.
123 |
124 | You may add Your own copyright statement to Your modifications and
125 | may provide additional or different license terms and conditions
126 | for use, reproduction, or distribution of Your modifications, or
127 | for any such Derivative Works as a whole, provided Your use,
128 | reproduction, and distribution of the Work otherwise complies with
129 | the conditions stated in this License.
130 |
131 | 5. Submission of Contributions. Unless You explicitly state otherwise,
132 | any Contribution intentionally submitted for inclusion in the Work
133 | by You to the Licensor shall be under the terms and conditions of
134 | this License, without any additional terms or conditions.
135 | Notwithstanding the above, nothing herein shall supersede or modify
136 | the terms of any separate license agreement you may have executed
137 | with Licensor regarding such Contributions.
138 |
139 | 6. Trademarks. This License does not grant permission to use the trade
140 | names, trademarks, service marks, or product names of the Licensor,
141 | except as required for reasonable and customary use in describing the
142 | origin of the Work and reproducing the content of the NOTICE file.
143 |
144 | 7. Disclaimer of Warranty. Unless required by applicable law or
145 | agreed to in writing, Licensor provides the Work (and each
146 | Contributor provides its Contributions) on an "AS IS" BASIS,
147 | WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or
148 | implied, including, without limitation, any warranties or conditions
149 | of TITLE, NON-INFRINGEMENT, MERCHANTABILITY, or FITNESS FOR A
150 | PARTICULAR PURPOSE. You are solely responsible for determining the
151 | appropriateness of using or redistributing the Work and assume any
152 | risks associated with Your exercise of permissions under this License.
153 |
154 | 8. Limitation of Liability. In no event and under no legal theory,
155 | whether in tort (including negligence), contract, or otherwise,
156 | unless required by applicable law (such as deliberate and grossly
157 | negligent acts) or agreed to in writing, shall any Contributor be
158 | liable to You for damages, including any direct, indirect, special,
159 | incidental, or consequential damages of any character arising as a
160 | result of this License or out of the use or inability to use the
161 | Work (including but not limited to damages for loss of goodwill,
162 | work stoppage, computer failure or malfunction, or any and all
163 | other commercial damages or losses), even if such Contributor
164 | has been advised of the possibility of such damages.
165 |
166 | 9. Accepting Warranty or Additional Liability. While redistributing
167 | the Work or Derivative Works thereof, You may choose to offer,
168 | and charge a fee for, acceptance of support, warranty, indemnity,
169 | or other liability obligations and/or rights consistent with this
170 | License. However, in accepting such obligations, You may act only
171 | on Your own behalf and on Your sole responsibility, not on behalf
172 | of any other Contributor, and only if You agree to indemnify,
173 | defend, and hold each Contributor harmless for any liability
174 | incurred by, or claims asserted against, such Contributor by reason
175 | of your accepting any such warranty or additional liability.
176 |
177 | END OF TERMS AND CONDITIONS
178 |
179 | APPENDIX: How to apply the Apache License to your work.
180 |
181 | To apply the Apache License to your work, attach the following
182 | boilerplate notice, with the fields enclosed by brackets "[]"
183 | replaced with your own identifying information. (Don't include
184 | the brackets!) The text should be enclosed in the appropriate
185 | comment syntax for the file format. We also recommend that a
186 | file or class name and description of purpose be included on the
187 | same "printed page" as the copyright notice for easier
188 | identification within third-party archives.
189 |
190 | Copyright [yyyy] [name of copyright owner]
191 |
192 | Licensed under the Apache License, Version 2.0 (the "License");
193 | you may not use this file except in compliance with the License.
194 | You may obtain a copy of the License at
195 |
196 | http://www.apache.org/licenses/LICENSE-2.0
197 |
198 | Unless required by applicable law or agreed to in writing, software
199 | distributed under the License is distributed on an "AS IS" BASIS,
200 | WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
201 | See the License for the specific language governing permissions and
202 | limitations under the License.
203 |
--------------------------------------------------------------------------------
/Makefile:
--------------------------------------------------------------------------------
1 |
2 | dumpfp: dumpfp.cc
3 | $(CXX) -o dumpfp dumpfp.cc -Wall -O3 -lm -lgmpxx -lgmp
4 |
--------------------------------------------------------------------------------
/README:
--------------------------------------------------------------------------------
1 |
2 | Tool for dumping floating-point numbers in their native
3 | representation.
4 |
5 | For more info see: http://blog.reverberate.org/2012/11/dumpfp-tool-to-inspect-floating-point.html
6 |
--------------------------------------------------------------------------------
/dumpfp.cc:
--------------------------------------------------------------------------------
1 |
2 | #include
3 | #include
4 | #include
5 | #include
6 | #include
7 | #include
8 | #include
9 | #include
10 | #include
11 |
12 | using namespace std;
13 |
14 | mpz_class bitat(int bit) { return mpz_class(1) << bit; }
15 | mpz_class mask(int bits) { return bitat(bits) - 1; }
16 |
17 | void dumpfp(mpz_class whole, mpz_class numer, int denom_lg2) {
18 | mpz_class denom = bitat(denom_lg2);
19 |
20 | // Reduce the fraction for nicer display.
21 | while ((numer & 0x1) == 0 && (denom & 0x1) == 0) {
22 | numer >>= 1;
23 | denom >>= 1;
24 | denom_lg2--;
25 | }
26 |
27 | // Pick the fewest number of digits that will precisely represent the value.
28 | int n = 1;
29 | mpz_class digits = 10;
30 | again:
31 | mpz_class product = numer * digits;
32 | mpz_class fraction = product / denom;
33 | mpz_class product2 = fraction * denom;
34 | if (product2 != product) {
35 | digits *= 10;
36 | n += 1;
37 | goto again;
38 | }
39 |
40 | if (whole > 0)
41 | cout << whole;
42 | if (whole > 0 && numer > 0)
43 | cout << " + ";
44 | if (numer > 0)
45 | cout << numer << "/2^" << denom_lg2;
46 |
47 | gmp_printf(" (%Zd.%.*Zd)", whole.get_mpz_t(), n, fraction.get_mpz_t());
48 | }
49 |
50 | void dumpfp2(mpz_class raw, int sig_bits, int exp_bits) {
51 | // In these calculations we rely on the general properties of IEEE floating
52 | // point: X bits of significant, followed by an implicit "1", followed by
53 | // Y bits of exponent and a sign bit. The exponent is expected to be biased
54 | // by half of the exponent's domain.
55 | //
56 | // Note that this will *not* work for "long double" because it does not use
57 | // an implicit 1 in its significand.
58 | mpz_class sig_scale = bitat(sig_bits),
59 | sig = raw & mask(sig_bits),
60 | sig2 = sig_scale + sig,
61 | expbias = mask(exp_bits - 1);
62 | int sign = mpz_class(raw >> (sig_bits + exp_bits)).get_si(),
63 | exp = mpz_class(raw >> sig_bits & mask(exp_bits)).get_si(),
64 | exp2 = mpz_class(exp - expbias).get_si();
65 |
66 | cout << " raw = 0x" << hex << raw << "\n";
67 | cout << " sign = 0x" << hex << sign << "\n";
68 |
69 | if (exp2 == (expbias + 1)) {
70 | cout << " exponent = 0x" << hex << exp << " (NaN or Infinity)\n";
71 | if (sig != 0)
72 | cout << " significand = " << hex << sig << " (non-zero indicates NaN)\n";
73 | else
74 | cout << " significand = " << hex << sig << " (zero indicates Infinity)\n";
75 | } else {
76 | cout << " exponent = 0x" << hex << exp << " (" << dec << exp2 << ")\n";
77 | cout << " significand = 0x" << hex << sig << dec << "\n";
78 | cout << "\n";
79 |
80 | cout << " VALUE CALCULATION =\n";
81 | cout << " significand (";
82 | dumpfp(1, sig, sig_bits);
83 | cout << ")\n";
84 | cout << " * 2^exponent (2^" << exp2 << ")\n";
85 | cout << " = VALUE (";
86 |
87 | if (exp2 >= sig_bits) {
88 | cout << (sig2 << (exp2 - sig_bits));
89 | } else {
90 | int split = sig_bits - exp2;
91 | mpz_class whole = sig2 >> split;
92 | mpz_class frac = sig2 - (whole << split);
93 | dumpfp(whole, frac, split);
94 | }
95 | cout << ")\n";
96 | }
97 |
98 | cout << "\n";
99 | }
100 |
101 | // Gets an integer value corresponding to the given raw bytes.
102 | // The currently will only work on a little-endian machine.
103 | mpz_class getraw(void *val, size_t size) {
104 | unsigned char bytes[128];
105 | memcpy(bytes, val, size);
106 | mpz_class raw(0);
107 | for (size_t i = 0; i < size; i++) {
108 | raw *= 256;
109 | raw += bytes[size-i-1];
110 | }
111 | return raw;
112 | }
113 |
114 | void dumpfloat(float x) {
115 | printf("Single Precision (IEEE 32-bit):\n");
116 | dumpfp2(getraw(&x, sizeof(float)), 23, 8);
117 | }
118 |
119 | void dumpdouble(double x) {
120 | printf("Double Precision (IEEE 64-bit):\n");
121 | dumpfp2(getraw(&x, sizeof(double)), 52, 11);
122 | }
123 |
124 | int main(int argc, char *argv[]) {
125 | if (argc < 2) {
126 | fprintf(stderr, "Usage: dumpfp 0.12345\n");
127 | return 1;
128 | }
129 |
130 | long double x = strtold(argv[1], NULL);
131 |
132 | dumpfloat(x);
133 | dumpdouble(x);
134 |
135 | return 0;
136 | }
137 |
--------------------------------------------------------------------------------