thfu
/
TailCorr
1
0
Fork 0
Code Issues Pull Requests Packages Projects Releases Wiki Activity

Compare commits

base: thfu:ddfb69f31b84b9966ec484752ae5c92feb1feb5b
Branches Tags
thfu:main
thfu:RealOnly
thfu:EZQ2p5
thfu:b2
..
compare: thfu:e297bf11e3f7f03845caeae39281f2d0e4918af3
Branches Tags
thfu:main
thfu:RealOnly
thfu:EZQ2p5
thfu:b2

No commits in common. "ddfb69f31b84b9966ec484752ae5c92feb1feb5b" and "e297bf11e3f7f03845caeae39281f2d0e4918af3" have entirely different histories.
ddfb69f31b ... e297bf11e3

11 changed files with 914 additions and 1419 deletions
Show Stats Expand all files Collapse all files

38
rtl/z_dsp/FixRound.v Normal file
View File

@ -0,0 +1,38 @@
module FixRound #(
parameter integer Data_width = 8
,parameter integer Fix_frac_coef_width = 31//division
)
(
input clk
,input rstn
,input en
,input signed [Data_width-1:0] din
,output signed [Data_width-1:0] dout
);
reg signed [Data_width-1:0] din_round;
always@(posedge clk or negedge rstn)
if(!rstn)
begin
din_round <= 'h0;
end
else if(en) begin
if(din[Data_width-1] == 1'b0)
begin
din_round <= din + {{1'b1},{(Fix_frac_coef_width-1){1'b0}}};
end
else if (din[Data_width-1] == 1'b1)
begin
din_round <= din + {{1'b1},{(Fix_frac_coef_width-1){1'b0}}} - 1'b1;
end
end
else begin
din_round <= din_round;
end
assign dout = din_round;
endmodule

132
rtl/z_dsp/IIR_Filter_p1.v
View File

@ -1,132 +0,0 @@
//+FHDR--------------------------------------------------------------------------------------------------------
// Company:
//-----------------------------------------------------------------------------------------------------------------
// File Name : IIR_Filter_p1.v
// Department :
// Author : hdzhang
// Author's Tel :
//-----------------------------------------------------------------------------------------------------------------
// Relese History
// Version Date Author Description
// 0.0 2025-03-09 hdzhang
//2024-05-28 10:22:49
//-----------------------------------------------------------------------------------------------------------------
// Keywords :
//
//-----------------------------------------------------------------------------------------------------------------
// Parameter
//
//-----------------------------------------------------------------------------------------------------------------
// Purpose :
//
//-----------------------------------------------------------------------------------------------------------------
// Target Device:
// Tool versions:
//-----------------------------------------------------------------------------------------------------------------
// Reuse Issues
// Reset Strategy:
// Clock Domains:
// Critical Timing:
// Asynchronous I/F:
// Synthesizable (y/n):
// Other:
//-FHDR--------------------------------------------------------------------------------------------------------
module IIR_Filter_p1 #(
parameter coef_width = 32
,parameter data_in_width = 16
,parameter cascade_in_width = 37
,parameter temp_var_width = cascade_in_width - 1
,parameter data_out_width = cascade_in_width - 2
)
//H(z) = a / (1 - b*z^-1)
(
input rstn
,input clk
,input en
,input signed [data_in_width-1 :0] din_re // Re(x(t))
,input signed [cascade_in_width-1:0] dout_r1_re // Re(y(t-1))
,input signed [cascade_in_width-1:0] dout_r1_im // Im(y(t-1))
,input signed [coef_width-1 :0] a_re
,input signed [coef_width-1 :0] a_im
,input signed [coef_width-1 :0] b_re
,input signed [coef_width-1 :0] b_im
,output signed [data_out_width-1:0] dout_re // Re(y(t-16))
,output signed [data_out_width-1:0] dout_im // Im(y(t-16))
);
wire signed [temp_var_width-1 :0] x1_re;
wire signed [temp_var_width-1 :0] x1_im;
wire signed [temp_var_width-1 :0] y1_re;
wire signed [temp_var_width-1 :0] y1_im;
wire signed [temp_var_width :0] y_re;
wire signed [temp_var_width :0] y_im;
wire signed [data_out_width-1:0] y_re_trunc;
wire signed [data_out_width-1:0] y_im_trunc;
// x1 = a * din delay M = a*x(t-8)
mult_x
#(
.A_width (data_in_width )
,.C_width (coef_width )
,.D_width (coef_width )
,.o_width (temp_var_width )
)
inst_c1 (
.clk (clk ),
.rstn (rstn ),
.en (en ),
.a (din_re ),
.c (a_re ),
.d (a_im ),
.Re (x1_re ),
.Im (x1_im )
);
// y1 = b * dout_r1 delay M = b*y(t-9)
// y = y1+x1 = a*x(t-8)+b*y(t-9) = y(t-8)
mult_C
#(
.A_width (cascade_in_width )
,.B_width (cascade_in_width )
,.C_width (coef_width )
,.D_width (coef_width )
,.o_width (temp_var_width )
)
inst_c3 (
.clk (clk ),
.rstn (rstn ),
.en (en ),
.a (dout_r1_re ),
.b (dout_r1_im ),
.c (b_re ),
.d (b_im ),
.Re (y1_re ),
.Im (y1_im )
);
assign y_re = x1_re + y1_re;
assign y_im = x1_im + y1_im;
// dout = round(y) delay M = round(y(t-16))
trunc #(
.diw (temp_var_width+1 )
,.msb (temp_var_width-1 )
,.lsb (temp_var_width-data_out_width )
) round_u1 (clk, rstn, en, y_re, y_re_trunc);
trunc #(
.diw (temp_var_width+1 )
,.msb (temp_var_width-1 )
,.lsb (temp_var_width-data_out_width )
) round_u2 (clk, rstn, en, y_im, y_im_trunc);
assign dout_re = y_re_trunc;
assign dout_im = y_im_trunc;
endmodule

215
rtl/z_dsp/IIR_Filter_p8.v
View File

