stats: Add statistics collector module and testbench

Signed-off-by: Alex Forencich <alex@alexforencich.com>

tb/stats/taxi_stats_collect/Makefile

@@ -0,0 +1,52 @@
+# SPDX-License-Identifier: CERN-OHL-S-2.0
+#
+# Copyright (c) 2021-2025 FPGA Ninja, LLC
+#
+# Authors:
+# - Alex Forencich
+
+TOPLEVEL_LANG = verilog
+
+SIM ?= verilator
+WAVES ?= 0
+
+COCOTB_HDL_TIMEUNIT = 1ns
+COCOTB_HDL_TIMEPRECISION = 1ps
+
+DUT      = taxi_stats_collect
+COCOTB_TEST_MODULES = test_$(DUT)
+COCOTB_TOPLEVEL     = test_$(DUT)
+MODULE   = $(COCOTB_TEST_MODULES)
+TOPLEVEL = $(COCOTB_TOPLEVEL)
+VERILOG_SOURCES += $(COCOTB_TOPLEVEL).sv
+VERILOG_SOURCES += ../../../rtl/stats/$(DUT).sv
+VERILOG_SOURCES += ../../../rtl/axis/taxi_axis_if.sv
+
+# handle file list files
+process_f_file = $(call process_f_files,$(addprefix $(dir $1),$(shell cat $1)))
+process_f_files = $(foreach f,$1,$(if $(filter %.f,$f),$(call process_f_file,$f),$f))
+uniq_base = $(if $1,$(call uniq_base,$(foreach f,$1,$(if $(filter-out $(notdir $(lastword $1)),$(notdir $f)),$f,))) $(lastword $1))
+VERILOG_SOURCES := $(call uniq_base,$(call process_f_files,$(VERILOG_SOURCES)))
+
+# module parameters
+export PARAM_CNT := 8
+export PARAM_INC_W := 8
+export PARAM_ID_BASE := 0
+export PARAM_UPDATE_PERIOD := 128
+export PARAM_STAT_INC_W := 16
+export PARAM_STAT_ID_W := $(shell python -c "print(($(PARAM_CNT)-1).bit_length())")
+
+ifeq ($(SIM), icarus)
+	PLUSARGS += -fst
+
+	COMPILE_ARGS += $(foreach v,$(filter PARAM_%,$(.VARIABLES)),-P $(COCOTB_TOPLEVEL).$(subst PARAM_,,$(v))=$($(v)))
+else ifeq ($(SIM), verilator)
+	COMPILE_ARGS += $(foreach v,$(filter PARAM_%,$(.VARIABLES)),-G$(subst PARAM_,,$(v))=$($(v)))
+
+	ifeq ($(WAVES), 1)
+		COMPILE_ARGS += --trace-fst
+		VERILATOR_TRACE = 1
+	endif
+endif
+
+include $(shell cocotb-config --makefiles)/Makefile.sim

tb/stats/taxi_stats_collect/test_taxi_stats_collect.py

