Efficient machine precision computation of impedance matrix elements in surface integral equation formulations