20030307, 11:31  #1 
"Mike"
Aug 2002
2^{2}·3·643 Posts 
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 selftest, <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 periteration error checking, <return> for no error checking >y p is prime...proceeding with LucasLehmer test... M( 9490001 ): using FFT length 480K = 491520 8byte 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 selftest, <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 periteration error checking, <return> for no error checking >y p is prime...proceeding with LucasLehmer test... M( 87700003 ): using FFT length 4608K = 4718592 8byte 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? 
20030307, 13:24  #2  
Mar 2003
2 Posts 
Re: Glucas/Mlucas errors...
Quote:
Glucas selftests use the same codepath as the regular tests, but use a different set of ini/quefiles. How did you setup Glucas for testing M9167489? Please post the contents of glucas.{ini,que}. 

20030307, 17:00  #3 
∂^{2}ω=0
Sep 2002
República de California
2^{2}·31·79 Posts 
Re: Glucas/Mlucas errors...
[quote="Xyzzy"]
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 selftest, <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 periteration error checking, <return> for no error checking >y p is prime...proceeding with LucasLehmer test... M( 9490001 ): using FFT length 480K = 491520 8byte 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] Apparently I was overly aggressive in setting the maximum exponent for the 480K FFT length  try using 512K instead. If you're planning a full test of an exponent in this range, you can override the program's default FFT length by starting the run in interactive mode as above and specifying the appropriate FFT length, then after the program gets to the first checkpoint, kill the run and restart from the worktodo.ini file  the run will continue with the manually set FFT length, but will write outputs to the pXXX.stat file, as it normally does for batch mode runs. As far as your 4608K run is concerned, I'm not sure how I settled on 87.7M as the max exponent for that FFT length, but it's much too large  my current version of the code lists 83.76M as pmax for 4608K. Cheers, Ernst 
20030307, 17:30  #4  
"Mike"
Aug 2002
2^{2}×3×643 Posts 
Re: Glucas/Mlucas errors...
Quote:
[code:1]9167489,1,0,0[/code:1] glucas.ini [code:1]alternative_output_flag=2 iteration_output=1000 last_error_flag=1 qa_interim_file=100000 roundoff_check=1 save_iterations=10000 time_flag=1 verbose_flag=1[/code:1] I also tried it with no glucas.ini and it had the same error... 

20030307, 17:43  #5  
"Mike"
Aug 2002
2^{2}×3×643 Posts 
Re: Glucas/Mlucas errors...
Quote:
[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,512 Enter 'y' to run a selftest, <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 periteration error checking, <return> for no error checking >y p is prime...proceeding with LucasLehmer test... M( 9490001 ): using FFT length 512K = 524288 8byte floats. this gives an average 18.1007404327393 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 8 32 32 32 100 iterations of M 9490001 with FFT length 524288 Res64: F9C0781482193C32. Program: E2.7b Clocks = 00:01:44.646 STOP[/code:1] Quote:


20030307, 18:07  #6  
∂^{2}ω=0
Sep 2002
República de California
2^{2}×31×79 Posts 
Re: Glucas/Mlucas errors...
Quote:


20030307, 18:40  #7  
Mar 2003
2 Posts 
Re: Glucas/Mlucas errors...
Quote:
Try if the following works: [code:1]9167489 1 0 0[/code:1] Looks like we need to add improved input validation here. As of Glucas 2.9 it's better to import PrimeNet assignments, instead of editing glucas.que directly. Example: DoubleCheck=9167489,60 Save the PrimeNet assignment in a file i.e. workotodo.ini Start Glucas with the filename as argument: ./Glucas wortodo.ini 

20030307, 19:29  #8 
"Mike"
Aug 2002
2^{2}×3×643 Posts 
It looks like I can get away with a 480K FFT... If my calculations are correct, it looks like it will take ~70 days... This is with error checking on, and the default radix... I'll try Glucas next to compare the speed...
[code:1][Xyzzy:~/Desktop/mlucas] mvang% nice 19 ./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) >9167489, 480 Enter 'y' to run a selftest, <return> for a full LL test > 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 periteration error checking, <return> for no error checking >y p is prime...proceeding with LucasLehmer test... M( 9167489 ): using FFT length 480K = 491520 8byte floats. this gives an average 18.6513041178385 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( 9167489 ) iteration = 2000 clocks = 00:39:09.309. Res64: 313B5721DFABE7CD WARN: radix set 1 not available  using defaults. Using complex FFT radices 15 16 32 32 M( 9167489 ) iteration = 4000 clocks = 00:21:57.501. Res64: D9783AAA0FF5EF2D WARN: radix set 1 not available  using defaults. Using complex FFT radices 15 16 32 32 ^C[/code:1] 21 [ENTER] 60 × 57.501 + 9167489 [ENTER] 2000 ÷ × 86400 ÷ 
20030307, 19:55  #9  
"Mike"
Aug 2002
2^{2}·3·643 Posts 
Re: Glucas/Mlucas errors...
Quote:
[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) >84000013,4608 Enter 'y' to run a selftest, <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 periteration error checking, <return> for no error checking >y p is prime...proceeding with LucasLehmer test... M( 84000013 ): using FFT length 4608K = 4718592 8byte floats. this gives an average 17.8019233279758 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 100 iterations of M84000013 with FFT length 4718592 Res64: 6EA52961F9CAD418. Program: E2.7b Clocks = 00:19:38.994 STOP[/code:1] 

20030307, 20:08  #10  
"Mike"
Aug 2002
2^{2}×3×643 Posts 
Re: Glucas/Mlucas errors...
Quote:
One other glitch I found is the contents of the glucas.ini file are case sensitive... I used all lowercase and none of the lines got recognized... I had to capitalize the first word in each line... :( I'll post timings in a few minutes... 

20030307, 20:44  #11 
"Mike"
Aug 2002
2^{2}×3×643 Posts 
Looks like ~32 days...
[code:1][Xyzzy:~/Desktop/Glucas2.9.0.ppc750.MacOSXdarwin_55] mvang% nice 19 ./glucas Appended ComputerID=408208707 in file glucas.ini. [Fri Mar 7 15:02:58 2003] Going to work with exponent 9167489 Starting from iteration 1. Exponent 9167489. Iter. 1000 ( 0.01%), Err= 0.074, 280.94 user 92% CPU (0.304 sec/iter). Iter. 2000 ( 0.02%), Err= 0.074, 281.95 user 93% CPU (0.302 sec/iter). Iter. 3000 ( 0.03%), Err= 0.066, 281.46 user 94% CPU (0.301 sec/iter). Iter. 4000 ( 0.04%), Err= 0.082, 281.49 user 93% CPU (0.302 sec/iter). Iter. 5000 ( 0.05%), Err= 0.062, 282.02 user 94% CPU (0.301 sec/iter). Iter. 6000 ( 0.07%), Err= 0.074, 281.92 user 94% CPU (0.301 sec/iter). Iter. 7000 ( 0.08%), Err= 0.070, 281.73 user 94% CPU (0.300 sec/iter).[/code:1] 9167489 [ENTER] .300 × 86400 ÷ 
Thread Tools  
Similar Threads  
Thread  Thread Starter  Forum  Replies  Last Post 
Glucas Source  nuggetprime  Software  13  20110114 19:51 
Mlucas/Glucas for LLR tests?  mdettweiler  Mlucas  4  20080312 22:25 
OS X Glucas build  rtharper  Software  3  20070613 23:28 
GLucas....  bayanne  Software  5  20030815 16:14 
Factoring with GLucas  bayanne  Software  10  20030207 18:25 