ECE 645: Lecture 5 Number Representation

Transcription

1 ECE 645: Lecture 5 Number Representation Part 2 Little-Endian vs. Big-Endian Representations Floating Point Representations Rounding Representation of the Galois Field elements

2 Required Reading Endianness, from Wikipedia, the free encyclopedia Behrooz Parhami, Computer Arithmetic: Algorithms and Hardware Design Chapter 17, Floating-Point Representations

3 Little-Endian vs. Big-Endian Representation of Integers

4 Little-Endian vs. Big-Endian Representation A0 B1 C2 D3 E4 F MSB Big-Endian LSB Little-Endian MSB = A0 B1 C2 D3 E4 F5 67 LSB = 89 0 address MAX LSB = F5 E4 D3 C2 B1 MSB = A0

5 Little-Endian vs. Big-Endian Camps MSB... LSB Big-Endian 0 address MAX LSB... MSB Little-Endian Motorola 68xx, 680x0 IBM Hewlett-Packard Sun SuperSPARC Internet TCP/IP Bi-Endian Motorola Power PC Silicon Graphics MIPS Intel AMD DEC VAX RS 232

6 Little-Endian vs. Big-Endian Origin of the terms Jonathan Swift, Gulliver s Travels A law requiring all citizens of Lilliput to break their soft-eggs at the little ends only A civil war breaking between the Little Endians and the Big-Endians, resulting in the Big Endians taking refuge on a nearby island, the kingdom of Blefuscu Satire over holy wars between Protestant Church of England and the Catholic Church of France

7

8 Little-Endian vs. Big-Endian Advantages and Disadvantages Big-Endian easier to determine a sign of the number Little-Endian easier addition and multiplication of multiprecision numbers easier to compare two numbers easier to divide two numbers easier to print

9 Pointers (1) Big-Endian Little-Endian 0 address MAX F5 E4 D3 C2 B1 A0 int * iptr; (* iptr) = 8967; (* iptr) = 6789; iptr+1

10 Pointers (2) Big-Endian Little-Endian 0 address MAX F5 E4 D3 C2 B1 A0 long int * lptr; (* lptr) = 8967F5E4; (* lptr) = E4F56789; lptr + 1

11 Floating Point Representations

12

13

14

15

16 The ANSI/IEEE standard floatingpoint number representation formats Short (32-bit) format Originally IEEE Superseded by IEEE Standard. 8 bits, bias = 127, ñ126 to 127 Sign Exponent 23 bits for fractional part (plus hidden 1 in integer part) Significand 11 bits, bias = 1023, ñ1022 to bits for fractional part (plus hidden 1 in integer part) Long (64-bit) format

17

18

19 Table 17.1 Some features of the ANSI/IEEE standard floatingpoint number representation formats

20 Exponent Encoding Exponent encoding in 8 bits for the single/short (32-bit) ANSI/IEEE format Decimal code Hex code E 7F 80 FE FF Exponent value f 2 e f = 0: Representation of ±0 f 0: Representation of denormals, 0.f f = 0: Representation of ± f 0: Representation of NaNs

21 Fig Denormals in the IEEE single-precision format.

22

23

24

25 New IEEE Standard Basic Formats

26 New IEEE Standard Binary Interchange Formats

27 Requirements for Arithmetic Results of the 4 basic arithmetic operations (+, -,, ) as well as squarerooting must match those obtained if all intermediate computations were infinitely precise That is, a floating-point arithmetic operation should introduce no more imprecision than the error attributable to the final rounding of a result that has no exact representation (this is the best possible) Example: ( ) ( ) Exact result Rounded result Error = ½ ulp

28 Rounding 101

29 Rounding Modes The IEEE standard includes five rounding modes: Default: Round to nearest, ties to even (rtne) Optional: Round to nearest, ties away from 0 (rtna) Round toward zero (inward) Round toward + Round toward (upward) (downward)

30 Required Reading Parhami, Chapter 17.5, Rounding schemes Rounding Algorithms

