[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