tolerance on Gauss-Seidel algorithm
1 Ansicht (letzte 30 Tage)
Ältere Kommentare anzeigen
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
0 Kommentare
Antworten (0)
Siehe auch
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!