/*----------------------------------------------------------------------------  ChucK Concurrent, On-the-fly Audio Programming Language    Compiler and Virtual Machine   Copyright (c) 2004 Ge Wang and Perry R. Cook.  All rights reserved.    http://chuck.stanford.edu/    http://chuck.cs.princeton.edu/   This program is free software; you can redistribute it and/or modify  it under the terms of the GNU General Public License as published by  the Free Software Foundation; either version 2 of the License, or  (at your option) any later version.   This program is distributed in the hope that it will be useful,  but WITHOUT ANY WARRANTY; without even the implied warranty of  MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the  GNU General Public License for more details.   You should have received a copy of the GNU General Public License  along with this program; if not, write to the Free Software  Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA 02111-1307  U.S.A.-----------------------------------------------------------------------------*/  //-----------------------------------------------------------------------------// name: util_xform.h// desc: transform utilities//// authors: Ge Wang (ge@ccrma.stanford.edu | gewang@cs.princeton.edu)//          Perry R. Cook (prc@cs.princeton.edu)//          FFT - from CARL distribution// date: 11.27.2003//-----------------------------------------------------------------------------#ifndef __UTIL_XFORM_H__#define __UTIL_XFORM_H__ #include <math.h>#include "chuck_def.h" // our float type#define FLOAT       SAMPLE // complex absolute value#define cmp_abs(x) ( sqrt( (x).re * (x).re + (x).im * (x).im ) )#define __modulus(x) cmp_abs(x)#define __phase(x) ( atan2((double)(x).im, (double)(x).re) ) #define FFT_FORWARD 1#define FFT_INVERSE 0#define make_window hanning // c linkage#if defined(_cplusplus) || defined(__cplusplus)extern "C" {#endif // make the windowvoid hanning( FLOAT * window, unsigned long length );void hamming( FLOAT * window, unsigned long length );void blackman( FLOAT * window, unsigned long length );void bartlett( FLOAT * window, unsigned long length );// apply the windowvoid apply_window( FLOAT * data, FLOAT * window, unsigned long length ); // real fft, N must be power of 2void rfft( FLOAT * x, long N, unsigned int forward );// complex fft, NC must be power of 2void cfft( FLOAT * x, long NC, unsigned int forward ); // type II dct, often referred to as "the dct"void the_dct( FLOAT * x, unsigned long N, FLOAT * out, unsigned long Nout );// generates NxN type II dct matrixvoid the_dct_matrix( FLOAT ** out, unsigned long N );// type III dct, often referred to as "the inverse_dct"void the_inverse_dct( FLOAT * x, unsigned long N, FLOAT * out, unsigned long Nout );// generates NxN type II dct matrixvoid the_inverse_dct_matrix( FLOAT ** out, unsigned long N ); // apply dct from matrixvoid the_dct_now( FLOAT * x, FLOAT ** matrix, unsigned long N, FLOAT * out, unsigned long Nout );// apply inverse dct from matrixvoid the_inverse_dct_now( FLOAT * x, FLOAT ** matrix, unsigned long N, FLOAT * out, unsigned long Nout ); // c linkage#if ( defined( __cplusplus ) || defined( _cplusplus ) )}#endif #endif