stream2 fails if coordinates are single precision?

Question

Chad Greene 2022-8-25

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1786600-stream2-fails-if-coordinates-are-single-precision

回答： Javed Mohd 2023-1-19

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?