Parallel cost-aware optimization of multidimensional black-box functions