Back to Forum | View unanswered posts | View active topics
Topic review - newbie stuck solving a parametrized system |
Author |
Message |
|
|
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));
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));[/code]
|
|
|
|
Posted: Thu Jul 19, 2007 6:58 pm |
|
|
|
|
|
It is currently Fri May 13, 2022 10:56 am
|
|