Reimplement the insertion of boundary lines in 3d grids. The old

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]]