Hi Denis, Thank you for your information.
Best, Shanglong 在 2016年7月25日星期一 UTC-4下午12:02:09，Denis Davydov写道： > > 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 dealii+unsubscr...@googlegroups.com. For more options, visit https://groups.google.com/d/optout.
