Add a time macro

Building the mesh currently takes a long time (12-14s for 32^3) because
of duplicate edge creation. If edges are only built once and then
maintained, this may not be a concern. Still, it may be better to go
back to a scheme which allowed for only unique edge creation.
Write and test the tree_size_atlevel function