20181030, 12:24  #1 
Jun 2003
Oxford, UK
13×149 Posts 
Primorial offsets
This thread will capture any work anyone is carrying out to generate large primes gaps using Chinese Remainder Theorem (CRT) offsets generated using ATH's prime interval software.
The thinking behind this is that, for any range of integers there are some ranges where there are relatively few members of the range without small prime factors. ATH's software looks, as a start point, to define a range of integers, (called an interval in the software) and then to look at possible modular values for the primes p in a primorial p# to determine those which generate low numbers of members with no factor <=p. Using CRT, those combinations of modular values can provide an integer, referred to as the offset, which if added or subtracted from the primorial, provides a starting point for an interval with these properties. If we multiply A, an integer, by the primorial, then add or subtract the offset, the properties of the interval demonstrate similar modular properties. Hence, for any given primorial and chosen offset, it is possible to check many gaps by altering A alone. The relative values of p that are the most efficient for any interval are yet to be determined. But if p is too small, then the chances of any A providing all members with composites (i.e. maximising the theoretical benefits of the interval) are very low. If p is too high, then the chances of meeting this criterion is quite high, but the absolute size of the members of the interval is also high, lessening the chance that the resulting prime gap is not a record. To ensure decent coordination, mersenneforum members may reserve an interval size here. Interval sizes can be any value but for convenience sake, the intervals that may be reserved are 50 apart, starting at 1000. Participating members should report any or all of three things: 1. results in aggregate coming from using the ATH software  for example the numbers of CRT offsets that provide the lowest numbers of composites for a given interval. For example, for the 2000 interval, carrying out a comprehensive test of all possible modular results for primes up to 37, there were 16 combinations providing solutions of 275, 166 providing 276, and 1806 at 277. Where a "solution" is the number of members of the interval with no factors smaller than or equal to 37 2. Any CRT offsets that are useful but which the member has no further interest in  typically this will be any offset that provides a solution within three of the best found. Three pieces of information are required  the interval, p, and the CRT solution. 3. Any records posted to Dr Nicely's site. 4. ATH output files can be posted if the member wants to release a range. The latest version of ATH's software, written in c is to be found here: https://www.mersenneforum.org/showpo...&postcount=114 A working piece of software written in perl, developed for another purpose by danaj, is to be found here: https://www.mersenneforum.org/showpo...&postcount=103 Reservations are in post #2 below 
20181030, 12:28  #2 
Jun 2003
Oxford, UK
13·149 Posts 
This post provides for reservations and results
A = Interval B = Reserved by C = Records sent to Dr Nicely Code:
A B C 1000 1050 1100 1150 1200 1250 1300 robert44444uk 1 1350 1400 1450 1500 1550 1600 1650 robert44444uk 1700 1750 1800 1850 1900 1950 2000 robert44444uk 2050 2100 2150 2200 robert44444uk 47 2250 2300 2350 2400 2450 2500 2550 2600 2650 2700 2750 2800 2850 2900 2950 3000 robert44444uk 3050 3100 3150 3200 3250 3300 3350 3400 3450 3500 3550 3600 3650 3700 3750 3800 3850 3900 3950 4000 robert44444uk 4050 4100 4150 4200 4250 4300 4350 4400 4450 4500 4550 4600 4650 4700 4750 4800 4850 4900 4950 5000 Last fiddled with by robert44444uk on 20181031 at 12:08 
20181031, 12:02  #3 
Jun 2003
Oxford, UK
13·149 Posts 
I have run 2000 and 4000 up to 79# and 101# respectively. I had 16 useful offsets that provide 197 positions within the 2000 interval where there are no prime factors smaller than 79, and 165 offsets that provide 405 positions within the 4000 interval where there are no prime factors less than 101.
I realise now that 4000 is probably too large an interval for ATH software to handle  I would need to run up to 137# to be competitive with the records. Reserving 3000 Last fiddled with by robert44444uk on 20181031 at 12:17 
20181031, 16:50  #4 
Einyen
Dec 2003
Denmark
7×449 Posts 
I can extend it to include higher primes than 101, but the problem is calculating CRT offset. It would exceed 128 bits so I would have to switch to GMP library to calculate it instead of the 128 bit variables I use now, and then you need to compile the GMP library first in order to compile the program
It should be no problem to compile GMP in MSYS2, I do it all the time, and you have all the needed packages installed, but I do not know if you feel it is worth it. The higher you go using only the best few solutions from the lower runs, the higher risk there is that you are no where near the real minimum. 
20181101, 08:02  #5  
Jun 2003
Oxford, UK
13×149 Posts 
Quote:
I agree that we are nowhere near the real minimum at higher primorials, but then again, I think a few tweaks to the program would allow for pretty decent offsets. For example if programmed to take results that are within 5 of the record at a given level, instead of 3, and then advanced those two primes at a time, then you could quickly get good results. To do that, a tweak would be required to eliminate from the results file those first few solutions that come out that are nowhere near the minimum. You could do that by calculating the theoretical expected result given an input file and only outputting those which are within 5. What do you think? Took one offset from 2000 interval forward last night and 4 records this morning between 2600 and 2750 :) Last fiddled with by robert44444uk on 20181101 at 08:04 

