Use petsc.h