In order to better control the parallelization, we add a grain size
parameter which controls the maximal size of a block that is treated
by a thread.