PrimeGrid
Please visit donation page to help the project cover running costs for this month

Toggle Menu

Join PrimeGrid

Returning Participants

Community

Leader Boards

Results

Other

drummers-lowrise

Advanced search

Message boards : Generalized Fermat Prime Search : Z-transform?

Author Message
Christopher Siegert
Send message
Joined: 5 Jul 09
Posts: 233
ID: 42986
Credit: 121,873,589
RAC: 0
321 LLR Amethyst: Earned 1,000,000 credits (1,041,657)Cullen LLR Amethyst: Earned 1,000,000 credits (1,366,558)PPS LLR Sapphire: Earned 20,000,000 credits (22,481,912)PSP LLR Amethyst: Earned 1,000,000 credits (1,521,444)SoB LLR Amethyst: Earned 1,000,000 credits (1,751,070)SR5 LLR Amethyst: Earned 1,000,000 credits (1,100,031)SGS LLR Amethyst: Earned 1,000,000 credits (1,014,953)TRP LLR Ruby: Earned 2,000,000 credits (3,647,371)Woodall LLR Amethyst: Earned 1,000,000 credits (1,479,507)Cullen/Woodall Sieve Bronze: Earned 10,000 credits (41,181)PPS Sieve Emerald: Earned 50,000,000 credits (74,260,244)TRP Sieve (suspended) Silver: Earned 100,000 credits (206,682)GFN Jade: Earned 10,000,000 credits (11,960,980)
Message 77210 - Posted: 14 Jun 2014 | 0:53:44 UTC

"all new genefer using a faster 'Z-transform'"

So what's the deal with this new Z-transform?

Is this new math or just a new programming algorithm?

What are the ramifications of this new Z-transform? Will it have an effect on the search for GFPs of higher N/B values? How much faster is it?

My curiosity is raging...

Profile rroonnaallddProject donor
Volunteer developer
Volunteer tester
Avatar
Send message
Joined: 3 Jul 09
Posts: 1213
ID: 42893
Credit: 34,634,263
RAC: 0
321 LLR Silver: Earned 100,000 credits (101,692)Cullen LLR Silver: Earned 100,000 credits (104,876)ESP LLR Silver: Earned 100,000 credits (101,979)PPS LLR Silver: Earned 100,000 credits (148,018)PSP LLR Silver: Earned 100,000 credits (140,441)SoB LLR Silver: Earned 100,000 credits (119,475)SR5 LLR Silver: Earned 100,000 credits (120,939)SGS LLR Silver: Earned 100,000 credits (122,783)TRP LLR Silver: Earned 100,000 credits (100,115)Woodall LLR Silver: Earned 100,000 credits (107,459)321 Sieve (suspended) Silver: Earned 100,000 credits (202,757)Cullen/Woodall Sieve Turquoise: Earned 5,000,000 credits (6,908,135)PPS Sieve Sapphire: Earned 20,000,000 credits (25,450,104)Sierpinski (ESP/PSP/SoB) Sieve (suspended) Silver: Earned 100,000 credits (130,966)TRP Sieve (suspended) Silver: Earned 100,000 credits (201,525)AP 26/27 Silver: Earned 100,000 credits (100,015)GFN Silver: Earned 100,000 credits (246,369)PSA Silver: Earned 100,000 credits (226,594)
Message 77217 - Posted: 14 Jun 2014 | 2:26:15 UTC - in response to Message 77210.

Please have a look in the thread Minimum CUDA requirement for GFN. Iain explained the "z-transform".
____________
Best wishes. Knowledge is power. by jjwhalen

Yves GallotProject donor
Volunteer developer
Project scientist
Send message
Joined: 19 Aug 12
Posts: 843
ID: 164101
Credit: 306,521,622
RAC: 5,385
GFN Double Silver: Earned 200,000,000 credits (306,521,622)
Message 77229 - Posted: 14 Jun 2014 | 10:23:47 UTC - in response to Message 77210.

I called this new algorithm the "z-Transform" because of Georg Bruun paper:
z-Transform DFT Filters and FFT's, IEEE Trans. ASSP, ASSP-26, No. 1, February 1978.
http://orbit.dtu.dk/fedora/objects/orbit:20643/datastreams/file_4658740/content

If you have a look at Figure 4, you see that we can use half a FFT to "filter" with the z-transfer function z^-N/2 - 1. Because (z^-N - 1) / (z^-N/2 - 1) = z^-N/2 + 1, we compute modulo a GF polynomial.

That is faster than the DWT that was used in the previous versions of genefer:
http://www.ams.org/journals/mcom/1994-62-205/S0025-5718-1994-1185244-1/S0025-5718-1994-1185244-1.pdf
http://www.ams.org/journals/mcom/2002-71-238/S0025-5718-01-01350-3/S0025-5718-01-01350-3.pdf

This is not really "new math" but it is a new algorithm.

Note that the chapter VI of Bruun's paper is known as Bruun's FFT
http://en.wikipedia.org/wiki/Bruun's_FFT_algorithm
It is twice as fast as the complex z-Transform (for real data) but rouding error increases quickly with FFT size. Because of this inaccurate numerical precision, it cannot be used with 64-bit floating point type for GF test.

Yves

Yves GallotProject donor
Volunteer developer
Project scientist
Send message
Joined: 19 Aug 12
Posts: 843
ID: 164101
Credit: 306,521,622
RAC: 5,385
GFN Double Silver: Earned 200,000,000 credits (306,521,622)
Message 77230 - Posted: 14 Jun 2014 | 11:19:22 UTC - in response to Message 77210.

Will it have an effect on the search for GFPs of higher N/B values? How much faster is it?

The rouding error is smaller then higher B can be tested.

How mush faster? It depends on hardware and on the number of instances:
10%-30% faster on GPU, 10%-60% faster on CPU.
The algorithm was completely rewritten on CPU with a different data ordering.
It was tuned on Ivy Bridge and Haswell processors and then is very fast on them.
On GPU, the "structure" is identical but the number of operations is smaller.

On Mike's Core i5-4670K, he measured:
N=20
3.2.0, AVX, 1 core: 8.77 ms
3.2.0, AVX, 4 cores: 14.6 ms.
3.2.1, FMA, 1 core: 5.27 ms
3.2.1, FMA, 4 cores: 6.46 ms

Yves

Message boards : Generalized Fermat Prime Search : Z-transform?

[Return to PrimeGrid main page]
DNS Powered by DNSEXIT.COM
Copyright © 2005 - 2023 Rytis Slatkevičius (contact) and PrimeGrid community. Server load 2.39, 2.05, 2.23
Generated 23 Sep 2023 | 14:35:39 UTC