On 2/28/26 05:28, Praveen C wrote:
We are exploring different linear solvers for a mixed H(div) method for
poisson equation.
We want to try both block solvers based on Schur form and also gmres, direct,
etc.
So we declare matrices and vectors of block form.
But then we are not able to use TrilinosWrappers::SolverGMRES with these
blocked variables.
Is this a fundamental limitation or something else ?
I have attached a MWE which shows compiler error when I call GMRES.
Hi Praveen,
yes, when you use TrilinosWrappers::SolverGMRES, then it wants a single
Trilinos matrix which a block matrix is not. But you can use
dealii::SolverGMRES with any kind of linear operator, which a block matrix is.
So just use the latter.
Best
W.
--
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 visit
https://groups.google.com/d/msgid/dealii/e753b8d7-df9b-4cd2-96b2-66c3651be488%40colostate.edu.
