Hi
I have made small change to circulant eig on the way it sets the
warning levels. It should now restore the warnings back to the
previous state when finished and does it in a unwind_protect block
(just in case).
Carnë
I added some circulant matrix functions to linear-algebra:
circulant_make_matrix Produce a full circulant matrix given the first column
circulant_matrix_vector_product Fast, compact calculation of the
product of a circulant matrix with a vector
circulant_eig Fast, compact calculation of eigenva
Yes, that's one reason why it's advantageous to treat them as special cases
On Wed, Mar 28, 2012 at 5:21 PM, Stephen Montgomery-Smith
wrote:
> Could all of these be implemented using fft and ifft? For example, to
> get the eigenvalues, just fft one of the rows of the matrix.
---
On 03/28/2012 03:06 PM, Nir Krakauer wrote:
> Are there any functions in Octave for working with circulant matrices
> ( finding inverses, eigenvalues, and the like --
> http://en.wikipedia.org/wiki/Circulant_matrix )? I've written a few --
> would a new category in the linear-algebra package be a g
Are there any functions in Octave for working with circulant matrices
( finding inverses, eigenvalues, and the like --
http://en.wikipedia.org/wiki/Circulant_matrix )? I've written a few --
would a new category in the linear-algebra package be a good place to
put them?
--Nir
-
