search

adipopescu / To-Do-List

My Singular Project
0 stars 0 forks source link

new segfault examples #45

Closed jakobkroeker closed 8 years ago

jakobkroeker commented 9 years ago 

ring rng = (integer),(x,y),(lp(1),lp(1),C);
short = 0 ;
ideal I = 19*x^2*y+2*x^2+14*x*y^2-33*x*y,-3*x^4-34*x^2*y^2,18*x^3-32*x^2*y^2-9*y-29;
ideal gI =  std(I);

//////////////////////////////////////////////////////////////////////////////////////////////////
ring rng = (integer),(x,y,z),(Wp(4,2,2),C);
short = 0 ;
ideal J = -42*x^2+25*x*y-3*z^3+10*x,29*z^4+18,-14*x^3*y+25*x*y*z^2+5*z;
ideal gJ =  std(J);

//////////////////////////////////////////////////////////////////////////////////////////////////
ring rng = (integer),(x,y,z),(lp(3),C);
short = 0 ;
ideal I = 29*x^2*z-27*x*y*z^2-16*x*y*z-7*y*z,-19*x*z+30*y*z^2-31*y*z+39,28*x^3*z-32*x*y*z^2-41*x*y;
ideal gI =  std(I);

//////////////////////////////////////////////////////////////////////////////////////////////////
// segfault
ring rng = (integer),(x,y,z),(lp(1),lp(1),lp(1),C);
short = 0 ;
intvec op = 100663296,16787537;
option(set, op) ;
option() ;
system("random", 2146987388);
ideal J = -17*x*y*z+18*z^2,-44*x*y^2*z-33*x*y*z+25*y+7,-x^3*z+31*x+19;
ideal gJ =  std(J);
quit;

// no segfault
ring rng = (integer),(x,y,z),(lp(1),lp(1),lp(1),C);
short = 0 ;
intvec op = 100663296,16787537;
option(set, op) ;
// option() ; 
system("random", 2146987388);
ideal J = -17*x*y*z+18*z^2,-44*x*y^2*z-33*x*y*z+25*y+7,-x^3*z+31*x+19;
ideal gJ =  std(J);
quit;

//////////////////////////////////////////////////////////////////////////////////////////////////
// segfault
ring rng = (integer),(x,y,z),(lp(1),dp(2),C);
short = 0 ;
intvec op = 100663298,16787537;
option(set, op) ;
option() ;
system("random", 1117796194);
ideal I = -2*x^2*y*z+41*y^2*z^2-16*y^2+40*z,20*x*y*z^2+43*x*z^2-19*z,-30*x^2*y^2-20*x^2*z^2-44;
ideal J = 21*x*y^2*z-2*x*y+24*y^2,34*x-26*y^3*z-5*y^2*z;
       ideal gI =  std(I);

//no  segfault      
ring rng = (integer),(x,y,z),(lp(1),dp(2),C);
short = 0 ;
intvec op = 100663298,16787537;
//option(set, op) ;
// option() ;
system("random", 1117796194);
ideal I = -2*x^2*y*z+41*y^2*z^2-16*y^2+40*z,20*x*y*z^2+43*x*z^2-19*z,-30*x^2*y^2-20*x^2*z^2-44;
ideal J = 21*x*y^2*z-2*x*y+24*y^2,34*x-26*y^3*z-5*y^2*z;
       ideal gI =  std(I);

//////////////////////////////////////////////////////////////////////////////////////////////////      

// segfault
ring rng = (integer),(x,y,z),(lp(1),dp(2),C);
short = 0 ;
intvec op = 100663296,16787537;
 option(set, op) ;
// option() ;
ideal I = 13*x^2*z^2-22*x*y+20,27*x^2+26*x*y*z^2+45*x*z^2,5*x^2*z+4*x*y^2*z+2*x*y*z-13*x*z^2;
ideal gI =  std(I);
OK 

// no segfault
ring rng = (integer),(x,y,z),(lp(1),dp(2),C);
short = 0 ;
intvec op = 100663296,16787537;
// option(set, op) ;
// option() ;
ideal I = 13*x^2*z^2-22*x*y+20,27*x^2+26*x*y*z^2+45*x*z^2,5*x^2*z+4*x*y^2*z+2*x*y*z-13*x*z^2;
ideal gI =  std(I);