31 Rounding Rounding occurs when we want to approximate a more precise number (i.e. more fractional bits L) with a less precise number (i.e. fewer fractional bits L') Example 1: old: (K=6, L=8) new: (K'=6, L'=2) Example 2: old: (K=6, L=8) new: (K'=6, L'=0) The following pages show rounding from L>0 fractional bits to L'=0 bits, but the mathematics hold true for any L' < L Usually, keep the number of integral bits the same K'=K 31

32 Rounding Equation Whole part Fractional part x k 1 x k 2... x 1 x 0. x 1 x 2... x Round l y k 1 y k 2... y 1 y 0 y = round(x) 32

33 Rounding Techniques There are different rounding techniques: 1) truncation results in round towards zero in signed magnitude results in round towards - in two's complement 2) round to nearest number 3) round to nearest even number (or odd number) 4) round towards + Other rounding techniques 5) jamming or von Neumann 6) ROM rounding Each of these techniques will differ in their error depending on representation of numbers i.e. signed magnitude versus two's complement Error = round(x) x 33

34 1) Truncation The simplest possible rounding scheme: chopping or truncation x k 1 x k 2... x 1 x 0. x 1 x 2... x trunc l x k 1 x k 2... x 1 x 0 ulp Truncation in signed-magnitude results in a number chop(x) that is always of smaller magnitude than x. This is called round towards zero or inward rounding (3.5) 10 à 011 (3) 10 Error = (-3.5) 10 à 111 (-3) 10 Error = +0.5 Truncation in two's complement results in a number chop(x) that is always smaller than x. This is called round towards - or downward-directed rounding (3.5) 10 à 011 (3) 10 Error = (-3.5) 10 à 100 (-4) 10 Error =

35 Truncation Function Graph: chop(x) chop( x ) Fig Truncation or chopping of a signed-magnitude number (same as round toward 0). x chop( x ) Fig Truncation or chopping of a 2 s-complement number (same as round to - ). x 35

36 Bias in two's complement truncation X (binary) X (decimal) chop(x) (binary) chop(x) (decimal) Error (decimal) Assuming all combinations of positive and negative values of x equally possible, average error is In general, average error = -(2 -L' -2 -L )/2, where L' = new number of fractional bits 36

37 Implementation truncation in hardware Easy, just ignore (i.e. truncate) the fractional digits from L to L'+1 x k-1 x k-2.. x 1 x 0. x -1 x -2.. x -L = y k-1 y k-2.. y 1 y 0. ignore (i.e. truncate the rest) 37

38 2) Round to nearest number Rounding to nearest number what we normally think of when say round (2.25) 10 à 010 (2) 10 Error = (2.75) 10 à 011 (3) 10 Error = (2.00) 10 à 010 (2) 10 Error = (2.5) 10 à 011 (3) 10 Error = +0.5 [round-half-up (arithmetic rounding)] (2.5) 10 à 010 (2) 10 Error = -0.5 [round-half-down] 38

39 Round-half-up: dealing with negative numbers Rounding to nearest number what we normally think of when say round (-2.25) 10 à 110 (-2) 10 Error = (-2.75) 10 à 101 (-3) 10 Error = (-2.00) 10 à 110 (-2) 10 Error = (-2.5) 10 à 110 (-2) 10 Error = +0.5 [asymmetric implementation] (-2.5) 10 à 101 (-3) 10 Error = -0.5 [symmetric implementation] 39

40 Round to Nearest Function Graph: rtn(x) Round-half-up version Asymmetric implementation Symmetric implementation rtn(x) ñ4 ñ3 ñ2 ñ ñ1 ñ2 ñ3 ñ4 x rtn(x) ñ4 ñ3 ñ2 ñ ñ1 ñ2 ñ3 ñ4 x 40

41 Bias in two's complement round to nearest Round-half-up asymmetric implementation X (binary) X (decimal) rtn(x) (binary) rtn(x) (decimal) Error (decimal) Assuming all combinations of positive and negative values of x equally possible, average error is Smaller average error than truncation, but still not symmetric error We have a problem with the midway value, i.e. exactly at 2.5 or -2.5 leads to positive error bias always Also have the problem that you can get overflow if only allocate K' = K integral bits Example: rtn(011.10) à overflow This overflow only occurs on positive numbers near the maximum positive value, not on negative numbers 41

