Whats up guys,
This is my first post here in the Fenics Project community.
I am trying to implement an adaptive meshing scheme solver for the Maxwell eigenvalue problem.
I already have the eigenvalue solver by following chapter 34 from Fenics book and I would like to extend this to an adaptive mesh solver.
I am trying to implement a posteriori residual based error indicator, which is defined in the following paper: “A posteriori Error Estimates for Maxwell’s Eigenvalue problem”.
The error indicator defined in this paper considers taking the divergence and the curl*curl of the computed solutions.
I would like to know if there is any direct way to evaluate these expressions from the computed solutions at each cell or if I need to implement this residual estimator expression by using another variational form from which I could form a matrix containing the residual contribution from each cell.
The residual expression that I want to use is the following:
I hope this is clear enough guys, thanks for your comments and help!