In preparation for experiments that start tomorrow, I added today 
real-array fft and power spectrum capabilities to my dataframe code.  It 
uses GSLL's fft functions (search for file fft.lisp).  The classes, generic 
functions, and methods are still a bit rough, and will be refined.

You can read more about it: 
https://github.com/mirkov/data-table/blob/master/user/example1/README.org

That same directory contains examples-w-fft.lisp that shows the code 
usage.  I annotated the code.

The png files show sunspot data and two power spectra.

The example shows how I create one data-table with raw data, and then 
another one with slightly massaged data.  I do the fft on the latter one.

The example also shows that the current syntax is `functional' vs 
`declarative'. I explicitly  build the second data table, instead of 
declaring its dependency on the first table.  I think that the latter 
approach is preferable from a user's point of view.  But I will refrain 
from trying to design and implement this declarative approach.  I need more 
real-world experience before trying to do this.

I may add windowing capability to the fft code later on.

Mirko
_______________________________________________
GSLL-devel mailing list
GSLL-devel@common-lisp.net
http://lists.common-lisp.net/cgi-bin/mailman/listinfo/gsll-devel

Reply via email to