High-accuracy Optimization by Parallel Iterative Discrete Approximation and Multi-GPU Computing