Hello everybody,
I still try to use a error estimator as presented by Bürg for an hp-refinement method for a Maxwell problem.
Therefore, I need to compute the term
[tex]
\begin{align*}
\Vert \nabla \times \left( \mu^{-1} \nabla \times u_h\right) - \varepsilon \omega^2 u_h \Vert_{L^2 (K)}^2
\end{align*}
[/tex]
on every cell K. I tried to set a new HCurl gridfunction for the curl of u_h, but since curl(u_h) returns the scalar curl
[tex]
\begin{align*}
\operatorname{curl}(u) = \left( \frac{\partial u_2}{\partial x_1} - \frac{\partial u_1}{\partial x_2} \right)
\end{align*}
[/tex]
this does not work (different dimensions).
Has anybody a suggestion how fix this issue?
Best regards,
Phil