Singular https://www.singular.uni-kl.de/forum/ |
|
newbie stuck solving a parametrized system https://www.singular.uni-kl.de/forum/viewtopic.php?f=10&t=1645 |
Page 1 of 1 |
Author: | guyarie [ Thu Jul 19, 2007 6:58 pm ] |
Post subject: | newbie stuck solving a parametrized system |
Hello, the subject says it all... i paste my code below, it ran all night on a AMD64 Turion dual core, used above 90% CPU almost all of the time, but still got stuck inside fglm. Thanks a lot Guy Code: LIB "solve.lib";
ring r4n = (0,n),(x,y,z),dp; poly p1 = z^2*x^2*y^2-y^2-x^2; poly p2 = -n+n*x+3*n*z^2*x^2-2*n*z^2*x^3+51*z^3*x^6-2*z^3*x^7; poly p3 = -n+n*y+3*n*z^2*y^2-2*n*z^2*y^3+51*z^3*y^6-2*z^3*y^7; ideal i1 = p1,p2,p3; ideal si1 = std(i1); printf("standard ideal should be zero, it is: %2s",dim(si1)); printf("number of complex solutions is: %2s",vdim(si1)); // it is 62 option(redSB); ideal maxI=maxideal(1); ideal jnumbers = sat(si1,maxI)[1]; // this takes about 10 seconds printf("number of non-zero solutions: %2s",vdim(jnumbers)); //still 62 // triangulate ring T4n = (0,n),(x,y,z),lp; print("calling fglm"); // the last line printed in the actual run ideal j4n = fglm(r4n,jnumbers); print("calling triangl"); list Sol4n = triangMH(j); print("DONE!, %s components",size(Sol4n)); |
Page 1 of 1 | All times are UTC + 1 hour [ DST ] |
Powered by phpBB © 2000, 2002, 2005, 2007 phpBB Group http://www.phpbb.com/ |