Adaptive Wilson element with at most one hangnode on one edge

[XuemeiZhou-1993]

Professor Chen,

Sorry to disturb you. There's some bug in my codes and it took a lot of effort and time to get nothing. I was coding the adaptive Wilson element to solve the Poisson equation. Since it is the rectangle element, I can rarely refer to you ifem package. I write the square bisect codes (bisectQuad.m). In the solving step (PoissonWilson.m), things go right when dealing with the uniform refinement (compSolFig_nohang.m). While with hangnodes, it goes wrong (compSolFig_hang.m).
I am hoping if you can have a look on my codes when you are spare. Thanks a lot!
The theoretical results refer to the following pdf.
Have a nice day!

Xuemei ZHou
A Ph.D. student of Chinese Academy of Sciences

hu2014 New a posteriori error estimate and quasi-optimal convergence of the adaptive nonconforming Wilson element.pdf

WilsonSolve.zip

One data m file is missing, which is just the Lshapedata.m in your ifem package.

CareyData.txt

[scaomath]

You have to recursively use adjacency matrix to find hanging nodes (assuming more than 1 level of irregular mesh). I have a quad sect code in iFEM data structure for bilinear element supporting multiple levels of hanging node. I will find time to clean up and encapsulate it and upload to iFEM.