![]() |
|
|
#100 |
|
Oct 2006
26010 Posts |
Great results Robert! Can you post the excel spreadsheet you used for reference?
Is there a way to sieve and test individual n values using the batch files that were posted here? For example, if I sieved and tested n-values to k=50M, and still had some left over, could I plug in the n-values I need into the batch files for it to sieve for 50M<k<100M? Thanks! |
|
|
|
|
|
#101 | |
|
Jun 2003
Suva, Fiji
2·1,021 Posts |
Quote:
To recreate: Download the data posted in this thread, paste into an excel sheet, use "data" "text to columns" to separate info into columns, using "space" as the separator. Construct formula using =median() over 50 values, construct another column next to it using the formula cell*2^cell. Create a graph using chart wizard, with the choice "xy scatter", set the font size to 2 point and you should be able to construct the graph easily I then screenprinted the graph and created a .jpg in Irfan jpg viewer before posting here Last fiddled with by robert44444uk on 2008-04-06 at 12:45 |
|
|
|
|
|
|
#102 | ||
|
"Gary"
May 2007
Overland Park, KS
23·29·53 Posts |
Quote:
1. On the main 'page', do everything like you would normally with a couple of exceptions: (a) The n-value will be the beginning of your n-range. (b) The output file name will be simply ".txt". (explanation as to why below) 2. Go to 'options' 'Sieve until' and do the following: (a) Select your 'Sieve up to p=' limit like you would normally. (Needs to be determined ahead of time.) (b) Go to the part that says 'When complete, automatically increment...', check the box on the left, and change the box on the right to 1. (c) On the next line for the stop point, simply tell it what n-value you want to stop at. You'll want to determine LLR times for the highest k in your range ahead of time to determine the optimal removal rate needed. When there is a jump in fftlen is when you'll want to increase your sieve depth. In other words, set 2(c) above to the highest n-value right before an fftlen change. One note: Because it is sieving a much smaller range of n, you will get to optimal depth at a much lower P-value because it is only removing candidates from that single-n value at a time. It can still only sieve one n at a time but this allows the process to be automated. This means that you will have to do more LLR tests than for a comparible single-n search BUT those LLR tests will be so much faster that you'll still be saving a lot of testing time overall. k=10G-100G takes a very long time to LLR, far longer than k=1M-3M. On the file names, Newpgen will create one file per n-value. It will give you a pop up that tells you what it will call them. NewPGen's default file name is: (n-value)_(the name you put in the output file name). So...by putting ".txt" in the output file name for a range starting at n=50K, you'll get file names of 50000_.txt, 50001_.txt, 50002_.txt, etc. I personally sieve 1000n at a time for my 'all twin prime search' here. I then go to the DOS prompt and copy all of the files into one big 1000n file. DOS is good about copying them into one big file in the correct numerical order as long as you allow Newpgen to name them numerically as specified above. Also, LLR ignores multiple instances of the header in the file. The first one is all that it pays attention to. That's all there is to it. It's extremely easy after you've done it once. Quote:
Also, I would greatly appreciate being the one to start testing above n=66K for a top-20 twin for k<1M. My goal with this effort is to find all of the twins up to the top-20 table and then continue on up to n=100K. Or if someone wants to coordinate with me on such an effort, that would be fine too and we could share any twin found. k>1M is fair game. I have no intention of testing any k's > 1M for twins other than what I already did with my n=100K twin. Gary |
||
|
|
|
|
|
#103 |
|
Oct 2006
22×5×13 Posts |
Here's a .xls chart of the data on Gary's website.
I notice of interest the swung s shape on the gap chart, and the sudden jumps in smallest twin k-values n-values of multiples of 1000. Perhaps there are smaller twins in some of these areas as a smooth curve makes a bit more sense to me than these jumps. (These are more evident when the graph is not in logarithmic mode) Is there a way to use the batch files posted earlier here for specific n-values, instead of a range? Thanks! |
|
|
|
|
|
#104 | |
|
Mar 2006
Germany
2·1,531 Posts |
Quote:
call the second batch (called 'do_one.bat') with the k you want to test, for example do_one 20000 to find the first twin for k=20000. make sure to set the kmax and pmax you need in this batch first. |
|
|
|
|
|
|
#105 | |
|
"Michael Kwok"
Mar 2006
22358 Posts |
Quote:
After we find a twin for n=333,333, a proposed n-range for finding our next twin is from n=450,001 to n=524,288. The k-range is from k=3 to k=2.8M, and there is a > 90% chance that a twin will be found in this range. Anyway, I found that the optimal sieve depth was about 3T. However, we won't sieve that high because 3T * 74,287 n-values = 222,861T to sieve, which is way too much. As a comparison, we've only sieved a total of 4772T for n=333,333. A more realistic sieve depth for that k and n range is 0.1T. That means that there will be a total of 7429T to sieve, which is reasonable compared with the previous sieve effort of 4772T. |
|
|
|
|
|
|
#106 | |
|
Oct 2006
4048 Posts |
Quote:
Thanks! EDIT: I think by k=20000 you mean n=20000. Last fiddled with by roger on 2008-04-09 at 23:54 |
|
|
|
|
|
|
#107 |
|
Mar 2006
Germany
1011111101102 Posts |
(yes: k<->n, i mixed them up)
to make specific values of k you can do this: create a file named 'do.bat' and put the following lines in: Code:
call do_one 9000 call do_one 9029 call do_one 9052 call do_one 9123 call do_one 9211 Last fiddled with by kar_bon on 2008-04-10 at 06:00 |
|
|
|
|
|
#108 | |
|
Jun 2003
Suva, Fiji
2×1,021 Posts |
Quote:
This is because the Gary's website grid has been incorrectly transposed, best instead to the take the first 1000 values that I posted in this thread. The jumpy pattern of the first 1,000 n then disappears I will be posting the 5,000 - 6,000 probably tomorrow. Then I will plan to do the test at 30-31000 n, ignoring the first 1,000,000 k - will keep timings Last fiddled with by robert44444uk on 2008-04-10 at 12:27 |
|
|
|
|
|
|
#109 | |
|
Oct 2006
22·5·13 Posts |
Quote:
![]() I'm looking to test specific n-values, not k-values in k*2^n+-1 |
|
|
|
|
|
|
#110 |
|
Mar 2006
Germany
57668 Posts |
|
|
|
|
![]() |
Similar Threads
|
||||
| Thread | Thread Starter | Forum | Replies | Last Post |
| Sieving with powers of small primes in the Small Prime variation of the Quadratic Sieve | mickfrancis | Factoring | 2 | 2016-05-06 08:13 |
| Relativistic Twins | davar55 | Science & Technology | 68 | 2015-01-20 21:01 |
| 3x*2^n-1 and 3x*2^n-1 possibly twins ? | science_man_88 | Riesel Prime Search | 10 | 2010-06-14 00:33 |
| The Twins | GP2 | Lounge | 1 | 2003-11-18 04:50 |
| NOT twins | graeme | Puzzles | 11 | 2003-09-04 00:41 |