#10952: better numerical accuracy testing
-----------------------------------------+----------------------------------
Reporter: robertwb | Owner: mvngu
Type: enhancement | Status: needs_review
Priority: critical | Milestone: sage-4.7.2
Component: doctest | Keywords:
Work_issues: | Upstream: N/A
Reviewer: Jason Grout, Mariah Lenox | Author: Robert Bradshaw
Merged: | Dependencies:
-----------------------------------------+----------------------------------
Changes (by rbeezer):
* status: needs_info => needs_review
Old description:
> If a line contains ``tol`` or ``tolerance, numerical results are only
> verified to the given tolerance. This may be prefixed by ``abs[olute``
> or ``rel[ative]`` to specify whether to measure absolute or relative
> error; defaults to relative error except when the expected value is
> exactly zero:
>
> {{{
> sage: RDF(pi) # abs tol 1e-5
> 3.14159
> sage: [10^n for n in [0.0 .. 4]] # rel tol 2e-4
> [0.9999, 10.001, 100.01, 999.9, 10001]
> }}}
>
> This can be useful when the exact output is subject to rounding error
> and/or processor floating point arithmetic variation.
New description:
If a line contains ``tol`` or ``tolerance, numerical results are only
verified to the given tolerance. This may be prefixed by ``abs[olute``
or ``rel[ative]`` to specify whether to measure absolute or relative
error; defaults to relative error except when the expected value is
exactly zero:
{{{
sage: RDF(pi) # abs tol 1e-5
3.14159
sage: [10^n for n in [0.0 .. 4]] # rel tol 2e-4
[0.9999, 10.001, 100.01, 999.9, 10001]
}}}
This can be useful when the exact output is subject to rounding error
and/or processor floating point arithmetic variation.
'''Apply''':
1. [attachment:10952-tol-bin.patch]
1. [attachment:10952-tol-doc.2.patch]
1. [attachment:trac_10952-reviewer-docs.patch]
--
Comment:
I have added an "Apply" list for this ticket.
My patch fixes one typo in the documentation (a missing ]).
I have added some non-trivial tests.
Jeroen:
1. abs() for a matrix computes the determinant. Not useful for testing.
2. Is there a class of polynomials with numerical coefficients? Does it
have abs() defined on it?
I am working on this as a high-priority ticket for Bug Days 32.
Rob
--
Ticket URL: <http://trac.sagemath.org/sage_trac/ticket/10952#comment:21>
Sage <http://www.sagemath.org>
Sage: Creating a Viable Open Source Alternative to Magma, Maple, Mathematica,
and MATLAB
--
You received this message because you are subscribed to the Google Groups
"sage-trac" group.
To post to this group, send email to [email protected].
To unsubscribe from this group, send email to
[email protected].
For more options, visit this group at
http://groups.google.com/group/sage-trac?hl=en.