@ -1,44 +1,42 @@
module IIR_Filter_p8 #( module IIR_Filter_p8 #(
parameter coef_width = 32 parameter data_in_width = 16
,parameter data_in_width = 16 ,parameter coef_width = 32
,parameter data_out_width = 37 ,parameter frac_data_out_width = 20//X for in,5
,parameter temp_var_width = data_out_width+5 ,parameter frac_coef_width = 31//division
) )
// H(z) = a(1 + b*z^-1 + b^2*z^-2 + b^3*z^-3 + b^4*z^-4 + b^5*z^-5 + b^6*z^-6 + b^7*z^-7) / (1 - b^8*z^-8)
( (
input rstn input rstn
,input clk ,input clk
,input en ,input en
,input signed [data_in_width-1 :0] dinp0 //x(8n+16) ,input signed [data_in_width-1:0] dinp0
,input signed [data_in_width-1 :0] dinp1 //x(8n+15) ,input signed [data_in_width-1:0] dinp1
,input signed [data_in_width-1 :0] dinp2 //x(8n+14) ,input signed [data_in_width-1:0] dinp2
,input signed [data_in_width-1 :0] dinp3 //x(8n+13) ,input signed [data_in_width-1:0] dinp3
,input signed [data_in_width-1 :0] dinp4 //x(8n+12) ,input signed [data_in_width-1:0] dinp4
,input signed [data_in_width-1 :0] dinp5 //x(8n+11) ,input signed [data_in_width-1:0] dinp5
,input signed [data_in_width-1 :0] dinp6 //x(8n+10) ,input signed [data_in_width-1:0] dinp6
,input signed [data_in_width-1 :0] dinp7 //x(8n+9) ,input signed [data_in_width-1:0] dinp7
,input signed [coef_width-1 :0] a_re ,input signed [coef_width-1 :0] a_re
,input signed [coef_width-1 :0] a_im ,input signed [coef_width-1 :0] a_im
,input signed [coef_width-1 :0] ab_re ,input signed [coef_width-1 :0] ab_re
,input signed [coef_width-1 :0] ab_im ,input signed [coef_width-1 :0] ab_im
,input signed [coef_width-1 :0] abb_re ,input signed [coef_width-1 :0] abb_re
,input signed [coef_width-1 :0] abb_im ,input signed [coef_width-1 :0] abb_im
,input signed [coef_width-1 :0] ab_pow3_re ,input signed [coef_width-1 :0] ab_pow3_re
,input signed [coef_width-1 :0] ab_pow3_im ,input signed [coef_width-1 :0] ab_pow3_im
,input signed [coef_width-1 :0] ab_pow4_re ,input signed [coef_width-1 :0] ab_pow4_re
,input signed [coef_width-1 :0] ab_pow4_im ,input signed [coef_width-1 :0] ab_pow4_im
,input signed [coef_width-1 :0] ab_pow5_re ,input signed [coef_width-1 :0] ab_pow5_re
,input signed [coef_width-1 :0] ab_pow5_im ,input signed [coef_width-1 :0] ab_pow5_im
,input signed [coef_width-1 :0] ab_pow6_re ,input signed [coef_width-1 :0] ab_pow6_re
,input signed [coef_width-1 :0] ab_pow6_im ,input signed [coef_width-1 :0] ab_pow6_im
,input signed [coef_width-1 :0] ab_pow7_re ,input signed [coef_width-1 :0] ab_pow7_re
,input signed [coef_width-1 :0] ab_pow7_im ,input signed [coef_width-1 :0] ab_pow7_im
,input signed [coef_width-1 :0] b_pow8_re ,input signed [coef_width-1 :0] b_pow8_re
,input signed [coef_width-1 :0] b_pow8_im ,input signed [coef_width-1 :0] b_pow8_im
,output signed [data_out_width-1:0] dout_re // Re(y(8n-8)) ,output signed [data_in_width-1:0] dout
,output signed [data_out_width-1:0] dout_im // Im(y(8n-8))
); );
wire signed [data_in_width-1 :0] dinp [7:0]; wire signed [data_in_width-1 :0] dinp [7:0];
@ -71,57 +69,39 @@ assign ab_pow_im[2] = abb_im;
assign ab_pow_im[1] = ab_im; assign ab_pow_im[1] = ab_im;
assign ab_pow_im[0] = a_im; assign ab_pow_im[0] = a_im;
wire signed [data_in_width+frac_data_out_width-1:0] x_re [0:7];
wire signed [data_in_width+frac_data_out_width-1:0] x_im [0:7];
wire signed [temp_var_width-1 :0] x_re [0:7];
wire signed [temp_var_width-1 :0] x_im [0:7];
wire signed [temp_var_width+3 :0] v_re;
wire signed [temp_var_width+3 :0] v_im;
reg signed [temp_var_width+3 :0] v1_re;
reg signed [temp_var_width+3 :0] v1_im;
wire signed [temp_var_width+3 :0] y_re;
wire signed [temp_var_width+3 :0] y_im;
wire signed [temp_var_width+3 :0] y1_re;
wire signed [temp_var_width+3 :0] y1_im;
wire signed [data_out_width-1:0] y_re_trunc;
wire signed [data_out_width-1:0] y_im_trunc;
// x[0] = (dinp0 * a_re) delay M = a*x(8n+8)
// x[1] = (dinp1 * ab_re) delay M = a*b*x(8n+7)
// x[2] = (dinp2 * abb_re) delay M = a*b^2*x(8n+6)
// x[3] = (dinp3 * ab_pow3_re) delay M = a*b^3*x(8n+5)
// x[4] = (dinp4 * ab_pow4_re) delay M = a*b^4*x(8n+4)
// x[5] = (dinp5 * ab_pow5_re) delay M = a*b^5*x(8n+3)
// x[6] = (dinp6 * ab_pow6_re) delay M = a*b^6*x(8n+2)
// x[7] = (dinp7 * ab_pow7_re) delay M = a*b^7*x(8n+1)
genvar i; genvar i;
generate generate
for (i = 0; i < 8; i = i + 1) begin: mult_c_inst for (i = 0; i < 8; i = i + 1) begin: mult_x_inst
mult_x #( mult_x #(
.A_width (data_in_width ), .A_width(data_in_width),
.C_width (coef_width ), .C_width(coef_width+frac_data_out_width),
.D_width (coef_width ), .D_width(coef_width+frac_data_out_width),
.o_width (temp_var_width ) .frac_coef_width(frac_coef_width)
) inst_c ( ) inst_mult_x (
.clk (clk ), .clk (clk),
.rstn (rstn ), .rstn (rstn),
.en (en ), .en (en),
.a (dinp[i] ), .a (dinp[i]),
.c (ab_pow_re[i] ), .c ({ab_pow_re[i],{frac_data_out_width{1'b0}}}),
.d (ab_pow_im[i] ), .d ({ab_pow_im[i],{frac_data_out_width{1'b0}}}),
.Re (x_re[i] ), .Re (x_re[i]),
.Im (x_im[i] ) .Im (x_im[i])
); );
end end
endgenerate endgenerate
wire signed [data_in_width+frac_data_out_width+3:0] v_re;
wire signed [data_in_width+frac_data_out_width+3:0] v_im;
// v1 = sum_{i=0,1,...,7}{x[i]} delay M = sum_{i=0,1,...,7}{a*b^i*x(8n-i)}
assign v_re = x_re[0] + x_re[1] +x_re[2] +x_re[3] +x_re[4] +x_re[5] +x_re[6] +x_re[7]; assign v_re = x_re[0] + x_re[1] +x_re[2] +x_re[3] +x_re[4] +x_re[5] +x_re[6] +x_re[7];
assign v_im = x_im[0] + x_im[1] +x_im[2] +x_im[3] +x_im[4] +x_im[5] +x_im[6] +x_im[7]; assign v_im = x_im[0] + x_im[1] +x_im[2] +x_im[3] +x_im[4] +x_im[5] +x_im[6] +x_im[7];
reg signed [data_in_width+frac_data_out_width+3:0] v1_re;
reg signed [data_in_width+frac_data_out_width+3:0] v1_im;
always @(posedge clk or negedge rstn) always @(posedge clk or negedge rstn)
if (!rstn) if (!rstn)
begin begin
@ -139,47 +119,76 @@ always @(posedge clk or negedge rstn)
v1_im <= v1_im; v1_im <= v1_im;
end end
wire signed [data_in_width+frac_data_out_width+3:0] y_re;
wire signed [data_in_width+frac_data_out_width+3:0] y_im;
wire signed [data_in_width+frac_data_out_width+3:0] y1_re;
wire signed [data_in_width+frac_data_out_width+3:0] y1_im;
reg signed [data_in_width-1:0] dout_re;
// y1 = (b^8 * y) delay M = b^8*y(8n-8)
// y = v1 + y1 = sum_{i=0,1,...,7}{a*b^i*x(8n-i)} + b^8*y(8n-8) = y(8n)
mult_C mult_C
#( #(
.A_width (temp_var_width+4 ) .A_width(data_in_width+frac_data_out_width+4)
,.B_width (temp_var_width+4 ) ,.B_width(data_in_width+frac_data_out_width+4)
,.C_width (coef_width ) ,.C_width(coef_width)
,.D_width (coef_width ) ,.D_width(coef_width)
,.o_width (temp_var_width+4 ) ,.frac_coef_width(frac_coef_width)
) )
inst_c9 ( inst_c9 (
.clk (clk ), .clk (clk ),
.rstn (rstn ), .rstn (rstn ),
.en (en ), .en (en ),
.a (y_re ), .a (y_re ),
.b (y_im ), .b (y_im ),
.c (b_pow8_re ), .c (b_pow8_re ),
.d (b_pow8_im ), .d (b_pow8_im ),
.Re (y1_re ), .Re (y1_re ),//b^8*y(n-1)
.Im (y1_im ) .Im (y1_im )
); );
assign y_re = v1_re + y1_re; assign y_re = v1_re + y1_re;
assign y_im = v1_im + y1_im; assign y_im = v1_im + y1_im;
// dout = round(y) delay M = round(y(8n-8)) wire signed [data_in_width+frac_data_out_width+3:0] dout_round;
trunc #(
.diw (temp_var_width+4 )
,.msb (temp_var_width-1 )
,.lsb (temp_var_width-data_out_width )
) round_u1 (clk, rstn, en, y_re, y_re_trunc);
trunc #(
.diw (temp_var_width+4 )
,.msb (temp_var_width-1 )
,.lsb (temp_var_width-data_out_width )
) round_u2 (clk, rstn, en, y_im, y_im_trunc);
FixRound #(data_in_width+frac_data_out_width+4,frac_data_out_width) u_round1 (clk, rstn, en, y_re, dout_round);
assign dout_re = y_re_trunc; always @(posedge clk or negedge rstn)
assign dout_im = y_im_trunc; if (!rstn)
begin
dout_re <= 'h0;
end
else if(en)
begin
dout_re <= dout_round[frac_data_out_width+15:frac_data_out_width];
end
else
begin
dout_re <= dout_re;
end
reg signed [data_in_width-1:0] dout_clip;
always @(posedge clk or negedge rstn)
if (!rstn)
begin
dout_clip <= 'h0;
end
else if(en)
begin
if(dout_round[frac_data_out_width+16:frac_data_out_width+15]==2'b01)
dout_clip <= 16'd32767;
else if(dout_round[frac_data_out_width+16:frac_data_out_width+15]==2'b10)
dout_clip <= -16'd32768;
else
dout_clip <= dout_re;
end
else
begin
dout_clip <= dout_clip;
end
assign dout = dout_clip;
endmodule endmodule

532
rtl/z_dsp/IIR_top.v
View File

@ -1,234 +1,322 @@
module IIR_top #( module IIR_top
parameter data_out_width = 23
,parameter temp_var_width = data_out_width + 14
)
( (
input rstn input rstn
,input clk ,input clk
,input en ,input en
,input signed [15 :0] IIRin_p0 // x(8n+9) ,input signed [15 :0] IIRin_p0
,input signed [15 :0] IIRin_p1 // x(8n+10) ,input signed [15 :0] IIRin_p1
,input signed [15 :0] IIRin_p2 // x(8n+11) ,input signed [15 :0] IIRin_p2
,input signed [15 :0] IIRin_p3 // x(8n+12) ,input signed [15 :0] IIRin_p3
,input signed [15 :0] IIRin_p4 // x(8n+13) ,input signed [15 :0] IIRin_p4
,input signed [15 :0] IIRin_p5 // x(8n+14) ,input signed [15 :0] IIRin_p5
,input signed [15 :0] IIRin_p6 // x(8n+15) ,input signed [15 :0] IIRin_p6
,input signed [15 :0] IIRin_p7 // x(8n+16) ,input signed [15 :0] IIRin_p7
,input signed [15 :0] IIRin_p0_r2 // x(8n+9) delay 2M -> x(8n- 7) ,input signed [31 :0] a_re
,input signed [15 :0] IIRin_p1_r4 // x(8n+10) delay 4M -> x(8n-22) ,input signed [31 :0] a_im
,input signed [15 :0] IIRin_p2_r6 // x(8n+11) delay 6M -> x(8n-37) ,input signed [31 :0] ab_re
,input signed [15 :0] IIRin_p3_r8 // x(8n+12) delay 8M -> x(8n-52) ,input signed [31 :0] ab_im
,input signed [15 :0] IIRin_p4_r10 // x(8n+13) delay 10M -> x(8n-67) ,input signed [31 :0] abb_re
,input signed [15 :0] IIRin_p5_r12 // x(8n+14) delay 12M -> x(8n-82) ,input signed [31 :0] abb_im
,input signed [15 :0] IIRin_p6_r14 // x(8n+15) delay 14M -> x(8n-97) ,input signed [31 :0] ab_pow3_re
,input signed [31 :0] a_re ,input signed [31 :0] ab_pow3_im
,input signed [31 :0] a_im ,input signed [31 :0] ab_pow4_re
,input signed [31 :0] b_re ,input signed [31 :0] ab_pow4_im
,input signed [31 :0] b_im ,input signed [31 :0] ab_pow5_re
,input signed [31 :0] ab_re ,input signed [31 :0] ab_pow5_im
,input signed [31 :0] ab_im ,input signed [31 :0] ab_pow6_re
,input signed [31 :0] abb_re ,input signed [31 :0] ab_pow6_im
,input signed [31 :0] abb_im ,input signed [31 :0] ab_pow7_re
,input signed [31 :0] ab_pow3_re ,input signed [31 :0] ab_pow7_im
,input signed [31 :0] ab_pow3_im ,input signed [31 :0] b_pow8_re
,input signed [31 :0] ab_pow4_re ,input signed [31 :0] b_pow8_im
,input signed [31 :0] ab_pow4_im
,input signed [31 :0] ab_pow5_re
,input signed [31 :0] ab_pow5_im
,input signed [31 :0] ab_pow6_re
,input signed [31 :0] ab_pow6_im
,input signed [31 :0] ab_pow7_re
,input signed [31 :0] ab_pow7_im
,input signed [31 :0] b_pow8_re
,input signed [31 :0] b_pow8_im
,output signed [data_out_width-1 :0] IIRout_p0 // y(8n-8) ,output signed [15 :0] IIRout_p0
,output signed [data_out_width-1 :0] IIRout_p1 // y(8n-23) ,output signed [15 :0] IIRout_p1
,output signed [data_out_width-1 :0] IIRout_p2 // y(8n-38) ,output signed [15 :0] IIRout_p2
,output signed [data_out_width-1 :0] IIRout_p3 // y(8n-53) ,output signed [15 :0] IIRout_p3
,output signed [data_out_width-1 :0] IIRout_p4 // y(8n-68) ,output signed [15 :0] IIRout_p4
,output signed [data_out_width-1 :0] IIRout_p5 // y(8n-83) ,output signed [15 :0] IIRout_p5
,output signed [data_out_width-1 :0] IIRout_p6 // y(8n-98) ,output signed [15 :0] IIRout_p6
,output signed [data_out_width-1 :0] IIRout_p7 // y(8n-113) ,output signed [15 :0] IIRout_p7
); );
wire signed [15:0] IIRin_p_r1 [7:1];
wire signed [15 : 0] IIRin_p [7:0];
assign IIRin_p[7] = IIRin_p7;
assign IIRin_p[6] = IIRin_p6;
assign IIRin_p[5] = IIRin_p5;
assign IIRin_p[4] = IIRin_p4;
assign IIRin_p[3] = IIRin_p3;
assign IIRin_p[2] = IIRin_p2;
assign IIRin_p[1] = IIRin_p1;
assign IIRin_p[0] = IIRin_p0;
wire signed [temp_var_width- 1:0] IIRout_p0_re; sirv_gnrl_dfflr #(16) dff_IIRin_p7_1(en,IIRin_p[7], IIRin_p_r1[7] ,clk,rstn);
wire signed [temp_var_width- 3:0] IIRout_p1_re; sirv_gnrl_dfflr #(16) dff_IIRin_p6_1(en,IIRin_p[6], IIRin_p_r1[6] ,clk,rstn);
wire signed [temp_var_width- 5:0] IIRout_p2_re; sirv_gnrl_dfflr #(16) dff_IIRin_p5_1(en,IIRin_p[5], IIRin_p_r1[5] ,clk,rstn);
wire signed [temp_var_width- 7:0] IIRout_p3_re; sirv_gnrl_dfflr #(16) dff_IIRin_p4_1(en,IIRin_p[4], IIRin_p_r1[4] ,clk,rstn);
wire signed [temp_var_width- 9:0] IIRout_p4_re; sirv_gnrl_dfflr #(16) dff_IIRin_p3_1(en,IIRin_p[3], IIRin_p_r1[3] ,clk,rstn);
wire signed [temp_var_width-11:0] IIRout_p5_re; sirv_gnrl_dfflr #(16) dff_IIRin_p2_1(en,IIRin_p[2], IIRin_p_r1[2] ,clk,rstn);
wire signed [temp_var_width-13:0] IIRout_p6_re; sirv_gnrl_dfflr #(16) dff_IIRin_p1_1(en,IIRin_p[1], IIRin_p_r1[1] ,clk,rstn);
wire signed [temp_var_width-15:0] IIRout_p7_re;
wire signed [temp_var_width- 1:0] IIRout_p0_im;
wire signed [temp_var_width- 3:0] IIRout_p1_im;
wire signed [temp_var_width- 5:0] IIRout_p2_im;
wire signed [temp_var_width- 7:0] IIRout_p3_im;
wire signed [temp_var_width- 9:0] IIRout_p4_im;
wire signed [temp_var_width-11:0] IIRout_p5_im;
wire signed [temp_var_width-13:0] IIRout_p6_im;
wire signed [temp_var_width-15:0] IIRout_p7_im;
IIR_Filter_p8 inst_iir_p0 (
.clk (clk ),
.rstn (rstn ),
.en (en ),
.dinp0 (IIRin_p[0] ),
.dinp1 (IIRin_p_r1[7] ),
.dinp2 (IIRin_p_r1[6] ),
.dinp3 (IIRin_p_r1[5] ),
.dinp4 (IIRin_p_r1[4] ),
.dinp5 (IIRin_p_r1[3] ),
.dinp6 (IIRin_p_r1[2] ),
.dinp7 (IIRin_p_r1[1] ),
.a_re (a_re ),
.a_im (a_im ),
.ab_re (ab_re ),
.ab_im (ab_im ),
.abb_re (abb_re ),
.abb_im (abb_im ),
.ab_pow3_re (ab_pow3_re ),
.ab_pow3_im (ab_pow3_im ),
.ab_pow4_re (ab_pow4_re ),
.ab_pow4_im (ab_pow4_im ),
.ab_pow5_re (ab_pow5_re ),
.ab_pow5_im (ab_pow5_im ),
.ab_pow6_re (ab_pow6_re ),
.ab_pow6_im (ab_pow6_im ),
.ab_pow7_re (ab_pow7_re ),
.ab_pow7_im (ab_pow7_im ),
.b_pow8_re (b_pow8_re ),
.b_pow8_im (b_pow8_im ),
.dout (IIRout_p0 )
);
IIR_Filter_p8 inst_iir_p1 (
IIR_Filter_p8 #( .clk (clk ),
.data_out_width (temp_var_width ) .rstn (rstn ),
) inst_iir_p0 ( .en (en ),
.clk (clk ), .dinp0 (IIRin_p[1] ),
.rstn (rstn ), .dinp1 (IIRin_p[0] ),
.en (en ), .dinp2 (IIRin_p_r1[7] ),
.dinp0 (IIRin_p7 ), // x(8n+16) .dinp3 (IIRin_p_r1[6] ),
.dinp1 (IIRin_p6 ), // x(8n+15) .dinp4 (IIRin_p_r1[5] ),
.dinp2 (IIRin_p5 ), // x(8n+14) .dinp5 (IIRin_p_r1[4] ),
.dinp3 (IIRin_p4 ), // x(8n+13) .dinp6 (IIRin_p_r1[3] ),
.dinp4 (IIRin_p3 ), // x(8n+12) .dinp7 (IIRin_p_r1[2] ),
.dinp5 (IIRin_p2 ), // x(8n+11) .a_re (a_re ),
.dinp6 (IIRin_p1 ), // x(8n+10) .a_im (a_im ),
.dinp7 (IIRin_p0 ), // x(8n+9) .ab_re (ab_re ),
.a_re (a_re ), .ab_im (ab_im ),
.a_im (a_im ), .abb_re (abb_re ),
.ab_re (ab_re ), .abb_im (abb_im ),
.ab_im (ab_im ), .ab_pow3_re (ab_pow3_re ),
.abb_re (abb_re ), .ab_pow3_im (ab_pow3_im ),
.abb_im (abb_im ), .ab_pow4_re (ab_pow4_re ),
.ab_pow3_re (ab_pow3_re ), .ab_pow4_im (ab_pow4_im ),
.ab_pow3_im (ab_pow3_im ), .ab_pow5_re (ab_pow5_re ),
.ab_pow4_re (ab_pow4_re ), .ab_pow5_im (ab_pow5_im ),
.ab_pow4_im (ab_pow4_im ), .ab_pow6_re (ab_pow6_re ),
.ab_pow5_re (ab_pow5_re ), .ab_pow6_im (ab_pow6_im ),
.ab_pow5_im (ab_pow5_im ), .ab_pow7_re (ab_pow7_re ),
.ab_pow6_re (ab_pow6_re ), .ab_pow7_im (ab_pow7_im ),
.ab_pow6_im (ab_pow6_im ), .b_pow8_re (b_pow8_re ),
.ab_pow7_re (ab_pow7_re ), .b_pow8_im (b_pow8_im ),
.ab_pow7_im (ab_pow7_im ), .dout (IIRout_p1 )
.b_pow8_re (b_pow8_re ), );
.b_pow8_im (b_pow8_im ), IIR_Filter_p8 inst_iir_p2 (
.dout_re (IIRout_p0_re ), // Re(y(8n-8)) .clk (clk ),
.dout_im (IIRout_p0_im ) // Im(y(8n-8)) .rstn (rstn ),
); .en (en ),
.dinp0 (IIRin_p[2] ),
IIR_Filter_p1 #( .dinp1 (IIRin_p[1] ),
.cascade_in_width (temp_var_width ) .dinp2 (IIRin_p[0] ),
) inst_iir_p1( .dinp3 (IIRin_p_r1[7] ),
.clk (clk ), .dinp4 (IIRin_p_r1[6] ),
.rstn (rstn ), .dinp5 (IIRin_p_r1[5] ),
.en (en ), .dinp6 (IIRin_p_r1[4] ),
.din_re (IIRin_p0_r2 ), // x(8n-7) .dinp7 (IIRin_p_r1[3] ),
.dout_r1_re (IIRout_p0_re ), // Re(y(8n-8)) .a_re (a_re ),
.dout_r1_im (IIRout_p0_im ), // Im(y(8n-8)) .a_im (a_im ),
.a_re (a_re ), .ab_re (ab_re ),
.a_im (a_im ), .ab_im (ab_im ),
.b_re (b_re ), .abb_re (abb_re ),
.b_im (b_im ), .abb_im (abb_im ),
.dout_re (IIRout_p1_re ), // Re(y(8n-23)) .ab_pow3_re (ab_pow3_re ),
.dout_im (IIRout_p1_im ) // Im(y(8n-23)) .ab_pow3_im (ab_pow3_im ),
); .ab_pow4_re (ab_pow4_re ),
IIR_Filter_p1 #( .ab_pow4_im (ab_pow4_im ),
.cascade_in_width (temp_var_width-2 ) .ab_pow5_re (ab_pow5_re ),
) inst_iir_p2 ( .ab_pow5_im (ab_pow5_im ),
.clk (clk ), .ab_pow6_re (ab_pow6_re ),
.rstn (rstn ), .ab_pow6_im (ab_pow6_im ),
.en (en ), .ab_pow7_re (ab_pow7_re ),
.din_re (IIRin_p1_r4 ), // x(8n-22) .ab_pow7_im (ab_pow7_im ),
.dout_r1_re (IIRout_p1_re ), // Re(y(8n-23)) .b_pow8_re (b_pow8_re ),
.dout_r1_im (IIRout_p1_im ), // Im(y(8n-23)) .b_pow8_im (b_pow8_im ),
.a_re (a_re ), .dout (IIRout_p2 )
.a_im (a_im ), );
.b_re (b_re ), IIR_Filter_p8 inst_iir_p3 (
.b_im (b_im ), .clk (clk ),
.dout_re (IIRout_p2_re ), // Re(y(8n-38)) .rstn (rstn ),
.dout_im (IIRout_p2_im ) // Im(y(8n-38)) .en (en ),
); .dinp0 (IIRin_p[3] ),
IIR_Filter_p1 #( .dinp1 (IIRin_p[2] ),
.cascade_in_width (temp_var_width-4 ) .dinp2 (IIRin_p[1] ),
) inst_iir_p3 ( .dinp3 (IIRin_p[0] ),
.clk (clk ), .dinp4 (IIRin_p_r1[7] ),
.rstn (rstn ), .dinp5 (IIRin_p_r1[6] ),
.en (en ), .dinp6 (IIRin_p_r1[5] ),
.din_re (IIRin_p2_r6 ), // x(8n-37) .dinp7 (IIRin_p_r1[4] ),
.dout_r1_re (IIRout_p2_re ), // Re(y(8n-38)) .a_re (a_re ),
.dout_r1_im (IIRout_p2_im ), // Im(y(8n-38)) .a_im (a_im ),
.a_re (a_re ), .ab_re (ab_re ),
.a_im (a_im ), .ab_im (ab_im ),
.b_re (b_re ), .abb_re (abb_re ),
.b_im (b_im ), .abb_im (abb_im ),
.dout_re (IIRout_p3_re ), // Re(y(8n-53)) .ab_pow3_re (ab_pow3_re ),
.dout_im (IIRout_p3_im ) // Im(y(8n-53)) .ab_pow3_im (ab_pow3_im ),
); .ab_pow4_re (ab_pow4_re ),
IIR_Filter_p1 #( .ab_pow4_im (ab_pow4_im ),
.cascade_in_width (temp_var_width-6 ) .ab_pow5_re (ab_pow5_re ),
) inst_iir_p4 ( .ab_pow5_im (ab_pow5_im ),
.clk (clk ), .ab_pow6_re (ab_pow6_re ),
.rstn (rstn ), .ab_pow6_im (ab_pow6_im ),
.en (en ), .ab_pow7_re (ab_pow7_re ),
.din_re (IIRin_p3_r8 ), // x(8n-52) .ab_pow7_im (ab_pow7_im ),
.dout_r1_re (IIRout_p3_re ), // Re(y(8n-53)) .b_pow8_re (b_pow8_re ),
.dout_r1_im (IIRout_p3_im ), // Im(y(8n-53)) .b_pow8_im (b_pow8_im ),
.a_re (a_re ), .dout (IIRout_p3 )
.a_im (a_im ), );
.b_re (b_re ), IIR_Filter_p8 inst_iir_p4 (
.b_im (b_im ), .clk (clk ),
.dout_re (IIRout_p4_re ), // Re(y(8n-68)) .rstn (rstn ),
.dout_im (IIRout_p4_im ) // Im(y(8n-68)) .en (en ),
); .dinp0 (IIRin_p[4] ),
IIR_Filter_p1 #( .dinp1 (IIRin_p[3] ),
.cascade_in_width (temp_var_width-8 ) .dinp2 (IIRin_p[2] ),
) inst_iir_p5 ( .dinp3 (IIRin_p[1] ),
.clk (clk ), .dinp4 (IIRin_p[0] ),
.rstn (rstn ), .dinp5 (IIRin_p_r1[7] ),
.en (en ), .dinp6 (IIRin_p_r1[6] ),
.din_re (IIRin_p4_r10 ), // x(8n-67) .dinp7 (IIRin_p_r1[5] ),
.dout_r1_re (IIRout_p4_re ), // Re(y(8n-68)) .a_re (a_re ),
.dout_r1_im (IIRout_p4_im ), // Im(y(8n-68)) .a_im (a_im ),
.a_re (a_re ), .ab_re (ab_re ),
.a_im (a_im ), .ab_im (ab_im ),
.b_re (b_re ), .abb_re (abb_re ),
.b_im (b_im ), .abb_im (abb_im ),
.dout_re (IIRout_p5_re ), // Re(y(8n-83)) .ab_pow3_re (ab_pow3_re ),
.dout_im (IIRout_p5_im ) // Im(y(8n-83)) .ab_pow3_im (ab_pow3_im ),
); .ab_pow4_re (ab_pow4_re ),
IIR_Filter_p1 #( .ab_pow4_im (ab_pow4_im ),
.cascade_in_width (temp_var_width-10 ) .ab_pow5_re (ab_pow5_re ),
) inst_iir_p6 ( .ab_pow5_im (ab_pow5_im ),
.clk (clk ), .ab_pow6_re (ab_pow6_re ),
.rstn (rstn ), .ab_pow6_im (ab_pow6_im ),
.en (en ), .ab_pow7_re (ab_pow7_re ),
.din_re (IIRin_p5_r12 ), // x(8n-82) .ab_pow7_im (ab_pow7_im ),
.dout_r1_re (IIRout_p5_re ), // Re(y(8n-83)) .b_pow8_re (b_pow8_re ),
.dout_r1_im (IIRout_p5_im ), // Im(y(8n-83)) .b_pow8_im (b_pow8_im ),
.a_re (a_re ), .dout (IIRout_p4 )
.a_im (a_im ), );
.b_re (b_re ), IIR_Filter_p8 inst_iir_p5 (
.b_im (b_im ), .clk (clk ),
.dout_re (IIRout_p6_re ), // Re(y(8n-98)) .rstn (rstn ),
.dout_im (IIRout_p6_im ) // Im(y(8n-98)) .en (en ),
); .dinp0 (IIRin_p[5] ),
IIR_Filter_p1 #( .dinp1 (IIRin_p[4] ),
.cascade_in_width (temp_var_width-12 ) .dinp2 (IIRin_p[3] ),
) inst_iir_p7 ( .dinp3 (IIRin_p[2] ),
.clk (clk ), .dinp4 (IIRin_p[1] ),
.rstn (rstn ), .dinp5 (IIRin_p[0] ),
.en (en ), .dinp6 (IIRin_p_r1[7] ),
.din_re (IIRin_p6_r14 ), // x(8n-97) .dinp7 (IIRin_p_r1[6] ),
.dout_r1_re (IIRout_p6_re ), // Re(y(8n-98)) .a_re (a_re ),
.dout_r1_im (IIRout_p6_im ), // Im(y(8n-98)) .a_im (a_im ),
.a_re (a_re ), .ab_re (ab_re ),
.a_im (a_im ), .ab_im (ab_im ),
.b_re (b_re ), .abb_re (abb_re ),
.b_im (b_im ), .abb_im (abb_im ),
.dout_re (IIRout_p7_re ), // Re(y(8n-113)) .ab_pow3_re (ab_pow3_re ),
.dout_im (IIRout_p7_im ) // Im(y(8n-113)) .ab_pow3_im (ab_pow3_im ),
); .ab_pow4_re (ab_pow4_re ),
.ab_pow4_im (ab_pow4_im ),
assign IIRout_p0 = IIRout_p0_re[temp_var_width- 0-1 : temp_var_width- 0-data_out_width]; // y(8n-8) .ab_pow5_re (ab_pow5_re ),
assign IIRout_p1 = IIRout_p1_re[temp_var_width- 2-1 : temp_var_width- 2-data_out_width]; // y(8n-23) .ab_pow5_im (ab_pow5_im ),
assign IIRout_p2 = IIRout_p2_re[temp_var_width- 4-1 : temp_var_width- 4-data_out_width]; // y(8n-38) .ab_pow6_re (ab_pow6_re ),
assign IIRout_p3 = IIRout_p3_re[temp_var_width- 6-1 : temp_var_width- 6-data_out_width]; // y(8n-53) .ab_pow6_im (ab_pow6_im ),
assign IIRout_p4 = IIRout_p4_re[temp_var_width- 8-1 : temp_var_width- 8-data_out_width]; // y(8n-68) .ab_pow7_re (ab_pow7_re ),
assign IIRout_p5 = IIRout_p5_re[temp_var_width-10-1 : temp_var_width-10-data_out_width]; // y(8n-83) .ab_pow7_im (ab_pow7_im ),
assign IIRout_p6 = IIRout_p6_re[temp_var_width-12-1 : temp_var_width-12-data_out_width]; // y(8n-98) .b_pow8_re (b_pow8_re ),
assign IIRout_p7 = IIRout_p7_re[temp_var_width-14-1 : temp_var_width-14-data_out_width]; // y(8n-113) .b_pow8_im (b_pow8_im ),
.dout (IIRout_p5 )
);
IIR_Filter_p8 inst_iir_p6 (
.clk (clk ),
.rstn (rstn ),
.en (en ),
.dinp0 (IIRin_p[6] ),
.dinp1 (IIRin_p[5] ),
.dinp2 (IIRin_p[4] ),
.dinp3 (IIRin_p[3] ),
.dinp4 (IIRin_p[2] ),
.dinp5 (IIRin_p[1] ),
.dinp6 (IIRin_p[0] ),
.dinp7 (IIRin_p_r1[7] ),
.a_re (a_re ),
.a_im (a_im ),
.ab_re (ab_re ),
.ab_im (ab_im ),
.abb_re (abb_re ),
.abb_im (abb_im ),
.ab_pow3_re (ab_pow3_re ),
.ab_pow3_im (ab_pow3_im ),
.ab_pow4_re (ab_pow4_re ),
.ab_pow4_im (ab_pow4_im ),
.ab_pow5_re (ab_pow5_re ),
.ab_pow5_im (ab_pow5_im ),
.ab_pow6_re (ab_pow6_re ),
.ab_pow6_im (ab_pow6_im ),
.ab_pow7_re (ab_pow7_re ),
.ab_pow7_im (ab_pow7_im ),
.b_pow8_re (b_pow8_re ),
.b_pow8_im (b_pow8_im ),
.dout (IIRout_p6 )
);
IIR_Filter_p8 inst_iir_p7 (
.clk (clk ),
.rstn (rstn ),
.en (en ),
.dinp0 (IIRin_p[7] ),
.dinp1 (IIRin_p[6] ),
.dinp2 (IIRin_p[5] ),
.dinp3 (IIRin_p[4] ),
.dinp4 (IIRin_p[3] ),
.dinp5 (IIRin_p[2] ),
.dinp6 (IIRin_p[1] ),
.dinp7 (IIRin_p[0] ),
.a_re (a_re ),
.a_im (a_im ),
.ab_re (ab_re ),
.ab_im (ab_im ),
.abb_re (abb_re ),
.abb_im (abb_im ),
.ab_pow3_re (ab_pow3_re ),
.ab_pow3_im (ab_pow3_im ),
.ab_pow4_re (ab_pow4_re ),
.ab_pow4_im (ab_pow4_im ),
.ab_pow5_re (ab_pow5_re ),
.ab_pow5_im (ab_pow5_im ),
.ab_pow6_re (ab_pow6_re ),
.ab_pow6_im (ab_pow6_im ),
.ab_pow7_re (ab_pow7_re ),
.ab_pow7_im (ab_pow7_im ),
.b_pow8_re (b_pow8_re ),
.b_pow8_im (b_pow8_im ),
.dout (IIRout_p7 )
);
endmodule endmodule

