def METHOD_NAME(benchmark, nqubits):
    benchmark.group = "SparseMatrix4Q"
    sparse_matrix = csr_matrix(([1], ([0], [0])), shape=(16, 16), dtype=complex)
    bench_gate(benchmark, nqubits, SparseMatrix([0, 1, 2, 3], sparse_matrix))