//////////////////////////////////////////////////////////////////////////////////////////////////
LIB "primdec.lib";
ring rng = (integer),(x,y,z),(dp(3),C);
short = 0 ;
ideal I = 9*x*y-1;
ideal J = 13*x^2*y-11*x*y-10*z,12*x^2*z+15*x*z^2-2*y+6,-6*x^2*z-z^3+4*x*y+4;
       ideal gI =  std(I);
       ideal ggI =  std(gI);
       ideal gJ =  std(J);
OK
jakobkroeker commented 9 years ago

update: the last example

LIB "primdec.lib";
ring rng = (integer),(x,y,z),(dp(3),C);
short = 0 ;
ideal I = 9*x*y-1;
ideal J = 13*x^2*y-11*x*y-10*z,12*x^2*z+15*x*z^2-2*y+6,-6*x^2*z-z^3+4*x*y+4;
       ideal gI =  std(I);
       ideal ggI =  std(gI);
       ideal gJ =  std(J);

finishes now, but seems big and not correct - compare with M2 result!

adipopescu commented 9 years ago

Wir würden gerne das Ergebnis vergleichen. Hast Du lead(gJ) aus m2 gerade greifbar?In Singular ist es_[1]=15_xz^2[2]=219700_y^2z[3]=66_x_yz[4]=6_x^2z[5]=114244y^3[6]=2_xy^2[7]=x^2y[8]=3276x^3[9]=83z^4[10]=yz^3[11]=xz^3[12]=26_y^2z^2[13]=x_yz^2[14]=3_x^2z^2[15]=338_y^3z[16]=1092_x^4  

 jakobkroeker <notifications@github.com> schrieb am 16:20 Mittwoch, 18.März 2015:

update: the last exampleLIB "primdec.lib"; ring rng = (integer),(x,y,z),(dp(3),C); short = 0 ; ideal I = 9_x_y-1; ideal J = 13_x^2_y-11_x_y-10_z,12_x^2_z+15_x_z^2-2_y+6,-6_x^2_z-z^3+4_x_y+4; ideal gI = std(I); ideal ggI = std(gI); ideal gJ = std(J); finishes now, but seems big and not correct - compare with M2 result!— Reply to this email directly or view it on GitHub.

jakobkroeker commented 9 years ago 
R = ZZ[x,y,z]
J = ideal(13*x^2*y-11*x*y-10*z,12*x^2*z+15*x*z^2-2*y+6,-6*x^2*z-z^3+4*x*y+4)

gJ = gb J;

toString entries leadTerm ideal  gens gJ
{{15*x*z^2,
 219700*y^2*z, 
66*x*y*z,
 6*x^2*z, 
114244*y^3, 
2*x*y^2, 
x^2*y, 
3276*x^3, 
83*z^4,
 y*z^3, 
x*z^3, 
26*y^2*z^2, 
x*y*z^2,
 3*x^2*z^2,
 338*y^3*z,
      1092*x^4}}
adipopescu commented 9 years ago

Hmm.. sieht genau so aus, oder?

jakobkroeker commented 9 years ago

Hmm.. sieht genau so aus, oder?

ja, du hast recht. Aber warum sind die Koeffizienten im tail so gross?

jakobkroeker commented 9 years ago

it seems hard to find segfault examples, which crash in debug mode, too. here are some:

ok now:

ring rng = (integer),(x,y,z),(Wp(3),Wp(5),Wp(3),C);
option(prot);
short = 0 ;
ideal J = -4*x^2+y+8,-13*x*z-2*z^2+12,7*x*y-20*x+15*y;
ideal gJ =  std(J);

not ok ?

ring rng = (integer),(x,y,z),(Dp(1),Dp(2),C);
option(prot);
short = 0 ;
ideal I = 15*x^3*y+34*y+15*z,-29*x^3*y+30*y*z+28*z,-3*x^2+27*x*y*z^2;
ideal gI =  std(I);
adipopescu commented 9 years ago

I have tested some of them. The last one seems not to terminate...

adipopescu commented 8 years ago

I don't see any more segfaults