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

Create successful ePaper yourself

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

106 CHAPTER 4. ALGORITHMS FOR STABLE IDEALSprint(Unquoted, "<strong>of</strong> all ideals to the given Hilbert polynomial");print(Unquoted, "equals the whole ring");return(M);end_if;if deg = 1 and args(0) = 3 and args(3) = hold(All) thencnt:= nops(M);Sat:= M;end_if;/* Successively compute all saturated ideals */for i from d+1 downto 2 do/* Consider the elements <strong>of</strong> M as ideals in the polynomialring with one more variable (lifting-process) */M:= map(M, Ideal -> map(Ideal, ind -> [op(ind), 0]));/* Compute the Hilbert polynomials and perform expansionsif necessary. */p1:= pList[i-1];tmp:= {};for Ideal in M dop2:= compute_Hilbert_polynomial(Ideal, n-i+2);steps:= p1 - p2;if degree(steps) = 0 and op(steps, 1) > 0 thenM:= M minus {Ideal};J:= Ideal;for j from 1 to op(steps, 1) doJ:= compute_expansion(J, J[nops(J)], n-i+2);end_for;tmp:= tmp union {J};elif degree(steps) > 0 or op(steps, 1) < 0 thenM:= M minus {Ideal};if i = 2 and deg > 1 thenId:= map(Ideal, ind -> [ind[k] $ k = 1..n]);Sat:= Sat minus {Id};end_if;end_if;end_for;

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

Saved successfully!

Ooh no, something went wrong!