From d88715489690e1d77558bb2d89adce92ecabdc84 Mon Sep 17 00:00:00 2001 From: "Pavel V. Shatov (Meister)" Date: Tue, 18 Jul 2017 02:26:18 +0300 Subject: Started adding exponentiator module w/ testbench. --- src/tb/modexp_fpga_model_vectors.v | 22 ++ src/tb/tb_exponentiator.v | 424 +++++++++++++++++++++++++++++++++++++ 2 files changed, 446 insertions(+) create mode 100644 src/tb/tb_exponentiator.v (limited to 'src/tb') 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