Fast integral equation method for metasurface antennas