991
rtl/z_dsp/TailCorr_top.v
View File

File diff suppressed because it is too large Load Diff

56
rtl/z_dsp/Trunc.v
View File

@ -1,56 +0,0 @@
module trunc #(
parameter integer diw = 8
//,parameter integer dow = msb - (lsb -1)
,parameter integer msb = 7
,parameter integer lsb = 1
,parameter integer half_precision = 0
)
(
input clk
,input rstn
,input en
,input signed [diw - 1 :0] din
,output signed [msb - lsb:0] dout
);
reg signed [msb - lsb : 0] d_tmp;
generate
if(lsb!=0 && half_precision != 0) begin
always @(posedge clk or negedge rstn) begin
if (!rstn) begin
d_tmp <= 'h0;
end
else if(en) begin
if(din[diw-1 : msb] != {(diw-msb){din[diw-1]}})
d_tmp <= {{din[diw-1]}, {(msb-lsb){!din[diw-1]}}};
else
d_tmp <= din[msb:lsb] + {{(msb-lsb){1'b0}},din[lsb-1]};
end
else begin
d_tmp <= d_tmp;
end
end
end
else begin
always @(posedge clk or negedge rstn) begin
if (!rstn) begin
d_tmp <= 'h0;
end
else if(en) begin
if(din[diw-1 : msb] != {(diw-msb){din[diw-1]}})
d_tmp <= {{din[diw-1]}, {(msb-lsb){!din[diw-1]}}};
else
d_tmp <= din[msb:lsb];
end
else begin
d_tmp <= d_tmp;
end
end
end
endgenerate
assign dout = d_tmp;
endmodule

96
rtl/z_dsp/mult_C.v
View File

@ -36,7 +36,7 @@ module mult_C #(
,parameter integer B_width = 8 ,parameter integer B_width = 8
,parameter integer C_width = 8 ,parameter integer C_width = 8
,parameter integer D_width = 8 ,parameter integer D_width = 8
,parameter integer o_width = 31//division ,parameter integer frac_coef_width = 31//division
) )
@ -46,72 +46,66 @@ module mult_C #(
en, en,
a, a,
b, b,
c, c,
d, d,
Re, Re,
Im Im
); );
input rstn; input rstn;
input clk; input clk;
input en; input en;
input signed [A_width-1 :0] a; input signed [A_width-1:0] a;
input signed [B_width-1 :0] b; input signed [B_width-1:0] b;
input signed [C_width-1 :0] c; input signed [C_width-1:0] c;
input signed [D_width-1 :0] d; input signed [D_width-1:0] d;
output signed [o_width-1 :0] Re; output signed [A_width+C_width-frac_coef_width-2:0] Re;
output signed [o_width-1 :0] Im; output signed [A_width+D_width-frac_coef_width-2:0] Im;
wire signed [A_width+C_width-1:0] ac; wire signed [A_width+C_width-1:0] ac;
wire signed [B_width+D_width-1:0] bd; wire signed [B_width+D_width-1:0] bd;
wire signed [A_width+D_width-1:0] ad; wire signed [A_width+D_width-1:0] ad;
wire signed [B_width+C_width-1:0] bc; wire signed [B_width+C_width-1:0] bc;
wire signed [A_width+C_width :0] Re_tmp;
wire signed [A_width+D_width :0] Im_tmp;
wire signed [o_width-1 :0] Re_trunc;
wire signed [o_width-1 :0] Im_trunc;
DW02_mult #(A_width,C_width) inst_c1( .A (a ),
.B (c ),
.TC (1'b1 ),
.PRODUCT (ac )
);
DW02_mult #(B_width,D_width) inst_c2( .A (b ), DW02_mult #(A_width,C_width) inst_c1( .A (a ),
.B (d ), .B (c ),
.TC (1'b1 ), .TC (1'b1 ),
.PRODUCT (bd ) .PRODUCT (ac )
); );
DW02_mult #(A_width,D_width) inst_c3( .A (a ), DW02_mult #(B_width,D_width) inst_c2( .A (b ),
.B (d ), .B (d ),
.TC (1'b1 ), .TC (1'b1 ),
.PRODUCT (ad ) .PRODUCT (bd )
); );
DW02_mult #(B_width,C_width) inst_c4( .A (b ),
.B (c ), DW02_mult #(A_width,D_width) inst_c3( .A (a ),
.TC (1'b1 ), .B (d ),
.PRODUCT (bc ) .TC (1'b1 ),
); .PRODUCT (ad )
);
DW02_mult #(B_width,C_width) inst_c4( .A (b ),
.B (c ),
.TC (1'b1 ),
.PRODUCT (bc )
);
wire signed [A_width+C_width:0] Re_tmp;
wire signed [A_width+D_width:0] Im_tmp;
assign Re_tmp = ac - bd; assign Re_tmp = ac - bd;
assign Im_tmp = ad + bc; assign Im_tmp = ad + bc;
trunc #( wire signed [A_width+C_width:0] Re_round;
.diw (A_width+C_width+1 ) wire signed [A_width+D_width:0] Im_round;
,.msb (A_width+C_width-2 )
,.lsb (A_width+C_width-o_width-1 ) FixRound #(A_width+C_width+1,frac_coef_width) u_round1 (clk, rstn, en, Re_tmp, Re_round);
) u_round1 (clk, rstn, en, Re_tmp, Re_trunc); FixRound #(A_width+C_width+1,frac_coef_width) u_round2 (clk, rstn, en, Im_tmp, Im_round);
trunc #(
.diw (A_width+D_width+1 )
,.msb (A_width+D_width-2 )
,.lsb (A_width+C_width-o_width-1 )
) u_round2 (clk, rstn, en, Im_tmp, Im_trunc);
// Since this is complex multiplication, the output bit width needs to be increased by one compared to the input. // Since this is complex multiplication, the output bit width needs to be increased by one compared to the input.
assign Re = Re_trunc; assign Re = Re_round[A_width+D_width-2:frac_coef_width];
assign Im = Im_trunc; assign Im = Im_round[A_width+D_width-2:frac_coef_width];
endmodule endmodule

