Implementation of the algorithm propused in "A PDE based approach to multi-domain partitioning and quadrilateral meshing" (see here). Although the implementation is general it has only been tested on a semicircle. The inputs are the radius of the circle, the center and the number of subdivisions of the centre square.
The following is the main result of the implementation. The resulting mesh can be used for any purpose.
I think the main contribution of the project is its pedagogical value. Through the code you can inspect all the main parts of the complex algorithm. These parts are shown below.
