The adaptive wavelet packet transform is applied to sparsify moment matrices for the fast solution of electromagnetic integral equations. An additive cost function is employed to adaptively select the optimal wavelet expansion/testing functions. It is found that the sparsified matrix has above-threshold elements that grow as O(N1.4) for typical scatterers.