|
|
|
last edited 16 years ago by japp |
Edit detail for #405 ContinuedFraction and Polynomial Domains revision 5 of 7
| 1 2 3 4 5 6 7 | ||
|
Editor: page
Time: 2007/12/12 14:40:06 GMT-8 |
||
| Note: chech change to version | ||
added:
From page Wed Dec 12 14:40:06 -0800 2007
From: page
Date: Wed, 12 Dec 2007 14:40:06 -0800
Subject: chech change to version
Message-ID: <20071212144006-0800@axiom-wiki.newsynthesis.org>
Axiom Version: /usr/local/lib/axiom/target/x86_64-unknown-linux => OpenAxiom-20071209
Consider
axiomd := continuedFraction(0,[1 for i in 1..], [1 for i in 0..])$CONTFRAC INT
 | (1) |
Type: ContinuedFraction? Integer
axiom2*d
 | (2) |
Type: ContinuedFraction? Integer
so far, so good. But changing the domain to CONTFRAC UP(x, FRAC INT) makes axiom crash.
The problem appears to be, that INT is ordered, while UP(x, FRAC INT) is not. In contfrac.spad we find:
eucWhole(a: Q): R == numer a quo denom a
eucWhole0(a: Q): R ==
isOrdered =>
n := numer a
d := denom a
q := n quo d
r := n - q*d
if r < 0 then q := q - 1
q
eucWhole a
If R is INT, eucWhole0(4/3) yields 1, in the other case it returns 4/3. I haven't been able to investigate further yet, though.
Martin
Specify version of Axiom --page, Wed, 12 Dec 2007 14:17:18 -0800 reply
Axiom Version: => friCAS-20071109
Axiom Version: /usr/local/lib/axiom/target/x86_64-unknown-linux => OpenAxiom?-20071209