[Devid2014]

Using something else as float or double would decrease performance significantly. Using array sorting would not only be solve too but actually do not improve precision that much at all.

Using double to sum floats has comparable precession as Kahan summation. But of course it is not possible to do the same to summing doubles. Or if there is no proper double precision available at all like on GPU.