algorithm was quadratic in the number of vertices, while the new
one is linear.  Not surprisingly, this speeds up reading large grids
a lot.  It appears to have been the last remaining quadratic algorithm
used when setting up a std_domain, and reading of very large grids
is now fast.

[[Imported from SVN: r8362]]