[yt-dev] Issue #1175: fcoords and tcoords not the same shape for rays from particle datasets (yt_analysis/yt)

Britton Smith issues-reply at bitbucket.org
Thu Feb 25 06:32:33 PST 2016


New issue 1175: fcoords and tcoords not the same shape for rays from particle datasets
https://bitbucket.org/yt_analysis/yt/issues/1175/fcoords-and-tcoords-not-the-same-shape-for

Britton Smith:

This script: http://paste.yt-project.org/show/6268/

gives the following error:

```
#!python

  File "test2.py", line 12, in <module>
    print ray["x"]
  File "/Users/britton/Documents/work/yt/yt-hg/yt/data_objects/data_containers.py", line 268, in __getitem__
    self.get_data(f)
  File "/Users/britton/Documents/work/yt/yt-hg/yt/data_objects/data_containers.py", line 1181, in get_data
    self._generate_fields(fields_to_generate)
  File "/Users/britton/Documents/work/yt/yt-hg/yt/data_objects/data_containers.py", line 1201, in _generate_fields
    fd = self._generate_field(field)
  File "/Users/britton/Documents/work/yt/yt-hg/yt/data_objects/data_containers.py", line 305, in _generate_field
    tr = self._generate_fluid_field(field)
  File "/Users/britton/Documents/work/yt/yt-hg/yt/data_objects/data_containers.py", line 325, in _generate_fluid_field
    rv = finfo(gen_obj)
  File "/Users/britton/Documents/work/yt/yt-hg/yt/fields/derived_field.py", line 184, in __call__
    dd = self._function(self, data)
  File "/Users/britton/Documents/work/yt/yt-hg/yt/geometry/coordinates/coordinate_handler.py", line 40, in _coords
    rv = data.ds.arr(data.fcoords[...,axi].copy(), units)
  File "/Users/britton/Documents/work/yt/yt-hg/yt/data_objects/data_containers.py", line 1282, in fcoords
    return self._current_chunk.fcoords
  File "/Users/britton/Documents/work/yt/yt-hg/yt/geometry/geometry_handler.py", line 271, in cached_func
    tr = func(self)
  File "/Users/britton/Documents/work/yt/yt-hg/yt/geometry/geometry_handler.py", line 322, in fcoords
    ci[ind:ind+c.shape[0], :] = c
ValueError: could not broadcast input array from shape (81,3) into shape (40,3)
```

This seems to be because the fcoords and tcoords for the ray are coming back with different shapes.





More information about the yt-dev mailing list