Maxwell solver example does not work after 0.5 update

dolfinx
2

Could you add what eigenvalues you expect to get printed?
When I run the code below in v0.4.1 i get:

eigenvalue: 2.467453845009
eigenvalue: 9.860312654523
eigenvalue: 9.870441650323
eigenvalue: 12.331992693807
eigenvalue: 19.747631128070
eigenvalue: 22.210852229644
eigenvalue: 32.109120054580
eigenvalue: 39.330212890661
eigenvalue: 39.491630029626
eigenvalue: 41.814442297205
eigenvalue: 49.267973337875
eigenvalue: 49.419419448067
eigenvalue: 61.692726728336
eigenvalue: 61.716767616632
eigenvalue: 71.681859411201
eigenvalue: 79.092163297278
eigenvalue: 88.079821259095

while if i run it with v0.5.0 I get

eigenvalue: 2.467453845010
eigenvalue: 9.860312654523
eigenvalue: 9.870441650323
eigenvalue: 12.331992693807
eigenvalue: 19.747631128070
eigenvalue: 22.210852229644
eigenvalue: 32.109120054579
eigenvalue: 39.330212890658
eigenvalue: 39.491630029626
eigenvalue: 41.814442297208
eigenvalue: 49.267973337872
eigenvalue: 49.419419448067
eigenvalue: 61.692726728332
eigenvalue: 61.716767616630
eigenvalue: 71.681859411199
eigenvalue: 79.092163297276

where the only difference (up to the 6-7 digit is that v0.4.1 finds one more eigenvalue).

Note that you could use

bc_facets = mesh.exterior_facet_indices(msh.topology)

in v0.5.0, producing the same result as v0.5.0 above.

Please note that you should never use XDMFFile for post-processing Nedelec-kind elements (as it interpolates into CG-1, which is not well defined).

You should use VTXWriter after interpolating the solution into DG-1.