Dear Bruno and Daniel, Thanks for your replies. I used linear_operator and then solve the problem like step-70. I guess there is no "Direct solver" for this situation, the one like we had for dealii BlockSparseMatrix.
Best regards, Masoud On Tuesday, 23 August 2022 at 16:11:12 UTC+1 [email protected] wrote: > Adding to what Bruno is saying, > > If you are not taking advantage of the block structure, you should be able > to treat the block matrix as "regular", i.e., non-block, matrix possibly by > copying it to a corresponding non-block type. > > Best, > Daniel > > On Mon, Aug 22, 2022 at 10:51 AM Bruno Turcksin <[email protected]> > wrote: > >> Hi, >> >> If you search for "block solver" here >> https://dealii.org/developer/doxygen/deal.II/Tutorial.html, you will see >> all the tutorials that use block solvers. I think that only deal.II's own >> solvers support BlockSparseMatrix directly. >> >> Best, >> >> Bruno >> >> On Monday, August 22, 2022 at 9:02:28 AM UTC-4 [email protected] wrote: >> >>> Dear All, >>> >>> The following system of equations: >>> KQ=R >>> where, >>> [image: Screenshot from 2022-08-22 13-45-45.png] >>> were solved using BlockSparseMatrix to form tangent matrix K. It was >>> solved by: >>> >>> SparseDirectUMFPACK A_direct; >>> A_direct.initialize(K); >>> A_direct.vmult(Q_stp, R); >>> >>> Now, I'm trying to run my code with MPI using PETSc/Trilinos, but the >>> solver does not except PETSc/Trilinos BlockSparseMatrix. How do we solve >>> such a general system of equations? >>> >> -- >> > 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/9c70b7e9-08ed-4cd6-bfd7-7eefe5eb2159n%40googlegroups.com >> >> <https://groups.google.com/d/msgid/dealii/9c70b7e9-08ed-4cd6-bfd7-7eefe5eb2159n%40googlegroups.com?utm_medium=email&utm_source=footer> >> . >> > -- 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/1fe4c5db-0e4d-4662-996c-c4363d6326adn%40googlegroups.com.
