Add norm23 to 3d.[ch], to compute the sqwuared norm.

This is convenient when, for exmaple, you just care about relative
magnitudes, not actual values, and want to save computing unnecessary
(for that purpose) square roots.