lcm and gcd of polynomials over the irrationals

asked 2013-11-26 06:10:45 +0100

gundamlh
167 ●15 ●21 ●29

Dear all,

I want to calculate the least common multiple of a list of polynomials with coefficients defined over RR or RDF in SAGE. Is it possible?

sage: R.<s> = PolynomialRing(RDF); R
Univariate Polynomial Ring in s over Real Double Field
sage: H = matrix(1,2, [(s+5)**2, s]); H
[s^2 + 10.0*s + 25.0                   s]
sage: lcm(H[0,0], H[0,1])
TypeError: Singular error:
   ? not implemented
   ? error occurred in or before poly.lib::lcm line 852: `           q=gcd(p,i[k]);`
   ? leaving poly.lib::lcm
   skipping text from `;` error at token `)`

Thanks in advance!

answered 2013-11-26 06:41:35 +0100

tmonteil
27413 ●31 ●204 ●517 http://wiki.sagemath.o...

updated 2013-11-26 06:44:24 +0100

lcm seems to work correctly on RR and QQ:

sage: R.<s> = PolynomialRing(RR); R
Univariate Polynomial Ring in s over Real Field with 53 bits of precision
sage: lcm([(s+5)**2, s])
s^3 + 10.0000000000000*s^2 + 25.0000000000000*s

sage: R.<s> = PolynomialRing(QQ); R
Univariate Polynomial Ring in s over Rational Field
sage: lcm([(s+5)**2, s])
s^3 + 10*s^2 + 25*s

But gcd works on RDF but not on RR:

sage: R.<s> = PolynomialRing(RDF); R
Univariate Polynomial Ring in s over Real Double Field
sage: gcd([(s+5)**2, s])
1.0

sage: R.<s> = PolynomialRing(RR); R
Univariate Polynomial Ring in s over Real Field with 53 bits of precision
sage: gcd([(s+5)**2, s])
TypeError: 'MinusInfinity' object cannot be interpreted as an index

But:

sage: ((s+5)**2).gcd(s)
1.00000000000000

This looks weird. If your polynomials have rational (or integers) coefficients, you can use QQ instead of RDF.

edit flag offensive delete link

Comments

Thanks! But if some coef. is sqrt(2) or pi ...

gundamlh ( 2013-11-26 07:20:29 +0100 )edit

for the gcd, working in the symbolic ring `SR` works, not for the lcm (lack of conversion to the singular lib).

tmonteil ( 2013-11-26 07:31:05 +0100 )edit

OK, perhaps it is not implemented in those cases. Oh, can we call NumPy or Maxima to do it instead?

gundamlh ( 2013-11-26 08:18:04 +0100 )edit

add a comment

Your Answer

Please start posting anonymously - your entry will be published after you log in or create a new account.

lcm and gcd of polynomials over the irrationals

1 Answer

Comments

Your Answer

Question Tools

Stats

Related questions

lcm and gcd of polynomials over the irrationals edit

1 Answer

Comments

Your Answer

Question Tools

Stats

Related questions

lcm and gcd of polynomials over the irrationals