Statistics for Fast and accurate analysis of optical metamaterials using surface integral equations and the parallel multilevel fast multipole algorithm