Chen, I am not sure if TrilinosWrappers::SolverDirect is supported but for your second point, you can use the linear_operator as a preconditioner. That's exactly what's done here <https://github.com/dealii/dealii/blob/master/tests/lac/step-40-linear_operator_06.cc#L220-L251> for example. The preconditioner is wrapped in a LinearOperator and this LinearOperator is then used as preconditioner.
Best, Bruno On Tuesday, April 19, 2022 at 10:39:17 AM UTC-4 [email protected] wrote: > Hi everyone, > > More questions: > 1. where can I find the supported linear solver of Trilinos of > linear_operator? TrilinosWrappers::SolverDirect seems not to be supported > by the linear_operator > 2. Does linear_operator can be a preconditioner? for example, > SparseMatrix A, B, C, D, E, F; Matrix G=A-B*C^{-1}*D-E^{-1}*F; So Matrix > G can only represent by a linear_operator, not its entity. I want to get > its inverse by inverse_operator using some linear solvers as: PreconditionILU > prec_G(op_G); const auto op_G_inv = inverse_operator(G, Solver, prec_G); But > it seems that the precondition of linear_operator is not supported. Is > there any solution? > > best > Chen > > 在2022年4月18日星期一 UTC+8 22:07:36<[email protected]> 写道： > >> Bruno Turcksin >> Thank you! it really helpful! >> >> best >> chen >> 在2022年4月18日星期一 UTC+8 21:07:45<[email protected]> 写道： >> >>> Chen, >>> >>> Yes, LinearOperator does support Trilinos and PETSc matrices. >>> LinearOperator even supports your own matrix type as long as you define >>> vmult and Tvmult (see here >>> <https://dealii.org/current/doxygen/deal.II/group__LAOperators.html#ga6b74b7ed76e4aeea2c67c0bbb03e3e24> >>> ) >>> >>> Best, >>> >>> Bruno >>> >>> On Monday, April 18, 2022 at 4:19:42 AM UTC-4 [email protected] wrote: >>> >>>> Hi everyone, >>>> >>>> When solving a block matrix system, the LinearOperator is a useful >>>> tool. My question is whether the LinearOpertor supports Matrix type (the >>>> same as the Vector type) like PETSCWrapper::MPI::SparseMatrix or >>>> TrilinosWrapper::MPI::SparseMatrix? >>>> >>>> Best >>>> Chen >>>> >>> -- 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/6ee40751-54ae-402d-903e-ae745545387bn%40googlegroups.com.
