The model order reduction techniques based on the multipoint projection Krylov methods have become the methods of choice to generate large macromodels of multiport RLC circuits. A well-known difficulty with such methods lies in the need for clever point selection to attain model compactness and accuracy. In this paper, we present an automatic methodology for optimizing sample point selection. This method is general, suitable for circuits obtained by the partial element equivalent circuit method coupled with an adaptive multilevel fast multipole method. Our algorithm has been validated on an industrial example to demonstrate the accuracy and robustness of the approach. Index Terms-Adaptive multilevel fast multipole method (AMLFMM), model order reduction (MOR), moment-matching method, partial element equivalent circuit (PEEC method), RLC. Trung-Son Nguyen received the M.Sc. and Ph.D. degrees in electrical engineering from Grenoble Institute of Technology, Grenoble, France, in 2009 and 2012, respectively. He held Postdoctoral position with Grenoble Electrical Engineering Laboratory from 2012 to 2013. In 2015, he moved to Vietnam Fire and Rescue Police Department, Vietnam, where he is currently a Researcher of Science-Technology Research & Apparatus and Equipment Inspection Division. His current research interests include model order reduction, circuit realization method, and development of numerical tools for modeling of electromagnetic phenomena (integral methods).