View Single Post
Old 2010-08-31, 19:38   #25
pacionet
 
pacionet's Avatar
 
Oct 2005
Italy

3×113 Posts
Default

Quote:
Originally Posted by rogue View Post
The -D values are specific to glucas and are unlikely to cause problems. -O3 could be a problem. Try rebuilding with -O1, then rebuild phrot and see what happens. If -O1 works, then try -O2 and do it again.
I was not lucky.
Rebuilt yeafft with -O1

Code:
gcc -DHAVE_CONFIG_H -I. -I. -I.     -O1  -DY_MEM_THRESHOLD=8192 -DY_KILL_BRANCHES -DY_VECTORIZE -c `test -f 'ynorm_16.c' || echo './'`ynorm_16.c
source='yeafft.c' object='yeafft.o' libtool=no \
depfile='.deps/yeafft.Po' tmpdepfile='.deps/yeafft.TPo' \
depmode=gcc3 /bin/sh ./depcomp \
gcc -DHAVE_CONFIG_H -I. -I. -I.     -O1  -DY_MEM_THRESHOLD=8192 -DY_KILL_BRANCHES -DY_VECTORIZE -c `test -f 'yeafft.c' || echo './'`yeafft.c
source='yeafft1.c' object='yeafft1.o' libtool=no \
depfile='.deps/yeafft1.Po' tmpdepfile='.deps/yeafft1.TPo' \
depmode=gcc3 /bin/sh ./depcomp \
gcc -DHAVE_CONFIG_H -I. -I. -I.     -O1  -DY_MEM_THRESHOLD=8192 -DY_KILL_BRANCHES -DY_VECTORIZE -c `test -f 'yeafft1.c' || echo './'`yeafft1.c
source='yeainit.c' object='yeainit.o' libtool=no \
depfile='.deps/yeainit.Po' tmpdepfile='.deps/yeainit.TPo' \
depmode=gcc3 /bin/sh ./depcomp \
gcc -DHAVE_CONFIG_H -I. -I. -I.     -O1  -DY_MEM_THRESHOLD=8192 -DY_KILL_BRANCHES -DY_VECTORIZE -c `test -f 'yeainit.c' || echo './'`yeainit.c
source='radix_2.c' object='radix_2.o' libtool=no \
depfile='.deps/radix_2.Po' tmpdepfile='.deps/radix_2.TPo' \
and phrot

Code:
gcc -maix64 -lm -lpthread  -Wall \
 \
-o phrot.g5 \
phrot.c \
-O1  -DY_MEM_THRESHOLD=8192 -DY_KILL_BRANCHES -DY_VECTORIZE \
-I../glucas ../glucas/libyeafft.a
Same wrong result ... with residues with a lot of zeros:

Code:
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=173.92s)
Code:
bash-3.2$ phrot.g5 -q 56093*36^58139-1
sizeof(int64)=8
Phil Carmody's Phrot (0.72)
Input 56093*36^58139-1 :  Actually testing 72696528*46656^19379-1 (witness=3 19381/40960 limbs)
56093*36^58139-1 is composite LLR64=0000000000000000. (e=0.00047 (0.00021126~1.53918e-15@0.000) t=267.87s)

On the other hand glucas seems to work

Code:
[Thu Aug 26 19:53:09 2010]
bash-3.2$ glucas 23209
Going to work with exponent 23209
Starting from iteration 1. Exponent 23209.
M23209. Saved file at iteration 4096. Res64: 60DD379FEFC017FE.
M23209. Saved file at iteration 8192. Res64: 57B014FAC6BAED3A.
Iter. 10000 ( 43.09%), Err= 0.086, 0.35 user  52% CPU (0.000 sec/iter).
M23209. Saved file at iteration 12288. Res64: 92B6CB2BE3CD9ECE.
M23209. Saved file at iteration 16384. Res64: 5322CA8405A26B7B.
Iter. 20000 ( 86.18%), Err= 0.078, 0.35 user  52% CPU (0.000 sec/iter).
M23209. Saved file at iteration 20480. Res64: 4EAA7F09C57C7B75.
[Tue Aug 31 21:24:56 2010]
M23209 is a known Mersenne prime!
Code:
bash-3.2$ glucas 216091
M216091. Saved file at iteration 167936. Res64: 65DAC3959C92014B.
Iter. 170000 ( 78.67%), Err= 0.000, 6.25 user 100% CPU (0.001 sec/iter).
M216091. Saved file at iteration 172032. Res64: 82AC0AEF60948C69.
M216091. Saved file at iteration 176128. Res64: 95582DC8C101C4E9.
Iter. 180000 ( 83.30%), Err= 0.000, 6.25 user 100% CPU (0.001 sec/iter).
M216091. Saved file at iteration 180224. Res64: DA6BE8F46193EDFC.
M216091. Saved file at iteration 184320. Res64: 4C0C54F6B3BB5100.
M216091. Saved file at iteration 188416. Res64: 406A52E8ECC1101E.
Iter. 190000 ( 87.93%), Err= 0.000, 6.26 user 100% CPU (0.001 sec/iter).
M216091. Saved file at iteration 192512. Res64: DD896660D12E8357.
M216091. Saved file at iteration 196608. Res64: BCFE65FA0F9F4778.
Iter. 200000 ( 92.55%), Err= 0.000, 6.25 user 100% CPU (0.001 sec/iter).
M216091. Saved file at iteration 200704. Res64: F4E92A3E32F77E9F.
M216091. Saved file at iteration 204800. Res64: 78897F6C7A3E982E.
M216091. Saved file at iteration 208896. Res64: 1007948E422A60B9.
Iter. 210000 ( 97.18%), Err= 0.000, 6.26 user 100% CPU (0.001 sec/iter).
M216091. Saved file at iteration 212992. Res64: 240FD3A839CB978D.
[Tue Aug 31 21:35:15 2010]
M216091 is a known Mersenne prime!
pacionet is offline   Reply With Quote