To reduce the human intervention in the gate design process, a new optimization program is designed to place multiple molding gates complex compoenent surface while minimizing the shrinkage and warpage. The main procedure of the analysis include: 1) using Poisson disk sampling (PDS) to preselect candidate gate locations automatically within the suitable gating region specified by designers; 2) using a space-filling initialization strategy and efficient global optimization to find the optimal gate locations. In the present setting, the molding gate design problem is formalized as a discrete optimization problem. The PDS is employed to construct the discrete solution space and EGO is used to efficiently search through a large solution space for the best design.