A mixed face-edge finite element formulation for 3D magnetostatic problems