Started adding exponentiator module w/ testbench.

2 files changed, 446 insertions, 0 deletions

diff --git a/src/tb/modexp_fpga_model_vectors.v b/src/tb/modexp_fpga_model_vectors.v
index 7a2b8e9..d5284c9 100644
--- a/src/tb/modexp_fpga_model_vectors.v
+++ b/src/tb/modexp_fpga_model_vectors.v
@@ -30,6 +30,16 @@ localparam [383:0] M_FACTOR_384 =
 	 32'h91e92683, 32'hc483bb6c, 32'h0ee1571d, 32'h6e28c2f5, 
 	 32'hff5e6b61, 32'h65fb6164, 32'hd3651e5a, 32'h746b8ca0};
 
+localparam [383:0] D_384 =
+	{32'had24a30c, 32'h766d8dc3, 32'he2100b02, 32'h24d1c4b0, 
+	 32'hbb6a6342, 32'h577df9be, 32'h89bb1ec3, 32'hdc3259f0, 
+	 32'h1a343f93, 32'h57a12599, 32'ha328ae2f, 32'hf85ef401};
+
+localparam [383:0] S_384 =
+	{32'h65752d0f, 32'h9a017293, 32'h36bfa115, 32'h4a7a81fc, 
+	 32'ha76b945b, 32'h49a3f645, 32'h76801499, 32'hb98e6a16, 
+	 32'hd2467b6a, 32'h75b7d614, 32'h0fff0fde, 32'hb31d1819};
+
 localparam [511:0] M_512 =
 	{32'h005536b6, 32'h43ea651f, 32'h2fd3c70a, 32'ha83659cb, 
 	 32'hd0c1f47b, 32'ha8033730, 32'h29c6b082, 32'h6db48613, 
@@ -66,3 +76,15 @@ localparam [511:0] M_FACTOR_512 =
 	 32'h663032a3, 32'h70734b62, 32'h2d30c132, 32'hefa75cc6, 
 	 32'h9f18b32a, 32'h97d6ddf8, 32'h2f6df2d0, 32'he9098874};
 
+localparam [511:0] D_512 =
+	{32'hc9686c43, 32'hbbe28d66, 32'h758ef8bc, 32'h9b7828e5, 
+	 32'h2ec2804a, 32'hb76745de, 32'h83fcbba0, 32'h2d9eba78, 
+	 32'h215f4cc2, 32'hf49387b3, 32'h8ed0b9dc, 32'h6c129231, 
+	 32'h944368be, 32'hdbf2db79, 32'h16323c49, 32'h34cdf801};
+
+localparam [511:0] S_512 =
+	{32'hcc2fc6b6, 32'he4849987, 32'h75773499, 32'hcb0792b0, 
+	 32'he79f4600, 32'hb2d739c5, 32'h1a661ac6, 32'hd3bf2db5, 
+	 32'hfd1e029d, 32'hfe887387, 32'h4312635f, 32'hb2b54b8d, 
+	 32'h5d3b379e, 32'h161eaa4f, 32'hedfd932b, 32'h780f0203};
+
diff --git a/src/tb/tb_exponentiator.v b/src/tb/tb_exponentiator.v
new file mode 100644
index 0000000..8ca9444
--- /dev/null
+++ b/src/tb/tb_exponentiator.v
@@ -0,0 +1,424 @@
+//======================================================================

+//

+// tb_expoentiator.v

+// -----------------------------------------------------------------------------

+// Testbench for Montgomery modular exponentiation block.

+//

+// Authors: Pavel Shatov

+//

+// Copyright (c) 2017, NORDUnet A/S All rights reserved.

+//

+// Redistribution and use in source and binary forms, with or without

+// modification, are permitted provided that the following conditions

+// are met:

+// - Redistributions of source code must retain the above copyright

+//   notice, this list of conditions and the following disclaimer.

+//

+// - Redistributions in binary form must reproduce the above copyright

+//   notice, this list of conditions and the following disclaimer in the

+//   documentation and/or other materials provided with the distribution.

+//

+// - Neither the name of the NORDUnet nor the names of its contributors may

+//   be used to endorse or promote products derived from this software

+//   without specific prior written permission.

+//

+// THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS "AS

+// IS" AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED

+// TO, THE IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A

+// PARTICULAR PURPOSE ARE DISCLAIMED. IN NO EVENT SHALL THE COPYRIGHT

+// HOLDER OR CONTRIBUTORS BE LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL,

+// SPECIAL, EXEMPLARY, OR CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT LIMITED

+// TO, PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES; LOSS OF USE, DATA, OR

