Added again the newMatrix (no race conditions now), added memory bandwith calculation (approx).
Changed the computation of the maximum absolute difference outside the kernel, using Thrust.
Added a kernel to  compute the analytical solution (not only the difference).
Functions for calculating index inlined, and remove unuseful indexes
Analytical computation removed (temporarily, doesn't work in release). Prepared for optimisation.