University of Paderborn Department of Mathematics Diploma Thesis ...
University of Paderborn Department of Mathematics Diploma Thesis ...
University of Paderborn Department of Mathematics Diploma Thesis ...
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;