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 :(
Will speak to my project supervisor on Monday...
< 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 :(
Will speak to my project supervisor on Monday...
Oh, Entropy, I've not been for aaages... Quite possibly.
There's a paper 'what every computer scientist should know about floating point' which you should probably read.