42
rtl/z_dsp/mult_x.v
View File

@ -35,7 +35,7 @@ module mult_x #(
parameter integer A_width = 8 parameter integer A_width = 8
,parameter integer C_width = 8 ,parameter integer C_width = 8
,parameter integer D_width = 8 ,parameter integer D_width = 8
,parameter integer o_width = 31//division ,parameter integer frac_coef_width = 31//division
) )
@ -53,17 +53,15 @@ module mult_x #(
input rstn; input rstn;
input clk; input clk;
input en; input en;
input signed [A_width-1 :0] a; input signed [A_width-1:0] a;
input signed [C_width-1 :0] c; input signed [C_width-1:0] c;
input signed [D_width-1 :0] d; input signed [D_width-1:0] d;
output signed [o_width-1 :0] Re; output signed [A_width+C_width-frac_coef_width-2:0] Re;
output signed [o_width-1 :0] Im; output signed [A_width+D_width-frac_coef_width-2:0] Im;
wire signed [A_width+C_width-1:0] ac; wire signed [A_width+C_width-1:0] ac;
wire signed [A_width+D_width-1:0] ad; wire signed [A_width+D_width-1:0] ad;
wire signed [o_width-1 :0] Re_trunc;
wire signed [o_width-1 :0] Im_trunc;
@ -73,27 +71,27 @@ DW02_mult #(A_width,C_width) inst_c1( .A (a ),
.PRODUCT (ac ) .PRODUCT (ac )
); );
DW02_mult #(A_width,D_width) inst_c3( .A (a ), DW02_mult #(A_width,D_width) inst_c3( .A (a ),
.B (d ), .B (d ),
.TC (1'b1 ), .TC (1'b1 ),
.PRODUCT (ad ) .PRODUCT (ad )
); );
wire signed [A_width+C_width:0] Re_tmp;
wire signed [A_width+D_width:0] Im_tmp;
assign Re_tmp = ac;
assign Im_tmp = ad;
trunc #( wire signed [A_width+C_width:0] Re_round;
.diw (A_width+C_width ) wire signed [A_width+D_width:0] Im_round;
,.msb (A_width+C_width-2 )
,.lsb (A_width+C_width-o_width-1 ) FixRound #(A_width+C_width+1,frac_coef_width) u_round1 (clk, rstn, en, Re_tmp, Re_round);
) u_round1 (clk, rstn, en, ac, Re_trunc); FixRound #(A_width+C_width+1,frac_coef_width) u_round2 (clk, rstn, en, Im_tmp, Im_round);
trunc #(
.diw (A_width+D_width )
,.msb (A_width+D_width-2 )
,.lsb (A_width+D_width-o_width-1 )
) u_round2 (clk, rstn, en, ad, Im_trunc);
// Since this is complex multiplication, the output bit width needs to be increased by one compared to the input. // Since this is complex multiplication, the output bit width needs to be increased by one compared to the input.
assign Re = Re_trunc; assign Re = Re_round[A_width+D_width-2:frac_coef_width];
assign Im = Im_trunc; assign Im = Im_round[A_width+D_width-2:frac_coef_width];
endmodule endmodule

