ENH Faster distance_multi

Profiling was showing operator + as a major contributor to the runtime,
so replace it by += in the inner loop.

Speeds up is about 40%.