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/f2bc8f82-6d9e-456e-b15f-c9a32f35721cn%40googlegroups.com.
