On 2013-01-04, Barton Willis <willisb at unk.edu> wrote:
> Using the current Maxima + CCL, rtest_eigen (in the share testsuite) doesn't finish. The problem seems to be
>
> similaritytransform(matrix ([3, 2, 1, -5], [-1, 2, 1, -3], [5, -1, 2, 4], [3, 1, -1, 5]));
Looks like this is a bug in the new 'invert' which punts to
invert_by_lu. When ratmx=true, invert_by_lu succeeds (with arithmetic in
crering) but otherwise it seems to get stuck (with arithmetic in
generalring). (Previously invertmx was called, which always bound ratmx
to true, and succeeded.)
I've appended the matrix which invert_by_lu gets stuck on. When
ratmx=true, it does run to completion.
Not sure where to go from here.
best,
Robert Dodier
PS.
matrix([sqrt(235*sqrt(337)+4291)/sqrt(951*sqrt(337)+18159),
sqrt(235*sqrt(337)+4291)/sqrt(951*sqrt(337)+18159),
43*2^(3/2)/sqrt(sqrt(sqrt(337)-17)
*(369*2^(7/2)*sqrt(337)+6385*2^(7/2))
+sqrt(2*sqrt(337)-34)*(1023*sqrt(337)+18393)
+15758*sqrt(337)+273170),
43*2^(3/2)/sqrt((-1023*sqrt(337)-18393)*sqrt(2*sqrt(337)-34)
+sqrt(sqrt(337)-17)
*(-369*2^(7/2)*sqrt(337)-6385*2^(7/2))
+15758*sqrt(337)+273170)],
[-sqrt(235*sqrt(337)+4291)*(sqrt(-2*sqrt(337)-34)*(3*sqrt(337)-21)
-6*sqrt(337)-246)
/(sqrt(951*sqrt(337)+18159)*(sqrt(-2*sqrt(337)-34)*(2*sqrt(337)+2)
+24*sqrt(337)+408)),
-sqrt(235*sqrt(337)+4291)*(sqrt(-2*sqrt(337)-34)*(3*sqrt(337)-21)
+6*sqrt(337)+246)
/(sqrt(951*sqrt(337)+18159)*(sqrt(-2*sqrt(337)-34)*(2*sqrt(337)+2)
-24*sqrt(337)-408)),
-(sqrt(2*sqrt(337)-34)*(9*sqrt(337)+183)+14*sqrt(337)+170)
/(2^(3/2)*sqrt(sqrt(sqrt(337)-17)
*(369*2^(7/2)*sqrt(337)+6385*2^(7/2))
+sqrt(2*sqrt(337)-34)*(1023*sqrt(337)+18393)
+15758*sqrt(337)+273170)),
(sqrt(2*sqrt(337)-34)*(9*sqrt(337)+183)-14*sqrt(337)-170)
/(2^(3/2)*sqrt((-1023*sqrt(337)-18393)*sqrt(2*sqrt(337)-34)
+sqrt(sqrt(337)-17)
*(-369*2^(7/2)*sqrt(337)-6385*2^(7/2))
+15758*sqrt(337)+273170))],
[sqrt(235*sqrt(337)+4291)*(sqrt(sqrt(337)+17)
*(2^(9/2)*sqrt(337)-67*2^(5/2))*%i
+sqrt(337)-135)
/(172*sqrt(951*sqrt(337)+18159)),
-sqrt(235*sqrt(337)+4291)*(sqrt(sqrt(337)+17)
*(2^(9/2)*sqrt(337)-67*2^(5/2))*%i
-sqrt(337)+135)
/(172*sqrt(951*sqrt(337)+18159)),
-(sqrt(sqrt(337)-17)*(2^(9/2)*sqrt(337)+67*2^(5/2))+sqrt(337)+135)
/(sqrt(2)*sqrt(sqrt(sqrt(337)-17)
*(369*2^(7/2)*sqrt(337)+6385*2^(7/2))
+sqrt(2*sqrt(337)-34)*(1023*sqrt(337)+18393)
+15758*sqrt(337)+273170)),
(sqrt(sqrt(337)-17)*(2^(9/2)*sqrt(337)+67*2^(5/2))-sqrt(337)-135)
/(sqrt(2)*sqrt((-1023*sqrt(337)-18393)*sqrt(2*sqrt(337)-34)
+sqrt(sqrt(337)-17)
*(-369*2^(7/2)*sqrt(337)-6385*2^(7/2))
+15758*sqrt(337)+273170))],
[sqrt(235*sqrt(337)+4291)*(sqrt(sqrt(337)+17)
*(5*sqrt(2)*sqrt(337)-73*sqrt(2))*%i
+3*sqrt(337)-61)
/(172*sqrt(951*sqrt(337)+18159)),
-sqrt(235*sqrt(337)+4291)*(sqrt(sqrt(337)+17)
*(5*sqrt(2)*sqrt(337)-73*sqrt(2))*%i
-3*sqrt(337)+61)
/(172*sqrt(951*sqrt(337)+18159)),
-(sqrt(sqrt(337)-17)*(5*sqrt(2)*sqrt(337)+73*sqrt(2))+3*sqrt(337)+61)
/(sqrt(2)*sqrt(sqrt(sqrt(337)-17)
*(369*2^(7/2)*sqrt(337)+6385*2^(7/2))
+sqrt(2*sqrt(337)-34)*(1023*sqrt(337)+18393)
+15758*sqrt(337)+273170)),
(sqrt(sqrt(337)-17)*(5*sqrt(2)*sqrt(337)+73*sqrt(2))-3*sqrt(337)-61)
/(sqrt(2)*sqrt((-1023*sqrt(337)-18393)*sqrt(2*sqrt(337)-34)
+sqrt(sqrt(337)-17)
*(-369*2^(7/2)*sqrt(337)-6385*2^(7/2))
+15758*sqrt(337)+273170))])