#14965: FloatingPoint exception when multiplying matrix over QQ by vector over
SR
-------------------------------------+-------------------------------------
Reporter: dkrenn | Owner:
Type: defect | Status: needs_review
Priority: major | Milestone: sage-
Component: linear algebra | duplicate/invalid/wontfix
Keywords: matrix, vector, | Resolution:
mutliplication, floatingpoint | Merged in:
exception | Reviewers: Jeroen Demeyer
Authors: | Work issues:
Report Upstream: N/A | Commit:
Branch: | Stopgaps:
Dependencies: |
-------------------------------------+-------------------------------------
Changes (by jdemeyer):
* status: new => needs_review
* reviewer: => Jeroen Demeyer
* milestone: => sage-duplicate/invalid/wontfix
Comment:
Works for me:
{{{
sage: Z=1/3
sage: def matrix_entry(s,r):
....: result = 0
....: result += Z^(s-r)
....: return result
....:
sage: A = matrix([[matrix_entry(s,r) for r in range(1,22)] for s in
range(1,22)])
sage: v = vector([(2/3)^s+x*(2/3+4/9*(2/3)^30)^s for s in range(1,22)])
sage: A*v
[...huge output...]
}}}
--
Ticket URL: <http://trac.sagemath.org/ticket/14965#comment:1>
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.
