Programming with PETSc/TAOΒΆ