10
script_m/TailCorr_Test_Verdi.m
View File

@ -2,11 +2,11 @@
clc;clear;close all clc;clear;close all
% addpath("/data/work/thfu/TailCorr/script_m"); % addpath("/data/work/thfu/TailCorr/script_m");
data_source = 'matlab'; data_source = 'matlab';
file_path = "/home/thfu/work/TailCorr/sim/"; file_path = "/home/thfu/work/TailCorr/sim/z_dsp/";
rng('shuffle'); rng('shuffle');
if strcmp(data_source, 'matlab') if strcmp(data_source, 'matlab')
in = floor(cat(1,0,30000*ones(4*2579+4,1))); in = floor(cat(1,0,3000*randn(4*2579+4,1)));
for i = 0:3 for i = 0:3
filename = strcat(file_path, "in", num2str(i), "_matlab.dat"); filename = strcat(file_path, "in", num2str(i), "_matlab.dat");
subset = in(i+1:4:end); subset = in(i+1:4:end);
@ -33,14 +33,14 @@ end
cs_wave = []; cs_wave = [];
for i = 0:7 for i = 0:3
filename = strcat(file_path, "dout", num2str(i), ".dat"); filename = strcat(file_path, "dout", num2str(i), ".dat");
dout_data = importdata(filename); dout_data = importdata(filename);
if isempty(cs_wave) if isempty(cs_wave)
N = length(dout_data); N = length(dout_data);
cs_wave = zeros(8*N, 1); cs_wave = zeros(4*N, 1);
end end
cs_wave(i+1:8:end) = dout_data; cs_wave(i+1:4:end) = dout_data;
end end
A = [0.025 0.015*1 0.0002*1 0]; A = [0.025 0.015*1 0.0002*1 0];

