mersenneforum.org  

Go Back   mersenneforum.org > Search Forums

Showing results 1 to 8 of 8
Search took 0.00 seconds.
Search: Posts Made By: ewmayer
Forum: Programming 2020-01-19, 19:47
Replies: 43
Views: 25,062
Posted By ewmayer
Found small bug in code I posted in #29 - when I...

Found small bug in code I posted in #29 - when I added the second computation (now both cos and sin approximations get computed) I forgot to reinit maxerr = 0 prior to part 2, which in the code is...
Forum: Programming 2020-01-19, 02:49
Replies: 43
Views: 25,062
Posted By ewmayer
I forgot to add, C++ defines...

I forgot to add, C++ defines reinterpret_cast<pointer> (www.cplusplus.com/doc/tutorial/typecasting/) as part of the language standard, thus in a sense enshrining the necessary-ugliness of...
Forum: Programming 2020-01-18, 23:26
Replies: 43
Views: 25,062
Posted By ewmayer
@Paul: I consider type-punning aliasing a...

@Paul: I consider type-punning aliasing a necessary evil, e.g. my quad-float code relies on it. There's also the "when you need to be 100% sure that double-const you are using is translated...
Forum: Programming 2020-01-17, 19:24
Replies: 43
Views: 25,062
Posted By ewmayer
Bzipped code attached. Note in the meanwhile I...

Bzipped code attached. Note in the meanwhile I got a PM form DrSardonicus about mapping a target interval [0,pi/2] to [-pi/4,+pi/4]. The code George posted in #14 of the thread was for sin and cos,...
Forum: Programming 2020-01-16, 21:42
Replies: 43
Views: 25,062
Posted By ewmayer
I've been playing with the use of a...

I've been playing with the use of a Chebyshev-polynomial basis expansion for approximating e.g. sin and cos - this gives Pade-similar accuracy but is much more computationally convenient. Longer...
Forum: Programming 2020-01-12, 19:36
Replies: 43
Views: 25,062
Posted By ewmayer
One could, say, loop over a bunch of...

One could, say, loop over a bunch of evenly-spaced points in the interval of interest, use the various formulae to compute sincos in DP, compute errors using expensive extended-precision sincos as a...
Forum: Programming 2020-01-11, 19:36
Replies: 43
Views: 25,062
Posted By ewmayer
That is more suitable for non-FMA architectures...

That is more suitable for non-FMA architectures ... for FMA-supporting ones, George's code needs 8 FMA (2 of which are pure-MUL), compared to 10 separate operations (mix of ADD, MUL, and FMA) for the...
Forum: Programming 2020-01-10, 22:44
Replies: 43
Views: 25,062
Posted By ewmayer
George, presumably you need this functionality...

George, presumably you need this functionality for on-the-fly computation of FFT twiddles, i.e. for complex roots of unity of form exp(I*2*Pi/n), where n is some highly composite integer, likely of...
Showing results 1 to 8 of 8

 
All times are UTC. The time now is 23:46.


Fri Jan 27 23:46:23 UTC 2023 up 162 days, 21:14, 0 users, load averages: 1.14, 1.07, 1.07

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.

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