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.

88 CHAPTER 4. ALGORITHMS FOR STABLE IDEALSOutput <strong>of</strong> the procedure compute expansion.◦ M — a list <strong>of</strong> lists, representing the set <strong>of</strong> minimal generators <strong>of</strong> the input ideal inlexicographic order after the expansion <strong>of</strong> the monomial mcompute_expansion:= proc(M, m, n)local newMonomials, last, mon, Ind, del, v, i, ind, j, monomial;begin/* Compute the set <strong>of</strong> new monomials to be included intothe set M <strong>of</strong> minimal generators. Hence, first computethe index ’last’ <strong>of</strong> the last variable occurring in m */for i from n+1 downto 0 doif m[i] 0 thenlast:= i;break;end_if;end_for;newMonomials:= [0 $ n-last+1];for i from last to n domon:= m;mon[i]:= mon[i] + 1;newMonomials[i-last+1]:= mon;end_for;/* Testing for minimality, i.e. removing all monomials,which are no longer minimal generators <strong>of</strong> the ideal */Ind:= {};for monomial in newMonomials d<strong>of</strong>or j from contains(M,m) + 1 to nops(M) dodel:= {};for i from 1 to n+1 doif monomial[i] - M[j][i] > 0 thendel:= del union {1};elif monomial[i] - M[j][i] < 0 thendel:= del union {-1};end_if;end_for;if contains(del, 1) = FALSE orcontains(del, -1) = FALSE then

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

Saved successfully!

Ooh no, something went wrong!