Yes, it automatically uses the same factorization.
> On Aug 4, 2020, at 6:55 PM, Anthony Paul Haas <a...@email.arizona.edu> wrote: > > Hello, > > When using Mumps to solve a linear system of equations (see below), can I > reuse the factorization to solve for multiple RHS, ie, can I use KSPSolve > multiple times while only building a different RHS in between the calls to > KSPSolve? > > Thanks, > > Anthony > > call KSPSetType(self%ksp,KSPPREONLY,self%ierr_ps) > call KSPGetPC(self%ksp,pc,self%ierr_ps) > > > > call PCSetType(pc,PCLU,self%ierr_ps) > > > > call PCFactorSetMatSolverPackage(pc,MATSOLVERMUMPS,self%ierr_ps) > call PCFactorSetUpMatSolverPackage(pc,self%ierr_ps)