let numpy handle the broadcasting in Coordinates

[ahms5]

same problem as in #664, but in the coordiantes class

we could think to shift the match_shape into the utils, so that we can use it anywhere