Comparing large sparse matrices for differences close to machine precision
6 views (last 30 days)
Basically what I am doing right now is creating the same matrix k in two very different ways. As it turned out simply comparing them by the find command results in almost no matches as the single elements are somehow very slightly different. as in if I subtract one element from the other the result is either close to or at machine precision. I am now looking for a way to compare those two matrices in a way, so I get an Idea of how bad the differences are.
Walter Roberson on 21 Sep 2022
[r1, c1, s1] = find(k1);
[r2, c2, s2] = find(k2);
rc1 = [r1, c1];
rc2 = [r2, c2];
[~, ia, ic] = intersect(rc1, rc2, 'rows');
common1 = s1(ia); common2 = s2(ic);
At this point, common1 and common2 are corresponding non-zero values that exist in both arrays, at locations k1(r1(ia(J)), c1(ia(J))) and k2(r2(ic(J)), c2(ic(J))) .. well except you should expect them to be at the same location in both.
You know these will be non-zero, so you can calculate relative change however is appropriate, such as
common2./common1 - 1 %also known as (common2 - common1)/common1