412073b367938792ec2e9b6a49da385077a8e5b1
1 # IEEE Floating Point Adder (Single Precision)
2 # Copyright (C) Jonathan P Dawson 2013
5 from nmigen
import Module
, Signal
, Cat
6 from nmigen
.cli
import main
10 def __init__(self
, width
):
13 self
.in_a
= Signal(width
)
14 self
.in_a_stb
= Signal()
15 self
.in_a_ack
= Signal()
17 self
.in_b
= Signal(width
)
18 self
.in_b_stb
= Signal()
19 self
.in_b_ack
= Signal()
21 self
.out_z
= Signal(width
)
22 self
.out_z_stb
= Signal()
23 self
.out_z_ack
= Signal()
25 s_out_z_stb
= Signal()
26 s_out_z
= Signal(width
)
30 def get_fragment(self
, platform
):
33 a
= Signal(self
.width
)
34 b
= Signal(self
.width
)
35 z
= Signal(self
.width
)
52 with m
.State("get_a"):
53 with m
.If((self
.in_a_ack
) & (self
.in_a_stb
)):
60 m
.d
.sync
+= self
.in_a_ack
.eq(1)
62 with m
.State("get_b"):
63 with m
.If((self
.in_b_ack
) & (self
.in_b_stb
)):
70 m
.d
.sync
+= self
.in_b_ack
.eq(1)
83 if (s_in_a_ack && in_a_stb) begin
93 if (s_in_b_ack && in_b_stb) begin
102 a_m <= {a[22 : 0], 3'd0};
103 b_m <= {b[22 : 0], 3'd0};
104 a_e <= a[30 : 23] - 127;
105 b_e <= b[30 : 23] - 127;
108 state <= special_cases;
113 //if a is NaN or b is NaN return NaN
114 if ((a_e == 128 && a_m != 0) || (b_e == 128 && b_m != 0)) begin
120 //if a is inf return inf
121 end else if (a_e == 128) begin
125 //if a is inf and signs don't match return nan
126 if ((b_e == 128) && (a_s != b_s)) begin
133 //if b is inf return inf
134 end else if (b_e == 128) begin
139 //if a is zero return b
140 end else if ((($signed(a_e) == -127) && (a_m == 0)) && (($signed(b_e) == -127) && (b_m == 0))) begin
142 z[30:23] <= b_e[7:0] + 127;
143 z[22:0] <= b_m[26:3];
145 //if a is zero return b
146 end else if (($signed(a_e) == -127) && (a_m == 0)) begin
148 z[30:23] <= b_e[7:0] + 127;
149 z[22:0] <= b_m[26:3];
151 //if b is zero return a
152 end else if (($signed(b_e) == -127) && (b_m == 0)) begin
154 z[30:23] <= a_e[7:0] + 127;
155 z[22:0] <= a_m[26:3];
158 //Denormalised Number
159 if ($signed(a_e) == -127) begin
164 //Denormalised Number
165 if ($signed(b_e) == -127) begin
176 if ($signed(a_e) > $signed(b_e)) begin
179 b_m[0] <= b_m[0] | b_m[1];
180 end else if ($signed(a_e) < $signed(b_e)) begin
183 a_m[0] <= a_m[0] | a_m[1];
192 if (a_s == b_s) begin
196 if (a_m >= b_m) begin
213 sticky <= tot[1] | tot[0];
221 state <= normalise_1;
226 if (z_m[23] == 0 && $signed(z_e) > -126) begin
233 state <= normalise_2;
239 if ($signed(z_e) < -126) begin
244 sticky <= sticky | round_bit;
252 if (guard && (round_bit | sticky | z_m[0])) begin
254 if (z_m == 24'hffffff) begin
263 z[22 : 0] <= z_m[22:0];
264 z[30 : 23] <= z_e[7:0] + 127;
266 if ($signed(z_e) == -126 && z_m[23] == 0) begin
269 if ($signed(z_e) == -126 && z_m[23:0] == 24'h0) begin
270 z[31] <= 1'b0; // FIX SIGN BUG: -a + a = +0.
272 //if overflow occurs, return inf
273 if ($signed(z_e) > 127) begin
285 if (s_out_z_stb && out_z_ack) begin
301 assign in_a_ack = s_in_a_ack;
302 assign in_b_ack = s_in_b_ack;
303 assign out_z_stb = s_out_z_stb;
304 assign out_z = s_out_z;
309 if __name__
== "__main__":
310 alu
= FPADD(width
=32)
312 alu
.in_a
, alu
.in_a_stb
, alu
.in_a_ack
,
313 alu
.in_b
, alu
.in_b_stb
, alu
.in_b_ack
,
314 alu
.out_z
, alu
.out_z_stb
, alu
.out_z_ack
,