On Thu, 8 Dec 2016, Kathrin Smetana wrote:
Unfortunately my matrix has a dimension of 1M \times 1M…
I was picturing smaller...
I would not try 1M dofs with a direct solver. Saving an incomplete
Cholesky factorization would probably still be a big win.
There's a LinearSolver::reuse_preconditioner() method which can tell
libMesh whether it should reuse or recompute a preconditioner on
subsequent solves, and the PetscLinearSolver respects that flag. I
don't think we have decent test coverage on this though so I'd not
trust it completely; run with all the PETSc info-printing options on
at first.
---
Roy
------------------------------------------------------------------------------
Developer Access Program for Intel Xeon Phi Processors
Access to Intel Xeon Phi processor-based developer platforms.
With one year of Intel Parallel Studio XE.
Training and support from Colfax.
Order your platform today.http://sdm.link/xeonphi
_______________________________________________
Libmesh-users mailing list
[email protected]
https://lists.sourceforge.net/lists/listinfo/libmesh-users
