Residual edge error indicator

1 week 3 days ago - 5 days 21 hours ago #1327 by THaubold

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.
I tackle the edge part by
+ Defining a bilinear form with the .other()-function
+ Using Apply
+ Integrate element_wise

Is this the correct way of doing this? (Without using the lifting approach)
And furthermore how can i incorporate the p_e for an individual edge?

Best regards,

Please Log in or Create an account to join the conversation.

© 2018 Netgen/NGSolve