Repository navigation
sc.midpoints does not operate on vector-valued data #3751
Description
sc.midpoint does not operate on vector-valued data, but a simple alternative works fine (probably missing the corner cases mentioned in the sc.midpoint doc string)
import scipp as sc
dim='x'
vecs = sc.vectors(values=[[0, 0, 0], [1, 1, 1], [2, 2, 2]], dims=[dim], unit='m')
midpoint = sc.vectors(values=[[0.5, 0.5, 0.5], [1.5, 1.5, 1.5]], dims=[dim], unit='m')
try:
calculated = sc.midpoints(vecs, dim)
except sc.DTypeError as er:
print(f'sc.midpoints failed with error `{er}`')
calculated = (vecs[dim,:-1] + vecs[dim,1:]) / 2
if not sc.allclose(midpoint, calculated):
print('Failed calculating the midpoint')raises an exception in the try statement which causes the following output on my machine:
sc.midpoints failed with error `'midpoints' does not support dtypes 'vector3', 'vector3', `Metadata
Metadata
Assignees
Labels
No labels
Activity