mersenneforum.org Leyland Primes (x^y+y^x primes)
 Register FAQ Search Today's Posts Mark Forums Read

2022-06-28, 16:50   #518
storm5510
Random Account

Aug 2009
Not U. + S.A.

3×751 Posts

Quote:
 Originally Posted by rogue He has a special build of xyyxsieve because the way he is testing ranges is based upon decimal length of the candidates. The "out of the box" xyyxsieve is more efficient with a relatively square workspace, i.e. the range of x and y are similar in size. At one time I was working on changing xyyxsieve to be provide the "best of both worlds", but I couldn't get it to work and gave up. I might return to it someday. The biggest challenge with xyyxsieve is that it is most efficient to sieve a very large search space once than to break that up into smaller chunks. But with a very large search space one needs a lot of memory.
On the surface, x^y + y^x does not seem all that complex. It is what surrounds it that could be very complex. Very large numbers, for one. What would qualify as a PRP candidate and what would not seems really fuzzy, to me, as each iteration would have some sort of result. I tried a few small ones with Windows Calculator. What I got was even numbers or those evenly divisible by 5. It did not take much to get into numbers larger than 1e9 either.

2022-06-28, 17:48   #519
xilman
Bamboozled!

"𒉺𒌌𒇷𒆷𒀭"
May 2003
Down not across

2×5,711 Posts

Quote:
 Originally Posted by storm5510 On the surface, x^y + y^x does not seem all that complex. It is what surrounds it that could be very complex. Very large numbers, for one.
That is why I drew them to people's attention. Simple representation, few clearly usable algebraic properties, observed high(-ish) density of (pseudo-)primes and a wide range of sizes.

2022-06-28, 18:47   #520
rogue

"Mark"
Apr 2003
Between here and the

3·17·131 Posts

Quote:
 Originally Posted by storm5510 On the surface, x^y + y^x does not seem all that complex. It is what surrounds it that could be very complex. Very large numbers, for one. What would qualify as a PRP candidate and what would not seems really fuzzy, to me, as each iteration would have some sort of result. I tried a few small ones with Windows Calculator. What I got was even numbers or those evenly divisible by 5. It did not take much to get into numbers larger than 1e9 either.
Removing terms with small divisors is trivial. One is basically removing terms where x^y = - y^x. The optimizations come from computing x^y from miny to maxy and y^x from minx to maxy. How much memory is needed is based upon the number of terms * number of primes tested per execution of an unrolled loop.

When I last tested a range for all 20000 < x < 30000 and 5000 < y < 30000, I still had over 3 million terms left at 15e9 after many days of sieving.

So if you consider that a GPU can do hundreds of primes at a time, one quickly sees that memory is going to be an issue. In fact the movement of data between GPU and CPU is the bottleneck. Still faster than CPU alone, but a noticeable bottleneck. I would like to try this on an Apple M1 since the CPU and GPU share memory. I have not finished the work to port to Metal as I haven't had time to work on it.

2022-06-28, 21:59   #521
pxp

Sep 2010
Weston, Ontario

22×3×19 Posts

Quote:
 Originally Posted by rogue So if you consider that a GPU can do hundreds of primes at a time, one quickly sees that memory is going to be an issue. In fact the movement of data between GPU and CPU is the bottleneck. Still faster than CPU alone, but a noticeable bottleneck. I would like to try this on an Apple M1 since the CPU and GPU share memory. I have not finished the work to port to Metal as I haven't had time to work on it.
