#18593: Relative tolerance in French Sage book
-------------------------------------+-------------------------------------
Reporter: rbeezer | Owner:
Type: defect | Status: positive_review
Priority: trivial | Milestone: sage-6.8
Component: documentation | Resolution:
Keywords: french book | Merged in:
doctest tolerance | Reviewers: Marc Mezzarobba
Authors: Rob Beezer | Work issues:
Report Upstream: N/A | Commit:
Branch: u/rbeezer | bf5b521bcbba59895cfae220284a3bc6ec05d12e
/tolerance-french-book | Stopgaps:
Dependencies: |
-------------------------------------+-------------------------------------
Changes (by zimmerma):
* cc: cpernet (added)
Comment:
just curious, I'd like to investigate more. Here is what I get with Sage
6.8.beta1 on my workstation:
{{{
sage: A = matrix(RDF, [[1,3,2],[1,4,2],[0,5,2],[1,3,2]]); A
[1.0 3.0 2.0]
[1.0 4.0 2.0]
[0.0 5.0 2.0]
[1.0 3.0 2.0]
sage: b = vector(RDF, [1,2,3,4]); b
(1.0, 2.0, 3.0, 4.0)
sage: transpose(A)
[1.0 1.0 0.0 1.0]
[3.0 4.0 5.0 3.0]
[2.0 2.0 2.0 2.0]
sage: R = transpose(A)*b; R
(7.0, 38.0, 20.0)
sage: Z = transpose(A)*A; Z
[ 3.0 10.0 6.0]
[10.0 59.0 30.0]
[ 6.0 30.0 16.0]
sage: Z.solve_right(R)
(-1.5000000000000049, -0.5000000000000011, 2.750000000000004)
sage: from sage.misc.citation import get_systems
sage: get_systems('Z.solve_right(R)')
['numpy', 'scipy']
}}}
What do you get on the machine where the test did fail?
Paul
--
Ticket URL: <http://trac.sagemath.org/ticket/18593#comment:4>
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 unsubscribe from this group and stop receiving emails from it, send an email
to [email protected].
To post to this group, send email to [email protected].
Visit this group at http://groups.google.com/group/sage-trac.
For more options, visit https://groups.google.com/d/optout.
