Output of "example Atkin;" from "atkins.lib"; Code: >example Atkin; // proc Atkin from lib atkins.lib EXAMPLE: ring R = 0,x,dp; Atkin(7691,100,5); -1 Atkin(3473,10,2); -1 printlevel=1; Atkin(10000079,100,2); Set i = 0, n = 0 and N(i) = N(0)= 10000079. pause>
List H of possibly suitable discriminants will be calculated. H = -3,-4,-7,-8,-11,-12,-15,-16,-19,-20,-23,-24,-27,-28,-31,-32,-35,-36,-39,-40,-43,-44,-47,-48,-51,-52,-55,-56,-59,-60,-63,-64,-67,-68,-71,-72,-75,-76,-79,-80,-83,-84,-87,-88,-91,-92,-95,-96,-99,-100,-103,-104,-107,-108,-111,-112,-115,-116,-119,-120,-123,-124,-127,-128,-131,-132,-135,-136,-139,-140,-143,-144,-147,-148,-151,-152,-155,-156,-159,-160,-163,-164,-167,-168,-171,-172,-175,-176,-179,-180,-183,-184,-187,-188,-191,-192,-195,-196,-199,-200 pause>
N(0) = 10000079 is divisible by 5. pause>
N(0) = N = 10000079 and therefore N is not prime. pause>
-1
>
In particular 10000079 should be a prime number, the message "10000079 is divisible by 5" is noncence. The reason for the problem is the usage of "a mod b" with number type. Code: > ring r =0, x, dp; > 10000079 mod 5; 4 > number(10000079) mod 5; 0 >
I suppose the problem with the library exists since years.
Output of "example Atkin;" from "atkins.lib"; [code] >example Atkin; // proc Atkin from lib atkins.lib EXAMPLE: ring R = 0,x,dp; Atkin(7691,100,5); -1 Atkin(3473,10,2); -1 printlevel=1; Atkin(10000079,100,2); Set i = 0, n = 0 and N(i) = N(0)= 10000079. pause>
List H of possibly suitable discriminants will be calculated. H = -3,-4,-7,-8,-11,-12,-15,-16,-19,-20,-23,-24,-27,-28,-31,-32,-35,-36,-39,-40,-43,-44,-47,-48,-51,-52,-55,-56,-59,-60,-63,-64,-67,-68,-71,-72,-75,-76,-79,-80,-83,-84,-87,-88,-91,-92,-95,-96,-99,-100,-103,-104,-107,-108,-111,-112,-115,-116,-119,-120,-123,-124,-127,-128,-131,-132,-135,-136,-139,-140,-143,-144,-147,-148,-151,-152,-155,-156,-159,-160,-163,-164,-167,-168,-171,-172,-175,-176,-179,-180,-183,-184,-187,-188,-191,-192,-195,-196,-199,-200 pause>
N(0) = 10000079 is divisible by 5. pause>
N(0) = N = 10000079 and therefore N is not prime. pause>
-1
> [/code]
In particular 10000079 should be a prime number, the message "10000079 is divisible by 5" is noncence.
The reason for the problem is the usage of "a mod b" with number type. [code] > ring r =0, x, dp; > 10000079 mod 5; 4 > number(10000079) mod 5; 0 > [/code] I suppose the problem with the library exists since years. :(
|