20120420, 18:46  #45  
Jun 2009
1010101011_{2} Posts 
Quote:
Mx experience is that if sieving to 1e9, the max. krange is 10T. It can take more, but it is slow. Try splitting the range in half and you should get much faster progress. 

20120420, 18:54  #46 
Apr 2010
Over the rainbow
2^{2}×5×127 Posts 
thank you, that did the trick... now the rate is at 60 k per sec, much quicker... ( test primality on several machine)
Last fiddled with by firejuggler on 20120420 at 18:55 
20120420, 20:43  #47  
Just call me Henry
"David"
Sep 2007
Cambridge (GMT/BST)
5858_{10} Posts 
Quote:
The 30/210 thing is just a speedup. The list axn posted assumed 5 1 +1 +5 +7 +11. So this cannot be used for the 7tuple combination I used. You can add to this list but cannot remove from it without missing candidates. Sometimes you can tighten it after adding more which is why we moved from 30 to 210. I have been trying to use candidates that include all the previous types as well. I also test them in this order so if I find 5 primes when searching for 6tuples then I have a 5tuple. This isn't always possible. 6 and 7 were just too different. I managed to find one for 7 that includes 5 and below though. I would suggest we search for something like 8 at the size for a record at 6. This would mean even if we don't find a 8 we find something useful. I think running separate searches at 6 and 8 would be less effficient. Maybe people at tps would like to search for larger tuples. Some of these records seem to be sitting here for the taken by a project/person with a few cpu years as far as I can see. @ puzzlepeter The site said there are no known 24tuples. This basically means there aren't any small optimal tuples which have been found. The trivial search space was exhausted before they came across the first one. Relying on the law of small numbers didn't work in this case. 

20120421, 06:03  #48 
Jun 2009
683 Posts 
Hm. I'm starting to realize just how much of the basics I'm missing.
A sixtuplet is a quadruplet with additions on each end, another 4 and +4. This is a pattern I don't see in any of the longer tuplets (well, I found one contained in a 19 *lol*). So sieving for 8 won't get any 6 as bycatch, right? But there is an 8pattern that includes a 7pattern, which is 8: 0 2 6 8 12 18 20 26 7: 0 2 6 8 12 18 20 The quadruplet pattern appears on most of the larger tuplets, which is why the quad code can be reused by adding to this pattern, correct? I'm just trying to make sure I am able to make use of everything I read here without messing things up. Thanks to everybody for your help and explanations! 
20120421, 11:47  #49 
Just call me Henry
"David"
Sep 2007
Cambridge (GMT/BST)
2·29·101 Posts 
I added a new constant tuple_size to all the source files. It basically changes the dimensions of ndx and the two loop counts. This means it is easier to do a change and is less error prone. All are attached. The only files with functional different are six and seven. I have changed seven to work mod 210 which should multiple the speed by 7 I believe.
Has anyone ever experimented with adding a prime to the presieve? Has anyone fiddled with SieveSize? 
20120421, 13:27  #50  
Jun 2009
683 Posts 
Quote:
I also tried different sievesizes, again with ambigous results. On one core it seemed to be faster when sievesize was larger than L2 cache size. Using all cores at the same time it seemed to be slower. I'm not sure if the constant disk access may have warped my timings as there were other processes needing a lot of bandwidth. I saw the iteration counter went up in irregular time intervals so I probably have to do the timings again. I fiddled around with the number of smallprimes also and saw it can be quite large and still be lucrative. Depending mainly of N of course. What's up with the 10000T limit? AFAIK qword can go higher by a factor 100. EDIT: in Quin, Six and Seven K_START is being size checked as given in G while K_END is handled as given in T. Search ranges grow rapidly with long tuples so I think T is fine. It's really fast anyway. Last fiddled with by PuzzlePeter on 20120421 at 13:39 

20120421, 15:32  #51 
Just call me Henry
"David"
Sep 2007
Cambridge (GMT/BST)
2·29·101 Posts 
k limit now increased to 10000P (10^19). There were a few arbitrary limits left in there from it being a twin prime sieve.
Both k_start and k_end are being checked in G now. I prefer G because it allows me to do small ranges for testing. I might try and work on a 5e9 sort of notation at some point. There are more changes needed for changing the presieve. 
20120421, 17:51  #52 
Jun 2009
683 Posts 

20120422, 15:40  #53 
Jun 2009
1253_{8} Posts 
Next question: the main part is do_sieve_iter. Its outer loop counts to sprime_ctr, so for each subrange sieved the runtime time should be linear with the number of small primes, right?
Decreasing SmallPrimes to a ridiculously low value like 1000 does not give very much of a speedup. Could it be that looping through the array to find the candidates left after sieving and write them to the outfile takes quite a lot of time? 
20120422, 17:32  #54  
Just call me Henry
"David"
Sep 2007
Cambridge (GMT/BST)
2·29·101 Posts 
Quote:
I am going to attempt running it with gprof. Win64 with gpof support isn't supported by fpc so trying linux. Unfortunately my ubuntu virtual machine is 9.04 and can't upgrade easily. Downloading the 12.04 Beta 2 iso now. 

20120422, 20:49  #55 
"Robert Gerbicz"
Oct 2005
Hungary
1,459 Posts 
I have written a sieve in c language for this problem, see https://sites.google.com/site/robertgerbicz/gsieve.
It should be faster than axn's version. It uses different type of sieves depending on the number of candidates and p value. Set bound_small_primes as 13, if 2*3*5*7*11*13*SmallPrimes<=end_kstart_k. You can try to use 13 earlier also. Or even higher. Play also with sieve_len. From the printed timing information and ratio of the progress you can get quickly an estimation for the completion time. Note that this sieve doesn't care about tiny ranges, it is really slow if you give say a 1 billion range. The current setup is fast if end_kstart_k>>2310*10^9. Due to the wheelsieve it doesn't print out the remaining k candidates in increasing order. Last fiddled with by R. Gerbicz on 20120422 at 20:50 
Thread Tools  
Similar Threads  
Thread  Thread Starter  Forum  Replies  Last Post 
How/Where to get Jens Kruse Andersen's prime constellation sieve?  Stargate38  And now for something completely different  2  20170428 00:08 
Efficiently finding a linear progression in data  fivemack  Math  27  20151212 18:42 
GPU Prime Sieve  tapion64  GPU Computing  7  20140410 06:15 
Sieve depth vs. prime probability  Unregistered  Information & Answers  2  20100525 20:51 
Prime in Riesel Sieve Project  Sloth  Prime Sierpinski Project  1  20060510 02:02 