Keyword Arguments in Calculating Distance Matrix

One of the challenges I had when writing up my RREML package in R was how to handle the coordinates as a keyword argument to my rreml function.  Well, I've pretty much abandoned my R package at this point in time, and would eventually like to implement the approach in Python.  I just stumbled across this very brief Reddit discussion:

Euclidian distance calculation with 2 matrices

And I found the scipy documentation referred to in the response to be a great reference for how to setup my new keyword argument.

scipy.spatial.distance.cdist