logo       

Re: MacOS X 10.3.9 build but load-save.m test failure: msg#00067

gnu.octave.bugs

Subject: Re: MacOS X 10.3.9 build but load-save.m test failure

On 28-Apr-2005, Dmitri A. Sergatskov wrote:

| John W. Eaton wrote:
| > On 28-Apr-2005, Dmitri A. Sergatskov wrote:
| >
| > | EXEC: ../src/octave -f -q -p .:../src//:./../scripts//
./octave.test/linalg/expm-2.m
| > |
| > | Octave Output:
| > | ans = 0
| > | FAIL: octave.test/linalg/expm-2.m
| >
| > The test is:
| >
| > arg = [1, 1; 0, 1];
| > result = [2.718281828459045, 2.718281828459045;
| > 0.000000000000000, 2.718281828459045];
| > all (all (abs (expm (arg) - result) < 2*eps))
| >
|
| But the error I get is :
|
| Octave Output:
| error: eig: argument must be a square matrix
| error: near line 1 of file `./octave.test/linalg/eig-5.m'
| child process exited abnormally
|
| ???

I think you are misreading the log file (it is easy to do). The text
you quote above applies to eig-5.m, which is supposed to produce an
error (so it should be marked with a PASS). The text I quoted above
is all that applies to the expm-2.m test (you want to look for
everything between the EXEC and either a PASS or FAIL line).

| > My guess is it is not computing a result to the expected precision.
| >
| > What is the max difference on your system? If it is reasonably small
| > (say < 256*eps or so) then probably we should just increase the
| > tolerance in the test file.
|
| The "make check" used to be OK up to 2.1.67 (the latest I run). Was the test
changed lately?

The CVS log shows that this test was added March 18 2005.

| Anyway:
|
|
| octave:6> format long
| octave:7> abs (expm (arg) - result)
| ans =
|
| 4.44089209850063e-16 4.44089209850063e-16
| 0.00000000000000e+00 4.44089209850063e-16
|
| octave:8> 2*eps
| ans = 4.44089209850063e-16

I suppose 4*eps would do it then. I made the change.

Thanks,

jwe



-------------------------------------------------------------
Octave is freely available under the terms of the GNU GPL.

Octave's home on the web: http://www.octave.org
How to fund new projects: http://www.octave.org/funding.html
Subscription information: http://www.octave.org/archive.html
-------------------------------------------------------------




<Prev in Thread] Current Thread [Next in Thread>
Google Custom Search

News | FAQ | advertise