Opened 13 years ago
Closed 8 years ago
#267 closed bug (fixed)
preimage versus eliminate
Reported by: | Wolfram Decker | Owned by: | |
---|---|---|---|
Priority: | major | Milestone: | 3-1-2 and higher |
Component: | dontKnow | Version: | 3-1-1 |
Keywords: | Cc: |
Description
I am worried about the timings when comparing preimage and eliminate:
ring SS = 0, (u,v,z), dp;
poly f = -2*u*v4*z4+u4*v5+12*u4*v3*z2+12*u2*v4*z3-u3*v*z5+11*u3*v2*z4-21*u3*v3*z3-4*u4*v*z4+2*u4*v2*z3-6*u4*v4*z+u5*z4-3*u5*v2*z2+u5*v3*z-3*u*v5*z3-2*u2*v3*z4+u3*v4*z2+v5*z4; f = subst(f,z,u+v+z);
ideal IC = f; ideal AI; AI[1]=v7+4v6z+6v5z2+4v4z3+v3z4; AI[2]=uv6+3uv5z+v6z+3uv4z2+3v5z2+uv3z3+3v4z3+v3z4; AI[3]=u2v5+3u2v4z+uv5z+3u2v3z2+3uv4z2+u2v2z3+3uv3z3+uv2z4; AI[4]=u3v4+2u3v3z+2u2v4z+u3v2z2+4u2v3z2+uv4z2+2u2v2z3+2uv3z3+uv2z4; AI[5]=u4v3+2u4v2z+2u3v3z+u4vz2+4u3v2z2+u2v3z2+2u3vz3+2u2v2z3+u2vz4; AI[6]=u5v2+u5vz+3u4v2z+3u4vz2+3u3v2z2+3u3vz3+u2v2z3+u2vz4; AI[7]=u6v+u6z+3u5vz+3u5z2+3u4vz2+3u4z3+u3vz3+u3z4; AI[8]=u7+4u6z+6u5z2+4u4z3+u3z4;
int aa = timer;
ring Rtarget = 0, (y(1..8)), dp; setring Rtarget; ideal RNC = preimage(SS, AI, IC);
timer-aa;
68
aa = timer;
def Rbig = Rtarget + SS; setring Rbig; ideal IC = imap(SS,IC); ideal AI = imap(SS,AI);
ideal J; J = IC; int k; for(k=1;k<=8;k++)
{
J=J,var(k)-AI[k];
}
ideal RNC = eliminate(J,uvz);
1
Change History (2)
comment:1 Changed 12 years ago by
Owner: | changed from somebody to decker@… |
---|
comment:2 Changed 8 years ago by
Resolution: | → fixed |
---|---|
Status: | new → closed |
preimage and eliminate are internally the nearly same: preimage has a little bit more freedom in chosing an eliminataion ordering, yielding in general a little faster results