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.

Reply via email to