Hello, I want to implement general procedures to generate discretized differential operators (like grad,div,curl,laplacian) using PETSc and staggered structured grids. Before I start, I am wondering if there is something similar in PETSc already? Regards, Alexander