Hi, I can't comment on the first part of your question, but as for the second part -- you are right, currently you can't reuse factorization. I you want it, it should be relatively easy to add. See Direct solver in TrilinosWrappers which was recently modified to support this: https://www.dealii.org/developer/doxygen/deal.II/classTrilinosWrappers_1_1SolverDirect.html#a6c0fecc921bc93f2c005e265cfbdd882
Regards, Denis. On Monday, July 25, 2016 at 5:52:06 PM UTC+2, shanglong zhang wrote: > > Dear all, > > I modify the step-17 so that it uses PETScWrappers::SparseDirectMUMPS as > solver. However, when I active the option of set_symmetric_mode, the > solution has discrepancy from CG or MUMPS without the symmetry option. (The > solutions from CG and MUMPS without symmetry option are identical though.) > Any suggestion? > > And it seems to me that PETScWrappers::SparseDirectMUMPS doesn't have the > functionality of reusing the factorization of system matrix for different > right-hand side vector. Do I miss anything? > > Thank you very much. > > Best, > Shanglong > -- The deal.II project is located at http://www.dealii.org/ For mailing list/forum options, see https://groups.google.com/d/forum/dealii?hl=en --- You received this message because you are subscribed to the Google Groups "deal.II User Group" group. To unsubscribe from this group and stop receiving emails from it, send an email to [email protected]. For more options, visit https://groups.google.com/d/optout.
