The prepare activation map function in unemap now uses vectorised
equations to carry out Hardy Interpolation and delaunay triangulation to
mesh non-rectangular domains.