[petsc-users] reuse ksp in time marching loop

Shriram Srinivasan shriram at ualberta.ca
Tue Jun 10 15:52:59 CDT 2014


Hi all,
I am implementing a time marching loop (My problem is not amenable to 
usage of the TS module) where the coefficient matrix A wont change with 
time. I was wondering, in the sequence of commands

/KSPCreate(PETSC_COMM_SELF, &ksp); //
//KSPSetOperators(ksp, A,  A, SAME_PRECONDITIONER); //
//KSPGetPC(ksp, &pc);//
//PCSetType(pc, PCLU); //
//PCFactorSetMatSolverPackage(pc, MATSOLVERMUMPS);//
//KSPSetFromOptions(ksp);//
//KSPSolve(ksp, b, x);//
//KSPDestroy(&ksp);/

would it be  efficient/possible to move some of them out of the loop.
If so, which ones can I move out.

Thanks,
Shriram

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20140610/ce505781/attachment.html>


More information about the petsc-users mailing list