Computing star discrepancies with numerical black-box optimization algorithms