Log in

My Journal

Procrastination, yup, that's right

Bah!

I finally get my code working, and then run a diff between the output of a known "working" version, and the output of my code. I get about 30 lines that show up as different.

< 426 8.05078e-09---> 426 8.05077e-09

Now, the trouble is that I'm working with floats. Given that up until one of the very last iterations all the variables match, I've a feeling that the difference may be that of the rounding errors from a straight sum of a matrix, as apposed to splitting the matrix, then summing the sums of the split matrices. But I have no way of knowing :(

I think they are rounding errors, especially since it's the least sig fig that's varying. Amusingly my code may actually be more accurate, since it adds small numbers to small numbers then adds the resultant bigger numbers, what an amusing side effect of parallelisation :)