Bruno, Thank you very much. I would check the time for the solve() function with different right hand side. If the factorization only be done the first time it is called, the time for solving the next right hand side should be lesser.
I looked through the documentation in Deal.II before. I think it could be done, too. It would be easier if there is an example to use the PETSc parameters for MATSOLVERMUMPS package in deal.II. Currently, I don't know where to put it in my codes. Best regards, Kien Vào Th 5, 27 thg 6, 2019 vào lúc 21:42 Bruno Turcksin < [email protected]> đã viết: > Kien, > > On Wednesday, June 26, 2019 at 10:20:38 PM UTC-4, Phạm Ngọc Kiên wrote: >> >> Hi colleagues, >> I am trying to build codes using PETScWrappers::SparseDirectMUMPS >> solver. >> Step-62 shows a very good example for using it. >> However, I really want to know if we can choose the factorization method >> in MUMPS by deal.II codes. >> > According to the documentation (see here > <https://www.dealii.org/developer/doxygen/deal.II/classPETScWrappers_1_1SparseDirectMUMPS.html>), > you can. I don't use PETSc but I think this done using the command line > > >> >> The second question is: do we have a way to get the factorized system >> matrix when using this direct solver? >> Because I want to reuse the factorized system matrix to solve the problem >> again with other right hand side. >> > Just call solve() with a different right-hand-side. The factorization > should only be done the first time the function is called. > > Best, > > Bruno > > -- > 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]. > To view this discussion on the web visit > https://groups.google.com/d/msgid/dealii/24c660f9-1450-4c1f-b0f7-d9c0a629fde8%40googlegroups.com > <https://groups.google.com/d/msgid/dealii/24c660f9-1450-4c1f-b0f7-d9c0a629fde8%40googlegroups.com?utm_medium=email&utm_source=footer> > . > For more options, visit https://groups.google.com/d/optout. > -- 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]. To view this discussion on the web visit https://groups.google.com/d/msgid/dealii/CAAo%2BSZep%2BHFiY_5St3U-MV%3DRVYK741xJ0QwJWOzG5xy3d-iULQ%40mail.gmail.com. For more options, visit https://groups.google.com/d/optout.
