Singular

 Page 1 of 1 [ 1 post ]
 Print view | E-mail friend Previous topic | Next topic
Author Message
 Post subject: newbie stuck solving a parametrized systemPosted: 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));

Top

 Display posts from previous: All posts1 day7 days2 weeks1 month3 months6 months1 year Sort by AuthorPost timeSubject AscendingDescending
 Page 1 of 1 [ 1 post ]

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

 It is currently Sat Oct 20, 2018 5:03 pm
Powered by phpBB © 2000, 2002, 2005, 2007 phpBB Group