This paper proposes an acceleration of the multi-level fast multipole algorithm(MLFMA) by using a double interpolation method. The MLFMA has been primarily used to conduct scattering analysis of electrically large targets, e.g. stealth aircraft. In the MLFMA, radiation functions of each basis functions are first precomputed, and then aggregated. After transfer calculations for the aggregations, each interaction is disaggregated, and then received in the testing function. The key idea of the proposed method is to decrease the sampling rates of the radiation and receiving functions. The computational complexity of the unit sphere integration in terms of the testing functions is thus highly alleviated. The remaining insufficient sampling rate is then complemented by using additional interpolation. We demonstrate the performance of the proposed method through radar cross-section(RCS) calculations for realistic aircraft.