pdkMeans function

K-means clustering for HPD matrices