[cig-commits] [commit] rajesh-petsc-schur: minor changes to get solve_Ahat_p_fhat_PETSc_Schur to work (b684666)
cig_noreply at geodynamics.org
cig_noreply at geodynamics.org
Wed Nov 5 19:04:39 PST 2014
Repository : https://github.com/geodynamics/citcoms
On branch : rajesh-petsc-schur
Link : https://github.com/geodynamics/citcoms/compare/464e1b32299b15819f93efd98d969cddb84dfe51...f97ae655a50bdbd6dac1923a3471ee4dae178fbd
>---------------------------------------------------------------
commit b6846661681f8ccc79943d187e2f57712910bb4b
Author: Rajesh Kommu <rajesh.kommu at gmail.com>
Date: Thu Sep 11 16:13:37 2014 -0700
minor changes to get solve_Ahat_p_fhat_PETSc_Schur to work
>---------------------------------------------------------------
b6846661681f8ccc79943d187e2f57712910bb4b
lib/Stokes_flow_Incomp.c | 8 ++++----
1 file changed, 4 insertions(+), 4 deletions(-)
diff --git a/lib/Stokes_flow_Incomp.c b/lib/Stokes_flow_Incomp.c
index 01f8763..cb6cfe2 100644
--- a/lib/Stokes_flow_Incomp.c
+++ b/lib/Stokes_flow_Incomp.c
@@ -289,10 +289,10 @@ static PetscErrorCode solve_Ahat_p_fhat_PETSc_Schur(struct All_variables *E,
/*-------------------------------------------*/
ierr = KSPCreate(PETSC_COMM_WORLD, &S_ksp); CHKERRQ(ierr);
ierr = KSPSetOperators(S_ksp, S, S); CHKERRQ(ierr);
- ierr = KSPGetPC(S_ksp, &S_pc); CHKERRQ(ierr);
+ //ierr = KSPGetPC(S_ksp, &S_pc); CHKERRQ(ierr);
ierr = KSPSetType(S_ksp, "cg"); CHKERRQ(ierr);
- ierr = PCSetType(S_pc, "none"); CHKERRQ(ierr);
- ierr = KSPSetInitialGuessNonzero(S_ksp, PETSC_TRUE); CHKERRQ(ierr);
+ //ierr = PCSetType(S_pc, "none"); CHKERRQ(ierr);
+ //ierr = KSPSetInitialGuessNonzero(S_ksp, PETSC_FALSE); CHKERRQ(ierr);
/*--------------------*/
/* Solve for pressure */
@@ -306,7 +306,7 @@ static PetscErrorCode solve_Ahat_p_fhat_PETSc_Schur(struct All_variables *E,
ierr = MatMult(E->G, PVec, fstar); CHKERRQ(ierr);
ierr = VecAYPX(fstar, -1.0, FF); CHKERRQ(ierr);
ierr = MatSchurComplementGetKSP(S, &inner_ksp); CHKERRQ(ierr);
- ierr = KSPSetInitialGuessNonzero(inner_ksp, PETSC_TRUE); CHKERRQ(ierr);
+ //ierr = KSPSetInitialGuessNonzero(inner_ksp, PETSC_FALSE); CHKERRQ(ierr);
ierr = KSPSolve(inner_ksp, fstar, VVec); CHKERRQ(ierr);
//strip_bcs_from_residual_PETSc( E, VVec, lev );
More information about the CIG-COMMITS
mailing list