42 Implementing round to nearest (rtn) in hardware Round-half-up asymmetric implementation Two methods Method 1: Add '1' in position one digit right of new LSB (i.e. digit L'+1) and keep only L' fractional bits x k-1 x k-2.. x 1 x 0. x -1 x -2.. x -L + 1 = y k-1 y k-2.. y 1 y 0. y -1 Method 2: Add the value of the digit one position to right of new LSB (i.e. digit L'+1) into the new LSB digit (i.e. digit L) and keep only L' fractional bits x k-1 x k-2.. x 1 x 0. x -1 x -2.. x -L + x -1 y k-1 y k-2.. y 1 y 0. ignore (i.e. truncate the rest) ignore (i.e truncate the rest) 42

43 Round to Nearest Even Function Graph: rtne(x) To solve the problem with the midway value we implement round to nearest-even number (or can round to nearest odd number) rtne(x) R*(x) ñ4 ñ3 ñ2 ñ ñ1 ñ2 ñ3 ñ4 Fig Rounding to the nearest even number x ñ4 ñ3 ñ2 ñ ñ1 ñ2 ñ3 ñ4 Fig R* rounding or rounding to the nearest odd number. x 43

44 Bias in two's complement round to nearest even (rtne) X (binary) X (decimal ) rtne(x) (binary) rtne(x) (decimal) Error (decimal) (overfl) average error is now 0 (ignoring the overflow) cost: more hardware 44

45 4) Rounding towards infinity We may need computation errors to be in a known direction Example: in computing upper bounds, larger results are acceptable, but results that are smaller than correct values could invalidate upper bound Use upward-directed rounding (round toward + ) up(x) always larger than or equal to x Similarly for lower bounds, use downward-directed rounding (round toward - ) down(x) always smaller than or equal to x We have already seen that round toward - in two's complement can be implemented by truncation 45

46 Rounding Toward Infinity Function Graph: up(x) and down(x) up(x) down(x) down(x) can be implemented by chop(x) in two's complement 46

47 Two's Complement Round to Zero Two's complement round to zero (inward rounding) also exists inward( x ) x 47

48 Other Methods Note that in two's complement round to nearest (rtn) involves an addition which may have a carry propagation from LSB to MSB Rounding may take as long as an adder takes Can break the adder chain using the following two techniques: Jamming or von Neumann ROM-based 48

49 5) Jamming or von Neumann jam(x) ñ4 ñ3 ñ2 ñ ñ1 ñ2 ñ3 x Chop and force the LSB of the result to 1 Simplicity of chopping, with the near-symmetry or ordinary rounding Max error is comparable to chopping (double that of rounding) ñ4 49

50 6) ROM Rounding ROM(x) ñ4 ñ3 ñ2 ñ ñ1 ñ2 ñ3 ñ4 x Fig ROM rounding with an 8 2 table. Example: Rounding with a 32 4 table Rounding result is the same as that of the round to nearest scheme in 31 of the 32 possible cases, but a larger error is introduced when x 3 = x 2 = x 1 = x 0 = x 1 = 1 x k 1... x 4 x 3 x 2 x 1 x 0. x 1 x 2... x l ROM x k 1... x 4 y 3 y 2 y 1 y 0 ROM address ROM data 50

51 Representation of the Galois Field elements

52 Evariste Galois ( )

53 Evariste Galois ( ) Studied the problem of finding algebraic solutions for the general equations of the degree 5, e.g., f(x) = a 5 x 5 + a 4 x 4 + a 3 x 3 + a 2 x 2 + a 1 x+ a 0 = 0 Answered definitely the question which specific equations of a given degree have algebraic solutions. On the way, he developed group theory, one of the most important branches of modern mathematics.

54 Evariste Galois ( ) 1829 Galois submits his results for the first time to the French Academy of Sciences Reviewer 1 Augustin-Luis Cauchy forgot or lost the communication Galois submits the revised version of his manuscript, hoping to enter the competition for the Grand Prize in mathematics Reviewer 2 Joseph Fourier died shortly after receiving the manuscript Third submission to the French Academy of Sciences Reviewer 3 Simeon-Denis Poisson did not understand the manuscript and rejected it.

