Subversion Repositories xilinx_virtex_fp_library

`timescale 1ns / 1ps

//////////////////////////////////////////////////////////////////////////////////

// Company: 

// Engineer: 

// 

// Create Date:    17:53:05 10/15/2013 

// Design Name: 

// Module Name:    Multiply_AccumulateConversion

// Project Name: 

// Target Devices: 

// Tool versions: 

// Description: C ± A*B with mapped conversions, conversion applies to C number

//

// Dependencies: 

//

// Revision: 

// Revision 0.01 - File Created

// Additional Comments: 

//

//////////////////////////////////////////////////////////////////////////////////

module Multiply_AccumulateConversion #( parameter size_mantissa = 24,   //mantissa bits(1.M)

                                                                                parameter size_exponent = 8,    //exponent bits

                                                                                parameter size_counter  = 5,    //log2(size_mantissa) + 1 = 5

                                                                                parameter size_exception_field = 2,     // zero/normal numbers/infinity/NaN

                                                                                parameter zero                  = 00, //00

                                                                                parameter normal_number = 01, //01

                                                                                parameter infinity              = 10, //10

                                                                                parameter NaN                   = 11, //11

                                                                                parameter size_integer                  = 32,

                                                                                parameter counter_integer               = 6, //log2(size_integer) + 1 = 6)

                                                                                parameter [1 : 0] FP_operation   = 0, //00 

                                                                                parameter [1 : 0] FP_to_int              = 1, //01 

                                                                                parameter [1 : 0] int_to_FP              = 2, //10 

                                                                                parameter pipeline              = 0,

                                                                                parameter pipeline_pos  = 0,  //8 bits

                                                                                parameter size = size_exponent + size_mantissa + size_exception_field)

                                                                        (       input clk,

                                                                                input rst,

                                                                                input [1 : 0] conversion,

                                                                                input [size - 1:0] c_number_i,

                                                                                input [size - 1:0] a_number_i,

                                                                                input [size - 1:0] b_number_i,

                                                                                input sub,

                                                                                output[size - 1:0] resulting_number_o);

        parameter size_mul_mantissa = size_mantissa + size_mantissa;

        parameter size_mul_counter      = size_counter + 1;

        parameter max_size                      = (size_integer > size_mantissa)? size_integer : size_mantissa;

        parameter max_counter           = (counter_integer > size_counter)? counter_integer : size_counter;

        parameter size_diff_i_m         = (size_integer > size_mantissa)? (size_integer - size_mantissa) : (size_mantissa - size_integer);

        parameter bias                          = {1'b0,{(size_exponent-1){1'b1}}};

        parameter exp_biased            = bias + size_mantissa;

        parameter exponent                      = (size_mul_mantissa - max_size) + exp_biased;

        parameter subtr                         = max_size -2'd2;

        parameter bias_0_bits           = size_exponent - 1;

        parameter shift_mantissa_0_bits = size_mantissa-1'b1;

        wire [size_exception_field - 1 : 0] sp_case_a_number, sp_case_b_number, sp_case_c_number;

        wire [size_mantissa - 1 : 0] m_a_number, m_b_number, m_c_number;

        wire [size_exponent - 1 : 0] e_a_number, e_b_number, e_c_number;

        wire s_a_number, s_b_number, s_c_number;

        wire [size_exponent     : 0] ab_greater_exponent, c_greater_exponent;

        wire [size_exponent - 1 : 0] exp_difference;

        wire [size_exponent     : 0] exp_inter;

        wire [size_mul_mantissa - 1     : 0] m_ab_mantissa, c_mantissa;

        wire [size_exponent                     : 0] e_ab_number_inter, e_ab_number;

        wire [size_mul_counter - 1      : 0] lz_mul;

        wire zero_flag;

        wire sign_res;

        wire eff_op;

        wire [size_mantissa - 1         : 0] initial_rounding_bits, inter_rounding_bits, final_rounding_bits;

        wire [size_mul_mantissa + 1 : 0] normalized_mantissa, adder_mantissa;

        wire [size_mul_mantissa         : 0] unnormalized_mantissa;

        wire [size_mul_mantissa - 1 : 0] shifted_m_ab, convert_neg_mantissa, mantissa_to_shift;

        wire [size_mul_mantissa - 1 : 0] m_c, m_ab;

        wire [size_exception_field - 1 : 0] sp_case_o, sp_case_result_o;

        wire [size_mantissa - 2 : 0] final_mantissa;

        wire [size_exponent - 1 : 0] final_exponent;

        wire [size_mantissa : 0] rounded_mantissa;

        wire [size_mantissa - 1 : 0] resulted_mantissa;

        wire [size_exponent - 1 : 0] resulted_exponent;

        wire [size_exponent  : 0] subtracter;

        wire [size_mul_mantissa-max_size : 0] max_entityINT_FP_msb;

        wire [size_exponent     : 0] shift_value_when_positive_exponent, shift_value_when_negative_exponent;

        wire [size_exponent - 1 : 0] shift_value, shft_val;

        wire [size_exponent - 1 : 0] max_unadjusted_exponent, max_adjust_exponent, adjust;

        wire [size_exponent - 1 : 0] max_exp_selection;

        wire [size_exponent - 1 : 0] max_resulted_e_o;

        wire [max_size - 1 : 0] max_entityINT_FP, max_entityFP_INT;

        wire lsb_shft_bit;

        wire arith_shift;

        wire max_ovf;

        wire do_conversion;

        assign do_conversion = |conversion; //let me know if there is a conversion

        assign m_a_number                       = {1'b1, a_number_i[size_mantissa - 2 :0]};

        assign m_b_number                       = {1'b1, b_number_i[size_mantissa - 2 :0]};

        assign m_c_number                       = {1'b1, c_number_i[size_mantissa - 2 :0]};

        assign e_a_number                       = a_number_i[size_mantissa + size_exponent - 1 : size_mantissa - 1];

        assign e_b_number                       = b_number_i[size_mantissa + size_exponent - 1 : size_mantissa - 1];

        assign e_c_number                       = c_number_i[size_mantissa + size_exponent - 1 : size_mantissa - 1];

        assign s_a_number                       = a_number_i[size - size_exception_field - 1];

        assign s_b_number                       = b_number_i[size - size_exception_field - 1];

        assign s_c_number                       = c_number_i[size - size_exception_field - 1];

        assign sp_case_a_number = a_number_i[size - 1 : size - size_exception_field];

        assign sp_case_b_number = b_number_i[size - 1 : size - size_exception_field];

        assign sp_case_c_number = c_number_i[size - 1 : size - size_exception_field];

        //instantiate multiply component

        multiply #(     .size_mantissa(size_mantissa),

                                        .size_counter(size_counter),

                                        .size_mul_mantissa(size_mul_mantissa))

                multiply_instance (     .a_mantissa_i(m_a_number),

                                                                        .b_mantissa_i(m_b_number),

                                                                        .mul_mantissa(m_ab_mantissa));

        assign c_mantissa       = {1'b0,m_c_number, {(shift_mantissa_0_bits){1'b0}}};

        assign e_ab_number_inter = e_a_number + e_b_number;

        assign e_ab_number = e_ab_number_inter  - {(bias_0_bits){1'b1}};

        //find the greater exponent

        assign ab_greater_exponent = e_ab_number - e_c_number;

        assign c_greater_exponent = e_c_number - e_ab_number;

        //find the difference between exponents

        assign exp_difference   = (ab_greater_exponent[size_exponent])? c_greater_exponent[size_exponent - 1 : 0] : ab_greater_exponent[size_exponent - 1 : 0];

        assign exp_inter                = (c_greater_exponent[size_exponent])? {1'b0, e_ab_number} : {1'b0, e_c_number};

        //set shifter always on m_ab_number

        assign {m_c, m_ab} = (ab_greater_exponent[size_exponent])? {c_mantissa, m_ab_mantissa} :

                                                        {m_ab_mantissa, c_mantissa};

        assign subtracter =  e_c_number - bias;

        assign shift_value_when_positive_exponent = subtr - subtracter[size_exponent-1 : 0];

        assign shift_value_when_negative_exponent = max_size + (~subtracter[size_exponent-1 : 0]);

        assign shift_value = (subtracter[size_exponent])? shift_value_when_negative_exponent[size_exponent - 1 : 0] :

                             (shift_value_when_positive_exponent[size_exponent])? (~shift_value_when_positive_exponent[size_exponent - 1 : 0]):

                                                                                   shift_value_when_positive_exponent[size_exponent - 1 : 0];

        assign shft_val = do_conversion? shift_value : exp_difference;

        assign convert_neg_mantissa = {1'b0, ~c_number_i[size_mantissa-2 : 0]};

        assign mantissa_to_shift = conversion[0]?        (s_c_number? {{size_mantissa{1'b0}}, convert_neg_mantissa + 1'b1} :

                                                                                                {{size_mantissa{1'b0}}, 1'b1, c_number_i[size_mantissa-2 : 0]}) : m_ab;

        assign arith_shift = conversion[0]? s_c_number : 1'b0;

        //shift m_ab_number                             

        shifter #(      .INPUT_SIZE(size_mul_mantissa),

                                .SHIFT_SIZE(size_exponent),

                                .OUTPUT_SIZE(size_mul_mantissa + size_mantissa),

                                .DIRECTION(1'b0), //0=right, 1=left

                                .PIPELINE(pipeline),

                                .POSITION(pipeline_pos))

                m_b_shifter_instance(   .a(mantissa_to_shift),//mantissa

                                                                .arith(arith_shift),//logical shift

                                                                .shft(shft_val),

                                                                .shifted_a({shifted_m_ab, initial_rounding_bits}));

        assign max_entityFP_INT = {s_c_number, shifted_m_ab[max_size - size_diff_i_m - 1 : 0], initial_rounding_bits[size_mantissa - 1 : size_mantissa - size_diff_i_m + 1]};

        //instantiate effective_op component

        effective_op effective_op_instance(     .sign_a(s_a_number),

                                                                                .sign_b(s_b_number),

                                                                                .sign_c(s_c_number),

                                                                                .sub(sub),

                                                                                .eff_sub(eff_op));

        //instantiate accumulate component

        accumulate #(.size_mul_mantissa(size_mul_mantissa))

                accumulate_instance (   .m_a(m_c),

                                                                .m_b(shifted_m_ab),

                                                                .eff_op(eff_op),

                                                                .adder_mantissa(adder_mantissa));

        //compute unnormalized_mantissa

        assign  unnormalized_mantissa =

                                (adder_mantissa[size_mul_mantissa + 1])?        (~adder_mantissa[size_mul_mantissa : 0]) : adder_mantissa[size_mul_mantissa      : 0];

        assign  inter_rounding_bits = do_conversion?    (s_c_number? {size_mantissa{1'b1}} : {size_mantissa{1'b0}}) :

                                                                                                        ((adder_mantissa[size_mul_mantissa + 1])? ~initial_rounding_bits : initial_rounding_bits);

        assign max_entityINT_FP = do_conversion? (s_c_number? (~c_number_i[max_size-1 : 0]) :  c_number_i[max_size-1 : 0]) :

                                                                                                        unnormalized_mantissa[max_size-1 : 0];

        assign max_entityINT_FP_msb = do_conversion? {(size_mul_mantissa-max_size+1){1'b0}} : unnormalized_mantissa[size_mul_mantissa : max_size];

        assign lsb_shft_bit = (do_conversion)? s_c_number : max_entityINT_FP[0];

        assign max_ovf = do_conversion? 1'b0 : unnormalized_mantissa[size_mul_mantissa];

        //instantiate leading_zeros component

        leading_zeros #(.SIZE_INT(size_mul_mantissa + 1'b1),

                                        .SIZE_COUNTER(size_mul_counter),

                                        .PIPELINE(pipeline))

                leading_zeros_instance( .a({max_entityINT_FP_msb, max_entityINT_FP}),

                                                                .ovf(max_ovf),

                                                                .lz(lz_mul));

        //instantiate shifter component

        shifter #(      .INPUT_SIZE(size_mul_mantissa + size_mantissa + 1),

                                .SHIFT_SIZE(size_mul_counter),

                                .OUTPUT_SIZE(size_mul_mantissa + size_mantissa + 2),

                                .DIRECTION(1'b1),

                                .PIPELINE(pipeline),

                                .POSITION(pipeline_pos))

                shifter_instance(       .a({{max_entityINT_FP_msb, max_entityINT_FP}, inter_rounding_bits}),

                                                        .arith(lsb_shft_bit),

                                                        .shft(lz_mul),

                                                        .shifted_a({normalized_mantissa, final_rounding_bits}));

        //instantiate rounding_component

        rounding #(     .SIZE_MOST_S_MANTISSA(size_mantissa+1),

                                .SIZE_LEAST_S_MANTISSA(size_mul_mantissa+2))

                rounding_instance(      .unrounded_mantissa({1'b0, normalized_mantissa[size_mul_mantissa+1 : size_mantissa + 2]}),

                                    .dummy_bits({normalized_mantissa[size_mantissa + 1 : 0],final_rounding_bits}),

                                    .rounded_mantissa(rounded_mantissa));

        assign max_exp_selection = do_conversion? exponent : exp_inter;

        assign max_adjust_exponent = max_exp_selection - lz_mul;

        assign adjust = do_conversion? size_diff_i_m : 2'd2;

        assign max_unadjusted_exponent = max_adjust_exponent + adjust;

        assign max_resulted_e_o = (do_conversion & ~(|{max_entityINT_FP_msb, max_entityINT_FP}))? bias : max_unadjusted_exponent + rounded_mantissa[size_mantissa];

        assign resulted_exponent = conversion[0]?        max_entityFP_INT[size_mantissa+size_exponent-2 : size_mantissa-1] : max_resulted_e_o;

        assign resulted_mantissa = conversion[0]?        max_entityFP_INT[size_mantissa-1 : 0] :

                                                                                                (rounded_mantissa[size_mantissa])?      (rounded_mantissa[size_mantissa : 1]) :

                                                                                                                                                                                (rounded_mantissa[size_mantissa-1 : 0]);

        //instantiate special_cases_mul_acc component

        special_cases_mul_acc   #(      .size_exception_field(size_exception_field),

                                                                .zero(zero),

                                                                .normal_number(normal_number),

                                                                .infinity(infinity),

                                                                .NaN(NaN))

                special_cases_mul_acc_instance  (       .sp_case_a_number(sp_case_a_number),

                                                                                        .sp_case_b_number(sp_case_b_number),

                                                                                        .sp_case_c_number(sp_case_c_number),

                                                                                        .sp_case_result_o(sp_case_o));

        assign sp_case_result_o = do_conversion? sp_case_c_number : sp_case_o;

        //set zero_flag in case of equal numbers

        assign zero_flag = ~((|{resulted_mantissa,sp_case_o[1]}) & (|sp_case_o));

        //compute resulted_sign

        assign sign_res =       do_conversion? s_c_number : ((eff_op)?  (!c_greater_exponent[size_exponent]?

                                                                                (!ab_greater_exponent[size_exponent]? ~adder_mantissa[size_mul_mantissa+1] : s_c_number) : ~(s_b_number^s_a_number)) : s_c_number);

        assign final_mantissa = resulted_mantissa;

        assign final_exponent = resulted_exponent;

        assign resulting_number_o = (zero_flag)? {size{1'b0}} :{sp_case_result_o, sign_res, final_exponent, final_mantissa};

endmodule