|
D.5.3.3 esStratumProcedure from library
Example: LIB "equising.lib"; int p=printlevel; printlevel=1; ring r = 0,(a,b,c,d,e,f,g,x,y),ds; poly F = (x2+2xy+y2+x5)+ax+by+cx2+dxy+ey2+fx3+gx4; list M = esStratum(F); M[1]; → _[1]=g → _[2]=f → _[3]=b → _[4]=a → _[5]=-4c+4d-4e+d2-4ce printlevel=3; // displays additional information esStratum(F,2); // es stratum over Q[a,b,c,d,e,f,g] / <a,b,c,d,e,f,g>^2 → // → // Compute HN development → // ---------------------- → // finished → // → // Blowup Step 1 completed → // Blowup Step 2 completed → // Blowup Step 3 completed → // 1 branch finished → // → // Elimination starts: → // ------------------- → // finished → // → // output of 'esStratum' is list consisting of: → // _[1] = ideal defining equisingular stratum → // _[2] = 0 → [1]: → _[1]=b → _[2]=a → _[3]=c-d+e → _[4]=g → _[5]=f → [2]: → 0 ideal I = f-fa,e+b; qring q = std(I); poly F = imap(r,F); esStratum(F); → // → // Compute HN development → // ---------------------- → // finished → // → // Blowup Step 1 completed → // Blowup Step 2 completed → // Blowup Step 3 completed → // 1 branch finished → // → // Elimination starts: → // ------------------- → // finished → // → // output of 'esStratum' is list consisting of: → // _[1] = ideal defining equisingular stratum → // _[2] = 0 → [1]: → _[1]=e → _[2]=a → _[3]=-4c+4d+d2 → _[4]=g → [2]: → 0 printlevel=p; See also: esIdeal; isEquising. |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() |