[aspect-devel] Using MUMPS through Trillions instead of PETSc

cedric thieulot c.thieulot at uu.nl
Sat Feb 25 10:57:30 PST 2017


Hi, 

out of curiosity, for our type of applications/matrices, which solver (SuperLU or MUMPS) do you think would be better/faster/less memory consuming/etc … ? 

Ce/




> On 22 Feb 2017, at 15:08, Timo Heister <heister at clemson.edu> wrote:
> 
>> I saw deal.II has the capability to use MUMPS through Trillions through the
>> SolverDirect class:
>> https://urldefense.proofpoint.com/v2/url?u=https-3A__www.dealii.org_8.4.1_doxygen_deal.II_classTrilinosWrappers-5F1-5F1SolverDirect.html&d=DwIFaQ&c=Ngd-ta5yRYsqeUsEDgxhcqsYYY1Xs5ogLxWPA_2Wlc4&r=c08Btfq4m9QEScXN3ZQwLZzzWQE7S8CYq1IYuzKV_Zk&m=WqHvTxGttdmdeKTLKkUCQDnDtDhtDP5fdaIgVKdwdZE&s=BoLSc86Lhci2Z6jfiArFykODyz62pKP7QBinfcqErOs&e= 
> 
> We currently don't expose this, but if you go into solver.cc you will
> find the line
> 
>        TrilinosWrappers::SolverDirect solver(cn);
> 
> If you change that to
> 
> TrilinosWrappers::SolverDirect::AdditionalData data(false, "Amesos_Mumps");
> TrilinosWrappers::SolverDirect solver(cn, data);
> 
> and you configure Trilinos with mumps (not enabled by default), you
> should be good to go.
> 
>> However, the few discussions I’ve seen regarding using MUMPS through
>> Trillions are not entirely encouraging :)
> 
> The first link you posted is about using MUMPS directly in deal.II
> (which we removed) and has nothing to do with using it through
> Trilinos. But yes, MUMPS is a nightmare to install/use.
> 
>> In reality, my interest is not explicitly with MUMPS but rather a parallel
>> direct solver option (MUMPS, PARDISO, etc).
> 
> You might try superlu (also through Trilinos). But let me ask, what is
> the reason you are looking into direct solvers? While they might help
> for very ill-conditioned problems, I doubt that you will have good
> results on large 3d problems.
> 
> -- 
> Timo Heister
> http://www.math.clemson.edu/~heister/
> _______________________________________________
> Aspect-devel mailing list
> Aspect-devel at geodynamics.org
> http://lists.geodynamics.org/cgi-bin/mailman/listinfo/aspect-devel



More information about the Aspect-devel mailing list