stream2 fails if coordinates are single precision?

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?

Answer 1

[Javed Mohd]

I also faced the same problem but thanks to your querry, I am able to rectify by casting my grid coordinates into double.