[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]
[Octave-bug-tracker] [bug #41624] pcg fails to converge with preconditio
|
From: |
Eugenio Gianniti |
|
Subject: |
[Octave-bug-tracker] [bug #41624] pcg fails to converge with preconditioner luinc |
|
Date: |
Mon, 17 Feb 2014 21:45:52 +0000 |
|
User-agent: |
Mozilla/5.0 (Macintosh; Intel Mac OS X 10_9_1) AppleWebKit/537.73.11 (KHTML, like Gecko) Version/7.0.1 Safari/537.73.11 |
Follow-up Comment #1, bug #41624 (project octave):
I checked your code on Matlab R2013b for Mac OS X. I found out that luinc has
been removed from Matlab, anyway an equivalent script:
A=[...]
b = ones(size(A,1),1);
b = A'*b;
A = A'*A;
setup.droptol = 0;
[l, u] = ilu(A,setup);
pcg(A,b,1e-8,size(A,1),l,u,[]);
yields the following output:
pcg stopped at iteration 1599 without converging to the desired tolerance
1e-08
because the maximum number of iterations was reached.
The iterate returned (number 1219) has relative residual 0.26.
_______________________________________________________
Reply to this item at:
<http://savannah.gnu.org/bugs/?41624>
_______________________________________________
Message sent via/by Savannah
http://savannah.gnu.org/