Forums before death by AOL, social media and spammers... "We can't have nice things"
|    sci.math.symbolic    |    Symbolic algebra discussion    |    10,432 messages    |
[   << oldest   |   < older   |   list   |   newer >   |   newest >>   ]
|    Message 9,519 of 10,432    |
|    bursejan@gmail.com to All    |
|    =?UTF-8?Q?Re=3A_Test_Cases_for_Gr=C3=B6b    |
|    05 Jul 17 17:59:52    |
      Now I have found myself an example, where my algorithm       still runs wrong. This is not desired:              ?- F is X*Y^2-X, G is Y*X^2-Y, R is F/G.       R is (X*Y-X*Y^3)/((1-X^2)*Y^2)              Checking my Gröbner Bases algorithm, there is something       going wrong somewhere:              ?- F is X*Y^2-X, G is Y*X^2-Y, S is F*Z, T is G*(1-Z),        sys_poly_groeb([T,S],H), R is H.       R is [(X-X^3)*Y^2-(X-X^3)*Y^4,       - (1-X^2)*Y^2-(X^2-Y^2)*Z,       - (1-X^2)*Y+(1-X^2)*Y*Z,       (X-X^3)*Y^2-(X-X^3)*Z]              Compare this with the maxima result, for the same        problem, the last element is just F*G, as desired:              (%i57) poly_buchberger([(x*y^2-x)*t,(y*x^2-y)*(1-t)],[t,x,y]);       (%o57) [t*x*y^2-t*x,(-t*x^2*y)+x^2*y+t*y-y,        (-x^2*y^2)-t*y^2+y^2+t*x^2,        (-x^2*y^3)-t*y^3+y^3+x^2*y+t*y-y,        x^3*y^3-x*y^3-x^3*y+x*y]              So going back to the lab, and        lets see what goes wrong.              Am Mittwoch, 5. Juli 2017 20:50:01 UTC+2 schrieb burs...@gmail.com:       > Maybe an effect of the used ordering. Dunno yet, might       > investigate more closely in the future. Dunno whether       > its a bad sign that maxima has a bigger base, or a good       > sign. Could be also a good sign.       >        > A more interesting case (test case) would be if my       > algorithm yields multiple members more than 3, and        > another system only yields 3. BTW: There are also       > situations where the resulting GB        >        > is single or two membered, more trivially:       >        > ?- X is 0, Y is A*B, S is X*Z, T is Y*(1-Z), sys_poly_groeb([S,T],H), R is H.       > R is [A*B-A*B*Z]       >        > ?- X is 1, Y is A*B, S is X*Z, T is Y*(1-Z), sys_poly_groeb([S,T],H), R is H.       > R is [Z,A*B],       >        > Quiz: Which one is the LCM?              --- SoupGate-Win32 v1.05        * Origin: you cannot sedate... all the things you hate (1:229/2)    |
[   << oldest   |   < older   |   list   |   newer >   |   newest >>   ]
(c) 1994, bbs@darkrealms.ca