New issue 845: Subclass of YTArray throwing errors https://bitbucket.org/yt_analysis/yt/issue/845/subclass-of-ytarray-throwing-...

Matthew Turk:

This script works for float64, but fails for float32.

``` #!python from yt.units import cm from yt.data_objects.octree_subset import YTPositionArray import numpy as np

a = np.random.random(100).astype("float32") * cm a *= np.finfo(a.dtype).eps

b = YTPositionArray(np.random.random((100, 3)).astype("float32"), input_units="cm") print type(b) b *= np.finfo(b.dtype).eps ```

It's also at paste 4674.

Responsible: ngoldbaum