20030307, 21:58  #12  
∂^{2}ω=0
Sep 2002
República de California
2·13·443 Posts 
Re: Glucas/Mlucas errors...
Quote:
Using the C version of 2.7c, I have no trouble running the selftest for your 9.49M exponent at 480K FFT length: [code:1] M9490001: using FFT length 480K = 491520 8byte floats. this gives an average 19.307456461588540 bits per digit INFO: Using real*16 for FFT sincos and DWT weights tables inits. WARN: radix set 1 not available  using defaults. Using complex FFT radices 15 16 32 32 100 iterations of M9490001 with FFT length 491520 Res64: 01A4E738255C522B. AvgMaxErr = 0.215053013. Program: E2.7c Clocks = 00:00:05.533 [/code:1] Also, to compare timingtest results between the f90 and C code, you need to run the f90 code for one more iteration than the C, or the C code for one less iteration than the f90. When I run the C/v2.7c code for 99 iterations, I get the same as your "100iteration" result for v2.7b: [code:1] M84000013: using FFT length 4608K = 4718592 8byte floats. this gives an average 17.801923327975803 bits per digit INFO: Using real*16 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 99 iterations of M84000013 with FFT length 4718592 Res64: 6EA52961F9CAD418. AvgMaxErr = 0.107991536. Program: E2.7c Clocks = 00:01:09.213 [/code:1] Using the moreaccurate C version, I can also run all the way up to 87.7M at 4608K (though the RO warnings and the AvgMaxErr > 0.3 indicate this is too large an exponent for a fulllength run at this FFT length  p ~= 87.0M is probably the maximum one could safely use): [code:1] M87700003: using FFT length 4608K = 4718592 8byte floats. this gives an average 18.586053424411350 bits per digit INFO: Using real*16 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 M 87700003 Roundoff warning on iteration 77 maxerr = 0.437500000000 M 87700003 Roundoff warning on iteration 91 maxerr = 0.406250000000 100 iterations of M87700003 with FFT length 4718592 Res64: 3C41DA97B5B87409. AvgMaxErr = 0.323716518. Program: E2.7c Clocks = 00:01:09.997 [/code:1] 

20030307, 22:17  #13  
"Mike"
Aug 2002
7,691 Posts 
Re: Glucas/Mlucas errors...
Quote:
If you want a SSH account on my box to work on a Mac client let me know... :) 

20160327, 16:24  #14 
Aug 2010
Republic of Belarus
2·5·17 Posts 
Hello. Oups. What happened with Mlucas?
I got this on IBM Power8 and Red Hat 7.2. Mlucas installed without problem ... I hope. Because now not sure. Code:
[Mar 02 21:59:06] M43892801 Iter# = 200000 clocks = 00:16:24.308 [ 0.0984 sec/iter] Res64: A82FA0F649C8F274. AvgMaxErr = 0.170602632. MaxErr = 0.250000000 [Mar 02 22:20:42] M43892801 Iter# = 210000 clocks = 00:21:35.630 [ 0.1296 sec/iter] Res64: BAFC8160194E1EFA. AvgMaxErr = 0.000000000. MaxErr = 0.000000000 
20160327, 21:08  #15  
∂^{2}ω=0
Sep 2002
República de California
2·13·443 Posts 
Quote:
This feature was one of the first things I ditched in my current v16.0 codedev plan ... too dangerous, as ROE warnings can catch not just dangerously high errors but many forms of data corruption. Your results should be fine, but I suggest you manually go into your Mlucas.c source file, change line 342 from relax_err=TRUE; to relax_err=FALSE; recompile that file, relink your binary and halt/restart your run. 

20160328, 20:09  #16 
Aug 2010
Republic of Belarus
2×5×17 Posts 
ewmayer, Thank you very much!
Now it working as before. Very strange with speed. Hope i will verify LLtest. Code:
[Mar 28 11:29:42] M43892801 Iter# = 17350000 clocks = 00:16:25.361 [ 0.0985 sec/iter] Res64: 9B98340B8DDCA51F. AvgMaxErr = 0.170527673. MaxErr = 0.226562500 [Mar 28 11:46:07] M43892801 Iter# = 17360000 clocks = 00:16:24.359 [ 0.0984 sec/iter] Res64: C4197B369FA08B12. AvgMaxErr = 0.170579383. MaxErr = 0.250000000 [Mar 28 12:02:32] M43892801 Iter# = 17370000 clocks = 00:16:25.449 [ 0.0985 sec/iter] Res64: B61947BA0900E712. AvgMaxErr = 0.170759450. MaxErr = 0.250000000 [Mar 28 12:18:57] M43892801 Iter# = 17380000 clocks = 00:16:24.465 [ 0.0984 sec/iter] Res64: 964B7F7B7BE64D97. AvgMaxErr = 0.170774248. MaxErr = 0.234375000 [Mar 28 12:35:23] M43892801 Iter# = 17390000 clocks = 00:16:25.531 [ 0.0986 sec/iter] Res64: 285FDC8DAC1094B6. AvgMaxErr = 0.170833351. MaxErr = 0.218750000 [Mar 28 12:51:48] M43892801 Iter# = 17400000 clocks = 00:16:24.742 [ 0.0985 sec/iter] Res64: B1C09B8BBF18AF65. AvgMaxErr = 0.170672593. MaxErr = 0.250000000 
20160329, 01:45  #17  
∂^{2}ω=0
Sep 2002
República de California
2×13×443 Posts 
Quote:
Still can't believe removal of the extra operations used in the roundoff check caused such a huge slowdown  your compiler must be doing some really bizarre 'optimizations'. 

20160507, 08:07  #18 
Aug 2010
Republic of Belarus
2×5×17 Posts 
Code:
[May 06 00:24:49] M43892801 Iter# = 41170000 clocks = 00:02:29.151 [ 0.0149 sec/iter] Res64: 853F7CB9D2160D0B. AvgMaxErr = 0.199261261. MaxErr = 0.281250000 [May 06 00:27:19] M43892801 Iter# = 41180000 clocks = 00:02:29.457 [ 0.0149 sec/iter] Res64: 038A7FDF252DD6C4. AvgMaxErr = 0.199309877. MaxErr = 0.281250000 There are Ok? Or MaxErr shouldn't be greater than 0.25? 
20160507, 16:59  #19 
P90 years forever!
Aug 2002
Yeehaw, FL
7,151 Posts 
MaxErr cannot exceed 0.5. You are looking good.

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 