Anisotropic Mesh Adaptation

I apply local mesh modification (such as edge swapping, node removal, refinement and node movement), in the context of anisotropic mesh adaptation.

It is possible to apply colouring for a thread parallel algorithm or just vectorized code suitable for an APL language style implementation, i.e. all the 0 edges in the below image can be flipped without any conflicts.

Below you can see it applied in two and three dimensions for a random initial mesh that is first modified to become isotropic and then to having a shock.

You can check out my proptotype MATLAB implementation or Gerard Gorman's C++ PRAgMaTIc library.