Post new topic Reply to topic  [ 1 post ] 
Author Message
 Post subject: newbie stuck solving a parametrized system
PostPosted: Thu Jul 19, 2007 6:58 pm 

Joined: Wed Jul 18, 2007 6:16 pm
Posts: 1
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));


Report this post
Top
 Profile  
Reply with quote  
Display posts from previous:  Sort by  
Post new topic Reply to topic  [ 1 post ] 

You can post new topics in this forum
You can reply to topics in this forum
You cannot edit your posts in this forum
You cannot delete your posts in this forum
You cannot post attachments in this forum

It is currently Fri May 13, 2022 10:54 am
cron
Powered by phpBB © 2000, 2002, 2005, 2007 phpBB Group