Residual edge error indicator

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

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
ngsolve.org/showcases/ngsolve/9-implemen...-a-lifting-technique

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,
Tim
Attachments:

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

© 2018 Netgen/NGSolve