Hi, looks like you are using the wrong constructor of yaspgrid and therefore using 4 different sequential grids. Try YaspGrid(MPI_COMM_WOLRD, length, size, periodic,overlap) instead of YaspGrid(length, size, periodic,overlap) Markus