@skipif(not hasattr(np, 'float96'), 'cannot find float96 so skipping')
    def test_1D_float96(self):
        a = ma.array((1,2,3,4), mask=(0,0,0,1))
        actual_dt = mstats.gmean(a, dtype=np.float96)
        desired_dt = np.power(1 * 2 * 3, 1. / 3.).astype(np.float96)
        assert_almost_equal(actual_dt, desired_dt, decimal=14)
        assert_(actual_dt.dtype == desired_dt.dtype)
dataset/ETHPy150Open scipy/scipy/scipy/stats/tests/test_mstats_basic.py/TestGMean.test_1D_float96/original