- **Mlucas**
(*https://www.mersenneforum.org/forumdisplay.php?f=118*)

- - **Glucas/Mlucas errors...**
(*https://www.mersenneforum.org/showthread.php?t=420*)

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? |

Re: Glucas/Mlucas errors...[quote="Xyzzy"]My iBook can pass all 31 Glucas self tests, but it fails immediately (bus error) when I try to run a real exponent...[/quote]
That's strange. Glucas selftests use the same codepath as the regular tests, but use a different set of ini/que-files. How did you setup Glucas for testing M9167489? Please post the contents of glucas.{ini,que}. |

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 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] 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 |

Re: Glucas/Mlucas errors...[quote="kiste"]How did you setup Glucas for testing M9167489?
Please post the contents of glucas.{ini,que}.[/quote] glucas.que [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... |

Re: Glucas/Mlucas errors...[quote="ewmayer"]Apparently I was overly aggressive in setting the maximum exponent for the 480K FFT length - try using 512K instead.[/quote]
It looks like that fixed it... [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 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 512K = 524288 8-byte 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="ewmayer"]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.[/quote] Do you have a known self-test result I can use to test this range? |

Re: Glucas/Mlucas errors...[quote="Xyzzy"]
Do you have a known self-test result I can use to test this range?[/quote] Try 84000013, which should yield a 100-iteration Res64 = C7323B8F92DCB618. Some notes I mad to myself abot this self-test seem to indicate that (based on the max. average RO error I saw for this test) 84M may actually be overconservative for 4608K, but in any event it shouldn't cause you any problems running it. |

Re: Glucas/Mlucas errors...[quote="Xyzzy"]
glucas.que [code:1]9167489,1,0,0[/code:1] [/quote] Apparently the docs aren't clear enough. The fields in glucas.que are supposed to be seperated by whitespace. 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 |

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 self-test, <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 per-iteration error checking, <return> for no error checking >y p is prime...proceeding with Lucas-Lehmer test... M( 9167489 ): using FFT length 480K = 491520 8-byte 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 ÷ |

Re: Glucas/Mlucas errors...[quote="ewmayer"]Try 84000013, which should yield a 100-iteration Res64 = C7323B8F92DCB618.[/quote]
I get a different Res64... [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 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( 84000013 ): using FFT length 4608K = 4718592 8-byte 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] |

Re: Glucas/Mlucas errors...[quote="kiste"]Apparently the docs aren't clear enough. The fields in glucas.que are supposed to be seperated by whitespace.
Try if the following works: [code:1]9167489 1 0 0[/code:1] Looks like we need to add improved input validation here.[/quote] The corrected glucas.que works... 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... |

Looks like ~32 days...
[code:1][Xyzzy:~/Desktop/Glucas-2.9.0.ppc750.MacOSX-darwin_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 ÷ |

All times are UTC. The time now is 21:38. |

Powered by vBulletin® Version 3.8.11

Copyright ©2000 - 2020, Jelsoft Enterprises Ltd.