Abstract
We prove a discrete version of the first Weber inequality on three-dimensional hybrid spaces spanned by vectors of polynomials attached to the elements and faces of a polyhedral mesh. We then introduce two Hybrid High-Order methods for the approximation of the magnetostatics model, in both its (first-order) field and (second-order) vector potential formulations. These methods are applicable on general polyhedral meshes, and allow for arbitrary orders of approximation. Leveraging the previously established discrete Weber inequality, we perform a comprehensive analysis of the two methods. We finally validate them on a set of test-cases.