Ask Your Question
1

resultant and gcd over real ?

asked 2014-02-12 18:03:49 +0100

etienne gravatar image

updated 2014-02-12 20:16:21 +0100

slelievre gravatar image

The resultant and gcd work well over rational but not over real or complex!

R.<x,y> = PolynomialRing(RR,2)
a = x^2 + y
b = x - y^2
a.resultant(b)
a.gcd(b)

This does not work while with QQ, it works. Anyone has an idea?

edit retag flag offensive close merge delete

2 Answers

Sort by ยป oldest newest most voted
1

answered 2014-02-16 04:38:03 +0100

Jeroen Demeyer gravatar image

If you define your polynomial ring as follows, it does work:

sage: P.<y> = PolynomialRing(RR)
sage: R.<x> = PolynomialRing(P)
sage: a = x^2 + y
sage: b = x - y^2
sage: a.resultant(b)
edit flag offensive delete link more
1

answered 2014-02-13 06:23:20 +0100

lftabera gravatar image

The problem is that for these kind of inexact rings the methods are not implemented.

For the resultant, you can compute the resultant with

a.sylvester_matrix(b,x).det()

But I cannot promise performance nor numeric stability.

The gcd is more tricky. For instance, if you take univariate polynomials over an inexact rings and perform Euclides method, you will end with a gcd 1 due to numerical errors. You are asking for an approximate gcd which is not easy to compute.

edit flag offensive delete link more

Your Answer

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

Add Answer

Question Tools

1 follower

Stats

Asked: 2014-02-12 18:03:49 +0100

Seen: 810 times

Last updated: Feb 16 '14