On the implementation of a boundary element method for 3-D multiply connected EM field problems