make check failure in chol
Kai Habel
kai.habel at gmx.de
Tue Dec 16 14:08:04 CST 2008
With a very recent tip I see the following make check error:
>>>>> processing /home/kai/hg-octave/octave/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
34 values do not match>>>>> processing
/home/kai/hg-octave/octave/src/DLD-FUNCTIONS/conv2.cc
I am not sure, but shouldn't we test for values smaller 1e-10 ?
E.g.:
assert (norm(Ainv-Ainv1),0,1e-10)
Kai
More information about the Octave-maintainers
mailing list