Just a couple of observations from my Mac-centric perspective. First, for my iMacs I've always been very generous with my RAM, making sure of course that it is purchased from a seller other than Apple! User-installable RAM is not an option on a Mac mini so I have been more frugal there. My current two main (Intel) iMacs each have 128 GB installed. Second, when I ordered an M1 Mac mini (lured by an extra 2 cores) to supplement my existing Intel versions, I was very disappointed by its performance and didn't use it at all for many months. I forget if I was doing xyyxsieve or your pfgw64 for the test but it was likely xyyxsieve because I have been using all 8 cores of the M1 to test some million-digit Leyland candidates for primality and its performance is quite good (in spite of my recent finding that perhaps I shouldn't be using all the cores available). More problematic (for me) is that 3 of my most-recently-acquired Intel Mac minis (running 10.15.7), which were running as fast (or faster) than the 9 older ones (running 10.14.6) when I was doing the 300000-digits search, are running poorly in the million-digit search. I think the only change in those newer minis is a security update.

2022-06-28, 23:50   #522
storm5510
Random Account

Aug 2009
Not U. + S.A.

1000110011012 Posts

Quote:
 Originally Posted by paulunderwood mtsieve + OpenPFGW is my guess
I have seen ABCD and similar before. It was when I was running Riesel primes for a project four years ago. It has been so long that I do not remember how to run them. I could not get OpenPFGW to accept anything I presented to it. Then, I tried the below:

Quote:
 ABC2 $a^$b+$b^$a a: from 2300 to 2320 b: from 2300 to 2320
It accepted this. It is based on a sample from a document file. It would overwrite each output line, something I prefer it not to do. I would have to study the documentation more.

Many thanks!

2022-06-29, 00:53   #523
rogue

"Mark"
Apr 2003
Between here and the

1A1916 Posts

Quote:
 Originally Posted by storm5510 I have seen ABCD and similar before. It was when I was running Riesel primes for a project four years ago. It has been so long that I do not remember how to run them. I could not get OpenPFGW to accept anything I presented to it. Then, I tried the below: It accepted this. It is based on a sample from a document file. It would overwrite each output line, something I prefer it not to do. I would have to study the documentation more. Many thanks!
Use -Cverbose with pfgw. But sieving, even to a small depth, saves significant time.

2022-06-29, 00:55   #524
rogue

"Mark"
Apr 2003
Between here and the

11010000110012 Posts

Quote:
 Originally Posted by pxp Just a couple of observations from my Mac-centric perspective. First, for my iMacs I've always been very generous with my RAM, making sure of course that it is purchased from a seller other than Apple! User-installable RAM is not an option on a Mac mini so I have been more frugal there. My current two main (Intel) iMacs each have 128 GB installed. Second, when I ordered an M1 Mac mini (lured by an extra 2 cores) to supplement my existing Intel versions, I was very disappointed by its performance and didn't use it at all for many months. I forget if I was doing xyyxsieve or your pfgw64 for the test but it was likely xyyxsieve because I have been using all 8 cores of the M1 to test some million-digit Leyland candidates for primality and its performance is quite good (in spite of my recent finding that perhaps I shouldn't be using all the cores available). More problematic (for me) is that 3 of my most-recently-acquired Intel Mac minis (running 10.15.7), which were running as fast (or faster) than the 9 older ones (running 10.14.6) when I was doing the 300000-digits search, are running poorly in the million-digit search. I think the only change in those newer minis is a security update.
From what I have read on the M1 Macs is that it isn't clear how to get applications to run on the performance cores. I want to get a PRP testing program running on it natively. Obviously gwnum is not an option for "natively". There are a couple of options, such as proth20, which I can convert to Metal (if necessary), but that is only for k*2^n+1 numbers right now.

Last fiddled with by rogue on 2022-06-29 at 00:57

2022-06-29, 05:59   #525
storm5510
Random Account

Aug 2009
Not U. + S.A.

43158 Posts

Quote:
 Originally Posted by rogue Use -Cverbose with pfgw. But sieving, even to a small depth, saves significant time.
Using -Cverbose, the output is nearly identical with what it was when I was running Riesel's.

It appears in many posts I looked at above, everyone was using five-digit numbers. I will just run something smaller for an overnight test. Do they have something they are choosing them from?

Thanks!

2022-06-29, 13:05   #526
rogue

"Mark"
Apr 2003
Between here and the

11010000110012 Posts

Quote:
 Originally Posted by storm5510 Do they have something they are choosing them from?
Go to their project page.

2022-06-29, 14:40   #527
storm5510
Random Account

Aug 2009
Not U. + S.A.

3·751 Posts

Quote:
 Originally Posted by rogue Go to their project page.
Very well.

Going back through the pages, I found a Windows build of xyyxsievecl. I studied the area around page 35. It has been a bit tough to get it to run.

2022-06-29, 16:27   #528
rogue

"Mark"
Apr 2003
Between here and the

3·17·131 Posts

Quote:
 Originally Posted by storm5510 Very well. Going back through the pages, I found a Windows build of xyyxsievecl. I studied the area around page 35. It has been a bit tough to get it to run.
The latest build of xyyxsievecl is in the mtsieve page. This assume you have an OpenCL compatible GPU. You might need to d/l drivers.

 Similar Threads Thread Thread Starter Forum Replies Last Post carpetpool Miscellaneous Math 4 2022-07-14 02:29 Batalov XYYXF Project 57 2022-06-30 17:24 emily Math 34 2017-07-16 18:44 davar55 Puzzles 9 2016-03-15 20:55 troels munkner Miscellaneous Math 4 2006-06-02 08:35

All times are UTC. The time now is 20:57.

Wed Aug 10 20:57:38 UTC 2022 up 34 days, 15:44, 3 users, load averages: 1.28, 1.27, 1.33

This forum has received and complied with 0 (zero) government requests for information.

Permission is granted to copy, distribute and/or modify this document under the terms of the GNU Free Documentation License, Version 1.2 or any later version published by the Free Software Foundation.
A copy of the license is included in the FAQ.

≠ ± ∓ ÷ × · − √ ‰ ⊗ ⊕ ⊖ ⊘ ⊙ ≤ ≥ ≦ ≧ ≨ ≩ ≺ ≻ ≼ ≽ ⊏ ⊐ ⊑ ⊒ ² ³ °
∠ ∟ ° ≅ ~ ‖ ⟂ ⫛
≡ ≜ ≈ ∝ ∞ ≪ ≫ ⌊⌋ ⌈⌉ ∘ ∏ ∐ ∑ ∧ ∨ ∩ ∪ ⨀ ⊕ ⊗ 𝖕 𝖖 𝖗 ⊲ ⊳
∅ ∖ ∁ ↦ ↣ ∩ ∪ ⊆ ⊂ ⊄ ⊊ ⊇ ⊃ ⊅ ⊋ ⊖ ∈ ∉ ∋ ∌ ℕ ℤ ℚ ℝ ℂ ℵ ℶ ℷ ℸ 𝓟
¬ ∨ ∧ ⊕ → ← ⇒ ⇐ ⇔ ∀ ∃ ∄ ∴ ∵ ⊤ ⊥ ⊢ ⊨ ⫤ ⊣ … ⋯ ⋮ ⋰ ⋱
∫ ∬ ∭ ∮ ∯ ∰ ∇ ∆ δ ∂ ℱ ℒ ℓ
𝛢𝛼 𝛣𝛽 𝛤𝛾 𝛥𝛿 𝛦𝜀𝜖 𝛧𝜁 𝛨𝜂 𝛩𝜃𝜗 𝛪𝜄 𝛫𝜅 𝛬𝜆 𝛭𝜇 𝛮𝜈 𝛯𝜉 𝛰𝜊 𝛱𝜋 𝛲𝜌 𝛴𝜎𝜍 𝛵𝜏 𝛶𝜐 𝛷𝜙𝜑 𝛸𝜒 𝛹𝜓 𝛺𝜔