In this paper, we deduced an iteration formula for the computation of central composite discrepancy. By using the iteration formula, the computational complexity of uniform design construction in flexible region can be greatly reduced. And we also made a refinement to threshold accepting algorithm to accelerate the algorithm’s convergence rate. Examples show that the refined algorithm can converge to the lower discrepancy design more stably. © 2020, T...