We present a class of non-standard numerical schemes which are modifications of the discrete gradient method. They preserve the energy integral exactly (up to the round-off error). The considered class contains locally exact discrete gradient schemes and integrators of arbitrary high order. In numerical experiments we compare our integrators with some other numerical schemes, including the stan...

Approximate, non-singular kernels are often used in moment-method formulations coping with thin-wire structures. Their use has important consequences, one of which is the appearance of oscillations in the computed currents when the number of sub-domain basis functions is sufficiently large. These oscillations are not due to round-off errors. In this paper, a smoothing procedure is used in conju...