Charles Greathouse on Tue, 10 Dec 2013 06:48:27 +0100


[Date Prev] [Date Next] [Thread Prev] [Thread Next] [Date Index] [Thread Index]

Failed bench linear


My compiled gp-sta comes up with a strange answer for the lindep bench:

? lindep(Mod(1,7)*[2,-1;1,3])
%1 = [Mod(6, 7), Mod(5, 7)]~

which does not match the expected [-3, 1]~ from src/test/64/linear or my dynamically-compiled version of gp. Any idea what might be wrong? Any other information I should provide? The whole contents of the dif file:

*** ../src/test/64/linear 2013-12-10 00:28:06.000000000 -0500
--- gp.out 2013-12-10 00:35:18.000000000 -0500
***************
*** 12,18 ****
  ? charpoly(Mod(1,8191)*[1,2;3,4],z,2)
  z^2 + Mod(8186, 8191)*z + Mod(8189, 8191)
  ? lindep(Mod(1,7)*[2,-1;1,3])
! [-3, 1]~
  ? lindep([(1-3*sqrt(2))/(3-2*sqrt(3)),1,sqrt(2),sqrt(3),sqrt(6)])
  [3, 3, -9, 2, -6]~
  ? lindep([(1-3*sqrt(2))/(3-2*sqrt(3)),1,sqrt(2),sqrt(3),sqrt(6)],14)
--- 12,18 ----
  ? charpoly(Mod(1,8191)*[1,2;3,4],z,2)
  z^2 + Mod(8186, 8191)*z + Mod(8189, 8191)
  ? lindep(Mod(1,7)*[2,-1;1,3])
! [Mod(6, 7), Mod(5, 7)]~
  ? lindep([(1-3*sqrt(2))/(3-2*sqrt(3)),1,sqrt(2),sqrt(3),sqrt(6)])
  [3, 3, -9, 2, -6]~
  ? lindep([(1-3*sqrt(2))/(3-2*sqrt(3)),1,sqrt(2),sqrt(3),sqrt(6)],14)
***************
*** 675,680 ****
  ? vector(10,x,1/x)
  [1, 1/2, 1/3, 1/4, 1/5, 1/6, 1/7, 1/8, 1/9, 1/10]
  ? getheap
! [96, 10148]
  ? print("Total time spent: ",gettime);
  Total time spent: 8
--- 675,680 ----
  ? vector(10,x,1/x)
  [1, 1/2, 1/3, 1/4, 1/5, 1/6, 1/7, 1/8, 1/9, 1/10]
  ? getheap
! [96, 10160]
  ? print("Total time spent: ",gettime);
  Total time spent: 8

I think this is my first experience with a seriously wrong bench, rather than just precision errors or the like.

Charles Greathouse
Analyst/Programmer
Case Western Reserve University