`axis` option for `stats.tmean` do not appear to be working as expected

My issue is about (see title)

Reproducing code example:

import scipy as sp
import numpy as np

r = np.array([[1, 1, 1], [2, 2, 2], [3, 3, 3]])

print(sp.stats.tmean(r, axis=1))
print(sp.stats.tstd(r, axis=1))

Output:

2.0
[0. 0. 0.]

Expectation:

[1.0, 2.0, 3.0]
[0., 0., 0.]

Scipy/Numpy/Python version information:

1.4.1 1.18.1 sys.version_info(major=3, minor=7, micro=6, releaselevel='final', serial=0)

[rgommers]

Thanks @siddhantwahal, please go for it!

[dburkhardt]

This is still a bug that makes tmean unusable for intended applications.