Now the program uses a 2-d torus blocked decomposition.