@romainhild here is a snippet of code for feelpp/po#21
it is untested/uncompiled though  (i will try to do that soon)

you will retrieve a data structure that allows the compute the non-zero entries
of the matrix C for the change of basis for the computation of the eigenmodes of
the curl operator.