@@ -0,0 +1,249 @@
+#!/usr/bin/env python
+# SPDX-License-Identifier: CERN-OHL-S-2.0
+"""
+
+Copyright (c) 2021-2025 FPGA Ninja, LLC
+
+Authors:
+- Alex Forencich
+
+"""
+
+import itertools
+import logging
+import os
+import random
+
+import cocotb_test.simulator
+
+import cocotb
+from cocotb.clock import Clock
+from cocotb.queue import Queue
+from cocotb.triggers import RisingEdge, Timer
+from cocotb.regression import TestFactory
+
+from cocotbext.axi import AxiStreamBus, AxiStreamSink
+
+
+class TB(object):
+    def __init__(self, dut):
+        self.dut = dut
+
+        self.log = logging.getLogger("cocotb.tb")
+        self.log.setLevel(logging.DEBUG)
+
+        cocotb.start_soon(Clock(dut.clk, 10, units="ns").start())
+
+        self.stat_sink = AxiStreamSink(AxiStreamBus.from_entity(dut.m_axis_stat), dut.clk, dut.rst)
+
+        for k in range(len(dut.stat_inc)):
+            dut.stat_inc[k].setimmediatevalue(0)
+            dut.stat_valid[k].setimmediatevalue(0)
+        dut.update.setimmediatevalue(0)
+
+    def set_backpressure_generator(self, generator=None):
+        if generator:
+            self.stat_sink.set_pause_generator(generator())
+
+    async def cycle_reset(self):
+        self.dut.rst.setimmediatevalue(0)
+        await RisingEdge(self.dut.clk)
+        await RisingEdge(self.dut.clk)
+        self.dut.rst.value = 1
+        await RisingEdge(self.dut.clk)
+        await RisingEdge(self.dut.clk)
+        self.dut.rst.value = 0
+        await RisingEdge(self.dut.clk)
+        await RisingEdge(self.dut.clk)
+
+
+async def run_test_acc(dut, backpressure_inserter=None):
+
+    tb = TB(dut)
+
+    stat_count = len(dut.stat_valid)
+
+    await tb.cycle_reset()
+
+    tb.set_backpressure_generator(backpressure_inserter)
+
+    for n in range(10):
+        await RisingEdge(dut.clk)
+        dut.stat_inc.value = [k for k in range(stat_count)]
+        dut.stat_valid.value = [1]*stat_count
+        await RisingEdge(dut.clk)
+        dut.stat_inc.value = [0]*stat_count
+        dut.stat_valid.value = [0]*stat_count
+
+        await Timer(1000, 'ns')
+
+    await Timer(1000, 'ns')
+
+    data = [0]*stat_count
+
+    while not tb.stat_sink.empty():
+        stat = await tb.stat_sink.recv()
+
+        assert stat.tdata[0] != 0
+
+        data[stat.tid] += stat.tdata[0]
+
+    print(data)
+
+    for n in range(stat_count):
+        assert data[n] == n*10
+
+    await RisingEdge(dut.clk)
+    await RisingEdge(dut.clk)
+
+
+async def run_stress_test(dut, backpressure_inserter=None):
+
+    tb = TB(dut)
+
+    stat_count = len(dut.stat_valid)
+    stat_inc_width = len(dut.stat_inc) // stat_count
+
+    await tb.cycle_reset()
+
+    tb.set_backpressure_generator(backpressure_inserter)
+
+    async def worker(num, queue_ref, queue_drive, count=1024):
+        for k in range(count):
+            count = random.randrange(1, 2**stat_inc_width)
+
+            await queue_drive.put(count)
+            await queue_ref.put((num, count))
+
+            await Timer(random.randint(1, 100), 'ns')
+
+    workers = []
+    queue_ref = Queue()
+    queue_drive = [Queue() for k in range(stat_count)]
+
+    for k in range(stat_count):
+        workers.append(cocotb.start_soon(worker(k, queue_ref, queue_drive[k], count=1024)))
+
+    async def driver(dut, queues):
+        while True:
+            await RisingEdge(dut.clk)
+
+            inc = [0]*stat_count
+            valid = [0]*stat_count
+            for num, queue in enumerate(queues):
+                if not queue.empty():
+                    count = await queue.get()
+                    inc[num] += count
+                    valid[num] = 1
+
+            dut.stat_inc.value = inc
+            dut.stat_valid.value = valid
+
+    driver = cocotb.start_soon(driver(dut, queue_drive))
+
+    while workers:
+        await workers.pop(0).join()
+
+    await Timer(1000, 'ns')
+
+    driver.kill()
+
+    await Timer(1000, 'ns')
+
+    data_ref = [0]*stat_count
+
+    while not queue_ref.empty():
+        num, count = await queue_ref.get()
+        data_ref[num] += count
+
+    print(data_ref)
+
+    data = [0]*stat_count
+
+    while not tb.stat_sink.empty():
+        stat = await tb.stat_sink.recv()
+
+        assert stat.tdata[0] != 0
+
+        data[stat.tid] += stat.tdata[0]
+
+    print(data)
+
+    assert data == data_ref
+
+    await RisingEdge(dut.clk)
+    await RisingEdge(dut.clk)
+
+
+def cycle_pause():
+    return itertools.cycle([1, 1, 1, 0])
+
+
+if cocotb.SIM_NAME:
+
+    for test in [run_test_acc]:
+
+        factory = TestFactory(test)
+        factory.add_option("backpressure_inserter", [None, cycle_pause])
+        factory.generate_tests()
+
+    factory = TestFactory(run_stress_test)
+    factory.generate_tests()
+
+
+# cocotb-test
+
+tests_dir = os.path.dirname(__file__)
+rtl_dir = os.path.abspath(os.path.join(tests_dir, '..', '..', '..', 'rtl'))
+
+
+def process_f_files(files):
+    lst = {}
+    for f in files:
+        if f[-2:].lower() == '.f':
+            with open(f, 'r') as fp:
+                l = fp.read().split()
+            for f in process_f_files([os.path.join(os.path.dirname(f), x) for x in l]):
+                lst[os.path.basename(f)] = f
+        else:
+            lst[os.path.basename(f)] = f
+    return list(lst.values())
+
+
+def test_taxi_stats_collect(request):
+    dut = "taxi_stats_collect"
+    module = os.path.splitext(os.path.basename(__file__))[0]
+    toplevel = module
+
+    verilog_sources = [
+        os.path.join(tests_dir, f"{toplevel}.sv"),
+        os.path.join(rtl_dir, "stats", f"{dut}.sv"),
+        os.path.join(rtl_dir, "axis", "taxi_axis_if.sv"),
+    ]
+
+    verilog_sources = process_f_files(verilog_sources)
+
+    parameters = {}
+
+    parameters['CNT'] = 8
+    parameters['INC_W'] = 8
+    parameters['ID_BASE'] = 0
+    parameters['UPDATE_PERIOD'] = 128
+    parameters['STAT_INC_W'] = 16
+    parameters['STAT_ID_W'] = (parameters['CNT']-1).bit_length()
+
+    extra_env = {f'PARAM_{k}': str(v) for k, v in parameters.items()}
+
+    sim_build = os.path.join(tests_dir, "sim_build",
+        request.node.name.replace('[', '-').replace(']', ''))
+
+    cocotb_test.simulator.run(
+        simulator="verilator",
+        python_search=[tests_dir],
+        verilog_sources=verilog_sources,
+        toplevel=toplevel,
+        module=module,
+        parameters=parameters,
+        sim_build=sim_build,
+        extra_env=extra_env,
+    )

