call VECCREATE(PETSC_COMM_WORLD, rhs, ierr)
call VECCREATE(PETSC_COMM_WORLD, solution, ierr)
call MATCREATEAIJ (PETSC_COMM_WORLD,
call KSPCREATE(PETSC_COMM_WORLD, ksp, ierr)