mersenneforum.org  

Go Back   mersenneforum.org > Search Forums

Showing results 1 to 11 of 11
Search took 0.01 seconds.
Search: Posts Made By: Prime95
Forum: Programming 2020-07-15, 20:14
Replies: 43
Views: 25,120
Posted By Prime95
No division is done. An example for the case...

No division is done. An example for the case where n (the FFT length) is known at compile-time to be 5*2^some_power.

double ksinpi(double k, const double n) {
const double multiplier = 935.0,...
Forum: Programming 2020-07-15, 15:00
Replies: 43
Views: 25,120
Posted By Prime95
. After lots of study, my final solution...

.

After lots of study, my final solution was based on your post. It differs in that I need sincos (2*pi*k/n) where n is known in advance and is highly composite. The advance knowledge of n lets...
Forum: Programming 2020-01-17, 02:06
Replies: 43
Views: 25,120
Posted By Prime95
Good work, Ernst. When you've polished up the...

Good work, Ernst. When you've polished up the code, could you upload it for all to enjoy?
Forum: Programming 2020-01-12, 20:50
Replies: 43
Views: 25,120
Posted By Prime95
Looks like that code does a sqrt to get the sin...

Looks like that code does a sqrt to get the sin value from the cos value. On my GPU that is more costly than the polynomial calculation.
Forum: Programming 2020-01-12, 05:47
Replies: 43
Views: 25,120
Posted By Prime95
I did not try your code. If I counted right, the...

I did not try your code. If I counted right, the op count is unchanged. A mul became an add.

I've not yet run across code that produces both sin and cos at a reduced cost.
Forum: Programming 2020-01-11, 20:20
Replies: 43
Views: 25,120
Posted By Prime95
One would also need to see an analysis of the...

One would also need to see an analysis of the absolute and relative error in this new code.
Forum: Programming 2020-01-10, 23:25
Replies: 43
Views: 25,120
Posted By Prime95
Yes, it is for FFT twiddles on a GPU. Take...

Yes, it is for FFT twiddles on a GPU.

Take the awesome Radeon VII. Using the memory bandwidth from the spec sheets, reading a sin/cos twiddle from memory is 53 clocks latency if memory is...
Forum: Programming 2020-01-10, 04:15
Replies: 43
Views: 25,120
Posted By Prime95
Here is my code. Any improvements are welcome! ...

Here is my code. Any improvements are welcome!

// This version of slowTrig assumes k is positive and k/n <= 0.5 which means we want cos and sin values in the range [0, pi/2]
// We found free...
Forum: Programming 2020-01-09, 22:10
Replies: 43
Views: 25,120
Posted By Prime95
From the stackexchage post, these 2 links are...

From the stackexchage post, these 2 links are just what the doctor ordered:





It is not too hard to convert my [0,pi/2] range to [-pi/4,pi/4]
Forum: Programming 2020-01-09, 05:28
Replies: 43
Views: 25,120
Posted By Prime95
Just found this stackoverflow page which has some...

Just found this stackoverflow page which has some promising leads for me to investigate:

https://stackoverflow.com/questions/2284860/how-does-c-compute-sin-and-other-math-functions
Forum: Programming 2020-01-09, 04:58
Replies: 43
Views: 25,120
Posted By Prime95
Algorithm wanted - sine and/or cosine

I cannot locate my Knuth book. I need the fastest possible algorithm to compute sine(pi * a/b) and/or cosine(pi * a/b) on a GPU where a and b are positive 32-bit integers and a is less than half of...
Showing results 1 to 11 of 11

 
All times are UTC. The time now is 21:47.


Sun Feb 5 21:47:10 UTC 2023 up 171 days, 19:15, 1 user, load averages: 0.72, 0.94, 1.02

Powered by vBulletin® Version 3.8.11
Copyright ©2000 - 2023, Jelsoft Enterprises Ltd.

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.

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