Fast and accurate analysis of large metamaterial structures using the multilevel fast multipole algorithm