stream2 fails if coordinates are single precision?

Question

0 个投票

I've encountered a strange behavior with stream2 which fails silently if any input coordinates are single precision. It looks like the velocity data can be single precision, but the grid coordinates and the query coordinates must both be double precision.

Here's an MWE:

x = 1:100; 
y = 1:100; 
vx = rand(100); 
vy = rand(100); 
% This works: 
st = cell2mat(stream2(x,y,vx,vy,50,50));
% This works: 
st = cell2mat(stream2(x,y,single(vx),single(vy),50,50));
% This fails (empty output): 
st = cell2mat(stream2(single(x),single(y),vx,vy,50,50));
% This fails: 
st = cell2mat(stream2(x,y,vx,vy,single(50),single(50)));
% This fails: 
st = cell2mat(stream2(single(x),single(y),single(vx),single(vy),50,50));
% This fails: 
st = cell2mat(stream2(single(x),single(y),single(vx),single(vy),single(50),single(50)));

Is there something I'm missing? Is there a way to use stream2 with single-precision coordinates? Shouldn't the function throw an error if inputs are single, instead of producing an empty output?