1 |
118 |
jguarin200 |
------------------------------------------------
|
2 |
119 |
jguarin200 |
--! @file fadd32.vhd
|
3 |
118 |
jguarin200 |
--! @brief RayTrac Floating Point Adder
|
4 |
|
|
--! @author Julián Andrés Guarín Reyes
|
5 |
|
|
--------------------------------------------------
|
6 |
|
|
|
7 |
|
|
|
8 |
|
|
-- RAYTRAC (FP BRANCH)
|
9 |
|
|
-- Author Julian Andres Guarin
|
10 |
119 |
jguarin200 |
-- fadd32.vhd
|
11 |
118 |
jguarin200 |
-- This file is part of raytrac.
|
12 |
|
|
--
|
13 |
|
|
-- raytrac is free software: you can redistribute it and/or modify
|
14 |
|
|
-- it under the terms of the GNU General Public License as published by
|
15 |
|
|
-- the Free Software Foundation, either version 3 of the License, or
|
16 |
|
|
-- (at your option) any later version.
|
17 |
|
|
--
|
18 |
|
|
-- raytrac is distributed in the hope that it will be useful,
|
19 |
|
|
-- but WITHOUT ANY WARRANTY; without even the implied warranty of
|
20 |
|
|
-- MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
|
21 |
|
|
-- GNU General Public License for more details.
|
22 |
|
|
--
|
23 |
|
|
-- You should have received a copy of the GNU General Public License
|
24 |
|
|
-- along with raytrac. If not, see <http://www.gnu.org/licenses/>
|
25 |
|
|
|
26 |
|
|
library ieee;
|
27 |
|
|
use ieee.std_logic_1164.all;
|
28 |
|
|
use ieee.std_logic_unsigned.all;
|
29 |
|
|
use ieee.std_logic_arith.all;
|
30 |
119 |
jguarin200 |
library lpm;
|
31 |
|
|
use lpm.all;
|
32 |
118 |
jguarin200 |
|
33 |
|
|
--! Esta entidad recibe dos números en formato punto flotante IEEE 754, de precision simple y devuelve las mantissas signadas y corridas, y el exponente correspondiente al resultado antes de normalizarlo al formato float.
|
34 |
|
|
--!\nLas 2 mantissas y el exponente entran despues a la entidad add2 que suma las mantissas y entrega el resultado en formato IEEE 754.
|
35 |
119 |
jguarin200 |
entity fadd32 is
|
36 |
118 |
jguarin200 |
port (
|
37 |
|
|
clk,dpc : in std_logic;
|
38 |
|
|
a32,b32 : in std_logic_vector (31 downto 0);
|
39 |
|
|
c32 : out std_logic_vector(31 downto 0)
|
40 |
|
|
);
|
41 |
119 |
jguarin200 |
end fadd32;
|
42 |
118 |
jguarin200 |
|
43 |
119 |
jguarin200 |
architecture fadd32_arch of fadd32 is
|
44 |
118 |
jguarin200 |
|
45 |
|
|
component lpm_mult
|
46 |
|
|
generic (
|
47 |
|
|
lpm_hint : string;
|
48 |
|
|
lpm_representation : string;
|
49 |
|
|
lpm_type : string;
|
50 |
|
|
lpm_widtha : natural;
|
51 |
|
|
lpm_widthb : natural;
|
52 |
|
|
lpm_widthp : natural
|
53 |
|
|
);
|
54 |
|
|
port (
|
55 |
|
|
dataa : in std_logic_vector ( lpm_widtha-1 downto 0 );
|
56 |
|
|
datab : in std_logic_vector ( lpm_widthb-1 downto 0 );
|
57 |
|
|
result : out std_logic_vector( lpm_widthp-1 downto 0 )
|
58 |
|
|
);
|
59 |
|
|
end component;
|
60 |
|
|
|
61 |
120 |
jguarin200 |
signal s1zero,s5tokena,s5tokenb,s5tokenc,s7sign : std_logic;
|
62 |
|
|
signal s5token : std_logic_vector(2 downto 0);
|
63 |
119 |
jguarin200 |
signal s1delta : std_logic_vector(5 downto 0);
|
64 |
|
|
signal s0delta,s1exp,s2exp,s3exp,s4exp,s5exp,s6exp,s5factor,s6factor,s7exp,s7factor : std_logic_vector(7 downto 0);
|
65 |
|
|
signal s1shifter,s5factorhot9,s6factorhot9 : std_logic_vector(8 downto 0);
|
66 |
|
|
signal s1pl,s6pl : std_logic_vector(17 downto 0);
|
67 |
|
|
signal s6postshift,s7postshift : std_logic_vector(22 downto 0);
|
68 |
|
|
signal s1umantshift,s1umantfixed,s1postshift,s1xorslab,s2xorslab : std_logic_vector(23 downto 0);
|
69 |
120 |
jguarin200 |
signal s5factorhot24 : std_logic_vector(23 downto 0);
|
70 |
119 |
jguarin200 |
signal s2umantshift,s2mantfixed,s3mantfixed,s3mantshift,s4xorslab : std_logic_vector(24 downto 0);
|
71 |
120 |
jguarin200 |
signal s4sresult,s5result,s6result : std_logic_vector(25 downto 0); -- Signed mantissa result
|
72 |
119 |
jguarin200 |
signal s1ph,s6ph : std_logic_vector(26 downto 0);
|
73 |
|
|
signal s0a,s0b : std_logic_vector(31 downto 0); -- Float 32 bit
|
74 |
118 |
jguarin200 |
|
75 |
|
|
begin
|
76 |
|
|
|
77 |
|
|
process (clk)
|
78 |
|
|
begin
|
79 |
|
|
if clk'event and clk='1' then
|
80 |
|
|
|
81 |
|
|
--!Registro de entrada
|
82 |
|
|
s0a <= a32;
|
83 |
|
|
s0b(31) <= dpc xor b32(31); --! Importante: Integrar el signo en el operando B
|
84 |
|
|
s0b(30 downto 0) <= b32(30 downto 0);
|
85 |
|
|
|
86 |
|
|
--!Etapa 0,Escoger el mayor exponente que sera el resultado desnormalizado, calcula cuanto debe ser el corrimiento de la mantissa con menor exponente y reorganiza los operandos, si el mayor es b, intercambia las posición si el mayor es a las posiciones la mantiene. Zero check.
|
87 |
|
|
--!signo,exponente,mantissa
|
88 |
|
|
if (s0b(30 downto 23)&s0a(30 downto 23))=x"0000" then
|
89 |
|
|
s1zero <= '0';
|
90 |
|
|
else
|
91 |
|
|
s1zero <= '1';
|
92 |
|
|
end if;
|
93 |
119 |
jguarin200 |
s1delta <= s0delta(7) & (s0delta(7) xor s0delta(4))&(s0delta(7) xor s0delta(3)) & s0delta(2 downto 0);
|
94 |
118 |
jguarin200 |
case s0delta(7) is
|
95 |
|
|
when '1' =>
|
96 |
|
|
s1exp <= s0b(30 downto 23);
|
97 |
|
|
s1umantshift <= s0a(31)&s0a(22 downto 0);
|
98 |
|
|
s1umantfixed <= s0b(31)&s0b(22 downto 0);
|
99 |
|
|
when others =>
|
100 |
|
|
s1exp <= s0a(30 downto 23);
|
101 |
|
|
s1umantshift <= s0b(31)&s0b(22 downto 0);
|
102 |
|
|
s1umantfixed <= s0a(31)&s0a(22 downto 0);
|
103 |
|
|
end case;
|
104 |
|
|
|
105 |
|
|
--! Etapa 1: Denormalización de la mantissas.
|
106 |
|
|
case s1delta(4 downto 3) is
|
107 |
|
|
when "00" => s2umantshift <= s1umantshift(23)&s1postshift(23 downto 0);
|
108 |
|
|
when "01" => s2umantshift <= s1umantshift(23)&x"00"&s1postshift(23 downto 8);
|
109 |
|
|
when "10" => s2umantshift <= s1umantshift(23)&x"0000"&s1postshift(23 downto 16);
|
110 |
|
|
when others => s2umantshift <= (others => '0');
|
111 |
|
|
end case;
|
112 |
|
|
s2mantfixed <= s1umantfixed(23) & ( ( ('1'&s1umantfixed(22 downto 0)) xor s1xorslab) + ( x"00000"&"000"&s1umantfixed(23) ) );
|
113 |
|
|
s2exp <= s1exp;
|
114 |
|
|
|
115 |
|
|
--! Etapa2: Signar la mantissa denormalizada.
|
116 |
|
|
s3mantfixed <= s2mantfixed;
|
117 |
|
|
s3mantshift <= s2umantshift(24)& ( ( s2umantshift(23 downto 0) xor s2xorslab) + ( x"00000"&"000"&s2umantshift(24) ) );
|
118 |
|
|
s3exp <= s2exp;
|
119 |
|
|
|
120 |
119 |
jguarin200 |
--! Etapa 3: Etapa 3 Realizar la suma, entre la mantissa corrida y la fija.
|
121 |
118 |
jguarin200 |
s4sresult <= (s3mantshift(24)&s3mantshift)+(s3mantfixed(24)&s3mantfixed);
|
122 |
|
|
s4exp <= s3exp;
|
123 |
|
|
|
124 |
|
|
--! Etapa 4: Quitar el signo a la mantissa resultante.
|
125 |
|
|
s5result <= s4sresult(25)&((s4sresult(24 downto 0) xor s4xorslab) +(x"000000"&s4sresult(25)));
|
126 |
|
|
s5exp <= s4exp;
|
127 |
|
|
|
128 |
|
|
|
129 |
|
|
--! Etapa 5: Codificar el corrimiento para la normalizacion de la mantissa resultante.
|
130 |
119 |
jguarin200 |
s6result <= s5result;
|
131 |
|
|
s6exp <= s5exp;
|
132 |
118 |
jguarin200 |
s6factor <= s5factor;
|
133 |
119 |
jguarin200 |
s6factorhot9 <= s5factorhot9;
|
134 |
118 |
jguarin200 |
|
135 |
119 |
jguarin200 |
--! Etapa 6: Ejecutar el corrimiento de la mantissa.
|
136 |
120 |
jguarin200 |
s7sign <= s6result(25);
|
137 |
119 |
jguarin200 |
s7exp <= s6exp;
|
138 |
120 |
jguarin200 |
s7factor <= s6factor;
|
139 |
119 |
jguarin200 |
s7postshift <= s6postshift;
|
140 |
|
|
|
141 |
|
|
--! Etapa 7: Entregar el resultado.
|
142 |
120 |
jguarin200 |
c32(31) <= s7sign;
|
143 |
119 |
jguarin200 |
c32(30 downto 23) <= s7exp+s7factor;
|
144 |
|
|
case s7factor(4 downto 3) is
|
145 |
|
|
when "01" => c32(22 downto 0) <= s7postshift(14 downto 00)&x"00";
|
146 |
|
|
when "10" => c32(22 downto 0) <= s7postshift(06 downto 00)&x"0000";
|
147 |
|
|
when others => c32(22 downto 0) <= s7postshift;
|
148 |
118 |
jguarin200 |
end case;
|
149 |
|
|
end if;
|
150 |
|
|
end process;
|
151 |
|
|
--! Combinatorial gremlin, Etapa 0 el corrimiento de la mantissa con menor exponente y reorganiza los operandos,\n
|
152 |
|
|
--! si el mayor es b, intercambia las posición si el mayor es a las posiciones la mantiene.
|
153 |
|
|
s0delta <= s0a(30 downto 23)-s0b(30 downto 23);
|
154 |
|
|
--! Combinatorial Gremlin, Etapa 1 Codificar el factor de corrimiento de denormalizacion y denormalizar la mantissa no fija. Signar la mantissa que se queda fija.
|
155 |
|
|
decodeshiftfactor:
|
156 |
|
|
process (s1delta(2 downto 0))
|
157 |
|
|
begin
|
158 |
|
|
case s1delta(2 downto 0) is
|
159 |
|
|
when "111" => s1shifter(8 downto 0) <= '0'&s1delta(5)&"00000"¬(s1delta(5))&'0';
|
160 |
|
|
when "110" => s1shifter(8 downto 0) <= "00"&s1delta(5)&"000"¬(s1delta(5))&"00";
|
161 |
|
|
when "101" => s1shifter(8 downto 0) <= "000"&s1delta(5)&'0'¬(s1delta(5))&"000";
|
162 |
|
|
when "100" => s1shifter(8 downto 0) <= '0'&x"10";
|
163 |
|
|
when "011" => s1shifter(8 downto 0) <= "000"¬(s1delta(5))&'0'&s1delta(5)&"000";
|
164 |
|
|
when "010" => s1shifter(8 downto 0) <= "00"¬(s1delta(5))&"000"&s1delta(5)&"00";
|
165 |
|
|
when "001" => s1shifter(8 downto 0) <= '0'¬(s1delta(5))&"00000"&s1delta(5)&'0';
|
166 |
|
|
when others => s1shifter(8 downto 0) <= not(s1delta(5))&"0000000"&s1delta(5);
|
167 |
|
|
end case;
|
168 |
|
|
end process;
|
169 |
|
|
denormhighshiftermult:lpm_mult
|
170 |
|
|
generic map ("DEDICATED_MULTIPLIER_CIRCUITRY=YES,MAXIMIZE_SPEED=9","UNSIGNED","LPM_MULT",9,18,27)
|
171 |
|
|
port map (s1shifter,s1zero&s1umantshift(22 downto 06),s1ph);
|
172 |
|
|
denormlowshiftermult:lpm_mult
|
173 |
|
|
generic map ("DEDICATED_MULTIPLIER_CIRCUITRY=YES,MAXIMIZE_SPEED=9","UNSIGNED","LPM_MULT",9,9,18)
|
174 |
|
|
port map (s1shifter,s1umantshift(5 downto 0)&"000",s1pl);
|
175 |
|
|
|
176 |
|
|
s1postshift(23 downto 7) <= s1ph(25 downto 9);
|
177 |
|
|
s1postshift(06 downto 0) <= s1ph(08 downto 2) or s1pl(17 downto 11);
|
178 |
|
|
s1xorslab(23 downto 0) <= (others => s1umantfixed(23));
|
179 |
|
|
|
180 |
|
|
--! Combinatorial Gremlin, Etapa 2: Signar la mantissa denormalizada.
|
181 |
|
|
s2xorslab <= (others => s2umantshift(24));
|
182 |
|
|
|
183 |
|
|
--! Combinatorial Gremlin, Etapa 4: Quitar el signo de la mantissa resultante.
|
184 |
|
|
s4xorslab <= (others => s4sresult(25));
|
185 |
|
|
|
186 |
|
|
--! Combinatorial Gremlin, Etapa 5: Codificar el factor de normalizacion de la mantissa resultante.
|
187 |
|
|
normalizerdecodeshift:
|
188 |
120 |
jguarin200 |
process (s5result,s5factorhot24,s5token,s5tokena,s5tokenb,s5tokenc,s5factorhot9)
|
189 |
118 |
jguarin200 |
begin
|
190 |
120 |
jguarin200 |
s5tokena <= not(s5result(24));
|
191 |
|
|
s5tokenb <= not(s5result(24));
|
192 |
|
|
s5tokenc <= not(s5result(24));
|
193 |
|
|
s5factor(7 downto 5) <= (others => s5result(24));
|
194 |
|
|
s5factorhot24 <= x"000000";
|
195 |
|
|
for i in 23 downto 16 loop
|
196 |
118 |
jguarin200 |
if s5result(i)='1' then
|
197 |
120 |
jguarin200 |
s5factorhot24(23-i) <= s5tokena;
|
198 |
|
|
s5tokenb <= '0';
|
199 |
|
|
s5tokenc <= '0';
|
200 |
118 |
jguarin200 |
exit;
|
201 |
|
|
end if;
|
202 |
|
|
end loop;
|
203 |
120 |
jguarin200 |
for i in 15 downto 8 loop
|
204 |
|
|
if s5result(i)='1' then
|
205 |
|
|
s5factorhot24(23-i) <= s5tokenb;
|
206 |
|
|
s5tokenc <= '0';
|
207 |
|
|
exit;
|
208 |
|
|
end if;
|
209 |
|
|
end loop;
|
210 |
|
|
for i in 7 downto 0 loop
|
211 |
|
|
if s5result(i)='1' then
|
212 |
|
|
s5factorhot24(23-i) <= s5tokenc;
|
213 |
|
|
exit;
|
214 |
|
|
end if;
|
215 |
|
|
end loop;
|
216 |
|
|
s5token <=s5tokena&s5tokenb&s5tokenc;
|
217 |
|
|
case (s5token) is
|
218 |
|
|
when "100" => s5factor(4 downto 3) <= "10";
|
219 |
|
|
when "110" => s5factor(4 downto 3) <= "01";
|
220 |
|
|
when "111" => s5factor(4 downto 3) <= "00";
|
221 |
|
|
when others => s5factor(4 downto 3) <= (others => s5result(24));
|
222 |
|
|
end case;
|
223 |
|
|
s5factorhot9 <= (s5factorhot24(7 downto 0)or s5factorhot24(15 downto 8)or s5factorhot24(23 downto 16)) & s5result(24);
|
224 |
|
|
case s5factorhot9 is
|
225 |
|
|
when "100000000" => s5factor(2 downto 0) <= "111";
|
226 |
|
|
when "010000000" => s5factor(2 downto 0) <= "110";
|
227 |
|
|
when "001000000" => s5factor(2 downto 0) <= "101";
|
228 |
|
|
when "000100000" => s5factor(2 downto 0) <= "100";
|
229 |
|
|
when "000010000" => s5factor(2 downto 0) <= "011";
|
230 |
|
|
when "000001000" => s5factor(2 downto 0) <= "010";
|
231 |
|
|
when "000000100" => s5factor(2 downto 0) <= "001";
|
232 |
|
|
when "000000010" => s5factor(2 downto 0) <= "000";
|
233 |
|
|
when others => s5factor (2 downto 0) <= (others => s5result(24));
|
234 |
|
|
end case;
|
235 |
|
|
|
236 |
118 |
jguarin200 |
end process;
|
237 |
119 |
jguarin200 |
|
238 |
|
|
--! Etapa 6: Ejecutar el corrimiento para normalizar la mantissa.
|
239 |
118 |
jguarin200 |
normhighshiftermult:lpm_mult
|
240 |
|
|
generic map ("DEDICATED_MULTIPLIER_CIRCUITRY=YES,MAXIMIZE_SPEED=9","UNSIGNED","LPM_MULT",9,18,27)
|
241 |
119 |
jguarin200 |
port map (s6factorhot9,s6result(24 downto 7),s6ph);
|
242 |
118 |
jguarin200 |
normlowshiftermult:lpm_mult
|
243 |
|
|
generic map ("DEDICATED_MULTIPLIER_CIRCUITRY=YES,MAXIMIZE_SPEED=9","UNSIGNED","LPM_MULT",9,9,18)
|
244 |
119 |
jguarin200 |
port map (s6factorhot9,s6result(06 downto 0)&"00",s6pl);
|
245 |
|
|
s6postshift(22 downto 15) <= s6ph(16 downto 09);
|
246 |
120 |
jguarin200 |
s6postshift(14 downto 06) <= s6ph(08 downto 00) + s6pl(17 downto 09);
|
247 |
119 |
jguarin200 |
s6postshift(05 downto 00) <= s6pl(08 downto 03);
|
248 |
118 |
jguarin200 |
|
249 |
|
|
|
250 |
|
|
|
251 |
|
|
|
252 |
|
|
|
253 |
119 |
jguarin200 |
end fadd32_arch;
|
254 |
118 |
jguarin200 |
|
255 |
|
|
|