Bug in single/double operations?

7 次查看（过去 30 天）
David Wilson
David Wilson 2011-6-27
Can anybody explain this behaviour?
If you run the following code, it behaves as expected. You finish up with approximately eps('single')=1e-7.
HOWEVER, if you comment out the disp(x) line, then you finish up with eps('double')=1e-16.
Why should a disp(x) change anything? Is it a sneaky invisible casting?
------------------------------------
format long
x = single(1) % single precision
while single(1+x) > 1 % keep going until no difference
x= x/2; % if you can tell any difference, half it.
disp(x) % comment this line out
end % while
---------------------------------------

请先登录，再回答此问题。

采纳的回答

James Tursa
James Tursa 2011-6-28
This exact issue has already been discussed on the newsgroup, with the JIT being the likely cause of the difference:
http://www.mathworks.com/matlabcentral/newsreader/view_thread/291577

更多回答（1 个）

Mike Hosea
Mike Hosea 2011-6-27
It appears to be something to do with the MATLAB JIT-Accelerator. I'm going to say it's probably unintentional. :)
I have to add, however, that while one can usually get away with an algorithm like that in MATLAB (except for code generation), we had to give up on it in numerical software in general after Intel came out with 80-bit extended-precision floating-point registers. Usually you want the eps for the data type, not for the register.

类别

在 Help Center 和 File Exchange 中查找有关 Creating and Concatenating Matrices 的更多信息

标签

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by