Removed unnecessary quadrature rules
Added reduction operators for arrays
Replaced manual calls to MPI_AllReduce() by reduce(), added an additional reduction for gamma_old in PipeCG for reasons of stability [Meurant 1987]