def METHOD_NAME():
    # if qset.num_cpus == 1:
    #    return qset.openmp_thresh
    jc_dims = np.unique(np.logspace(0.45, 1.78, 20, dtype=int))
    jc_result = system_bench(_jc_liouvillian, jc_dims)
    opto_dims = np.unique(np.logspace(0.4, 1.33, 12, dtype=int))
    opto_result = system_bench(_opto_liouvillian, opto_dims)
    spin_dims = np.unique(np.logspace(0.45, 1.17, 10, dtype=int))
    spin_result = system_bench(_spin_hamiltonian, spin_dims)
    # Double result to be conservative
    thresh = 2*int(max([jc_result, opto_result, spin_result]))
    if thresh < 0:
        thresh = np.iinfo(np.int32).max
    return thresh