26
sim/TailCorr_en/files.f
View File

@ -1,21 +1,11 @@
../rtl/z_dsp/mult_C.v ../../rtl/z_dsp/mult_C.v
../rtl/z_dsp/mult_x.v ../../rtl/z_dsp/FixRound.v
../rtl/z_dsp/Trunc.v ../../rtl/z_dsp/TailCorr_top.v
../rtl/z_dsp/TailCorr_top.v ../../rtl/z_dsp/IIR_top.v
../rtl/z_dsp/IIR_top.v ../../rtl/z_dsp/diff_p.v
../rtl/z_dsp/diff_p.v ../../rtl/z_dsp/s2p_2.v
../rtl/z_dsp/s2p_2.v ../../rtl/z_dsp/IIR_Filter_p8.v
../rtl/z_dsp/IIR_Filter_p8.v ../../rtl/model/DW02_mult.v
../rtl/z_dsp/IIR_Filter_p1.v
../rtl/z_dsp/sirv_gnrl_dffs.v
../rtl/ref/mult_C.v
../rtl/ref/FixRound.v
../rtl/ref/TailCorr_top.v
../rtl/ref/IIR_top.v
../rtl/ref/diff_p.v
../rtl/ref/s2p_2.v
../rtl/ref/IIR_Filter_p8.v
../rtl/model/DW02_mult.v
tb_TailCorr_en.v tb_TailCorr_en.v

