Add decent petsc formula