octave-maintainers
[Top][All Lists]
Advanced

[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]

Re: test fftfilt fail w/ current tip


From: Daniel J Sebald
Subject: Re: test fftfilt fail w/ current tip
Date: Wed, 24 Jul 2013 00:14:20 -0500
User-agent: Mozilla/5.0 (X11; U; Linux x86_64; en-US; rv:1.9.2.24) Gecko/20111108 Fedora/3.1.16-1.fc14 Thunderbird/3.1.16

On 07/23/2013 09:07 PM, Dmitri A. Sergatskov wrote:
On Tue, Jul 23, 2013 at 7:46 PM, Daniel J Sebald <address@hidden
<mailto:address@hidden>> wrote:


    I looks like the test that is failing is supposed to test the
    arithmetic of the case when the inputs are not integer.  When the
    inputs are integer (second column), the result is "cast" to integer
    by dropping the floating portion.  In the first column the inputs
    are not all integer so the epsilon should be left intact.  However,
    perhaps because of machine precision the epsilon is being lost and
    it just so happens the output is exactly 1.

    Dmitri, on your system, please test a larger discrepancy, say 2*eps
    or 3*eps:

    assert (fftfilt (b, [x.'+2*eps x.']) == [1 1 0 0 0 0 0 0 0 0].'*[1
    1], [false(10, 1) true(10, 1)]);



Test passes with 2*eps.

OK, thanks. If no one gets around to this by the weekend, I'll post a patch to Savannah then.

Dan


reply via email to

[Prev in Thread] Current Thread [Next in Thread]