|
From: | Gerard BURNSIDE |
Subject: | [Octave-bug-tracker] [bug #36568] wrong flag returned by gmres |
Date: | Thu, 31 May 2012 15:20:48 +0000 |
User-agent: | Mozilla/5.0 (Windows NT 5.1) AppleWebKit/536.11 (KHTML, like Gecko) Chrome/20.0.1132.21 Safari/536.11 |
URL: <http://savannah.gnu.org/bugs/?36568> Summary: wrong flag returned by gmres Project: GNU Octave Submitted by: gerardburnside Submitted on: jeu. 31 mai 2012 15:20:47 GMT Category: None Severity: 3 - Normal Priority: 5 - Normal Item Group: Inaccurate Result Status: None Assigned to: None Originator Name: Originator Email: Open/Closed: Open Discussion Lock: Any Release: 3.6.1 Operating System: Any _______________________________________________________ Details: While checking out the samples for gmres and bicgstab, I noticed that gmres returned a flag 1 (max num iterations reached) where I expected it to return a flag 0 (converged). Looking at gmres.m it seems that you just need to change the inequality after the while loop from: if (presn > rtol * prec_b_norm) flag = 0; to: if (presn <= rtol * prec_b_norm) flag = 0; (by default the flag is initializd to 1). _______________________________________________________ Reply to this item at: <http://savannah.gnu.org/bugs/?36568> _______________________________________________ Message posté via/par Savannah http://savannah.gnu.org/
[Prev in Thread] | Current Thread | [Next in Thread] |