дискретное преобразование фурье

[alexdadaev]

Никакк не могу решить проблему с дискретным преобразованием фурье. Всё бы ничего, только класс у меня - шаблон и поэтому всё время сталкиваюсь с проблемой несовпадения типов. Даже банальную функцию pow не реализовать, так как она возвращает double, а у меня явно то тип не указывается.
В общем подскажите пожалуйста

template<class type> class complex {

	type re, im;
	
	public:
	complex();
	complex(type _re, type _im);
	~complex();
	
	type get_re() { return this->re; };
	type get_im() { return this->im; };
	type set_re(type re_val) { re = re_val; }
	type set_im(type im_val) { im = im_val; }

	complex operator+(complex& a);
	complex operator-(complex& a);
	complex operator*(complex& a);
	complex operator/(complex& a);
	
	void show();
	complex fft(complex<type> *x,type N);
	
};

	template<class type> void complex<type>::show()
	{ std::cout << re << " + " << im << "i\n"; }
	
	template<class type> complex<type>::complex():
	re(0),im(0)
	{ }
	
	template<class type> complex<type>::complex(type _re, type _im):
	re(_re),im(_im)
	{ }
	
	template<class type> complex<type>::~complex()
	{ }

	template<class type> complex<type> complex<type>::operator+(complex<type> &a)
	{
		type re_, im_;
		re_ = re + a.re;
		im_ = im + a.im;
		return complex(re_,im_);
	}
	
	template<class type> complex<type> complex<type>::operator-(complex<type> &a)
	{
		type re_, im_;
		re = re - a.re;
		im = im - a.im;
		return complex(re_,im_);
	}
	
	template<class type> complex<type> complex<type>::operator/(complex<type> &a)
	{
		type re_,im_;
		re_ = (re * a.re + im * a.im) / (a.re * a.re + im * im);  
	    im_ = (im * a.re - re * a.im) / (a.re * a.re + im * im);  
   
		return complex(re_,im_); 
	}
	
	template<class type> complex<type> complex<type>::operator*(complex<type> &a)
	{
		type re_,im_;
		re = re*a.re;
		im = im*a.im;
		return complex(re_,im_);
	}

[like-nix]

Не удивительно вы в каждом операторе определяете локальные переменные для оператора класса type re_,im_; и далее передаете их в return complex(re_,im_); а с неинициализированные переменными базового типа нельзя выполнять арефметические операции. уберите в операторах строку type re_,im_; и замените return complex(re_,im_); на return complex(re,im); и оператор умножения комплексных чисел у вас не правильный

[Obey-Kun]

Вместо type используйте T, это красивее и общепринято. Вопрос не читал, извините.