tb/stats/taxi_stats_collect/test_taxi_stats_collect.sv

@@ -0,0 +1,76 @@
+// SPDX-License-Identifier: CERN-OHL-S-2.0
+/*
+
+Copyright (c) 2025 FPGA Ninja, LLC
+
+Authors:
+- Alex Forencich
+
+*/
+
+`resetall
+`timescale 1ns / 1ps
+`default_nettype none
+
+/*
+ * Statistics collector testbench
+ */
+module test_taxi_stats_collect #
+(
+    /* verilator lint_off WIDTHTRUNC */
+    parameter CNT = 8,
+    parameter INC_W = 8,
+    parameter ID_BASE = 0,
+    parameter UPDATE_PERIOD = 128,
+    parameter STAT_INC_W = 16,
+    parameter STAT_ID_W = $clog2(CNT)
+    /* verilator lint_on WIDTHTRUNC */
+)
+();
+
+logic clk;
+logic rst;
+
+logic [INC_W-1:0] stat_inc[CNT];
+logic [0:0] stat_valid[CNT];
+
+taxi_axis_if #(
+    .DATA_W(STAT_INC_W),
+    .KEEP_EN(0),
+    .KEEP_W(1),
+    .ID_EN(1),
+    .ID_W(STAT_ID_W)
+) m_axis_stat();
+
+logic update;
+
+taxi_stats_collect #(
+    .CNT(CNT),
+    .INC_W(INC_W),
+    .ID_BASE(ID_BASE),
+    .UPDATE_PERIOD(UPDATE_PERIOD)
+)
+uut (
+    .clk(clk),
+    .rst(rst),
+
+    /*
+     * Increment inputs
+     */
+    .stat_inc(stat_inc),
+    .stat_valid(stat_valid),
+
+    /*
+     * Statistics increment output
+     */
+    .m_axis_stat(m_axis_stat),
+
+    /*
+     * Control inputs
+     */
+    .update(update)
+);
+
+endmodule
+
+`resetall