Singular https://www.singular.uni-kl.de/forum/ |
|
Bug in minAssGTZ with real coefficients? https://www.singular.uni-kl.de/forum/viewtopic.php?f=10&t=2777 |
Page 1 of 1 |
Author: | uffish [ Fri Aug 17, 2018 7:29 pm ] |
Post subject: | Bug in minAssGTZ with real coefficients? |
Hi all, The following program causes Singular to receive a segmentation fault and exit. (background: if you replace the coefficient ring to Q, the program works, but is very slow: the std in the beginning takes about an hour. I guessed that the reason was the appearnace of very large coefficients and tried to switch to real coefficients. It really made the program faster, but...) Thanks, ring R=(real,200),(a,b,c,d,e,f),dp; poly pxb=a+b*d-a*b*d+c*e-a*c*e-b*c*d*e+a*b*c*d*e+c*d*f-a*c*d*f-b*c*d*f+a*b*c*d*f+b*e*f-a*b*e*f-b*c*e*f+a*b*c*e*f-b*d*e*f+a*b*d*e*f-c*d*e*f+a*c*d*e*f+2b*c*d*e*f-2a*b*c*d*e*f; poly pxA=a+b-a*b+c-a*c-b*c+a*b*c; poly pab=a*b+d-a*b*d+b*c*e-a*b*c*e-b*c*d*e+a*b*c*d*e+a*c*f-a*b*c*f-a*c*d*f+a*b*c*d*f+e*f-a*b*e*f-a*c*e*f-b*c*e*f+2a*b*c*e*f-d*e*f+a*b*d*e*f+a*c*d*e*f+b*c*d*e*f-2a*b*c*d*e*f; poly pa2b=a*c+b*c*d-a*b*c*d+e-a*c*e-b*c*d*e+a*b*c*d*e+a*b*f-a*b*c*f+d*f-a*b*d*f-a*c*d*f-b*c*d*f+2a*b*c*d*f-a*b*e*f+a*b*c*e*f-d*e*f+a*b*d*e*f+a*c*d*e*f+b*c*d*e*f-2a*b*c*d*e*f; poly dpxb0=1-b*d-c*e+b*c*d*e-c*d*f+b*c*d*f-b*e*f+b*c*e*f+b*d*e*f+c*d*e*f-2b*c*d*e*f; poly dpxA0=1-b-c+b*c; poly dpab0=b-b*d-b*c*e+b*c*d*e+c*f-b*c*f-c*d*f+b*c*d*f-b*e*f-c*e*f+2b*c*e*f+b*d*e*f+c*d*e*f-2b*c*d*e*f; poly dpxb1=d-a*d-c*d*e+a*c*d*e-c*d*f+a*c*d*f+e*f-a*e*f-c*e*f+a*c*e*f-d*e*f+a*d*e*f+2c*d*e*f-2a*c*d*e*f; poly dpxA1=1-a-c+a*c; poly dpab1=a-a*d+c*e-a*c*e-c*d*e+a*c*d*e-a*c*f+a*c*d*f-a*e*f-c*e*f+2a*c*e*f+a*d*e*f+c*d*e*f-2a*c*d*e*f; poly dpxb2=e-a*e-b*d*e+a*b*d*e+d*f-a*d*f-b*d*f+a*b*d*f-b*e*f+a*b*e*f-d*e*f+a*d*e*f+2b*d*e*f-2a*b*d*e*f; poly dpxA2=1-a-b+a*b; poly dpab2=b*e-a*b*e-b*d*e+a*b*d*e+a*f-a*b*f-a*d*f+a*b*d*f-a*e*f-b*e*f+2a*b*e*f+a*d*e*f+b*d*e*f-2a*b*d*e*f; poly dpxb3=b-a*b-b*c*e+a*b*c*e+c*f-a*c*f-b*c*f+a*b*c*f-b*e*f+a*b*e*f-c*e*f+a*c*e*f+2b*c*e*f-2a*b*c*e*f; poly dpxA3=0; poly dpab3=1-a*b-b*c*e+a*b*c*e-a*c*f+a*b*c*f-e*f+a*b*e*f+a*c*e*f+b*c*e*f-2a*b*c*e*f; poly dpxb4=c-a*c-b*c*d+a*b*c*d+b*f-a*b*f-b*c*f+a*b*c*f-b*d*f+a*b*d*f-c*d*f+a*c*d*f+2b*c*d*f-2a*b*c*d*f; poly dpxA4=0; poly dpab4=b*c-a*b*c-b*c*d+a*b*c*d+f-a*b*f-a*c*f-b*c*f+2a*b*c*f-d*f+a*b*d*f+a*c*d*f+b*c*d*f-2a*b*c*d*f; poly dpxb5=c*d-a*c*d-b*c*d+a*b*c*d+b*e-a*b*e-b*c*e+a*b*c*e-b*d*e+a*b*d*e-c*d*e+a*c*d*e+2b*c*d*e-2a*b*c*d*e; poly dpxA5=0; poly dpab5=a*c-a*b*c-a*c*d+a*b*c*d+e-a*b*e-a*c*e-b*c*e+2a*b*c*e-d*e+a*b*d*e+a*c*d*e+b*c*d*e-2a*b*c*d*e; ideal I=pab-pa2b,dpxb0-pxA*dpab0-dpxA0*pab,dpxb1-pxA*dpab1-dpxA1*pab,dpxb2-pxA*dpab2-dpxA2*pab,dpxb3-pxA*dpab3-dpxA3*pab,dpxb4-pxA*dpab4-dpxA4*pab,dpxb5-pxA*dpab5-dpxA5*pab; option(prot); LIB "primdec.lib"; list l=minAssGTZ(I); |
Author: | hannes [ Sat Aug 25, 2018 4:52 pm ] |
Post subject: | Re: Bug in minAssGTZ with real coefficients? |
main incredience of minAssGTZ are Groebner bases, which are not useful for inexact coeffcients, and factorization, which is not implemented for inexact coeffcients: in short: that cannot work. |
Page 1 of 1 | All times are UTC + 1 hour [ DST ] |
Powered by phpBB © 2000, 2002, 2005, 2007 phpBB Group http://www.phpbb.com/ |