Full wave electromagnetic optimizations using surface integral equations and the multilevel fast multipole algorithm