roots_helper< CC, UU, SS > Struct Template Reference

#include <fft_naive.hpp>

List of all members.

Classes

struct fft_mul_sc_op

Public Types

typedef CC C
typedef UU U
typedef SS S

Static Public Member Functions

static U * create_roots (nat n, const format< U > &fm)
static void destroy_roots (U *u, nat n)
static void fft_cross (C *c1, C *c2)
static void dfft_cross (C *c1, C *c2, const U *u)
static void ifft_cross (C *c1, C *c2, const U *u)
static void dtft_cross (C *c1, C *c2)
static void dtft_cross (C *c1, C *c2, const U *u)
static void itft_flip (C *c1, C *c2, const U *u)
static void itft_flip (C *c1, C *c2)

Detailed Description

template<typename CC, typename UU, typename SS>
struct mmx::roots_helper< CC, UU, SS >

Definition at line 24 of file fft_naive.hpp.

Member Typedef Documentation

typedef CC C

Definition at line 25 of file fft_naive.hpp.

typedef SS S

Definition at line 27 of file fft_naive.hpp.

typedef UU U

Definition at line 26 of file fft_naive.hpp.

Member Function Documentation

static U* create_roots	(	nat	n,
		const format< U > &	fm
	)			`[inline, static]`

Definition at line 30 of file fft_naive.hpp.

References mmx::bit_mirror().

00030                                             {
00031     nat k= primitive_root_max_order<U> (2); (void) k;
00032     VERIFY (k == 0 || n <= k, "maximum order exceeded");
00033     VERIFY (n >= 2, "size must be at least two");
00034     U* roots= mmx_new<U> (n);
00035     for (nat i=0; i<n; i+=2) {
00036       roots[i]  = primitive_root<U> (n, bit_mirror (i, n), fm);
00037       roots[i+1]= primitive_root<U> (n, i==0? 0: n - bit_mirror (i, n), fm);
00038     }
00039     return roots; }

static void destroy_roots	(	U *	u,
		nat	n
	)			`[inline, static]`

Definition at line 42 of file fft_naive.hpp.

00042                               {
00043     mmx_delete<U> (u, n); }

static void dfft_cross	(	C *	c1,
		C *	c2,
		const U *	u
	)			`[inline, static]`

Definition at line 52 of file fft_naive.hpp.

Referenced by roots_helper< CC, UU, SS >::dtft_cross().

00052                                         {
00053     C temp= (*u) * (*c2);
00054     *c2   = (*c1) - temp;
00055     *c1   = (*c1) + temp; }

static void dtft_cross	(	C *	c1,
		C *	c2,
		const U *	u
	)			`[inline, static]`

Definition at line 71 of file fft_naive.hpp.

References roots_helper< CC, UU, SS >::dfft_cross(), and mmx::invert().

00071                                         {
00072     static S h= invert (S (2));
00073     dfft_cross (c1, c2, u);
00074     *c1 *= h;
00075     *c2 *= h; }

static void dtft_cross	(	C *	c1,
		C *	c2
	)			`[inline, static]`

Definition at line 64 of file fft_naive.hpp.

References roots_helper< CC, UU, SS >::fft_cross(), and mmx::invert().

00064                             {
00065     static S h= invert (S (2));
00066     fft_cross (c1, c2);
00067     *c1 *= h;
00068     *c2 *= h; }

static void fft_cross	(	C *	c1,
		C *	c2
	)			`[inline, static]`

Definition at line 46 of file fft_naive.hpp.

Referenced by roots_helper< CC, UU, SS >::dtft_cross().

00046                            {
00047     C temp= (*c2);
00048     *c2   = (*c1) - temp;
00049     *c1   = (*c1) + temp; }

static void ifft_cross	(	C *	c1,
		C *	c2,
		const U *	u
	)			`[inline, static]`

Definition at line 58 of file fft_naive.hpp.

00058                                         {
00059     C temp= *c2;
00060     *c2   = (*u ) * ((*c1) - temp);
00061     *c1   = (*c1) + temp; }

static void itft_flip	(	C *	c1,
		C *	c2
	)			`[inline, static]`

Definition at line 85 of file fft_naive.hpp.

References mmx::invert().

00085                            {
00086     static S h= invert (S(2));
00087     *c1  += (*c1) - (*c2);
00088     *c2   = h * ((*c1) - (*c2)); }

static void itft_flip	(	C *	c1,
		C *	c2,
		const U *	u
	)			`[inline, static]`

Definition at line 78 of file fft_naive.hpp.