From f4e533c64e0c005e567b3fa76bc3456facbe5484 Mon Sep 17 00:00:00 2001 From: Jukka Ojanen Date: Wed, 16 Sep 2015 17:52:55 +0300 Subject: Add double-double arithmetic to generate "exact" double precision cosine and sine tables. Correct rounding verified using MPFR upto 2^28. SSE2 optimized ffts_generate_cosine_sine_pow2_64f takes twice as long as ffts_generate_cosine_sine_pow2_32f. --- src/ffts_trig.h | 3 +++ 1 file changed, 3 insertions(+) (limited to 'src/ffts_trig.h') diff --git a/src/ffts_trig.h b/src/ffts_trig.h index cfed2fb..0b22738 100644 --- a/src/ffts_trig.h +++ b/src/ffts_trig.h @@ -46,6 +46,9 @@ int ffts_generate_cosine_sine_pow2_32f(ffts_cpx_32f *const table, int table_size); int +ffts_generate_cosine_sine_pow2_64f(ffts_cpx_64f *const table, int table_size); + +int ffts_generate_table_1d_real_32f(struct _ffts_plan_t *const p, int sign, int invert); -- cgit v1.1