The harmonic mean function was added (to the `develop` branch) by merge commit 8b55806. There is also a weighted harmonic mean -- see [Wikipedia](https://en.wikipedia.org/wiki/Harmonic_mean#Weighted_harmonic_mean), which should also be implemented.