+// PROFITS; OR BUSINESS INTERRUPTION) HOWEVER CAUSED AND ON ANY THEORY OF

+// LIABILITY, WHETHER IN CONTRACT, STRICT LIABILITY, OR TORT (INCLUDING

+// NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY OUT OF THE USE OF THIS

+// SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF SUCH DAMAGE.

+//

+//======================================================================

+

+`timescale 1ns / 1ps

+

+module tb_exponentiator;

+

+		//

+		// Test Vectors

+		//

+	`include "modexp_fpga_model_vectors.v";

+	

+		//

+		// Parameters

+		//

+	localparam NUM_WORDS_384 = 384 / 32;

+	localparam NUM_WORDS_512 = 512 / 32;

+		

+		//

+		// Clock (100 MHz)

+		//

+	reg clk = 1'b0;

+	always #5 clk = ~clk;

+			

+		//

+		// Inputs

+		//

+	reg				rst_n;

+	reg				ena;

+	

+	reg	[ 3: 0]	n_num_words;

+	reg	[ 8: 0]	d_num_bits;

+

+		//

+		// Outputs

+		//

+	wire				rdy;

+

+		//

+		// Integers

+		//

+	integer w;

+	

+		//

+		// BRAM Interfaces

+		//

+	wire	[ 3: 0]	core_m_addr;

+	wire	[ 3: 0]	core_d_addr;

+	wire	[ 3: 0]	core_n1_addr;

+	wire	[ 3: 0]	core_n2_addr;

+	wire	[ 3: 0]	core_n_coeff1_addr;

+	wire	[ 3: 0]	core_n_coeff2_addr;

+	wire	[ 3: 0]	core_r_addr;

+	

+	wire	[31: 0]	core_m_data;

+	wire	[31: 0]	core_d_data;

+	wire	[31: 0]	core_n1_data;

+	wire	[31: 0]	core_n2_data;

+	wire	[31: 0]	core_n_coeff1_data;

+	wire	[31: 0]	core_n_coeff2_data;

+	wire	[31: 0]	core_r_data_in;

+

+	wire				core_r_wren;

+

+	reg	[ 3: 0]	tb_mdn_addr;
+	reg	[ 3: 0]	tb_r_addr;
+

+	reg	[31:0]	tb_m_data;
+	reg	[31:0]	tb_d_data;
+	reg	[31:0]	tb_n_data;

+	reg	[31:0]	tb_n_coeff_data;
+	wire	[31:0]	tb_r_data;
+	
+	reg				tb_mdn_wren;
+	

+		//

+		// BRAMs

+		//

+	bram_1rw_1ro_readfirst #(.MEM_WIDTH(32), .MEM_ADDR_BITS(4))
+	bram_m (.clk(clk),

+		.a_addr(tb_mdn_addr), .a_wr(tb_mdn_wren), .a_in(tb_m_data), .a_out(),
+		.b_addr(core_m_addr), .b_out(core_m_data));
+

+	bram_1rw_1ro_readfirst #(.MEM_WIDTH(32), .MEM_ADDR_BITS(4))
+	bram_d (.clk(clk),

+		.a_addr(tb_mdn_addr), .a_wr(tb_mdn_wren), .a_in(tb_d_data), .a_out(),
+		.b_addr(core_d_addr), .b_out(core_d_data));
+

+	bram_1rw_1ro_readfirst #(.MEM_WIDTH(32), .MEM_ADDR_BITS(4))
+	bram_n1 (.clk(clk),

+		.a_addr(tb_mdn_addr), .a_wr(tb_mdn_wren), .a_in(tb_n_data), .a_out(),
+		.b_addr(core_n1_addr), .b_out(core_n1_data));
+

+	bram_1rw_1ro_readfirst #(.MEM_WIDTH(32), .MEM_ADDR_BITS(4))
+	bram_n2 (.clk(clk),

+		.a_addr(tb_mdn_addr), .a_wr(tb_mdn_wren), .a_in(tb_n_data), .a_out(),
+		.b_addr(core_n2_addr), .b_out(core_n2_data));
+

+	bram_1rw_1ro_readfirst #(.MEM_WIDTH(32), .MEM_ADDR_BITS(4))
+	bram_n_coeff1 (.clk(clk),

+		.a_addr(tb_mdn_addr), .a_wr(tb_mdn_wren), .a_in(tb_n_coeff_data), .a_out(),
+		.b_addr(core_n_coeff1_addr), .b_out(core_n_coeff1_data));

+

