mersenneforum.org  

Go Back   mersenneforum.org > Great Internet Mersenne Prime Search > Software > Mlucas

Reply
 
Thread Tools
Old 2003-03-07, 11:31   #1
Xyzzy
 
Xyzzy's Avatar
 
"Mike"
Aug 2002

22·3·643 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
Old 2003-03-07, 13:24   #2
kiste
 
Mar 2003

2 Posts
Default Re: Glucas/Mlucas errors...

Quote:
Originally Posted by Xyzzy
My iBook can pass all 31 Glucas self tests, but it fails immediately (bus error) when I try to run a real exponent...
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}.
kiste is offline   Reply With Quote
Old 2003-03-07, 17:00   #3
ewmayer
2ω=0
 
ewmayer's Avatar
 
Sep 2002
República de California

22·31·79 Posts
Default 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
ewmayer is offline   Reply With Quote
Old 2003-03-07, 17:30   #4
Xyzzy
 
Xyzzy's Avatar
 
"Mike"
Aug 2002

22×3×643 Posts
Default Re: Glucas/Mlucas errors...

Quote:
Originally Posted by kiste
How did you setup Glucas for testing M9167489?
Please post the contents of glucas.{ini,que}.
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...
Xyzzy is offline   Reply With Quote
Old 2003-03-07, 17:43   #5
Xyzzy
 
Xyzzy's Avatar
 
"Mike"
Aug 2002

22×3×643 Posts
Default Re: Glucas/Mlucas errors...

Quote:
Originally Posted by ewmayer
Apparently I was overly aggressive in setting the maximum exponent for the 480K FFT length - try using 512K instead.
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:
Originally Posted by 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.
Do you have a known self-test result I can use to test this range?
Xyzzy is offline   Reply With Quote
Old 2003-03-07, 18:07   #6
ewmayer
2ω=0
 
ewmayer's Avatar
 
Sep 2002
República de California

22×31×79 Posts
Default Re: Glucas/Mlucas errors...

Quote:
Originally Posted by Xyzzy
Do you have a known self-test result I can use to test this range?
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.
ewmayer is offline   Reply With Quote
Old 2003-03-07, 18:40   #7
kiste
 
Mar 2003

2 Posts
Default Re: Glucas/Mlucas errors...

Quote:
Originally Posted by Xyzzy
glucas.que
[code:1]9167489,1,0,0[/code:1]
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
kiste is offline   Reply With Quote
Old 2003-03-07, 19:29   #8
Xyzzy
 
Xyzzy's Avatar
 
"Mike"
Aug 2002

22×3×643 Posts
Default

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 ÷
Xyzzy is offline   Reply With Quote
Old 2003-03-07, 19:55   #9
Xyzzy
 
Xyzzy's Avatar
 
"Mike"
Aug 2002

22·3·643 Posts
Default Re: Glucas/Mlucas errors...

Quote:
Originally Posted by ewmayer
Try 84000013, which should yield a 100-iteration Res64 = C7323B8F92DCB618.
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]
Xyzzy is offline   Reply With Quote
Old 2003-03-07, 20:08   #10
Xyzzy
 
Xyzzy's Avatar
 
"Mike"
Aug 2002

22×3×643 Posts
Default Re: Glucas/Mlucas errors...

Quote:
Originally Posted by 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.
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...
Xyzzy is offline   Reply With Quote
Old 2003-03-07, 20:44   #11
Xyzzy
 
Xyzzy's Avatar
 
"Mike"
Aug 2002

22×3×643 Posts
Default

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 ÷
Xyzzy is offline   Reply With Quote
Reply

Thread Tools


Similar Threads
Thread Thread Starter Forum Replies Last Post
Glucas Source nuggetprime Software 13 2011-01-14 19:51
Mlucas/Glucas for LLR tests? mdettweiler Mlucas 4 2008-03-12 22:25
OS X Glucas build rtharper Software 3 2007-06-13 23:28
GLucas.... bayanne Software 5 2003-08-15 16:14
Factoring with GLucas bayanne Software 10 2003-02-07 18:25

All times are UTC. The time now is 23:36.

Fri Oct 30 23:36:52 UTC 2020 up 50 days, 20:47, 2 users, load averages: 2.14, 2.14, 2.08

Powered by vBulletin® Version 3.8.11
Copyright ©2000 - 2020, Jelsoft Enterprises Ltd.

This forum has received and complied with 0 (zero) government requests for information.

Permission is granted to copy, distribute and/or modify this document under the terms of the GNU Free Documentation License, Version 1.2 or any later version published by the Free Software Foundation.
A copy of the license is included in the FAQ.