55 Evariste Galois ( ) May 1832 Galois provoked into a duel The night before the duel he wrote a letter to his friend containing the summary of his discoveries. The letter ended with a plea: Eventually there will be, I hope, some people who will find it profitable to decipher this mess. May 30, 1832 Galois was grievously wounded in the duel and died in the hospital the following day Galois manuscript rediscovered by Joseph Liouville 1846 Galois manuscript published for the first time in a mathematical journal.

56

57 Field Set F, and two operations typically denoted by (but not necessarily equivalent to) + and * Set F, and definitions of these two operations must fulfill special conditions.

58 Infinite fields Finite fields Examples of fields { R= set of real numbers, + addition of real numbers * multiplication of real numbers } { set Zp={0, 1, 2,, p-1}, + (mod p): addition modulo p, * (mod p): multiplication modulo p }

59 Finite Fields = Galois Fields GF(p m ) p prime p m number of elements in the field Arithmetic operations present in many libraries GF(p) GF(2 m ) Polynomial basis representation Most significant special cases Normal basis representation Fast in hardware Fast squaring

60 Quotient and remainder Given integers a and n, n>0! q, r Z such that a = q n + r and 0 r < n q quotient r remainder (of a divided by n) q = a n = a div n r = a - q n = a = a mod n a n n =

61 32 mod 5 = -32 mod 5 =

62 Integers coungruent modulo n Two integers a and b are congruent modulo n (equivalent modulo n) written a b iff a mod n = b mod n or a = b + kn, k Z or n a - b

63 Laws of modular arithmetic

64 Rules of addition, subtraction and multiplication modulo n a + b mod n = ((a mod n) + (b mod n)) mod n a - b mod n = ((a mod n) - (b mod n)) mod n a b mod n = ((a mod n) (b mod n)) mod n

65 9 13 mod 5 = mod 26 =

66 Laws of modular arithmetic Regular addition a+b = a+c iff b=c Regular multiplication If a b = a c and a 0 then b = c Modular addition a+b a+c (mod n) iff b c (mod n) Modular multiplication If a b a c (mod n) and gcd (a, n) = 1 then b c (mod n)

67 Modular Multiplication: Example (mod 8) (mod 8) 3 7 (mod 8) x 6 x mod x 5 x mod

68 Finite Fields = Galois Fields GF(p m ) p prime p m number of elements in the field Arithmetic operations present in many libraries GF(p) GF(2 m ) Polynomial basis representation Most significant special cases Normal basis representation Fast in hardware Fast squaring

69 Elements of the Galois Field GF(2 m ) Binary representation (used for storing and processing in computer systems): A = (a m-1, a m-2,, a 2, a 1, a 0 ) a i {0, 1} Polynomial representation (used for the definition of basic arithmetic operations): m-1 A(x) = a i x i = a m-1 x m-1 + a m-2 x m a 2 x 2 + a 1 x+a 0 i=0 multiplication + addition modulo 2 (XOR)

70 Addition and Multiplication in the Galois Field GF(2 m ) Inputs A = (a m-1, a m-2,, a 2, a 1, a 0 ) B = (b m-1, b m-2,, b 2, b 1, b 0 ) a i, b i {0, 1} Output C = (c m-1, c m-2,, c 2, c 1, c 0 ) c i {0, 1}

71 Addition Addition in the Galois Field GF(2 m ) A A(x) B B(x) C C(x) = A(x) + B(x) = = (a m-1 +b m-1 ) x m-1 + (a m-2 +b m-2 ) x m (a 2 +b 2 ) x 2 + (a 1 +b 1 ) x + (a 0 +b 0 ) = = c m-1 x m-1 + c m-2 x m c 2 x 2 + c 1 x+c 0 multiplication + addition modulo 2 (XOR) c i = a i + b i = a i XOR b i C = A XOR B

72 Multiplication in the Galois Field GF(2 m ) Multiplication A A(x) B B(x) C C(x) = A(x) B(x) mod P(X) = c m-1 x m-1 + c m-2 x m c 2 x 2 + c 1 x+c 0 P(x) - irreducible polynomial of the degree m P(x) = p m x m + p m-1 x m p 2 x 2 + p 1 x+p 0

73