+	bram_1rw_1ro_readfirst #(.MEM_WIDTH(32), .MEM_ADDR_BITS(4))
+	bram_n_coeff2 (.clk(clk),

+		.a_addr(tb_mdn_addr), .a_wr(tb_mdn_wren), .a_in(tb_n_coeff_data), .a_out(),
+		.b_addr(core_n_coeff2_addr), .b_out(core_n_coeff2_data));

+		
+	bram_1rw_1ro_readfirst #(.MEM_WIDTH(32), .MEM_ADDR_BITS(4))
+	bram_r (.clk(clk),

+		.a_addr(core_r_addr), .a_wr(core_r_wren), .a_in(core_r_data_in), .a_out(),
+		.b_addr(tb_r_addr), .b_out(tb_r_data));
+

+		//

+		// UUT

+		//

+	modexpa7_exponentiator #

+	(

+		.OPERAND_ADDR_WIDTH		(4),	// 32 * (2**4) = 512-bit operands

+		.SYSTOLIC_ARRAY_POWER	(2)	// 2 ** 2 = 4-tap systolic array

+	)

+	uut

+	(

+		.clk						(clk), 

+		.rst_n					(rst_n), 

+		

+		.ena						(ena), 

+		.rdy						(rdy), 

+		

+		.m_bram_addr			(core_m_addr),

+		.d_bram_addr			(core_d_addr),

+		.n1_bram_addr			(core_n1_addr),

+		.n2_bram_addr			(core_n2_addr),

+		.n_coeff1_bram_addr	(core_n_coeff1_addr),

+		.n_coeff2_bram_addr	(core_n_coeff2_addr),

+		.r_bram_addr			(core_r_addr),

+

+		.m_bram_out				(core_m_data), 

+		.d_bram_out				(core_d_data), 

+		.n1_bram_out			(core_n1_data), 

+		.n2_bram_out			(core_n2_data), 

+		.n_coeff1_bram_out	(core_n_coeff1_data), 

+		.n_coeff2_bram_out	(core_n_coeff1_data), 

+		

+		.r_bram_in				(core_r_data_in), 

+		.r_bram_wr				(core_r_wren), 

+		

+		.n_num_words			(n_num_words),

+		.d_num_bits				(d_num_bits)

+	);

+

+

+		//

+		// Script

+		//

+	initial begin

+

+		rst_n = 1'b0;

+		ena = 1'b0;

+		

+		#200;		

+		rst_n = 1'b1;

+		#100;

+		

+		test_exponent_384(M_FACTOR_384, D_384, N_384, N_COEFF_384, S_384);

+		//test_exponent_512(M_512);

+		

+	end

+      

+		

+		//

+		// Test Tasks

+		//

+		

+	task test_exponent_384;

+		//
+		input	[383:0] m;

+		input	[383:0] d;

+		input	[383:0] n;

+		input	[383:0] n_coeff;

+		input	[383:0] s;

+		reg   [383:0] r;

+		//

+		integer i;

+		//
+		begin

+			//						

+			n_num_words = 4'd11;								// set number of words

+			d_num_bits = 9'd383;								// set number of bits

+			//

+			write_memory_384(m, d, n, n_coeff);			// fill memory

+			
+			ena = 1;												// start operation
+			#10;													//

+			ena = 0;												// clear flag

+			
+			while (!rdy) #10;									// wait for operation to complete

+			read_memory_384(r);								// get result from memory
+						

+			$display("    calculated: %x", r);			// display result

+			$display("    expected:   %x", s);			//

+							

+				// check calculated value

+			if (r === s) begin
+				$display("        OK");

+				$display("SUCCESS: Test passed.");
+			end else begin

+				$display("        ERROR");

+				$display("FAILURE: Test not passed.");

+			end

+			//

+		end
+		//

+	endtask
+	/*

+	task test_factor_512;

+		//
+		input	[511:0] n;

+		reg	[511:0] f;

+		reg	[511:0] factor;

+		integer i;

+		//
+		begin

+			//			

+			calc_factor_512(n, f);	// calculate factor on-the-fly

+			

+				// make sure, that the value matches the one saved in the include file

+			if (f !== FACTOR_512) begin

+				$display("ERROR: Calculated factor value differs from the one in the test vector!");

+				$finish;

+			end

+			

+			

+			n_num_words = 4'd15;								// set number of words

+			write_memory_512(n);								// fill memory

+			
+			ena = 1;												// start operation
+			#10;													//

+			ena = 0;												// clear flag

+			
+			while (!rdy) #10;									// wait for operation to complete

+			read_memory_512(factor);						// get result from memory
+						

+			$display("    calculated: %x", factor);	// display result

+			$display("    expected:   %x", f);			//

+							

+				// check calculated value

+			if (f === factor) begin
+				$display("        OK");

+				$display("SUCCESS: Test passed.");
+			end else begin

+				$display("        ERROR");

+				$display("FAILURE: Test not passed.");

+			end

+			//

+		end
+		//

+	endtask
+	*/

