#13662: matrix() fails for polynomial quoutient rings over inexact rings
-------------------------------------+-------------------------------------
Reporter: saraedum | Owner: roed
Type: defect | Status: needs_review
Priority: minor | Milestone: sage-6.3
Component: padics | Resolution:
Keywords: sd59 | Merged in:
Authors: Julian Rueth | Reviewers:
Report Upstream: N/A | Work issues:
Branch: | Commit:
u/saraedum/ticket/13662 | 50132187410c11a5c4a5407a28453d14446b87d2
Dependencies: #14482 | Stopgaps:
-------------------------------------+-------------------------------------
Comment (by saraedum):
Replying to [comment:1 nbruin]:
> However, the alternative is to create the companion matrix of the
generator (that's just a matter of putting coefficients into a matrix; no
arithmetic required), take its powers, and take the appropriate linear
combination of that. The big advantage there is that you completely avoid
doing polynomial divisions/remainders with inexact coefficients, so the
numerical stability of your computations is much easier to
analyse/guarantee.
I decided to just reduce leading zero coefficient mod `f` when creating
the quotient ring elements. This appears to be the right thing to do to
me. A different algorithm to compute the matrix would be an alternative of
course. I just did what appeared to be easier to implement.
--
Ticket URL: <http://trac.sagemath.org/ticket/13662#comment:13>
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.
