10.07.2015 Views

University of Paderborn Department of Mathematics Diploma Thesis ...

University of Paderborn Department of Mathematics Diploma Thesis ...

University of Paderborn Department of Mathematics Diploma Thesis ...

SHOW MORE
SHOW LESS

You also want an ePaper? Increase the reach of your titles

YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.

104 CHAPTER 4. ALGORITHMS FOR STABLE IDEALSMuPAD>> L_p:= compute_L_p(poly(2*z^2 + 2*z + 1, [z]), 4):>> nops(compute_ideals(L_p, 4));Output77In the next section, we take advantage <strong>of</strong> procedure 4.17 to compute all Hilbert seriesassociated to a given Hilbert polynomial.The last part <strong>of</strong> this section is dedicated to an implementation <strong>of</strong> Algorithm 3.45, i.e. thenext procedure can be used to compute all saturated stable ideals to a given Hilbertpolynomial.MuPAD Source Code 4.20. Let p be the Hilbert polynomial <strong>of</strong> R/I for I ⊂ R a saturatedstable ideal and n the index <strong>of</strong> the last variable <strong>of</strong> R. The procedure below outputsall saturated stable ideals to the Hilbert polynomial p in R. If one uses the option ’All’as third argument for the procedure, it additionally returns the set and the number <strong>of</strong> alldouble saturations <strong>of</strong> the saturated stable ideals associated to the Hilbert polynomial p.Input for the procedure compute all ideals.◦ p — the Hilbert polynomial <strong>of</strong> R/I◦ n — the index <strong>of</strong> the last variable <strong>of</strong> the polynomial ring R◦ All — an optional, third input parameter (for details, see description <strong>of</strong> the output<strong>of</strong> the procedure below)Output <strong>of</strong> the procedure compute all ideals.◦ M — a set <strong>of</strong> lists <strong>of</strong> lists, encoding all saturated stable ideals J ⊂ R with p R/J (z) =p(z).◦ Sat — a set <strong>of</strong> lists <strong>of</strong> lists, encoding all possible double saturations (as ideals in ¯R :=K[x 0 , . . . , x n−1 ]). Note that Sat will only be returned, if the third input parameterAll is used.◦ cnt — the number <strong>of</strong> elements <strong>of</strong> the set Sat. Note that cnt will only be returned,if the third input parameter All is used.compute_all_ideals:= proc(p, n /*, All */)local d, pList, i, L, j, M, MM, Ideal, p1, p2, steps, J, tmp,Sat, deg, Id, k;begin

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

Saved successfully!

Ooh no, something went wrong!