1. Technical Field
The present disclosure relates to digital arithmetic circuits in general, and, in particular, to fused multiply-adder circuits.
2. Description of Related Art
The institute of Electrical and Electronics Engineers (IEEE) standard for floating-point arithmetic defines specific formats for representing floating-point numbers. According to the IEEE standard, a floating-point number includes a sign bit, an exponent, and a fraction. The value of a floating-point number X is represented by:
X=(−1)Xs*(1.Xf)*2(Xe-bias)
where Xs is a sign bit, 1.Xf is intended to represent the binary number created by prefixing Xf with an implicit leading 1 and a binary point, and Xe is the unsigned binary exponent. If Xe=0, X is considered a denormalized number, and its value is instead represented by:
X=(−1)Xs*(0.Xf)*2(1-bias)
The IEEE standard also defines floating-point numbers in multiple levels of precisions. For example, a single-precision floating-point number has an 8-bit exponent and a 23-bit fraction, a double-precision floating-point number has an 11-bit exponent and a 52-bit fraction, and a quadruple-precision floating-point number has a 15-bit exponent and a 112-bit fraction.
Modern computer processors typically include a floating-point unit to perform mathematical operations on floating-point numbers according to the IEEE standard. One important floating-point instruction is the multiply-add instruction that implements the operation
T=A*B+C
in one step with only one rounding error (instead of two that would result from executing a multiply instruction followed by an add instruction). Two different approaches have been used for implementing the fused multiply-add instruction to support multiple precisions of floating-point numbers. The first approach uses separate data paths for each number precision, so instructions using single- and double-precision numbers, for example, can be executed at the same time, but at the expense of a larger chip area. The second approach uses only one data path capable of handling both single- and double-precision numbers, but only half the operand bandwidth is utilized when handling single-precision numbers.
Consequently, it would be desirable to provide an improved fused multiply-adder for performing multiply-add instructions.
In accordance with a preferred embodiment of the present disclosure, a fused multiply-adder includes a Booth encoder, a fraction multiplier, a carry corrector, and an adder. Operands for the fused multiply-adder can be full-precision numbers and/or lesser precision numbers having at least one zero guard bit separating individual numbers. The Booth encoder initially encodes a first operand. The fraction multiplier multiplies the Booth-encoded first operand by a second operand to produce partial products, and then reduces the partial products into a set of redundant sum and carry vectors. The carry corrector then generates a carry correction factor for correcting the carry vectors. The adder adds the redundant sum and carry vectors and the carry correction factor to a third operand to yield a final result of a multiply-add operation.
All features and advantages of the present disclosure will become apparent in the following detailed written description.
The disclosure itself, as well as a preferred mode of use, further objects, and advantages thereof, will best be understood by reference to the following detailed description of an illustrative embodiment when read in conjunction with the accompanying drawings, wherein:
A multiply-add instruction implements the equation
T=A*B+C
where A, B and C are three separate operands. For ease of illustration, the method of the present invention will be illustrated using two single-precision operations computed simultaneously in double-precision hardware. However, it is understood by those skilled in the art that the method of the present invention can readily be extended to k-way operations, each computed in n/k bit precision, where n is the maximum precision of the underlying hardware. Thus, for example, four single-precision or two double-precision operations can be simultaneously supported by quadruple-precision hardware.
Referring to the drawings and in particular to
Fused multiply-adder 10 can operate in full-precision mode or two-way mode. With full-precision mode, each of A register 11, B register 12, and C register 13 contains an operand occupying the full width of the register. With two-way mode, each of A register 11, B register 12, and C register 13 contains two operands that are separated by sufficient guard zeros to allow the partial products P of operand A * operand B to retain at least one guard zero between its constituent parts. With full-precision mode, fused multiply-adder 10 outputs one full-precision number. With two-way mode, fused multiply-adder 10 outputs two numbers of the same precision as their respective operands.
When operating in full-precision mode, the Booth-encoding of one set of input multiplicands will not lead to any errors in the final result. However, when operating in two-way mode, the Booth-encoding of one set of input multiplicands may cause error in the final result unless a set of carry correction factors is applied to correct the carry vectors within adder 17. Carry corrector 16 provides the computations in order to generate the necessary carry correction factors within adder 17.
Adder 17 then adds the carry and sum vectors and the carry correction factors to an operand within C register 13, which has been correctly aligned via aligner 19. After normalizing and rounding the result from adder 17 via normalizer/rounder 18, fused multiply-adder 10 outputs a floating-point result of A*B+C.
Under two-way mode, each of carry (C) and sum (S) vectors is split into lower halves, CLO, SLO, and higher halves, CHI, SHI, with at least one guard zero in between. Thus, the half-width partial products may be defined as
P
HI
=A
HI
*B
HI
=C
HI
+S
HI
P
LO
=A
LO
*B
LO
=C
LO
+S
LO
resulting in the full-width partial products as:
P=P
HI
+P
LO
=C
HI
+S
HI
+C
LO
+S
LO
where the HI and LO terms do not overlap. The boundary between the HI and LO terms is the same for all terms, and is chosen to ensure that there is at least one guard zero between the HI and LO terms.
When the sum and carry vectors are added, the result becomes
S
LO
+C
LO
=P
LO
+r (1)
S
HI
+C
HI
=P
HI
−r (2)
where r is a carry bit produced from the addition of the lower halves, as depicted in
However, r is an artifact of the Booth encoding of an operand within B register 12 and the addition of the subsequent partial products. Thus, r should not be included in the final calculation result, and must be corrected for. Also, the correction cannot be delayed until r becomes known as it will ripple from the lower half to the higher half during the addition operation. Since r must be added and subtracted before it is known, it must be predicted.
M can be defined as the most significant bit position in SLO and CLO, so that r has the significance M−1 (for the present embodiment, 0 is the most significant bit). The fact that the boundary between the higher half and lower half was chosen such that Pm=0 can be used to predict r in all cases, without waiting for r to be computed.
Case 1: If SM+CM=0, a possible carry generated in the lower order bits of SLO+CLO can never ripple beyond bit position M. Hence r=0.
Case 2: If SM+CM=1, PM must be 0 because the boundary was chosen that way, and r=1.
Case 3: If SM+CM=2, bit position M is generating the carry out from SLO+CLO, so r=1. Any possible carry generated in the lower order bits of SLO CLO can never ripple beyond position M.
The above-mentioned three cases can be reduced to
r=S
M OR CM.
Hence, r can be predicted solely by examining SM and CM without waiting for actual calculation results.
Equations (1) and (2) can be rearranged into
P
LO
=S
LO
+C
LO
−r
P
HI
=S
HI
+C
HI
+r
In a fused multiply-add pipeline, partial products P are not computed explicitly. Instead the aligned addend fraction is added (or subtracted) before a full carry propagating addition begins. So, in full-precision mode, adder 17 (from
sum=S+C±addend
In two-way mode, taking r into account, adder 17 instead computes
sumHI=SHI+CHI+r±addend (for the higher half)
sumLO=SLO+CLO−r±addend (for the lower half)
A carry-save adder is used to compress S+C+addend=ss+cc, where ss and cc are redundant sum and carry vectors, respectively. In this format, the equations for each individual bit i are
ss
i
=S
i
XOR C
i
XOR addendi
cc
i-1=(Si AND Ci) OR (Si AND addendi) OR (Ci AND addend)
where bit position 0 is the most significant bit position.
In two-way mode, M is the most significant bit position of the lower half and M−1 is the least significant bit position of the higher half. Since data in the lower half does not influence the data in the higher half, ccM-1 is set to “0” in two-way mode; instead the predicted r can be inserted in this bit position.
In full-precision mode, from the general equation above,
cc
M-1=(SM AND CM) OR (SM AND addendM) OR (CM AND addendM).
But in two-way mode, the desired result is
cc
M-1
=r=S
M OR CM.
This can be rewritten in the same form as the full-precision equation
cc
M-1=(SM AND CM) OR (SM AND 1) OR (CM AND 1)
Both forms can be generalized into
cc
M-1=(SM AND CM) OR (SM AND (addendM OR 2-way))
OR (CM AND (addendM OR 2-way)) (3)
where 2-way is a control signal set to 1 if this is a two-way mode operation and 0 otherwise. Equation (3) is valid for two-way mode as well as full-precision mode.
Hence, by substituting the ccM-1 output from the carry-save adder with a circuit generating
cc
M-1=(SM AND CM) OR (SM AND (addendM OR 2-way))
OR (CM AND (addendM OR 2-way))
during addend alignment, off the critical multiplier-adder path, the higher-half result of the dual-mode multiply-add operation is now correct, as illustrated in
As previously described, bit r has a weight of M−1, which is the carry-out position of the lower half (or sumLO) addition. Without correction, this carry-out bit will be wrong. The correct carry-out bit is denoted as coLO.
The early carry-out can be computed using the general equation
cc
M-1LO=(SM AND CM) OR (SM AND addendM) OR (CM AND addendM)
while during the carry-propagating addition of ssLO+ccLO, a late carry-out co′LO is produced. All of these elements are related by the equation
cc
M-1LO
+co′
LO
=r+co
LO
which can be rearranged into
co
LO
=cc
M-1LO
+co′
LO
−r
Then, the following three cases must be accounted for:
Case 1: If addendM=1, by the general equation
cc
M-1LO=(SM AND CM) OR (SM AND 1) OR (CM AND 1)=SM OR CM=r
then
co
LO
=r+co′
LO
−r=co′
LO
which can also be written as
co
LO
=co′
LO AND addendM
since addendM=1 and the AND operation with 1 has no effect.
Case 2: If addendM=0 and SM+CM< >1, by the general equation
cc
M-1LO=(SM AND CM) OR (SM AND 0) OR (CM AND 0)=SM OR CM=r
then as before,
co
LO
=r+co′
LO
−r=co′
LO.
But here, co′LO will be 0, because as the guard bit, the sum of SM and CM must be 0, and addendM=0. Hence,
co
LO
=co′
LO=0=co′LO AND addendM.
Case 3: If addendM=0 and SM+CM=1, by the general equation
cc
M-1LO=(SM AND CM) OR (SM AND addendM) OR (CM AND addendM)=0
i but
r=S
M OR CM=1.
Hence,
cc
LO
=cc
M-1LO
+co′
LO
−r=0+cc′LO−1
which is only possible if co′LO=1 and coLO=0. As a result, once again
co
LO
=co′
LO AND addendM
is true.
Thus, in all three cases, the correct carry out coLO from the lower half of the dual-mode multiply-add operation can be computed by ignoring ccM-1LO and instead using the result of co′LO and addendM, as illustrated in
With reference now to
As has been described, the present invention provides an improved fused multiply-adder. By incorporating appropriate corrections during the addition phase, it is possible to implement a fused multiply-adder that supports multiple-precision floating-point operands with a Booth-encoded multiplication operand.
While the disclosure has been particularly shown and described with reference to a preferred embodiment, it will be understood by those skilled in the art that various changes in form and detail may be made therein without departing from the spirit and scope of the disclosure.