Virtual Element approximation of 2D magnetostatic problems

We consider the use of nodal and edge Virtual Element spaces for the discretization of magnetostatic problems in two dimensions, following the variational formulation of Kikuchi. In addition, we present a novel Serendipity variant of the same spaces that allow to save many internal degrees of freedom. These Virtual Element Spaces of different type can be useful in applications where an exact sequence is particularly convenient, together with commuting-diagram interpolation operators, as is definitely the case in electromagnetic problems. We prove stability and optimal error estimates, and we check the performance with some academic numerical experiments.