20181101, 09:03  #6 
Jun 2003
Oxford, UK
13·149 Posts 
If the CRT calculation is time intensive then maybe the program should only calculate it if the result is within 23 of the best solution as it is unlikely that people would look at solutions that are far way from the best.
Last fiddled with by robert44444uk on 20181101 at 09:05 
20181122, 08:28  #7 
Jun 2003
Oxford, UK
13·149 Posts 
I've been concentrating on two prime gap ranges (20003000 and 60008000) where it is clear there is relatively low hanging fruit in terms of records. It is surprising how weak the range 20003000 is. I think this partly due to gapcoin's efforts, which appear as a very clear line starting at 5000  clearly smaller gaps were not that interesting to them  and Spielaur, who has clearly invested a lot of computer resources in the 15002000 range but perhaps less after that.
Last fiddled with by robert44444uk on 20181122 at 08:33 
20181129, 08:40  #8 
Jun 2003
Oxford, UK
791_{16} Posts 
There are some fairly big gaps in terms of merit coming from this study, these are from the last two weeks and have already been submitted, The format is the simplified version of that used on the Nicely site, i,e. gap, merit, number of decimal digits of the lower of two primes defining the gap, and the detail of that prime, # representing the primorial. The value following the + is the offset developed by ATH's program.
The last gap was the smallest record found. Code:
2982 32.39 40 3012136000*79#+2903102709744756044297448892452178 2368 32.26 32 39771103997*59#+1710142320516174370411652 8646 32.09 118 899511684*269#+59573743018125804909853587106173441453569059671308281347688363825715194782198214520873698524917753936040155912 8392 31.23 117 418821175*269#+58966551729571283794666649805410544398324935677917064024228678474697709088842192419726944457327554156922052924 8150 30.78 115 8497670*269#+58966551729571283794666649805410544398324935677917064024228678474697709088842192419726944457327554156922054438 2800 30.01 41 10443985146*79#+29654909871668022980835590278211998 3708 29.82 54 3537439094*109#+2431542602822126822439136112423298092582872353046 7838 29.34 117 86297041*269#+59573743018125804909853587106173441453569059671308281347688363825715194782198214520873698524917753936040152692 7870 29.33 117 284758679*269#+74923557143241058066754483986704478908429625542794075459695816298291881932291735831576463152954478323396753506 2178 29.32 33 94343930124*59#+1067076873163923483832 2732 29.22 41 12351800046*79#+29654909871668022980835590278211184 7810 29.19 117 139673464*269#+52406140972143071826528924557389853249351158030171732299242680240432650872235779255887173651550099391722652912 7650 29.12 115 1070393*269#+74923557143241058066754483986704478908429625542794075459695816298291881932291735831576463152954478323396754364 3920 29.11 59 7708801521*127#+11111323287912892558856522182104492640353031648951378 …. 2012 27.10 33 90267121489*59#+10670768731639234831216 Last fiddled with by robert44444uk on 20181129 at 08:46 
Thread Tools  
Similar Threads  
Thread  Thread Starter  Forum  Replies  Last Post 
Largest Known Primorial  a1call  Miscellaneous Math  11  20161214 21:35 
Primorial calculation  FreakyPotato  Programming  7  20150206 10:33 
primorial primes  jasong  Math  1  20060812 01:38 
Primorial puzzle  Citrix  Puzzles  3  20060307 15:07 
Primorial question  Dougy  Math  2  20050728 13:13 