Decimal64 floating-point format: Difference between revisions

Revision as of 20:46, 1 January 2025

Latest revision as of 14:33, 25 August 2025

(9 intermediate revisions by 5 users not shown)

Line 1:

{{Short description|64-bit computer number format}}

{{lowercase title}}

{{Use dmy dates|date=July 2020|cs1-dates=y}}

{{floating-point}}

In [[computing]], '''decimal64''' is a [[decimal floating point|decimal floating-point]] [[computer number format]] that occupies 8 bytes (64 bits) in computer memory.

Decimal64 is a decimal floating-point format, formally introduced in the [[IEEE 754-2008 revision|2008 revision]]<ref name="IEEE-754_2008">{{cite book |title=IEEE Standard for Floating-Point Arithmetic |author=IEEE Computer Society |date=2008-08-29 |publisher=[[IEEE]] |id=IEEE Std 754-2008 |doi=10.1109/IEEESTD.2008.4610935 |ref=CITEREFIEEE_7542008 |isbn=978-0-7381-5753-5 }}</ref> of the [[IEEE 754]] standard, also known as ISO/IEC/IEEE 60559:2011.<ref name="ISO-60559_2011">{{Cite book |last=ISO/IEC JTC 1/SC 25|title=ISO/IEC/IEEE 60559:2011 — Information technology — Microprocessor Systems — Floating-Point arithmetic |url=https://www.iso.org/standard/57469.html |publisher=ISO |pages=1–58 |date=June 2011}}</ref>

== Purpose and use ==

decimal64 fits well to replace binary64 format in applications where 'small deviations' are unwanted and speed isn't extremely crucial.

== Format ==

In contrast to the '''binary'''xxx data formats the '''decimal'''xxx formats provide exact representation of decimal fractions, exact calculations with them and enable human common 'ties away from zero' rounding (in some range, to some precision, to some degree). In a trade-off for reduced performance. They are intended for applications where it's requested to come near to schoolhouse math, such as financial and tax computations. (In short they avoid plenty of problems like 0.2 + 0.1 -> 0.30000000000000004 which happen with binary64 datatypes.)

Decimal64 supports 'normal' values that can have 16 digit precision from {{gaps|±1.000|000|000|000|000|e=-383}} to {{gaps|±9.999|999|999|999|999|e=384}}, plus 'denormal' values with ramp-down relative precision down to ±1.×10⁻³⁹⁸, [[signed zero]]s, signed infinities and [[NaN]] (Not a Number). This format supports two different encodings.

The binary format of the same size supports a range from denormal-min {{gaps|±5|||||e=-324|}}, over normal-min with full 53-bit precision {{gaps|±2.225|073|858|507|201|e=-308|4}} to max {{gaps|±1.797|693|134|862|315|e=+308|7}}.

== Range and precision ==

Decimal64 supports 'normal' values that can have 16 digit precision from {{gaps|±1.000|000|000|000|000|e=-383}} to {{gaps|±9.999|999|999|999|999|e=384}}, plus 'denormal' values with ramp-down relative precision down to ±1 × 10⁻³⁹⁸ (only one digit left), [[signed zero]]s, signed infinities and [[NaN]] (Not a Number).

Because the significand for the [[IEEE 754]] decimal formats is not normalized, most values with less than 16 [[significant digits]] have multiple possible representations; 1000000 × 10⁻²=100000 × 10⁻¹=10000 × 10⁰=1000 × 10¹ all have the value 10000. These sets of representations for a same value are called ''[[Cohort (floating point)|cohorts]]'', the different members can be used to denote how many digits of the value are known precisely. Each signed zero has 768 possible representations (1536 for all zeros, in two different cohorts).

The binary format of the same bit-size supports a range from denormal-min {{gaps|±5|||||e=-324|}}, over normal-min with full 53-bit precision {{gaps|±2.225|073|858|507|201|e=-308|4}} to max {{gaps|±1.797|693|134|862|315|e=+308|7}}.

== Encoding of decimal64 values ==

== Performance ==

Performance comparison is not easy, not very accurate and lacks reproducibility on modern IT systems for various reasons. One can roughly say that in a current 64-bit Intel(r) / linux / gcc / libdfp / BID implementation, basic arithmetic operations with decimal64 values are between factor 2 and 15 slower than with binary64 data types, while 'higher' functions like powers ( ~600 ) and trigonometric functions like tangent ( {{gaps|~10|000}} ) suffer more performance penalties. To get an idea about performance on a specific system the code in 'Addendum - code' can be used. Perhaps the [https://gcc.gnu.org/contribute.html GNU gcc project] and the [https://github.com/libdfp/libdfp/issues?q=is%3Aissue 'libdfp' project on github] could like some help to improve.

== Representation / encoding of decimal64 values ==

decimal64 values are represented in a 'not normalized' near to 'scientific format', with combining some bits of the exponent with the leading bits of the significand in a 'combination field'.

{| class="wikitable"

|+ Generic encoding

|-

! Sign !! Combination !! Significand continuation

! <u>S</u>ign !! Co<u>m</u>bination !! <u>T</u>railing significand bits

|-

! 1 bit !! 13 bits !! 50 bits

|-

| {{mono|s}} || {{mono|mmmmmmmmmmmmm}} || tttttttttttttttttttttttttttttttttttttttttttttttttt{{mono|cccccccccccccccccccccccccccccccccccccccccccccccccc}}

|}

IEEE 754 allows two alternative encodings for decimal64 values. The standard does not specify how to signify which representation is used, for instance in a situation where decimal64 values are communicated between systems:

Besides the special cases infinities and NaNs there are four points relevant to understand the encoding of decimal64.

* In the [[#Binary integer significand field|binary encoding]], the 16-digit significand is represented as a binary coded positive integer, based on [[binary integer decimal]] (BID).

* BID vs. DPD encoding, '''B'''inary '''I'''nteger '''D'''ecimal using a positive [[binary integer decimal]] for the significand, software centric and designed by Intel(r), vs. '''D'''ensely '''P'''acked '''D'''ecimal based on [[densely packed decimal]] encoding for all except the first digit of the significand, hardware centric and promoted by IBM(r), differences see below. Both alternatives provide exactly the same range of representable numbers: 16 digits of significand and {{math|size=100%|1=3 × 2⁸ = 768}} possible exponent values. IEEE 754 allows these two different encodings, without a concept to denote which is used, for instance in a situation where decimal64 values are communicated between systems. CAUTION!: Be aware that transferring binary data between systems using different encodings will mostly produce valid decimal64 numbers, '''but with different value'''. Prefer data exchange in íntegral or ASCII 'triplets' for sign, exponent and significand.

* In the [[#Densely packed decimal significand field|decimal encoding]], the 16-digit significand is represented as a decimal coded positive integer, based on [[densely packed decimal]] (DPD) with 5 groups of 3 digits (except the most significant digit encoded specially) are each represented in declets (10-bit sequences). This is pretty efficient, because 2¹⁰ = 1024, is only little more than needed to still contain all numbers from 0 to 999.

Both alternatives provide exactly the same set of representable numbers: 16 digits of significand and {{math|size=100%|1=3 × 2⁸ = 768}} possible decimal exponent values. (All the possible decimal exponent values storable in a [[binary64]] number are representable in decimal64, and most bits of the significand of a binary64 are stored keeping roughly the same number of decimal digits in the significand.)

* The significands are not 'normalized' (the leading digit(s) are allowed to be "0"), and thus most values with less than 7 [[significant digits]] have multiple possible representations; 1000000 × 10^-2=100000 × 10^-1=10000 × 10⁰=1000 × 10¹ all have the value 10000. These sets of representations for a same value are called [[Cohort (floating point)|cohorts]]'','' the different members can be used to denote how many digits of the value are known precisely.

In both cases, the most significant 4 bits of the significand (which actually only have 10 possible values) are combined with two bits of the exponent (3 possible values) to use 30 of the 32 possible values of a 5-bit field. The remaining combinations encode [[infinity|infinities]] and [[NaN]]s. BID and DPD use different bits of the combination field for that.

* The encodings combine two bits of the exponent with the leading 3 to 4 bits of the significand in a 'combination field', different for 'big' vs. 'small' significands. That enables bigger precision and range, in trade-off that some simple functions like sort and compare, very frequently used in coding, do not work on the bit pattern but require computations to extract exponent and significand and then try to obtain an exponent aligned representation. This effort is partly balanced by saving the effort for normalization, but contributes to the slower performance of the decimal datatypes. Beware: BID and DPD use different bits of the combination field for that, see below.

In the cases of Infinity and NaN, all other bits of the encoding are ignored. Thus, it is possible to initialize an array to Infinities or NaNs by filling it with a single byte value.

* Different understanding of significand as integer or fraction, and acc. different bias to apply for the exponent (for decimal64 what is stored in bits can be decoded as base to the power of 'stored value for the exponent minus '''bias of 383'''<nowiki/>' times significand understood as d₀ '''.''' d₋₁ d₋₂ d₋₃ d₋₄ d₋₅ d₋₆ d₋₇ d₋₈ d₋₉ d₋₁₀ d₋₁₁ d₋₁₂ d₋₁₃ d₋₁₄ d₋₁₅ (note: radix dot after first digit, significand '''fractional'''), or base to the power of 'stored value for the exponent minus '''bias of 398'''<nowiki/>' times significand understood as d15 d14 d13 d12 d11 d10 d9 d8 d7 d₆ d₅ d₄ d₃ d₂ d₁ d₀ (note: no radix dot, significand '''integral'''), both produce the same result [2019 version of IEEE 754 in clause 3.3, page 18]. For decimal datatypes the second view is more common, while for binary datatypes the first, the biases are different for each datatype.)

=== Binary integer significand field ===

In all cases for decimal364, the value represented is

: (−1)^''sign'' × 10^{''exponent''−'''398'''} × ''significand'', with the ''significand'' understood as positive integer.

Alternatively it can be understood as (−1)^''sign'' × 10^{''exponent''−'''383'''} × ''significand'' with the ''significand'' digits understood as d₀ '''.''' d₋₁ d₋₂ d₋₃ d₋₄ d₋₅ d_{−6 ...}, note the radix dot making it a fraction.

decimal64 is superpower to binary64 in range, and head to head in precision, dec64 all normal values 16-digit, while bin64 alternates between 'not fully 16' and about 16.5 decimal digits ( 17 digits with gaps ).

=== BID encoding ===

This format uses a binary significand from 0 to {{math|size=100%|1=10¹⁶ − 1 = {{gaps|9|999|999|999|999|999}} = 2386F26FC0FFFF₁₆ = {{gaps|1000|1110000110|1111001001|1011111100|0000111111|1111111111₂}}.}}The encoding, completely stored on 64 bits, can represent binary significands up to {{math|size=100%|1=10 × 2⁵⁰ − 1 = {{gaps|11|258|999|068|426|239}} = 27FFFFFFFFFFFF₁₆,}} but values larger than {{math|size=100%|1=10¹⁶ − 1}} are illegal (and the standard requires implementations to treat them as 0, if encountered on input).

Line 56 ⟶ 44:

If the {{val|2|u=bits}} after the sign bit are "11", then the 10-bit exponent field is shifted {{val|2|u=bits}} to the right (after both the sign bit and the "11" bits thereafter), and the represented significand is in the remaining {{val|51|u=bits}}. In this case there is an implicit (that is, not stored) leading 3-bit sequence "100" for the MSB bits of the true significand (in the remaining lower bits ''ttt...ttt'' of the significand, not all possible values are used).

Be aware that the bit numbering used in the tables for e.g. m₁₂ … m₀  is in opposite direction than that used in the paper for the IEEE 754 standard G₀ … G₁₂.

{| class="wikitable" style="text-align:left; border-width:0;"

Line 67 ⟶ 53:

! rowspan="2" |Significand / Description

|-

! g12 !! g11 !! g10 !! g9 !! g8 !! g7 !! g6 !! g5 !! g4 !! g3 !! g2

! m₁₂!! m₁₁!! m₁₀!! m₉!! m₈!! m₇!! m₆!! m₅!! m₄!! m₃!! m₂

!g1

!m₁

!g0

!m₀

|-

| colspan="16" |combination field not! starting with '11', bits ab = 00, 01 or 10

|-

| style="font-family:monospace; background:#cedff2;" | '''a''' || style="font-family:monospace; background:#cedff2;" | '''b''' || style="font-family:monospace; background:#cedff2;" | '''c''' || style="font-family:monospace; background:#cedff2;" | '''d''' || style="font-family:monospace; background:#cedff2;" | '''m''' || style="font-family:monospace; background:#cedff2;" | '''m''' || style="font-family:monospace; background:#cedff2;" | '''m''' || style="font-family:monospace; background:#cedff2;" | '''m''' || style="font-family:monospace; background:#cedff2;" | '''m''' || style="font-family:monospace; background:#cedff2;" | '''m''' || style="font-family:monospace; background:#cef2e0;" | '''e''' || style="font-family:monospace; background:#cef2e0;" |'''f''' || style="font-family:monospace; background:#cef2e0;" |'''g'''

| || style="font-family:monospace; background:#cedff2;" | '''abcdmmmmmm''' || style="background:#cef2e0;" | {{mono|(0)'''efgtttttttttttttttttttttttttttttttttttttttttttttttttt''' }}

Finite number with 'small' significand,first beingdigit <of 9007199254740992, fits into 53significand bits(0&nbsp;..&nbsp;7).

|-

| colspan="16" |combination field starting with '11', but not 1111, bits ab = 11, bits cd = 00, 01 or 10

|-

| 1 || 1 || style="font-family:monospace; background:#cedff2;" | '''c'''|| style="font-family:monospace; background:#cedff2;" | '''d''' || style="font-family:monospace; background:#cedff2;" | '''m''' || style="font-family:monospace; background:#cedff2;" | '''m''' || style="font-family:monospace; background:#cedff2;" | '''m''' || style="font-family:monospace; background:#cedff2;" | '''m''' || style="font-family:monospace; background:#cedff2;" | '''m''' || style="font-family:monospace; background:#cedff2;" | '''m''' || style="font-family:monospace; background:#cedff2;" | '''e''' || style="font-family:monospace; background:#cedff2;" | '''f''' || style="font-family:monospace; background:#cef2e0;" | '''g'''

| || style="font-family:monospace; background:#cedff2;" | '''cdmmmmmmef''' || style="background:#cef2e0;" | {{mono|'''100gtttttttttttttttttttttttttttttttttttttttttttttttttt''' }}

Finite number with 'big' significand,first beingdigit >of 9007199254740991,significand needs(8 54or bits9).

|-

| colspan="16" |combination field starting with '1111', bits abcd = 1111

Line 103 ⟶ 89:

|signaling NaN (with payload in significand)

|}

In contrast to DPD format below theThe leading bits of the significand field do ''not'' encode the most significant decimal digit; they are simply part of a larger pure-binary number. For example, combineda significand of {{gaps|8|000|000|000|000|000}} is encoded as binary {{gaps|0111|0001101011|1111010100|1001100011|0100000000|0000000000}}₂, with the implicitleading prefix{{val|4|u=bits}} ofencoding 7; the first significand which requires a 54th bit is {{math|size=100%|1=2⁵³ for= big{{gaps|9|007|199|254|740|992}}.}} significands,The simplyhighest partvalid ofsignificant ais larger{{gaps|9|999|999|999|999|999}} whose pure-binary number.encoding is

{{gaps|(100)0|1110000110|1111001001|1011111100|0000111111|1111111111}}₂ (with the 3 most significant bits (100) not stored but implicit as shown above; and the next bit is always zero in valid encodings).

The resulting 'raw' exponent is a 10 bit binary integer where the leading bits are not '11', thus values 0 ... 1011111111_b = 0 ... 767_d, appr. bias is to be subtracted. The resulting significand could be a positive binary integer of 54 bits up to 1001 1111111111 1111111111 1111111111 1111111111 1111111111_b = 11258999068426239_d, but values above 10¹⁶ − 1 = 9999999999999999 = 2386F26FC0FFFF₁₆ = 100011100001101111001001101111110000001111111111111111₂ are 'illegal' and have to be treated as zeroes. To obtain the individual decimal digits the significand has to be divided by 10 repeatedly.

In the above cases, the value represented is

Line 111 ⟶ 96:

: {{math|1=(−1)^sign × 10^{exponent−398} × significand}} <!-- Remember, significand is defined as an integer: 0 <= significand < 10^16 -->

If the four bits after the sign bit are "1111" then the value is an infinity or a NaN, as described above:

=== DPD encoding ===

0 11110 xx...x +infinity

1 11110 xx...x -infinity

x 11111 0x...x a quiet NaN

x 11111 1x...x a signalling NaN

=== Densely packed decimal significand field ===

In this version, the significand is stored as a series of decimal digits. The leading digit is between 0 and 9 (3 or 4 binary bits), and the rest of the significand uses the [[densely packed decimal]] (DPD) encoding.

Line 122 ⟶ 114:

If the first two bits after the sign bit are "00", "01", or "10", then those are the leading bits of the exponent, and the three bits "cde" after that are interpreted as the leading decimal digit (0 to 7):

If the first two bits after the sign bit are "11", then the second 2-bits are the leading bits of the exponent, and the next bit "e" is prefixed with implicit bits "100" to form the leading decimal digit of the significand (8 or 9):

The remaining two combinations (11 110 and 11 111) of the 5-bit field after the sign bit are used to represent ±infinity and NaNs, respectively.

Line 134 ⟶ 126:

! rowspan="2" |Significand / Description

|-

! g12 !! g11 !! g10 !! g9 !! g8 !! g7 !! g6 !! g5 !! g4 !! g3 !! g2

! m₁₂!! m₁₁!! m₁₀!! m₉!! m₈!! m₇!! m₆!! m₅!! m₄!! m₃!! m₂

!g1

!m₁

!g0

!m₀

|-

| colspan="16" |combination field not! starting with '11', bits ab = 00, 01 or 10

|-

| style="font-family:monospace; background:#cedff2;" | '''a''' || style="font-family:monospace; background:#cedff2;" | '''b''' || style="font-family:monospace; background:#cef2e0;" | '''c''' || style="font-family:monospace; background:#cef2e0;" | '''d''' || style="font-family:monospace; background:#cef2e0;" | '''e''' || style="font-family:monospace; background:#cedff2;" | '''m''' || style="font-family:monospace; background:#cedff2;" | '''m''' || style="font-family:monospace; background:#cedff2;" | '''m''' || style="font-family:monospace; background:#cedff2;" | '''m''' || style="font-family:monospace; background:#cedff2;" | '''m''' || style="font-family:monospace; background:#cedff2;" | '''m''' || style="font-family:monospace; background:#cedff2;" | '''m''' || style="font-family:monospace; background:#cedff2;" | '''m'''

| || style="font-family:monospace; background:#cedff2;" | '''abmmmmmmmm'''|| style="background:#cef2e0;" | {{nowrap|{{mono|(0)'''cde tttttttttt tttttttttt tttttttttt tttttttttt tttttttttt''' }}}}

Finite number with small first digit of significand (0&nbsp;…&nbsp;7).

|-

| colspan="16" |combination field starting with '11', but not 1111, bits ab = 11, bits cd = 00, 01 or 10

|-

| 1 || 1 || style="font-family:monospace; background:#cedff2;" | '''c''' || style="font-family:monospace; background:#cedff2;" | '''d''' || style="font-family:monospace; background:#cef2e0;" | '''e''' || style="font-family:monospace; background:#cedff2;" | '''m''' || style="font-family:monospace; background:#cedff2;" | '''m''' || style="font-family:monospace; background:#cedff2;" | '''m''' || style="font-family:monospace; background:#cedff2;" | '''m''' || style="font-family:monospace; background:#cedff2;" | '''m''' || style="font-family:monospace; background:#cedff2;" | '''m''' || style="font-family:monospace; background:#cedff2;" | '''m''' || style="font-family:monospace; background:#cedff2;" | '''m'''

| || style="font-family:monospace; background:#cedff2;" | '''cdmmmmmmmm'''|| style="background:#cef2e0;" | {{nowrap|{{mono|'''100e tttttttttt tttttttttt tttttttttt tttttttttt tttttttttt''' }}}}

Finite number with big first digit of significand (8 or 9).

|-

Line 171 ⟶ 163:

|}

The DPD/3BCD transcoding for the declets is given by the following table. b9...b0 are the bits of the DPD, and d2...d0 are the three BCD digits.

The resulting 'raw' exponent is a 10 bit binary integer where the leading bits are not '11', thus values 0 ... 1011111111_b = 0 ... 767_d, appr. bias is to be subtracted. The significand's leading decimal digit forms from the '''(0)cde''' or '''100e''' bits as binary integer. The subsequent digits are encoded in the 10 bit 'declet' fields 'tttttttttt' according the DPD rules (see below). The full decimal significand is then obtained by concatenating the leading and trailing decimal digits.

The 10-bit DPD to 3-digit BCD transcoding for the declets is given by the following table. b₉ … b₀ are the bits of the DPD, and d₂ … d₀ are the three BCD digits. Be aware that the bit numbering used here for e.g. b₉ … b₀ is in opposite direction than that used in the paper for the IEEE 754 standard b₀ … b₉, add. the decimal digits are numbered 0-based here while in opposite direction and 1-based in the IEEE 754 paper. The bits on white background are not counting for the value, but signal how to understand / shift the other bits. The concept is to denote which digits are small (0 … 7) and encoded in three bits, and which are not, then calculated from a prefix of '100', and one bit specifying if 8 or 9.

{{Densely packed decimal}}

Line 179 ⟶ 169:

The 8 decimal values whose digits are all 8s or 9s have four codings each.

The bits marked x in the table above are [[don't care|ignored]] on input, but will always be 0 in computed results.

(The {{math|size=100%|1=8 × 3 = 24}} non-standard encodings fill in the unusedgap range frombetween {{math|size=100%|1=10³ = 1000 toand 2¹⁰ - 1 = 10231024.}})

In the above cases, with the ''true significand'' as the sequence of decimal digits decoded, the value represented is

:<math>(-1)^\text{signbit}\times 10^{\text{exponentbits}_2-398_{10}}\times \text{truesignificand}_{10}</math>

== History ==

decimal64 was formally introduced in the [[IEEE 754-2008 revision|2008 revision]]<ref name="IEEE-754_2008">{{cite book |author=IEEE Computer Society |url=https://ieeexplore.ieee.org/document/4610935 |title=IEEE Standard for Floating-Point Arithmetic |date=2008-08-29 |publisher=[[IEEE]] |isbn=978-0-7381-5753-5 |doi=10.1109/IEEESTD.2008.4610935 |id=IEEE Std 754-2008 |ref=CITEREFIEEE_7542008 |access-date=2016-02-08}}</ref> of the [[IEEE 754]] standard, which was taken over into the ISO/IEC/IEEE 60559:2011<ref name="ISO-60559_2011">{{Cite book |last=ISO/IEC JTC 1/SC 25 |url=https://www.iso.org/standard/57469.html |title=ISO/IEC/IEEE 60559:2011 — Information technology — Microprocessor Systems — Floating-Point arithmetic |date=June 2011 |publisher=ISO |pages=1–58}}</ref> standard.

== Less important information, side effects of the encoding ==

DPD encoding is pretty efficient, not wasting more than about 2.4 percent of space vs. BID, because the 2¹⁰ = 1024 possible values in 10 bit is only little more than what is used to encode all numbers from 0 to 999.

Zero has 768 possible representations (1536 accounting signed zeroes, in two different cohorts), (even many more if you account the 'illegal' significands which have to be treated as zeroes).

The gain in range and precision by the 'combination encoding' evolves because the taken 2 bits from the exponent only use three of four possible states, and the 4 MSBs of the significand stay within 0000&nbsp;…&nbsp;1001 (10 of 16 possible states). In total that is {{math|1=3&nbsp;×&nbsp;10&nbsp;=&nbsp;30}} possible values when combined in one encoding, which is representable in 5 instead of 6 bits ({{tmath|1=2^5=32}}).

The decimalxxx formats include denormal values, for a graceful degradation of precision near zero, but in contrast to the binaryxxx formats they are not marked / do not need a special exponent, in decimal64 they are just values too small to have full 16 digit precision even with the smallest exponent.

In the cases of Infinity and NaN, all other bits of the encoding are ignored. Thus, it is possible to initialize an array to Infinities or NaNs by filling it with a single byte value.

== Addendum - code ==

try this to see performance of addition and tangent for bin32, bin64, dec32 and dec64 datatypes, compile and run instructions see in header.

The output consists of: clock cycles taken; iterations; result; expression tested .

<pre>

// program to compare the performance of binary vs. decimal datatypes,

// WIP, covering elementary operations,

// requires 'libdfp' installed,

// compile with: 'gcc -I /usr/local/include/dfp -o decxx_perf_sample.c -ldfp -lm -lquadmath'

// or - optimized - with: 'gcc -O2 -I /usr/local/include/dfp -o decxx_perf_sample.c -ldfp -lm -lquadmath'

// run with: './dec32_perf_sample value_1 value_2 (count)

// e.g. './dec32_perf_sample 8.0 5.0E-7 1000'

#define __STDC_WANT_IEC_60559_TYPES_EXT__

#define __STDC_WANT_DEC_FP__

#define __STDC_WANT_IEC_60559_DFP_EXT__

#include <fenv.h>

#include <stdio.h> // reg. e.g. printf,

#include <float.h>

#include <limits.h>

#include <math.h> // reg. e.g. pow,

#include <stdlib.h> // reg. e.g. atof,

#include <time.h> // reg. e.g. clock(),

#include <locale.h> // reg. formatted print of integers, not yet sufficient,

// #include <decimal.h> // reg. ???

#include <string.h> // reg. e.g. strcat,

#include <quadmath.h> // reg. e.g. quadmath_snprintf,

clock_t start1, end1;

#define TIMEITcf( expr, N ) \

start1 = clock(); \

for( int i = 1; i < N; ++i ) \

{ \

expr; \

} \

end1 = clock(); \

printf( "%07d; %d; %.9E; %s \n", end1 - start1, N, expr, #expr )

#define TIMEITcd( expr, N ) \

start1 = clock(); \

for( int i = 1; i < N; ++i ) \

{ \

expr; \

} \

end1 = clock(); \

printf( "%07d; %d; %.18E; %s \n", end1 - start1, N, expr, #expr )

#define TIMEITcDF( expr, N ) \

start1 = clock(); \

for(int i = 1; i < N; ++i) \

{ \

expr; \

} \

end1 = clock(); \

printf( "%07d; %d; %.8HE; %s \n", end1 - start1, N, expr, #expr )

#define TIMEITcDD( expr, N ) \

start1 = clock(); \

for(int i = 1; i < N; ++i) \

{ \

expr; \

} \

end1 = clock(); \

printf( "%07d; %d; %.17DE; %s \n", end1 - start1, N, expr, #expr )

int main( int argc, char *argv[] )

{

fe_dec_setround( 4 ); // round ties away from zero for decimal datatypes,

// setlocale(LC_ALL, "en_US"); // or any other locale that supports thousands separators

volatile float x1f = 0.0, x2f = 0.0, x3f = 0.0;

volatile double x1d = 0.0, x2d = 0.0, x3d = 0.0;

volatile _Decimal32 x1DF = 0.0DF, x2DF = 0.0DF, x3DF = 0.0DF;

volatile _Decimal64 x1DD = 0.0DD, x2DD = 0.0DD, x3DD = 0.0DD;

volatile int count = 1000000; // how many times to run the loop,

if( argv[ 3 ] ) count = atoi( argv[ 3 ] );

printf("add two values from command line arguments \n" );

printf("no benefit from '-O2'? \n" );

x1f = strtof( argv[ 1 ], NULL );

x2f = strtof( argv[ 2 ], NULL );

TIMEITcf( x1f = x1f + x2f, count );

x1d = strtod( argv[ 1 ], NULL );

x2d = strtod( argv[ 2 ], NULL );

TIMEITcd( x1d = x1d + x2d, count );

x1DF = strtod32( argv[ 1 ], NULL );

x2DF = strtod32( argv[ 2 ], NULL );

TIMEITcDF( x1DF = x1DF + x2DF, count );

x1DD = strtod64( argv[ 1 ], NULL );

x2DD = strtod64( argv[ 2 ], NULL );

TIMEITcDD( x1DD = x1DD + x2DD, count );

printf(" \n" );

printf("tangent of value from command line argument \n" );

printf("no benefit from '-O2'? \n" );

x1f = strtof( argv[ 1 ], NULL );

TIMEITcf( x3f = tan( x1f ), count );

x1d = strtod( argv[ 1 ], NULL );

TIMEITcd( x3d = tan( x1d ), count );

x1DF = strtod32( argv[ 1 ], NULL );

TIMEITcDF( x3DF = tand32( x1DF ), count );

x1DD = strtod64( argv[ 1 ], NULL );

TIMEITcDD( x3DD = tand64( x1DD ), count );

printf(" \n" );

return 0;

}

</pre>

== See also ==

* [[ISO/IEC 10967]], Language Independent Arithmetic

* [[Primitive data type]]

* [[Q notation (scientific notation)|D (E) notation (scientific notation)]]

== References ==