http://people.csail.mit.edu/jaffer/SLIB  
The SLIB Portable Scheme Library  

SLIB supports Bigloo, Chez, ELK 3.0, Gambit 4.0, Gauche, Guile, JScheme, Kawa, Larceny, MacScheme, MIT/GNU Scheme, Pocket Scheme, RScheme, scheme>C, Scheme48, SCM, SCM Mac, scsh, sisc, Stk, T3.1, umbscheme, and VSCM.Documentation includes a manifest, installation instructions, and coding standards for the library. Each library package is documented.
SLIB, Guile, Kawa, MIT/GNU Scheme, and SCM are GNU packages.
SLIB's new multidimensional spacefilling package has algorithms which unify the encoding and decoding of Hilbert and Peano curves of any rank, their inverses, and novel variations. Read the paper at http://arxiv.org/abs/1402.1807
Otherwise, slib3b5 is a minor release. Details at http://cvs.savannah.gnu.org/viewvc/*checkout*/slib/slib/ChangeLog
Volunteer Opportunities
 Add a Scheme procedure to eliminate specified variables from a list of equations to the commutativering (symbolic algebra) package.
 Augment dft.scm to work efficiently on nonpoweroftwo sized arrays.
Rupinder Singh (for the GSoC) proposed a plan of work:
We start with a basic implementation of the CooleyTukey Algorithm. CooleyTukey is probably the most popular and efficient (subject to limitations) algorithm for calculating DFT. Quoting from it's Wikipedia article: "The Cooleyâ€“Tukey algorithm, named after J.W. Cooley and John Tukey, is the most common fast Fourier transform (FFT) algorithm". It has been successful in achiveing a complexity of O(N logN) for composite numbers. CooleyTukey functions by expressing, recursively, the DFT of size N as N = N1 x N2 (where N1, N2 ≥ 1).
There are a couple of improvements / pruning that can be made to this algorithm. For example, if we incorporate Chinese Remainder Theorem to choose N1 and N2, we can avoid the 'Twiddle Factors' arising as Phase terms. Also, if at any stage of the recursion N comes out to be a larger prime, the CooleyTukey algorithm deviates from its O(N logN) performance and moves towards a O(N^{2}) performance. Here comes the Rader's Algorithm.
Rader's Algorithm computes the DFT of prime sizes through a cyclic convolution. The performance of this Algorithm is O(N logN) (after zero padding N to a power of 2). If a DHT (Discrete Hartley Transform) is incorporated, using modified versions of the reindexing routine, there can be a gain factor of two (Johnson and Frigo, 2007).
Hence, the algorithm can be reasonably expected to achieve a O(1/2 N logN) performance, faster than original CooleyTukey Algorithm (O(N^{2}) in prime case), Bluestein's Algorithm, Winograd's Algorithm and also better than a more recent Montium TileProcessor (exp(N)).
Therefore, this calls for a restructuring of the existing DFT mechanism. The major changes would be introduction of a new module to implement Chinese Remainder Theorem, another for Rader's Algorithm, and then it's cyclic convolution, Optimizing it by DHT, and Finally, the encapsulating CooleyTukey algorithm.
I am a guest and not a member of the MIT Computer Science and Artificial Intelligence Laboratory.
My actions and comments do not reflect in any way on MIT.  
agj @ alum.mit.edu  Go Figure! 