Thread: Combined Sieve
View Single Post
Old 2020-12-14, 13:21   #26
henryzz
Just call me Henry
 
henryzz's Avatar
 
"David"
Sep 2007
Cambridge (GMT/BST)

133418 Posts
Default

I am getting an error when I try to get the plots:

Code:
david@DESKTOP-BVMNJE2:/mnt/e/Mersenne/primegaps/sethsieve/build/prime-gap$ ./gap_test.py --unknown-filename 9973_208110_1_1999999_s100000_l1000000M.txt --no-one-side-skip --num-plots 3 --prp-top-percent 100
K = 14260 bits, 4293 digits, log(K) = 9883.75
Min Gap ~= 118604 (for merit > 12.0)

        Loaded 80726 records (112558 to 345816) from 'gaps.db'
        Min merit for record 11.39, gapsize for likely record 15.22 150448
        Loading unknowns from '9973_208110_1_1999999_s100000_l1000000M.txt'

Found 4584 (1.0%) results
        expect 878 left, 0.878% of SL=100000 after 1000000M
        0.010% of 4293 digit numbers are prime
        0.498% of tests should be prime (49.2x speedup)
        ~2x200.9 = 401.7 PRP tests per m
        sieve_length=100000 is insufficient ~1.26% of time


Starting m(456682) 1 to 1999999

Testing 456682 m where prob(record) >= 0

        Thread 0 started

^CReceived first  Ctrl+C | Waiting for current work to finish
No more work pushing NONE
Waiting on 1 of 1 results
        Thread 0 stopped
         61  879 <- unknowns ->  887    2160 <- gap -> 42062
            tests      1         (67.9 secs/test, 64.6 secs/test)  68, 65 secs
            unknowns   1766      (avg: 1766.00), 99.12% composite  49.77% <- % -> 50.23%
            prp tests  259       (avg: 129.50/side, 259.00/m) (3.816 tests/sec)
            sum(prob_minmerit):  1.7154e-06, 0.00218/day        found: 0
            sum(prob_record):    0.010897, 13.9/day     found: 0
            merit      4.472     (at m=61)
Joining work_q (should be instant)
Joining 1 processes
Done!
experimental_gap size mismatch 456682 787

R^2 for expected gap: 0.001, gap = 156430.5 + -0.765 * expected

|P(gap)| = 6996, Sum(P(gap)) = 450916.2
|P(gap)| = 48724, Sum(P(gap)) = 36292.3
gap > min_merit(12)  | sum(P) = 3923.5862
gap = record         | sum(P) = 4997.0288
Traceback (most recent call last):
  File "./gap_test.py", line 520, in <module>
    prime_gap_test(args)
  File "./gap_test.py", line 506, in prime_gap_test
    gap_test_plotting.plot_stuff(
  File "/mnt/e/Mersenne/primegaps/sethsieve/build/prime-gap/gap_test_plotting.py", line 439, in plot_stuff
    stats_plots(
  File "/mnt/e/Mersenne/primegaps/sethsieve/build/prime-gap/gap_test_plotting.py", line 312, in stats_plots
    verify_no_trend(valid_m, d)
  File "/mnt/e/Mersenne/primegaps/sethsieve/build/prime-gap/gap_test_plotting.py", line 43, in verify_no_trend
    trend, _ = np.polyfit(m_values, d, 1)
  File "<__array_function__ internals>", line 5, in polyfit
  File "/home/david/.local/lib/python3.8/site-packages/numpy/lib/polynomial.py", line 603, in polyfit
    raise TypeError("expected x and y to have same length")
TypeError: expected x and y to have same length
Am I doing something wrong? I haven't been able to work out how to call the plots without running at least 1 m and stopping the script either?


edit: What happens if you get a record or first occurrence gap below min_merit?

Last fiddled with by henryzz on 2020-12-14 at 15:03
henryzz is online now   Reply With Quote