[eigen] Circulant matrix and FFT

["Philipp Wissmann" <phwissmann@xxxxxxxxx>]

Hi

 

I was thinking of adding a circulant matrix type and a way to use that type to solve linear systems by means of FFT.

 

There are two issues I came across

 

1) There is/was some discussions about special matrix types and how they should be implemented.

http://eigen.tuxfamily.org/index.php?title=SpecialMatrix#Toeplitz_matrix

 

Can anybody comment on this? If I do this, should I implement it in a more abstract way to allow for different "shapes"?

 

2) The FFT module is listed as unsupported; this would make any method for solving linear systems based on FFT unsupported as well.

Is this

http://eigen.tuxfamily.org/index.php?title=EigenFFT

and this

http://eigen.tuxfamily.org/index.php?title=Status_of_unsupported_modules

up to date considering the status of the FFT module?

What would be needed to move at least the 1d FFT kissfft implementation to 'supported'?

 

Best,

Philipp