[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]
test failures for chol.cc
From: |
John W. Eaton |
Subject: |
test failures for chol.cc |
Date: |
Wed, 24 Dec 2008 14:51:03 -0500 |
On 24-Dec-2008, Ben Abbott wrote:
| I pulled sources yesterday and built today. I encountered many new
| failures. Most result (151) because I'm missing ARPACK.
|
| I also encountered a chol.cc test failure (lines from fntests.log
| below).
|
| 15 >>>>> processing /Users/bpabbott/Development/mercurial/
| octave-3.1.51/src/DLD-FUNCTIONS/chol.cc
| 16 ***** test
| 17 A = [2,0.2;0.2,1];
| 18 issymmetric(A)
| 19 min(eig(A))
| 20 Ainv = inv(A);
| 21 Ainv1 = cholinv(A);
| 22 Ainv2 = inv(sparse(A));
| 23 Ainv3 = cholinv(sparse(A));
| 24 Ainv4 = spcholinv(sparse(A));
| 25 assert (norm(Ainv-Ainv1),1e-10)
| 26 assert (norm(Ainv-Ainv2),1e-10)
| 27 assert (norm(Ainv-Ainv3),1e-10)
| 28 assert (norm(Ainv-Ainv4),1e-10)
| 29 !!!!! test failed
| 30 assert (norm (Ainv - Ainv1),1e-10) expected
| 31 1.0000e-10
| 32 but got
| 33 0
|
| Perhaps, I'm missing the intent, but shouldn't the proper value be
| zero (David)?
|
| 25 assert (norm(Ainv-Ainv1),0,1e-10)
| 26 assert (norm(Ainv-Ainv2),0,1e-10)
| 27 assert (norm(Ainv-Ainv3),0,1e-10)
| 28 assert (norm(Ainv-Ainv4),0,1e-10)
|
| Assuming my understanding is correct, I've attached the trivial
| changeset.
I applied this changeset.
Thanks,
jwe