[Octave-bug-tracker] [bug #53258] test pcg.m fail

[Dmitri A. Sergatskov]

Follow-up Comment #1, bug #53258 (project octave):

I usually get few failures if I repeat test 1000 times. 
This is one of the failed cases:


b = [
   0.12278472
   0.10887554
   0.37156328
   0.05834160
   0.52981297
   0.68811238
   0.00012278
   0.05733196
   0.78572195
   0.06894106
];


N = 10;
A = diag ([1:N]);
X = A \ b;  # X is the true solution
[x, flag] = pcg (A, b, [], N+1);
norm (x - X) / norm (X)

# ans =  0.00000057028



This is with OpenBLAS.

Dmitri.
-- 




    _______________________________________________________

Reply to this item at:

  <http://savannah.gnu.org/bugs/?53258>

_______________________________________________
  Message sent via/by Savannah
  http://savannah.gnu.org/