Subversion Repositories dblclockfft

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

//

// Filename:    mpy_tb.cpp

//

// Project:     A Doubletime Pipelined FFT

//

// Purpose:     A test-bench for the shift and add shiftaddmpy.v subfile of

//              the double clocked FFT.  This file may be run autonomously. 

//              If so, the last line output will either read "SUCCESS" on

//              success, or some other failure message otherwise.

//

//              This file depends upon verilator to both compile, run, and

//              therefore test shiftaddmpy.v

//

// Creator:     Dan Gisselquist, Ph.D.

//              Gisselquist Tecnology, LLC

//

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

//

// Copyright (C) 2015, Gisselquist Technology, LLC

//

// This program is free software (firmware): you can redistribute it and/or

// modify it under the terms of  the GNU General Public License as published

// by the Free Software Foundation, either version 3 of the License, or (at

// your option) any later version.

//

// This program is distributed in the hope that it will be useful, but WITHOUT

// ANY WARRANTY; without even the implied warranty of MERCHANTIBILITY or

// FITNESS FOR A PARTICULAR PURPOSE.  See the GNU General Public License

// for more details.

//

// You should have received a copy of the GNU General Public License along

// with this program.  (It's in the $(ROOT)/doc directory, run make with no

// target there if the PDF file isn't present.)  If not, see

// <http://www.gnu.org/licenses/> for a copy.

//

// License:     GPL, v3, as defined and found on www.gnu.org,

//              http://www.gnu.org/licenses/gpl.html

//

//

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

#include "Vshiftaddmpy.h"

#include "verilated.h"

class   MPYTB {

public:

        Vshiftaddmpy    *mpy;

        long    vals[32];

        int     m_addr;

        MPYTB(void) {

                mpy = new Vshiftaddmpy;

                for(int i=0; i<32; i++)

                        vals[i] = 0;

                m_addr = 0;

        }

        ~MPYTB(void) {

                delete mpy;

        }

        void    tick(void) {

                mpy->i_clk = 0;

                mpy->eval();

                mpy->i_clk = 1;

                mpy->eval();

        }

        void    reset(void) {

                mpy->i_clk = 0;

                mpy->i_ce = 1;

                mpy->i_a = 0;

                mpy->i_b = 0;

                for(int k=0; k<20; k++)

                        tick();

        }

        bool    test(const int ia, const int ib) {

                bool    success;

                int     a, b;

                long    out;

                a = ia & 0x0ffff;

                b = ib & 0x0ffff;

                mpy->i_ce = 1;

                mpy->i_a = a & 0x0ffff;

                mpy->i_b = b & 0x0ffff;

                if (a&0x8000) a |= (-1 << 15);

                if (b&0x8000) b |= (-1 << 15);

                vals[m_addr&31] = (long)a * (long)b;

                tick();

                printf("k=%3d: A =%06x, B =%06x, ANS =%10lx, S=%3d,%3d,%3d,%3d, O = %8x\n",

                        m_addr, a & 0x0ffffff, b & 0x0ffffff,

                        vals[m_addr&31] & (~(-1l<<40)),

                        mpy->v__DOT__sgn,

                        mpy->v__DOT__r_s[0],

                        mpy->v__DOT__r_s[1],

                        mpy->v__DOT__r_s[2],

                        mpy->o_r);

                out = mpy->o_r;

                if (out & (1<<31)) out |= (-1 << 31);

                m_addr++;

                success = (m_addr < 20)||(out == vals[(m_addr-18)&31]);

                if (!success) {

                        fprintf(stderr, "WRONG ANSWER: %8lx != %8lx\n", vals[(m_addr-18)&0x01f], out);

                        exit(-1);

                }

                return success;

        }

};

int     main(int argc, char **argv, char **envp) {

        Verilated::commandArgs(argc, argv);

        MPYTB           *tb = new MPYTB;

        tb->reset();

        for(int k=0; k<15; k++) {

                int     a, b;

                a = (1<<k);

                b = 1;

                tb->test(a, b);

        }

        for(int k=0; k<15; k++) {

                int     a, b, out;

                a = (1<<15);

                b = (1<<k);

                tb->test(a, b);

        }

        for(int k=0; k<200; k++) {

                int     a, b, out;

                tb->test(rand(), rand());

        }

        delete  tb;

        printf("SUCCESS!\n");

        exit(0);

}