Residual edge error indicator


i’m trying to implement a weighted residual-based error estimator for a 2D Poisson equation. It consists of an inner and an edge part. Similar to the one Philip Lederer implemented in the showcase

My problem is, that the weight is given by (h_T/p_T)^2 for the elements and h_e/p_e for the individual edges.
The first part is no problem since Joachim added the GetOrder()-function to the python interface.

And furthermore how can i incorporate the p_e for an individual edge?

Best regards,