185
sim/TailCorr_en/tb_TailCorr_en.v
View File

@ -48,8 +48,6 @@ end
reg rstn; reg rstn;
reg [31:0] a_re0; reg [31:0] a_re0;
reg [31:0] a_im0; reg [31:0] a_im0;
reg [31:0] b_re0;
reg [31:0] b_im0;
reg [31:0] ab_re0; reg [31:0] ab_re0;
reg [31:0] ab_im0; reg [31:0] ab_im0;
reg [31:0] abb_re0; reg [31:0] abb_re0;
@ -68,8 +66,6 @@ reg [31:0] b_pow8_re0;
reg [31:0] b_pow8_im0; reg [31:0] b_pow8_im0;
reg [31:0] a_re1; reg [31:0] a_re1;
reg [31:0] a_im1; reg [31:0] a_im1;
reg [31:0] b_re1;
reg [31:0] b_im1;
reg [31:0] ab_re1; reg [31:0] ab_re1;
reg [31:0] ab_im1; reg [31:0] ab_im1;
reg [31:0] abb_re1; reg [31:0] abb_re1;
@ -88,8 +84,6 @@ reg [31:0] b_pow8_re1;
reg [31:0] b_pow8_im1; reg [31:0] b_pow8_im1;
reg [31:0] a_re2; reg [31:0] a_re2;
reg [31:0] a_im2; reg [31:0] a_im2;
reg [31:0] b_re2;
reg [31:0] b_im2;
reg [31:0] ab_re2; reg [31:0] ab_re2;
reg [31:0] ab_im2; reg [31:0] ab_im2;
reg [31:0] abb_re2; reg [31:0] abb_re2;
@ -108,8 +102,6 @@ reg [31:0] b_pow8_re2;
reg [31:0] b_pow8_im2; reg [31:0] b_pow8_im2;
reg [31:0] a_re3; reg [31:0] a_re3;
reg [31:0] a_im3; reg [31:0] a_im3;
reg [31:0] b_re3;
reg [31:0] b_im3;
reg [31:0] ab_re3; reg [31:0] ab_re3;
reg [31:0] ab_im3; reg [31:0] ab_im3;
reg [31:0] abb_re3; reg [31:0] abb_re3;
@ -128,8 +120,6 @@ reg [31:0] b_pow8_re3;
reg [31:0] b_pow8_im3; reg [31:0] b_pow8_im3;
reg [31:0] a_re4; reg [31:0] a_re4;
reg [31:0] a_im4; reg [31:0] a_im4;
reg [31:0] b_re4;
reg [31:0] b_im4;
reg [31:0] ab_re4; reg [31:0] ab_re4;
reg [31:0] ab_im4; reg [31:0] ab_im4;
reg [31:0] abb_re4; reg [31:0] abb_re4;
@ -148,8 +138,6 @@ reg [31:0] b_pow8_re4;
reg [31:0] b_pow8_im4; reg [31:0] b_pow8_im4;
reg [31:0] a_re5; reg [31:0] a_re5;
reg [31:0] a_im5; reg [31:0] a_im5;
reg [31:0] b_re5;
reg [31:0] b_im5;
reg [31:0] ab_re5; reg [31:0] ab_re5;
reg [31:0] ab_im5; reg [31:0] ab_im5;
reg [31:0] abb_re5; reg [31:0] abb_re5;
@ -192,18 +180,6 @@ begin
a_im3 = 32'd0; a_im3 = 32'd0;
a_im4 = 32'd0; a_im4 = 32'd0;
a_im5 = 32'd0; a_im5 = 32'd0;
b_re0 = 32'd2143083068;
b_re1 = 32'd2145807236;
b_re2 = 32'd2146812530;
b_re3 = 32'd2147483648;
b_re4 = 32'd0;
b_re5 = 32'd0;
b_im0 = 32'd0;
b_im1 = 32'd0;
b_im2 = 32'd0;
b_im3 = 32'd0;
b_im4 = 32'd0;
b_im5 = 32'd0;
ab_re0 = 32'd54894517; ab_re0 = 32'd54894517;
ab_re1 = 32'd32664510; ab_re1 = 32'd32664510;
ab_re2 = 32'd429381 ; ab_re2 = 32'd429381 ;
@ -292,7 +268,7 @@ begin
b_pow8_re0 = 32'd2112530470; b_pow8_re0 = 32'd2112530470;
b_pow8_re1 = 32'd2134108939; b_pow8_re1 = 32'd2134108939;
b_pow8_re2 = 32'd2142120573; b_pow8_re2 = 32'd2142120573;
b_pow8_re3 = 32'd2147483648; b_pow8_re3 = 32'd0;
b_pow8_re4 = 32'd0; b_pow8_re4 = 32'd0;
b_pow8_re5 = 32'd0; b_pow8_re5 = 32'd0;
b_pow8_im0 = 32'd0; b_pow8_im0 = 32'd0;
@ -436,7 +412,6 @@ assign dac_mode_sel = 2'b00;
wire tc_bypass; wire tc_bypass;
wire vldo; wire vldo;
wire vldo_ref;
assign tc_bypass = 1'b0; assign tc_bypass = 1'b0;
@ -450,7 +425,6 @@ always @(posedge clk or negedge rstn)begin
end end
end end
wire signed [15:0] dout_p[7:0]; wire signed [15:0] dout_p[7:0];
wire signed [15:0] dout_ref_p[7:0];
TailCorr_top inst_TailCorr_top TailCorr_top inst_TailCorr_top
@ -467,8 +441,6 @@ TailCorr_top inst_TailCorr_top
.din3 (iir_in[3]), .din3 (iir_in[3]),
.a_re0 (a_re0), .a_re0 (a_re0),
.a_im0 (a_im0), .a_im0 (a_im0),
.b_re0 (b_re0),
.b_im0 (b_im0),
.ab_re0 (ab_re0), .ab_re0 (ab_re0),
.ab_im0 (ab_im0), .ab_im0 (ab_im0),
.abb_re0 (abb_re0), .abb_re0 (abb_re0),
@ -487,8 +459,6 @@ TailCorr_top inst_TailCorr_top
.b_pow8_im0 (b_pow8_im0), .b_pow8_im0 (b_pow8_im0),
.a_re1 (a_re1), .a_re1 (a_re1),
.a_im1 (a_im1), .a_im1 (a_im1),
.b_re1 (b_re1),
.b_im1 (b_im1),
.ab_re1 (ab_re1), .ab_re1 (ab_re1),
.ab_im1 (ab_im1), .ab_im1 (ab_im1),
.abb_re1 (abb_re1), .abb_re1 (abb_re1),
@ -507,8 +477,6 @@ TailCorr_top inst_TailCorr_top
.b_pow8_im1 (b_pow8_im1), .b_pow8_im1 (b_pow8_im1),
.a_re2 (a_re2), .a_re2 (a_re2),
.a_im2 (a_im2), .a_im2 (a_im2),
.b_re2 (b_re2),
.b_im2 (b_im2),
.ab_re2 (ab_re2), .ab_re2 (ab_re2),
.ab_im2 (ab_im2), .ab_im2 (ab_im2),
.abb_re2 (abb_re2), .abb_re2 (abb_re2),
@ -527,8 +495,6 @@ TailCorr_top inst_TailCorr_top
.b_pow8_im2 (b_pow8_im2), .b_pow8_im2 (b_pow8_im2),
.a_re3 (a_re3), .a_re3 (a_re3),
.a_im3 (a_im3), .a_im3 (a_im3),
.b_re3 (b_re3),
.b_im3 (b_im3),
.ab_re3 (ab_re3), .ab_re3 (ab_re3),
.ab_im3 (ab_im3), .ab_im3 (ab_im3),
.abb_re3 (abb_re3), .abb_re3 (abb_re3),
@ -547,8 +513,6 @@ TailCorr_top inst_TailCorr_top
.b_pow8_im3 (b_pow8_im3), .b_pow8_im3 (b_pow8_im3),
.a_re4 (a_re4), .a_re4 (a_re4),
.a_im4 (a_im4), .a_im4 (a_im4),
.b_re4 (b_re4),
.b_im4 (b_im4),
.ab_re4 (ab_re4), .ab_re4 (ab_re4),
.ab_im4 (ab_im4), .ab_im4 (ab_im4),
.abb_re4 (abb_re4), .abb_re4 (abb_re4),
@ -567,8 +531,6 @@ TailCorr_top inst_TailCorr_top
.b_pow8_im4 (b_pow8_im4), .b_pow8_im4 (b_pow8_im4),
.a_re5 (a_re5), .a_re5 (a_re5),
.a_im5 (a_im5), .a_im5 (a_im5),
.b_re5 (b_re5),
.b_im5 (b_im5),
.ab_re5 (ab_re5), .ab_re5 (ab_re5),
.ab_im5 (ab_im5), .ab_im5 (ab_im5),
.abb_re5 (abb_re5), .abb_re5 (abb_re5),
@ -598,151 +560,6 @@ TailCorr_top inst_TailCorr_top
); );
TailCorr_top_ref refm_TailCorr_top
(
.clk (clk ),
.en (en ),
.rstn (rstn ),
.vldi (vldi_matlab[0] ),
// .dac_mode_sel (dac_mode_sel ),
// .intp_mode (intp_mode ),
.din0 (iir_in[0]),
.din1 (iir_in[1]),
.din2 (iir_in[2]),
.din3 (iir_in[3]),
.a_re0 (a_re0),
.a_im0 (a_im0),
//.b_re0 (b_re0),
//.b_im0 (b_im0),
.ab_re0 (ab_re0),
.ab_im0 (ab_im0),
.abb_re0 (abb_re0),
.abb_im0 (abb_im0),
.ab_pow3_re0 (ab_pow3_re0),
.ab_pow3_im0 (ab_pow3_im0),
.ab_pow4_re0 (ab_pow4_re0),
.ab_pow4_im0 (ab_pow4_im0),
.ab_pow5_re0 (ab_pow5_re0),
.ab_pow5_im0 (ab_pow5_im0),
.ab_pow6_re0 (ab_pow6_re0),
.ab_pow6_im0 (ab_pow6_im0),
.ab_pow7_re0 (ab_pow7_re0),
.ab_pow7_im0 (ab_pow7_im0),
.b_pow8_re0 (b_pow8_re0),
.b_pow8_im0 (b_pow8_im0),
.a_re1 (a_re1),
.a_im1 (a_im1),
//.b_re1 (b_re1),
//.b_im1 (b_im1),
.ab_re1 (ab_re1),
.ab_im1 (ab_im1),
.abb_re1 (abb_re1),
.abb_im1 (abb_im1),
.ab_pow3_re1 (ab_pow3_re1),
.ab_pow3_im1 (ab_pow3_im1),
.ab_pow4_re1 (ab_pow4_re1),
.ab_pow4_im1 (ab_pow4_im1),
.ab_pow5_re1 (ab_pow5_re1),
.ab_pow5_im1 (ab_pow5_im1),
.ab_pow6_re1 (ab_pow6_re1),
.ab_pow6_im1 (ab_pow6_im1),
.ab_pow7_re1 (ab_pow7_re1),
.ab_pow7_im1 (ab_pow7_im1),
.b_pow8_re1 (b_pow8_re1),
.b_pow8_im1 (b_pow8_im1),
.a_re2 (a_re2),
.a_im2 (a_im2),
//.b_re2 (b_re2),
//.b_im2 (b_im2),
.ab_re2 (ab_re2),
.ab_im2 (ab_im2),
.abb_re2 (abb_re2),
.abb_im2 (abb_im2),
.ab_pow3_re2 (ab_pow3_re2),
.ab_pow3_im2 (ab_pow3_im2),
.ab_pow4_re2 (ab_pow4_re2),
.ab_pow4_im2 (ab_pow4_im2),
.ab_pow5_re2 (ab_pow5_re2),
.ab_pow5_im2 (ab_pow5_im2),
.ab_pow6_re2 (ab_pow6_re2),
.ab_pow6_im2 (ab_pow6_im2),
.ab_pow7_re2 (ab_pow7_re2),
.ab_pow7_im2 (ab_pow7_im2),
.b_pow8_re2 (b_pow8_re2),
.b_pow8_im2 (b_pow8_im2),
.a_re3 (a_re3),
.a_im3 (a_im3),
//.b_re3 (b_re3),
//.b_im3 (b_im3),
.ab_re3 (ab_re3),
.ab_im3 (ab_im3),
.abb_re3 (abb_re3),
.abb_im3 (abb_im3),
.ab_pow3_re3 (ab_pow3_re3),
.ab_pow3_im3 (ab_pow3_im3),
.ab_pow4_re3 (ab_pow4_re3),
.ab_pow4_im3 (ab_pow4_im3),
.ab_pow5_re3 (ab_pow5_re3),
.ab_pow5_im3 (ab_pow5_im3),
.ab_pow6_re3 (ab_pow6_re3),
.ab_pow6_im3 (ab_pow6_im3),
.ab_pow7_re3 (ab_pow7_re3),
.ab_pow7_im3 (ab_pow7_im3),
.b_pow8_re3 (b_pow8_re3),
.b_pow8_im3 (b_pow8_im3),
.a_re4 (a_re4),
.a_im4 (a_im4),
//.b_re4 (b_re4),
//.b_im4 (b_im4),
.ab_re4 (ab_re4),
.ab_im4 (ab_im4),
.abb_re4 (abb_re4),
.abb_im4 (abb_im4),
.ab_pow3_re4 (ab_pow3_re4),
.ab_pow3_im4 (ab_pow3_im4),
.ab_pow4_re4 (ab_pow4_re4),
.ab_pow4_im4 (ab_pow4_im4),
.ab_pow5_re4 (ab_pow5_re4),
.ab_pow5_im4 (ab_pow5_im4),
.ab_pow6_re4 (ab_pow6_re4),
.ab_pow6_im4 (ab_pow6_im4),
.ab_pow7_re4 (ab_pow7_re4),
.ab_pow7_im4 (ab_pow7_im4),
.b_pow8_re4 (b_pow8_re4),
.b_pow8_im4 (b_pow8_im4),
.a_re5 (a_re5),
.a_im5 (a_im5),
//.b_re5 (b_re5),
//.b_im5 (b_im5),
.ab_re5 (ab_re5),
.ab_im5 (ab_im5),
.abb_re5 (abb_re5),
.abb_im5 (abb_im5),
.ab_pow3_re5 (ab_pow3_re5),
.ab_pow3_im5 (ab_pow3_im5),
.ab_pow4_re5 (ab_pow4_re5),
.ab_pow4_im5 (ab_pow4_im5),
.ab_pow5_re5 (ab_pow5_re5),
.ab_pow5_im5 (ab_pow5_im5),
.ab_pow6_re5 (ab_pow6_re5),
.ab_pow6_im5 (ab_pow6_im5),
.ab_pow7_re5 (ab_pow7_re5),
.ab_pow7_im5 (ab_pow7_im5),
.b_pow8_re5 (b_pow8_re5),
.b_pow8_im5 (b_pow8_im5),
.dout_p0 (dout_ref_p[0] ),
.dout_p1 (dout_ref_p[1] ),
.dout_p2 (dout_ref_p[2] ),
.dout_p3 (dout_ref_p[3] ),
.dout_p4 (dout_ref_p[4] ),
.dout_p5 (dout_ref_p[5] ),
.dout_p6 (dout_ref_p[6] ),
.dout_p7 (dout_ref_p[7] ),
.vldo (vldo_ref )
);
integer signed In_fid[0:3]; integer signed In_fid[0:3];
integer signed dout_fid[0:7]; integer signed dout_fid[0:7];