<div dir="ltr"><div class="gmail_extra"><div class="gmail_quote">On Thu, Jul 10, 2014 at 1:37 AM, Jin, Shuangshuang <span dir="ltr"><<a href="mailto:Shuangshuang.Jin@pnnl.gov" target="_blank">Shuangshuang.Jin@pnnl.gov</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">






<div>
<font face="Calibri"><span style="font-size:11pt">
<div>Hi, I’m using the KSP solver as following to solve a linear system Ax=b, where A is a 16767x16767 square matrix, b is a 16767 length vector. </div>
<div> </div>
<div>EXTERN_C_BEGIN</div>
<div>extern PetscErrorCode PCCreate_Jacobi(PC);</div>
<div>EXTERN_C_END</div>
<div> </div>
<div>static PetscErrorCode solvingAxb(Mat A, Vec b, PetscInt nbus, Vec &x, const int me)</div>
<div>{</div>
<div>  PetscErrorCode ierr;</div>
<div> </div>
<div>  KSP ksp; // linear solver context</div>
<div>  Vec u; // exact solution</div>
<div>  PetscInt its;</div>
<div>  PetscReal norm; // norm of solution error</div>
<div>  PetscLogDouble t1, t2;</div>
<div>  PetscViewer viewer;</div>
<div>  PC pc; // preconditioner context</div>
<div>  PetscInt Istart, Iend;</div>
<div> </div>
<div>  ierr = VecDuplicate(b, &u); CHKERRQ(ierr);</div>
<div>  ierr = VecDuplicate(b, &x); CHKERRQ(ierr);</div>
<div> </div>
<div>  ierr = KSPCreate(PETSC_COMM_WORLD, &ksp); CHKERRQ(ierr);</div>
<div> </div>
<div>  ierr = KSPSetOperators(ksp, A, A, DIFFERENT_NONZERO_PATTERN); CHKERRQ(ierr);</div>
<div> </div>
<div>  PCRegister("ourjacobi", PCCreate_Jacobi);</div>
<div> </div>
<div>  KSPGetPC(ksp, &pc);</div>
<div>  PCSetType(pc, "ourjacobi");</div>
<div> </div>
<div>  ierr = KSPSetFromOptions(ksp); CHKERRQ(ierr);</div>
<div> </div>
<div>  ierr = KSPSolve(ksp, b, x); CHKERRQ(ierr);</div>
<div> </div>
<div>  ierr = MatMult(A,x,u); CHKERRQ(ierr);</div>
<div>  ierr = VecAXPY(u, -1.0, b); CHKERRQ(ierr);</div>
<div>  ierr = VecNorm(u, NORM_2, &norm); CHKERRQ(ierr);</div>
<div>  ierr = KSPGetIterationNumber(ksp, &its); CHKERRQ(ierr);</div>
<div> </div>
<div>  ierr = VecDestroy(&b); CHKERRQ(ierr);</div>
<div>  ierr = VecDestroy(&u); CHKERRQ(ierr);</div>
<div> </div>
<div>  ierr = KSPDestroy(&ksp);CHKERRQ(ierr);</div>
<div> </div>
<div>  PetscFunctionReturn(0);</div>
<div>}</div>
<div><font face="Times New Roman"> </font></div>
<div>I would like to know how to set the solver options at runtime to make it run faster, such as ksp_type, pc_type, and etc? It takes very long time to solve the system if I use no options.</div></span></font></div></blockquote>
<div><br></div><div>The preconditioner always depends on the particular problem being solved. You will have to tell us more about it.</div><div><br></div><div>  Thanks,</div><div><br></div><div>   Matt</div><div> </div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<div><font face="Calibri"><span style="font-size:11pt">
<div>Thanks,</div>
<div>Shuangshuang</div>
<div><font face="Times New Roman"> </font></div>
</span></font>
</div>

</blockquote></div><br><br clear="all"><div><br></div>-- <br>What most experimenters take for granted before they begin their experiments is infinitely more interesting than any results to which their experiments lead.<br>
-- Norbert Wiener
</div></div>