Hei, is there a function to square/cube etc. each element of a matrix or vector in PETSc? It is rather useful for things like intensity calculations (I = |E|^2), and should work also for distributed matrices. Thanks! Regards, Roland Richter