tolerance on Gauss-Seidel algorithm
Afficher commentaires plus anciens
Hi all, I'm writing a function that solves a linear system with the Gauss-Seidel method. I'm going to use this type of arrest:
norm(residuo, Inf) < toll_r*b_norm + toll_a
where
residuo = A*x - b
and
b_ norm = norm(b, Inf)
toll_r is the relative tolerance and toll_a is the absolute tolerance.
I use this criterion because of is on my study book.
My problem is that I don't know how I have to choose the tolerance, for example:
toll_r = 2.22e-16 (Machine epsilon)
toll_a = 1e-18
I don't want write numbers at random
Réponses (0)
Catégories
En savoir plus sur Linear Algebra dans Centre d'aide et File Exchange
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!