I would suggest you have a look at: https://github.com/jpdean/maxwell/blob/master/solver.py
The Nedelec elements in DOLFINx is properly implemented wrt interpolation, meaning that the functionals described the basis functions are integral moments, not point evaluations, as in the case of legacy dolfin. The effect of this is optimal interpolation properties.
You can Also see: Magnetostatics 3D - wrong 3D magnetic vector potential - #5 by ftrillaudp