bugGNU Octave - Bugs: bug #33863, One test fails in quadl.m (MinGW)

 
 

bug #33863: One test fails in quadl.m (MinGW)

Submitted by:  Tatsuro MATSUOKA <tmacchant>
Submitted on:  Tue 26 Jul 2011 04:29:56 AM UTC  
 
Category: InterpreterSeverity: 3 - Normal
Priority: 5 - NormalItem Group: Incorrect Result
Status: FixedAssigned to: None
Originator Name: Open/Closed: Closed
Release: devOperating System: Microsoft Windows

Add a New Comment(Rich Markup)
   

You are not logged in

Please log in, so followups can be emailed to you.

 

(Jump to the original submission Jump to the original submission)

Mon 19 Dec 2011 09:54:49 AM UTC, comment #16:

I have confirmed the fix. Thanks!

Tatsuro MATSUOKA <tmacchant>
Mon 19 Dec 2011 05:18:44 AM UTC, comment #15:

I relaxed the tolerance to -1.1e-15 (http://hg.savannah.gnu.org/hgweb/octave/rev/6eeb9e8e63cf) and it should work now. This will be part of the 3.6.0 release.

Rik <rik5>
Project Administrator
Sun 18 Dec 2011 09:19:08 PM UTC, comment #14:

I have attached trace.txt.
Hope this result is useful.

(file #24626)

Tatsuro MATSUOKA <tmacchant>
Fri 16 Dec 2011 06:49:37 PM UTC, comment #13:

Tatsuro,

Maybe we debug this by running the "TRACE" feature of quadl to output the intermediate results.

Can you try the following commands and attach the output, "trace.txt", to this bug report.

Rik <rik5>
Project Administrator
Mon 03 Oct 2011 08:08:13 AM UTC, comment #12:

Rik,

you are right. I checked your result with maxima. I cannot reproduce my old "exact" value. Sorry for the noise.

Marco

Marco Caliari <caliari>
Project Member
Mon 03 Oct 2011 03:21:39 AM UTC, comment #11:

Marco suggested that the exact value for the integral test that is failing is incorrect in the last digit. I'm unable to verify this. I used the website wolframalpha.com which provides a front end to Mathematica and their arbitrary precision output. The value for the integral in question is attached as .gif. Mathematica believes the current value for the integral (2.588424538641647e+10) is correct.

This means that there will still be one test failing for quadl.m on the MinGW platform.

(file #24087)

Rik <rik5>
Project Administrator
Sun 02 Oct 2011 03:02:27 PM UTC, comment #10:

Re-titling bug report now that only one issue remains.

Rik <rik5>
Project Administrator
Fri 26 Aug 2011 03:36:28 AM UTC, comment #9:

Sorry I have misled

The result for quad (but not quadl)
octave:3> quad (@(x) sin (3x).cosh (x).*sinh (x),10,15)
ans = 25884245386.4172

is little differ from that by quadl.

However, the last digit may be corrected as Marco said

Tatsuro MATSUOKA <tmacchant>
Fri 26 Aug 2011 02:39:24 AM UTC, comment #8:

Follow-up Comment #7: (Rik)

The error with quadl and large tolerances has been fixed in this changeset
(http://hg.savannah.gnu.org/hgweb/octave/rev/85dac13a911b). If you can verify this then we can close the bug report.

The changeset solved the first error.

The second error is not related the above changeset.

As Marco pointed out, the change
should be attached.

Regards

Tatsuro

Tatsuro MATSUOKA <tmacchant>
Thu 25 Aug 2011 05:04:25 PM UTC, comment #7:

The error with quadl and large tolerances has been fixed in this changeset (http://hg.savannah.gnu.org/hgweb/octave/rev/85dac13a911b). If you can verify this then we can close the bug report.

Rik <rik5>
Project Administrator
Thu 18 Aug 2011 12:05:29 AM UTC, comment #6:

Sorry for my late reply.

The change you have shown

is OK for the MinGW build.

Regards

Tatsuro

Tatsuro MATSUOKA <tmacchant>
Tue 09 Aug 2011 06:50:02 AM UTC, comment #5:

Dear Tatsuro,

the last digit in the reference solution is wrong. I think you should change your changeset into

Do you confirm that with 2 instead of 7 as last digit, you get a relative error < 9e-15 on MinGW?

Marco

Marco Caliari <caliari>
Project Member
Tue 09 Aug 2011 05:44:42 AM UTC, comment #4:

I will propose a changeset the below to avoid error on the MinGW platform

Tatsuro MATSUOKA <tmacchant>
Thu 04 Aug 2011 09:28:44 PM UTC, comment #3:

Hello Macro

For the first erro, you are right. Can you make a changeset for the test?

Regards

Tatsuro

Tatsuro MATSUOKA <tmacchant>
Mon 01 Aug 2011 02:28:45 PM UTC, comment #2:

Dear Tatsuro,

regarding the first failure, do you know that quadl uses a default relative tolerance equal to eps? Then,

which is not eps, but still quite good (and less than 9e-15). About the value of exact, it seems that the value provided in the test is wrong in the last digit: in fact, I also get

but if I use the exact value

everything is fine. Regarding the second failure, I will say something in #33792.

Marco

Marco Caliari <caliari>
Project Member
Tue 26 Jul 2011 05:25:57 AM UTC, comment #1:

The second fail is the same as
http://savannah.gnu.org/bugs/?33792

Tatsuro MATSUOKA <tmacchant>
Tue 26 Jul 2011 04:29:56 AM UTC, original submission:

Hello

I have encountered two test fails in quadl.m

For the first fail,

> quadl (@(x) sin (3 * x) .* cosh (x) .* sinh (x), 10, 15) - 2.588424538641647e+10

ans = -2.70843505859375e-004

Error seems large.

For second fail, error is apparently large.
I tried

quadl does not work correctly in some cases.

Regards

Tatsuro

Tatsuro MATSUOKA <tmacchant>

 

(Note: upload size limit is set to 16384 kB, after insertion of the required escape characters.)

Attach File(s):
   
   
Comment:
   

Attached Files
file #24626:  trace.txt added by tmacchant (39KiB - text/plain)
file #24087:  integral.gif added by rik5 (8KiB - image/gif)

 

Depends on the following items: None found

Digest:
   bug dependencies.

 

Carbon-Copy List
  • -unavailable- added by rik5 (Posted a comment)
  • -unavailable- added by caliari (Posted a comment)
  • -unavailable- added by tmacchant (Submitted the item)
  •  

    Do you think this task is very important?
    If so, you can click here to add your encouragement to it.
    This task has 0 encouragements so far.

    Only project members can vote.

     

    Please enter the title of George Orwell's famous dystopian book (it's a date):

     

     

    Follow 7 latest changes.

    Date Changed By Updated Field Previous Value => Replaced By
    Mon 19 Dec 2011 05:18:44 AM UTCrik5StatusNeed Info=>Fixed
      Open/ClosedOpen=>Closed
    Sun 18 Dec 2011 09:19:08 PM UTCtmacchantAttached File-=>Added trace.txt, #24626
    Wed 30 Nov 2011 06:19:21 AM UTCrik5Dependencies-=>bugs #34942 is dependent
    Mon 03 Oct 2011 03:21:39 AM UTCrik5Attached File-=>Added integral.gif, #24087
    Sun 02 Oct 2011 03:02:27 PM UTCrik5SummaryTwo test fails in quadl.m (dev. source MinGW)=>One test fails in quadl.m (MinGW)
    Thu 25 Aug 2011 05:04:25 PM UTCrik5StatusNone=>Need Info

    Back to the top


    Powered by Savane 3.1-cleanup1