actually, we have to be super careful about floating point numbers here.
There is a recently added method ```CMath::fequals``` that can deal with all cases. Pls use that (or modify it to work), its already unit tested, so your unit test really only need to deal with matrix positions differing, the actual difference (NaN, inf, etc) do not need to be checked again.