Why does my subtraction yield different results?
2 views (last 30 days)
Show older comments
I have been working on an analysis and after simplifying my code for speed the results came out differently even though I think I should have obtained the same results.
L = eye(n * c) / (eye(n * c) - A);
x_test = L * F * ones(p , 1);
isequal(eps(x),eps(x_test))
% so x and x_test are equal
x_1 = L * F * ones(p , 1) - (eye(n * c) - A_tilde) \ F_tilde * ones(p , 1);
x_2 = x - (eye(n * c) - A_tilde) \ F_tilde * ones(p , 1);
isequal(eps(x_1),eps(x_2))
The first isequal return me a one while the second isequal returns me a zero.
This means that x_test is equal to x but if I subtract the same vector on the RHS of the minus sign, I get a different result. I am not sure how this is possible. Could anyone point me in the right direction?
0 Comments
Accepted Answer
More Answers (1)
William Rose
on 28 Aug 2023
Can you provide more of your code? Why do you use the eps() function when testing the equality of two numbers or vectors or matrices? The following example shows that two unequal numbers can have equal eps's.
x1=1; x2=x1+1e-8; isequal(x1,x2), isequal(eps(x1),eps(x2))
0 Comments
See Also
Products
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!