A PETSc-Based Parallel Implementation of Finite Element Method for Elasticity Problems

Algorithm 3

call KSPCreate(PETSC_COMM_WORLD,ksp,ierr)

call KSPSetOperators(ksp,AK,AK,DIFFERENT_NONZERO_PATTERN,ierr)

call KSPSetType(Ksp,KSPCG,ierr)

call KSPCGSetType(Ksp,KSP_CG_SYMMETRIC,ierr)

call KSPGetPC(ksp,pc,ierr)

call PCSetType(pc,PCJACOBI,ierr)

tol = 1.0d-7

call KSPSetTolerances(ksp,tol,PETSC_DEFAULT_DOUBLE_PRECISION,   &

  PETSC_DEFAULT_DOUBLE_PRECISION,PETSC_DEFAULT_INTEGER,ierr)

call KSPSolve(ksp,F,x,ierr)