+

+	//

+	// write_memory_384

+	//

+	task write_memory_384;

+		//

+		input	[383:0] m;

+		input	[383:0] d;

+		input	[383:0] n;

+		input	[383:0] n_coeff;

+		reg	[383:0] m_shreg;

+		reg	[383:0] d_shreg;

+		reg	[383:0] n_shreg;

+		reg	[383:0] n_coeff_shreg;

+		//

+		begin
+			//

+			tb_mdn_wren		= 1;			// start filling memories
+			m_shreg			= m;			// preload shift register

+			d_shreg			= d;			// preload shift register

+			n_shreg			= n;			// preload shift register

+			n_coeff_shreg	= n_coeff;	// preload shift register

+			//
+			for (w=0; w<NUM_WORDS_384; w=w+1) begin							// write all words
+				tb_mdn_addr			= w[3:0];											// set address
+				tb_m_data			= m_shreg[31:0];									// set data
+				tb_d_data			= d_shreg[31:0];									// set data
+				tb_n_data			= n_shreg[31:0];									// set data
+				tb_n_coeff_data	= n_coeff_shreg[31:0];							// set data
+				m_shreg				= {{32{1'bX}}, m_shreg[383:32]};				// update shift register
+				d_shreg				= {{32{1'bX}}, d_shreg[383:32]};				// update shift register
+				n_shreg				= {{32{1'bX}}, n_shreg[383:32]};				// update shift register
+				n_coeff_shreg		= {{32{1'bX}}, n_coeff_shreg[383:32]};		// update shift register
+				#10;																			// wait for 1 clock tick
+			end
+			//
+			tb_mdn_addr			= {4{1'bX}};	// wipe addresses
+			tb_m_data			= {32{1'bX}};	// wipe data
+			tb_d_data			= {32{1'bX}};	// wipe data
+			tb_n_data			= {32{1'bX}};	// wipe data
+			tb_n_coeff_data	= {32{1'bX}};	// wipe data
+			tb_mdn_wren	= 0;				// stop filling memory

+			//

+		end

+		//

+	endtask
+			

+	/*	

+	//

+	// write_memory_512

+	//

+	task write_memory_512;

+		//

+		input	[511:0] n;

+		reg	[511:0] n_shreg;

+		//

+		begin
+			//

+			tb_n_wren	= 1;	// start filling memories
+			n_shreg		= n;	// preload shift register

+			//
+			for (w=0; w<NUM_WORDS_512; w=w+1) begin				// write all words
+				tb_n_addr	= w[3:0];									// set address
+				tb_n_data	= n_shreg[31:0];							// set data
+				n_shreg		= {{32{1'bX}}, n_shreg[511:32]};		// update shift register
+				#10;															// wait for 1 clock tick
+			end
+			//
+			tb_n_addr	= {4{1'bX}};	// wipe addresses
+			tb_n_data	= {32{1'bX}};	// wipe data
+			tb_n_wren	= 0;				// stop filling memory

+			//

+		end

+		//

+	endtask
+	*/

+

+	//

+	// read_memory_384

+	//

+	task read_memory_384;

+		//

+		output	[383:0] r;

+		reg		[383:0] r_shreg;

+		//

+		begin

+			//

+			for (w=0; w<NUM_WORDS_384; w=w+1) begin		// read result word-by-word
+				tb_r_addr	= w[3:0];							// set address
+				#10;													// wait for 1 clock tick
+				r_shreg = {tb_r_data, r_shreg[383:32]};	// store data word
+			end

+			//

+			tb_r_addr = {4{1'bX}};								// wipe address

+			r = r_shreg;											// return

+			//

+		end

+		//

+	endtask
+

+	/*

+	//

+	// read_memory_512

+	//

+	task read_memory_512;

+		//

+		output	[511:0] f;

+		reg		[511:0] f_shreg;

+		//

+		begin

+			//

+			for (w=0; w<NUM_WORDS_512; w=w+1) begin		// read result word-by-word
+				tb_f_addr	= w[3:0];							// set address
+				#10;													// wait for 1 clock tick
+				f_shreg = {tb_f_data, f_shreg[511:32]};	// store data word
+			end

+			//

+			tb_f_addr = {4{1'bX}};								// wipe address

+			f = f_shreg;											// return

+			//

+		end

+		//

+	endtask

+	*/

+

+endmodule

+

+//======================================================================

+// End of file

+//======================================================================