501   KSP                ksp;  in DMSwarmProjectFromSwarmToCells()  local
527   // Setup KSP  in DMSwarmProjectFromSwarmToCells()
531     PetscCall(KSPCreate(comm, &ksp));  in DMSwarmProjectFromSwarmToCells()
532     PetscCall(KSPGetPC(ksp, &pc));  in DMSwarmProjectFromSwarmToCells()
535     PetscCall(KSPSetType(ksp, KSPCG));  in DMSwarmProjectFromSwarmToCells()
536     PetscCall(KSPSetNormType(ksp, KSP_NORM_NATURAL));  in DMSwarmProjectFromSwarmToCells()
537     PetscCall(KSPSetTolerances(ksp, 1e-10, PETSC_DEFAULT, PETSC_DEFAULT, PETSC_DEFAULT));  in DMSwarmProjectFromSwarmToCells()
538     PetscCall(KSPSetOperators(ksp, M, M));  in DMSwarmProjectFromSwarmToCells()
539     PetscCall(KSPSetFromOptions(ksp));  in DMSwarmProjectFromSwarmToCells()
540     PetscCall(PetscObjectSetName((PetscObject)ksp, "Swarm-to-Mesh Projection"));  in DMSwarmProjectFromSwarmToCells()
541     PetscCall(KSPViewFromOptions(ksp, NULL, "-ksp_projection_view"));  in DMSwarmProjectFromSwarmToCells()
549   PetscCall(KSPSolve(ksp, B_mesh, U_mesh));  in DMSwarmProjectFromSwarmToCells()
551   // KSP summary  in DMSwarmProjectFromSwarmToCells()
557     PetscCall(KSPGetType(ksp, &ksp_type));  in DMSwarmProjectFromSwarmToCells()
558     PetscCall(KSPGetConvergedReason(ksp, &reason));  in DMSwarmProjectFromSwarmToCells()
559     PetscCall(KSPGetIterationNumber(ksp, &its));  in DMSwarmProjectFromSwarmToCells()
560     PetscCall(KSPGetResidualNorm(ksp, &rnorm));  in DMSwarmProjectFromSwarmToCells()
564                             "Swarm-to-Mesh Projection KSP Solve:\n"  in DMSwarmProjectFromSwarmToCells()
565                             "  KSP type: %s\n"  in DMSwarmProjectFromSwarmToCells()
566                             "  KSP convergence: %s\n"  in DMSwarmProjectFromSwarmToCells()
567                             "  Total KSP iterations: %" PetscInt_FMT "\n"  in DMSwarmProjectFromSwarmToCells()
574   PetscCall(KSPViewFromOptions(ksp, NULL, "-ksp_view"));  in DMSwarmProjectFromSwarmToCells()
579   PetscCall(KSPDestroy(&ksp));  in DMSwarmProjectFromSwarmToCells()