Now it actually works pretty well. The time calculation was changed from
clock() to omp's clock function. This makes the time accurate on
multiple CPUs.

In this code 5/6 of the time is spent on added the already calculated
values into two global variables. This is far from optimal. We should do
something about that.