Bulk xy-coordinate calculation with np.matmul

[byewokko]

The vertices' coordinates are calculated one-by-one, which is very slow. It could be sped up by stacking the vertices' multigrid coordinates into a single 2D array and then performing the matrix multiplication.