19.07.2014 Views

Contents - Student subdomain for University of Bath

Contents - Student subdomain for University of Bath

Contents - Student subdomain for University of Bath

SHOW MORE
SHOW LESS

Create successful ePaper yourself

Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.

130 CHAPTER 4. MODULAR METHODS<br />

Figure 4.9: Diagrammatic illustration <strong>of</strong> g.c.d.s in Z[x, y] (1)<br />

Z[y][x] - - - - - - - - - - - - - - - - - - - - - - - - -> - Z[y][x]<br />

k×reduce ↓<br />

↑ interpret<br />

& check<br />

Z[y] y=v1 [x]<br />

.<br />

gcd<br />

−→<br />

.<br />

gcd<br />

Z[y] y=v1 [x]<br />

Z[y] y=vk [x] −→ Z[y] y=vk [x] ⎪⎭<br />

↑<br />

using Algorithm 14/15<br />

.<br />

⎫<br />

⎪⎬<br />

C.R.T.<br />

−→<br />

Z[y][x]<br />

∏ ′<br />

(y − v 1 ) · · · (y − v k )<br />

Figure 4.10: Diagrammatic illustration <strong>of</strong> g.c.d.s in Z[x, y] (2)<br />

Z[y][x] - - - - - - - - - - - - - - - - - - - - - - - - -> - Z[y][x]<br />

k×reduce ↓<br />

↑ interpret<br />

& check<br />

Z p1 [y][x]<br />

.<br />

gcd<br />

−→<br />

.<br />

gcd<br />

Z p1 [y][x]<br />

Z pk [y][x] −→ Z pk [y][x] ⎪⎭<br />

↑<br />

using Algorithm 17<br />

.<br />

⎫<br />

⎪⎬<br />

C.R.T.<br />

−→<br />

Z[y][x]<br />

∏ ′<br />

p 1 · · · p k<br />

Should this happen, we actually know that all reductions so far are wrong,<br />

because they are compatible in degree, but cannot be right. Hence we start<br />

again.<br />

Z p [y][x] Here we use an analogue <strong>of</strong> Figure 4.2 to reduce us to the case <strong>of</strong> Figure<br />

4.7. The overall structure is shown in Figure 4.10.<br />

Open Problem 6 (Which is the Better Route <strong>for</strong> Bivariate g.c.d.?) As<br />

far as the author can tell, the question <strong>of</strong> which route to follow has not been systematically<br />

explored. The initial literature [Bro71b, Algorithm M] and the more<br />

recent survey [Lau82] assume the route in Figure 4.10. [Bro71b] explicitly assumes<br />

that p is large enough that we never run out <strong>of</strong> values, i.e. that the<br />

algebraic extension at the end <strong>of</strong> Figure 4.8 is never needed.<br />

Implementations the author has seen tend to follow the route in Figure 4.9.<br />

This is probably <strong>for</strong> pragmatic reasons: as one is writing a system one first<br />

wishes <strong>for</strong> univariate g.c.d.s, so implements Algorithm 14 (or 15). Once one<br />

has this, Figure 4.9 is less work.<br />

There is a natural tendency to believe that Figure 4.10 is better, as all numbers<br />

involved are bounded by p untl the very last Chinese Remainder calculations.

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!