Dear PETSc-Users, I would like to use the inverted block diagonals of a a matrix. I have seen the function MatInvertBlockDiagonal() but I don’t know how to create a matrix out of them or an array of block matrizes. Does anyone have an example on how to use **values to create a PETSc matrix? Thanks Cyrill