View Single Post
Old 2003-03-07, 11:31   #1
Xyzzy
 
Xyzzy's Avatar
 
"Mike"
Aug 2002

1E8016 Posts
Default Glucas/Mlucas errors...

My iBook can pass all 31 Glucas self tests, but it fails immediately (bus error) when I try to run a real exponent... Mlucas fails on every self test I try... The real exponent I'm trying with Glucas is 9167489...

Glucas... http://sourceforge.net/project/showfiles.php?group_id=24518
Mlucas... http://gravity.psu.edu/~khanna/hpc.html

Here is the output from the Glucas self test... http://www.teamprimerib.com/txt/selftest.res

Here is the output from Mlucas...

[code:1][Xyzzy:~/Desktop/Mlucas] mvang% ./Mlucas
looking for worktodo.ini file...
no worktodo.ini file found...switching to interactive mode.
Enter exponent, FFT length in K (set K = 0 for default FFT length) >9490001,480
Enter 'y' to run a self-test, <return> for a full LL test >y
Enter number of iterations for timing test>
100
Enter index of radix set to be used for the FFT:
(See file fft_radix.txt for a list of available choices; enter -1 to get the default) >-1
Enter 'y' to enable per-iteration error checking, <return> for no error checking >y
p is prime...proceeding with Lucas-Lehmer test...
M( 9490001 ): using FFT length 480K = 491520 8-byte floats.
this gives an average 19.3074564615885 bits per digit
INFO: Using real* 8 for FFT sincos and DWT weights tables inits.
WARN: radix set -1 not available - using defaults.
Using complex FFT radices 15 16 32 32
M 9490001 Roundoff warning on iteration 24 maxerr = 0.406250000000
M 9490001 Roundoff warning on iteration 25 maxerr = 0.500000000000
FATAL ERROR...Halting execution.
STOP
[/code:1]
[code:1][Xyzzy:~/Desktop/Mlucas] mvang% ./Mlucas
looking for worktodo.ini file...
no worktodo.ini file found...switching to interactive mode.
Enter exponent, FFT length in K (set K = 0 for default FFT length) >87700003,4608
Enter 'y' to run a self-test, <return> for a full LL test >y
Enter number of iterations for timing test>
100
Enter index of radix set to be used for the FFT:
(See file fft_radix.txt for a list of available choices; enter -1 to get the default) >-1
Enter 'y' to enable per-iteration error checking, <return> for no error checking >y
p is prime...proceeding with Lucas-Lehmer test...
M( 87700003 ): using FFT length 4608K = 4718592 8-byte floats.
this gives an average 18.5860534244114 bits per digit
INFO: Using real* 8 for FFT sincos and DWT weights tables inits.
WARN: radix set -1 not available - using defaults.
Using complex FFT radices 9 16 16 32 32
M87700003 Roundoff warning on iteration 27 maxerr = 0.500000000000
FATAL ERROR...Halting execution.
STOP
[/code:1]
The hardware is a 800MHz PowerPC 750FX with 640MB of memory and 512KB L2...

What am I doing wrong?
Xyzzy is offline   Reply With Quote