media armónica en python

La función media armónica en Python (scipy.stats.hmean) requiere que la entrada sea números positivos.

Por ejemplo:

from scipy import stats
print stats.hmean([ -50.2 , 100.5 ])

resultados en:

ValueError: Harmonic mean only defined if all elements greater than zero

No entiendo matemáticamente por qué este debería ser el caso, excepto en el raro caso en que terminaría dividiendo por cero. En lugar de buscar una división por cero,hmean() luego arroja un error al ingresar cualquier número positivo, ya sea que se pueda encontrar una media armónica o no.

¿Me estoy perdiendo algo aquí en las matemáticas? ¿O es esto realmente una limitación enSciPy?

¿Cómo irías para encontrar la media armónica de un conjunto de números que podrían ser positivos o negativos en python?