generic initial ideals berechnen

[kirsten]

Ich hätte gerne eine Methode, die mir für ein graduiertes Ideal und eine gegebene Termordnung (mit großer Wahrscheinlichkeit) das generische Initialideal berechnet.

[decker]

Forwarded message from Viktor Levandovskyy <levandov@…> (die ging an die reply -Adresse, trot des "do not reply" im Subject.)

Hallo,

wir haben darueber mal im Forum diskutiert. Ich habe damals folgenden Kod geschrieben. Viele Gruesse, Viktor

LIB "random.lib";
LIB "poly.lib";
proc gin(ideal I)
{
 matrix M = randommat(1,nvars(basering), maxideal(1), 1000);
 ideal m = ideal(M); map F = basering, m;
 ideal J = F(I);
 J = groebner(J);
 J = normalize(lead(J));
 return(J);
}
example
{
// revlex gin
ring r = 0,(a,b,c,d),dp;
ideal I = a^3+c^2*d, b^3-a*d^2;
gin(I);
// lex gin, quite big
ring r = 0,(a,b,c,d),lp;
ideal I = a^3+c^2*d, b^3-a*d^2;
gin(I);
}