If it's of any help, our CG and GMRES solvers (and maybe some of the
others?) can return the estimated condition number during / after a call
to solve().
https://dealii.org/developer/doxygen/deal.II/classSolverCG.html#a52f1fd3d41765ec06d2d0add1a74a37c
<https://dealii.org/developer/doxygen/deal.II/classSolverCG.html#a52f1fd3d41765ec06d2d0add1a74a37c>
https://dealii.org/developer/doxygen/deal.II/classSolverGMRES.html#a3025d81c3d7db3e1a1b0114bef794202
<https://dealii.org/developer/doxygen/deal.II/classSolverGMRES.html#a3025d81c3d7db3e1a1b0114bef794202>
Jean-Paul
On 03.12.20 05:46, Wolfgang Bangerth wrote:
On 12/2/20 9:30 PM, Reza Rastak wrote:
I am trying to calculate the condition number of a Trilinos Block
Sparse matrix. Is there an example of that in the documentation that
you can point me to?
No, and that's actually a quite difficult task because you have to
estimate eigenvalues. The easiest way is probably to use the power
iteration and inverse power iteration to estimate largest and smallest
eigenvalues -- see the EigenPower and EigenInverse classes.
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 on the web visit
https://groups.google.com/d/msgid/dealii/85c9726c-da34-1979-e0c1-5ae555c1659d%40gmail.com.