Post a reply
Username:
Note:If not registered, provide any username. For more comfort, register here.
Subject:
Message body:
Enter your message here, it may contain no more than 60000 characters. 

Smilies
:D :) :( :o :shock: :? 8) :lol: :x :P :oops: :cry: :evil: :twisted: :roll: :wink: :!: :?: :idea: :arrow: :| :mrgreen:
Font size:
Font colour
Options:
BBCode is ON
[img] is ON
[flash] is OFF
[url] is ON
Smilies are ON
Disable BBCode
Disable smilies
Do not automatically parse URLs
Confirmation of post
To prevent automated posts the board requires you to enter a confirmation code. The code is displayed in the image you should see below. If you are visually impaired or cannot otherwise read this code please contact the %sBoard Administrator%s.
Confirmation code:
Enter the code exactly as it appears. All letters are case insensitive, there is no zero.
   

Topic review - newbie stuck solving a parametrized system
Author Message
  Post subject:  newbie stuck solving a parametrized system  Reply with quote
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));
Post Posted: Thu Jul 19, 2007 6:58 pm


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