Hi all, 

Thank you all for all the helpful comments and suggestions! 

I will let you know how it went.

Thanks again,
Kathrin 

> Am 08.12.2016 um 14:29 schrieb Roy Stogner <[email protected]>:
> 
> 
> 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

Reply via email to