Thread: Software for IBM AIX View Single Post
 2010-08-31, 16:47 #20 pacionet     Oct 2005 Italy 15316 Posts Without -ffast-math and with -O1: Code: gcc -maix64 -lm -lpthread -Wall \ -O1 \ -o phrot.g5 \ phrot.c \ -O1 -DY_MEM_THRESHOLD=8192 -DY_KILL_BRANCHES -DY_VECTORIZE \ -I../glucas ../glucas/libyeafft.a The testm seems work: Input 103914*4001^4001-1 : Actually testing 103914*4001^4001-1 (witness=3 4003/8192 limbs) 103914*4001^4001-1 is PRP. (e=0.00000 (6.28077e-07~3.53103e-16@0.000) t=15.16s) Input 103951848*10^2555-1 : Actually testing 103951848*100000^511-1 (witness=3 513/1152 limbs) 103951848*10^2555-1 is PRP. (e=0.00005 (0.000150342~2.35792e-16@0.000) t=0.35s) Input 1041845805*10^2562-1 : Actually testing 1041845805*1000^854-1 (witness=3 858/1792 limbs) 1041845805*10^2562-1 is PRP. (e=0.00000 (1.53901e-08~2.54359e-16@0.000) t=0.54s) Input 104297193*10^5069-1 : Actually testing 1042971930*10000^1267-1 (witness=3 1270/2560 limbs) 104297193*10^5069-1 is PRP. (e=0.00000 (1.88708e-06~3.51235e-16@0.000) t=1.52s) [...] but this test fails: Code: bash-3.2\$ phrot.g5 -q 5028*10^83982+1 sizeof(int64)=8 Phil Carmody's Phrot (0.72) Input 5028*10^83982+1 : Actually testing 5028*1000000^13997+1 (witness=3 13998/28672 limbs) 5028*10^83982+1 is composite LLR64=0000000000000000. (e=0.00098 (0.0824785~8.27893e-18@0.000) t=164.23s) How can I test if the yeafft is correctly compiled ? Other tests with composite and known residues (failed) (with the same options, ffast-math removed and O1) 56093*36^58139-1 is composite: RES64: [00BFD059739B458C] My result: LLR64=0000000000000000. 56093*36^59413-1 is composite: RES64: [3A3314C49046512B] My result: LLR64=0000000000000000 Last fiddled with by pacionet on 2010-08-31 at 17:00