"Brian J. Beesley" wrote:
> I've posted on my ftp server a file which may be of interest to
> anyone developing or verifying the operation of LL testing programs.
>
> The URL of the file is
> ftp://lettuce.edsc.ulst.ac.uk/gimps/PrimeQA/QADATA.TXT
>
> The format is
> exponent,# iterations,residual,program id,max error
>
> The residual listed is the hexadecimal representation of the low
> orfer 64 bits of the complete residual. The program id is always
> "lucdwt"; the max error can be ignored (basically it's there as a
> self-check to make sure that a sensible FFT size was used to generate
> the data).
>
> The range of exponents is selected exponents between 8,000 and
> 80,000,000. The low 64 bits of the residual after 400 iterations is
> listed for every exponent in the file. (After one iteration the
> residual should be 000000000000000E for every exponent). For
> exponents less than 4,000,000, the residual after 1,000 iterations is
> also listed. For exponents greater than 20,000,000, the residual
> after 100 iterations is also listed.
>
> The results were generated on a 533 MHz Alpha 21164LX system; the
> program used is a version of lucdwt (included in Richard Crandall's
> giantint package), modified to produce output in a suitable format.
>
> I hope that this data assists anyone developing or modifying LL
> testing software to ensure that their code is operating correctly.
I've got some patches to Will's mers package to print out results that
can be compared with this. For example, here's the first few lines
of QADATA.TXT:
64511,100,EBA1B280B7C00FF7,lucdwt,0.000107
64511,200,96F428A92412ADA8,lucdwt,0.000122
64511,300,38498B59F3A614B7,lucdwt,0.000122
64511,400,3389BC878321980A,lucdwt,0.000122
64513,100,D24F0966C2E9B3C4,lucdwt,0.000092
64513,200,2376508C726DFB47,lucdwt,0.000092
64513,300,D4BA71A0879202D5,lucdwt,0.000092
64513,400,FF4069B4F7F0BDA4,lucdwt,0.000095
and the start of a test run of MacLucasUNIX running on an UltraSPARC:
M( 64511 )C, 0xeba1b280b7c00ff7, iter = 100, n = 4096, MacLucasUNIX v6.20
Sweeney
M( 64511 )C, 0x96f428a92412ada8, iter = 200, n = 4096, MacLucasUNIX v6.20
Sweeney
M( 64511 )C, 0x38498b59f3a614b7, iter = 300, n = 4096, MacLucasUNIX v6.20
Sweeney
M( 64511 )C, 0x3389bc878321980a, iter = 400, n = 4096, MacLucasUNIX v6.20
Sweeney
M( 64513 )C, 0xd24f0966c2e9b3c4, iter = 100, n = 4096, MacLucasUNIX v6.20
Sweeney
M( 64513 )C, 0x2376508c726dfb47, iter = 200, n = 4096, MacLucasUNIX v6.20
Sweeney
M( 64513 )C, 0xd4ba71a0879202d5, iter = 300, n = 4096, MacLucasUNIX v6.20
Sweeney
M( 64513 )C, 0xff4069b4f7f0bda4, iter = 400, n = 4096, MacLucasUNIX v6.20
Sweeney
The tests are still running - when it is finished I'll write a little
perl script that compares the two files and prints out any error
messages.
Simon.
_________________________________________________________________
Unsubscribe & list info -- http://www.scruz.net/~luke/signup.htm
Mersenne Prime FAQ -- http://www.tasam.com/~lrwiman/FAQ-mers
