Mon 19 Dec 2011 09:54:49 AM UTC, comment #16:
I have confirmed the fix. Thanks!
|
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.
|
Sun 18 Dec 2011 09:19:08 PM UTC, comment #14:
I have attached trace.txt.
Hope this result is useful.
(file #24626)
|
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.
|
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
|
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)
|
Sun 02 Oct 2011 03:02:27 PM UTC, comment #10:
Re-titling bug report now that only one issue remains.
|
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
|
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
|
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.
|
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
|
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
|
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
|
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
|
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
|
Tue 26 Jul 2011 05:25:57 AM UTC, comment #1:
The second fail is the same as
http://savannah.gnu.org/bugs/?33792
|
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
|