def METHOD_NAME(self):
    """pauli cob multi qubit test"""
    # Helpful prints for debugging commented out below.
    multis = [Y ^ X ^ I ^ I, I ^ Z ^ Y ^ X, X ^ Y ^ I ^ Z, I ^ I ^ I ^ X, X ^ X ^ X ^ X]
    for pauli, dest in itertools.product(multis, reversed(multis)):
        # print(pauli)
        # print(dest)
        converter = PauliBasisChange(destination_basis=dest)
        inst, dest = converter.get_cob_circuit(pauli.primitive)
        cob = converter.convert(pauli)
        # print(inst)
        # print(pauli.to_matrix())
        # print(np.round(inst.adjoint().to_matrix() @ cob.to_matrix()))
        np.testing.assert_array_almost_equal(
            pauli.to_matrix(), inst.adjoint().to_matrix() @ dest.to_matrix() @ inst.to_matrix()
        )
        np.testing.assert_array_almost_equal(pauli.to_matrix(), cob.to_matrix())
        np.testing.assert_array_almost_equal(
            inst.compose(pauli).compose(inst.adjoint()).to_matrix(), dest.to_matrix()
        )