A Boundary Integral Equation Scheme for Simulating the Nonlocal Hydrodynamic Response of Metallic Antennas at Deep-Nanometer Scales