Update of /cvsroot/pure-data/externals/iem/iem_bin_ambi/src
In directory sc8-pr-cvs1.sourceforge.net:/tmp/cvs-serv443/src
Added Files:
bin_ambi_calc_HRTF.c bin_ambi_reduced_decode.c iem_bin_ambi.c
iem_bin_ambi.dsp iem_bin_ambi.dsw iem_bin_ambi.h iemlib.h
Log Message:
added the binaural ambisonic decoder
--- NEW FILE: iem_bin_ambi.c ---
/* For information on usage and redistribution, and for a DISCLAIMER OF ALL
* WARRANTIES, see the file, "LICENSE.txt," in this distribution.
iem_bin_ambi written by Thomas Musil, Copyright (c) IEM KUG Graz Austria 2000 - 2003 */
#ifdef NT
#pragma warning( disable : 4244 )
#pragma warning( disable : 4305 )
#endif
#include "m_pd.h"
#include "iemlib.h"
static t_class *iem_bin_ambi_class;
static void *iem_bin_ambi_new(void)
{
t_object *x = (t_object *)pd_new(iem_bin_ambi_class);
return (x);
}
void bin_ambi_calc_HRTF_setup(void);
void bin_ambi_reduced_decode_setup(void);
/* ------------------------ setup routine ------------------------- */
void iem_bin_ambi_setup(void)
{
bin_ambi_calc_HRTF_setup();
bin_ambi_reduced_decode_setup();
post("iem_bin_ambi (R-1.15) library loaded!");
}
--- NEW FILE: iem_bin_ambi.h ---
/* For information on usage and redistribution, and for a DISCLAIMER OF ALL
* WARRANTIES, see the file, "LICENSE.txt," in this distribution.
iem_bin_ambi written by Thomas Musil, Copyright (c) IEM KUG Graz Austria 2000 - 2003 */
#ifndef __IEMBINAMBI_H__
#define __IEMBINAMBI_H__
typedef struct
{
float real;
float imag;
}
BIN_AMBI_COMPLEX;
#endif
--- NEW FILE: iem_bin_ambi.dsw ---
(This appears to be a binary file; contents omitted.)
--- NEW FILE: iemlib.h ---
/* For information on usage and redistribution, and for a DISCLAIMER OF ALL
* WARRANTIES, see the file, "LICENSE.txt," in this distribution.
iem_bin_ambi written by Thomas Musil, Copyright (c) IEM KUG Graz Austria 2000 - 2003 */
#ifndef __IEMLIB_H__
#define __IEMLIB_H__
#define IS_A_POINTER(atom,index) ((atom+index)->a_type == A_POINTER)
#define IS_A_FLOAT(atom,index) ((atom+index)->a_type == A_FLOAT)
#define IS_A_SYMBOL(atom,index) ((atom+index)->a_type == A_SYMBOL)
#define IS_A_DOLLAR(atom,index) ((atom+index)->a_type == A_DOLLAR)
#define IS_A_DOLLSYM(atom,index) ((atom+index)->a_type == A_DOLLSYM)
#define IS_A_SEMI(atom,index) ((atom+index)->a_type == A_SEMI)
#define IS_A_COMMA(atom,index) ((atom+index)->a_type == A_COMMA)
#ifdef NT
int sys_noloadbang;
//t_symbol *iemgui_key_sym=0;
#include <io.h>
#else
extern int sys_noloadbang;
//extern t_symbol *iemgui_key_sym;
#include <unistd.h>
#endif
#define DEFDELVS 64
#define XTRASAMPS 4
#define SAMPBLK 4
#define UNITBIT32 1572864. /* 3*2^19; bit 32 has place value 1 */
/* machine-dependent definitions. These ifdefs really
should have been by CPU type and not by operating system! */
#ifdef IRIX
/* big-endian. Most significant byte is at low address in memory */
#define HIOFFSET 0 /* word offset to find MSB */
#define LOWOFFSET 1 /* word offset to find LSB */
#define int32 long /* a data type that has 32 bits */
#else
#ifdef NT
/* little-endian; most significant byte is at highest address */
#define HIOFFSET 1
#define LOWOFFSET 0
#define int32 long
#else
#ifdef __linux__
#include <endian.h>
#if !defined(__BYTE_ORDER) || !defined(__LITTLE_ENDIAN)
#error No byte order defined
#endif
#if __BYTE_ORDER == __LITTLE_ENDIAN
#define HIOFFSET 1
#define LOWOFFSET 0
#else
#define HIOFFSET 0 /* word offset to find MSB */
#define LOWOFFSET 1 /* word offset to find LSB */
#endif /* __BYTE_ORDER */
#include <sys/types.h>
#define int32 int32_t
#endif /* __linux__ */
#endif /* NT */
#endif /* SGI */
union tabfudge
{
double tf_d;
int32 tf_i[2];
};
#endif
--- NEW FILE: iem_bin_ambi.dsp ---
# Microsoft Developer Studio Project File - Name="iem_bin_ambi" - Package Owner=<4>
# Microsoft Developer Studio Generated Build File, Format Version 6.00
# ** NICHT BEARBEITEN **
# TARGTYPE "Win32 (x86) External Target" 0x0106
CFG=iem_bin_ambi - Win32 Debug
!MESSAGE Dies ist kein gültiges Makefile. Zum Erstellen dieses Projekts mit NMAKE
!MESSAGE verwenden Sie den Befehl "Makefile exportieren" und führen Sie den Befehl
!MESSAGE
!MESSAGE NMAKE /f "iem_bin_ambi.mak".
!MESSAGE
!MESSAGE Sie können beim Ausführen von NMAKE eine Konfiguration angeben
!MESSAGE durch Definieren des Makros CFG in der Befehlszeile. Zum Beispiel:
!MESSAGE
!MESSAGE NMAKE /f "iem_bin_ambi.mak" CFG="iem_bin_ambi - Win32 Debug"
!MESSAGE
!MESSAGE Für die Konfiguration stehen zur Auswahl:
!MESSAGE
!MESSAGE "iem_bin_ambi - Win32 Release" (basierend auf "Win32 (x86) External Target")
!MESSAGE "iem_bin_ambi - Win32 Debug" (basierend auf "Win32 (x86) External Target")
!MESSAGE
# Begin Project
# PROP AllowPerConfigDependencies 0
# PROP Scc_ProjName ""
# PROP Scc_LocalPath ""
!IF "$(CFG)" == "iem_bin_ambi - Win32 Release"
# PROP BASE Use_Debug_Libraries 0
# PROP BASE Output_Dir "Release"
# PROP BASE Intermediate_Dir "Release"
# PROP BASE Cmd_Line "NMAKE /f makefile_win"
# PROP BASE Rebuild_Opt "/a"
# PROP BASE Target_File "makefile_win.exe"
# PROP BASE Bsc_Name "makefile_win.bsc"
# PROP BASE Target_Dir ""
# PROP Use_Debug_Libraries 0
# PROP Output_Dir "Release"
# PROP Intermediate_Dir "Release"
# PROP Cmd_Line "NMAKE /f makefile_win"
# PROP Rebuild_Opt "/a"
# PROP Target_File "iem_bin_ambi.exe"
# PROP Bsc_Name "iem_bin_ambi.bsc"
# PROP Target_Dir ""
!ELSEIF "$(CFG)" == "iem_bin_ambi - Win32 Debug"
# PROP BASE Use_Debug_Libraries 1
# PROP BASE Output_Dir "Debug"
# PROP BASE Intermediate_Dir "Debug"
# PROP BASE Cmd_Line "NMAKE /f makefile_win"
# PROP BASE Rebuild_Opt "/a"
# PROP BASE Target_File "makefile_win.exe"
# PROP BASE Bsc_Name "makefile_win.bsc"
# PROP BASE Target_Dir ""
# PROP Use_Debug_Libraries 1
# PROP Output_Dir "Debug"
# PROP Intermediate_Dir "Debug"
# PROP Cmd_Line "NMAKE /f makefile_win"
# PROP Rebuild_Opt "/a"
# PROP Target_File "iem_bin_ambi.exe"
# PROP Bsc_Name "iem_bin_ambi.bsc"
# PROP Target_Dir ""
!ENDIF
# Begin Target
# Name "iem_bin_ambi - Win32 Release"
# Name "iem_bin_ambi - Win32 Debug"
!IF "$(CFG)" == "iem_bin_ambi - Win32 Release"
!ELSEIF "$(CFG)" == "iem_bin_ambi - Win32 Debug"
!ENDIF
# Begin Source File
SOURCE=.\makefile_win
# End Source File
# End Target
# End Project
--- NEW FILE: bin_ambi_reduced_decode.c ---
/* For information on usage and redistribution, and for a DISCLAIMER OF ALL
* WARRANTIES, see the file, "LICENSE.txt," in this distribution.
iem_bin_ambi written by Thomas Musil, Copyright (c) IEM KUG Graz Austria 2000 - 2003 */
#ifdef NT
#pragma warning( disable : 4244 )
#pragma warning( disable : 4305 )
#endif
#include "m_pd.h"
#include "iemlib.h"
#include "iem_bin_ambi.h"
#include <math.h>
#include <stdio.h>
#include <string.h>
[...1209 lines suppressed...]
class_addmethod(bin_ambi_reduced_decode_class, (t_method)bin_ambi_reduced_decode_decoding, gensym("decoding"), 0);
class_addmethod(bin_ambi_reduced_decode_class, (t_method)bin_ambi_reduced_decode_ambi_weight, gensym("ambi_weight"), A_GIMME, 0);
class_addmethod(bin_ambi_reduced_decode_class, (t_method)bin_ambi_reduced_decode_sing_range, gensym("sing_range"), A_DEFFLOAT, 0);
class_sethelpsymbol(bin_ambi_reduced_decode_class, gensym("iemhelp2/help-bin_ambi_reduced_decode"));
}
/*
Reihenfolge:
n_ambi x bin_ambi_reduced_decode_ls
bin_ambi_reduced_decode_calc_pinv
n_ambi x bin_ambi_reduced_decode_load_HRIR
n_ambi x bin_ambi_reduced_decode_check_arrays
n_ambi x bin_ambi_reduced_decode_calc_reduced
bin_ambi_reduced_decode_decoding
*/
--- NEW FILE: bin_ambi_calc_HRTF.c ---
/* For information on usage and redistribution, and for a DISCLAIMER OF ALL
* WARRANTIES, see the file, "LICENSE.txt," in this distribution.
iem_bin_ambi written by Thomas Musil, Copyright (c) IEM KUG Graz Austria 2000 - 2003 */
#ifdef NT
#pragma warning( disable : 4244 )
#pragma warning( disable : 4305 )
#endif
#include "m_pd.h"
#include "iemlib.h"
#include "iem_bin_ambi.h"
#include <math.h>
#include <stdio.h>
#include <string.h>
/* -------------------------- bin_ambi_calc_HRTF ------------------------------ */
/*
** berechnet ein reduziertes Ambisonic-Decoder-Set in die HRTF-Spektren **
** Inputs: ls + Liste von 3 floats: Index [1 .. 16] + Elevation [-90 .. +90 degree] + Azimut [0 .. 360 degree] **
** Inputs: calc_inv **
** Inputs: load_HRIR + float index1..16 **
** Outputs: List of 2 symbols: left-HRIR-File-name + HRIR-table-name **
** Inputs: calc_reduced **
** "output" ... writes the HRTF into tables **
** **
** **
** setzt voraus , dass die HRIR-tabele-names von LS1_L_HRIR .. LS16_L_HRIR heissen und existieren **
** setzt voraus , dass die HRTF-tabele-names von LS1_HRTF_re .. LS16_HRTF_re heissen und existieren **
** setzt voraus , dass die HRTF-tabele-names von LS1_HRTF_im .. LS16_HRTF_im heissen und existieren **
*/
typedef struct _bin_ambi_calc_HRTF
{
t_object x_obj;
t_atom x_at[2];
int x_n_ls;
int x_fftsize;
int *x_delta;
int *x_phi;
BIN_AMBI_COMPLEX *x_spec;
BIN_AMBI_COMPLEX *x_sin_cos;
float *x_beg_fade_out_hrir;
float *x_beg_hrir;
float **x_beg_hrtf_re;
float **x_beg_hrtf_im;
t_symbol **x_hrir_filename;
t_symbol **x_s_hrir;
t_symbol **x_s_hrtf_re;
t_symbol **x_s_hrtf_im;
t_symbol *x_s_fade_out_hrir;
double x_pi_over_180;
} t_bin_ambi_calc_HRTF;
static t_class *bin_ambi_calc_HRTF_class;
static void bin_ambi_calc_HRTF_init_cos(t_bin_ambi_calc_HRTF *x)
{
int i, fftsize = x->x_fftsize;
float f, g;
BIN_AMBI_COMPLEX *sincos = x->x_sin_cos;
g = 2.0f * 3.1415926538f / (float)fftsize;
for(i=0; i<fftsize; i++)
{
f = g * (float)i;
(*sincos).real = cos(f);
(*sincos).imag = sin(f);/*change*/
sincos++;
}
}
static void bin_ambi_calc_HRTF_quant(t_bin_ambi_calc_HRTF *x, double *delta_deg2rad, double *phi_deg2rad, int index)
{
double q = 1.0;
double d = *delta_deg2rad;
double p = *phi_deg2rad;
int i;
if(d < -40.0)
d = -40.0;
if(d > 90.0)
d = 90.0;
while(p < 0.0)
p += 360.0;
while(p >= 360.0)
p -= 360.0;
if(d < -35.0)
{
*delta_deg2rad = -40.0;
q = 360.0 / 56.0;
}
else if(d < -25.0)
{
*delta_deg2rad = -30.0;
q = 6.0;
}
else if(d < -15.0)
{
*delta_deg2rad = -20.0;
q = 5.0;
}
else if(d < -5.0)
{
*delta_deg2rad = -10.0;
q = 5.0;
}
else if(d < 5.0)
{
*delta_deg2rad = 0.0;
q = 5.0;
}
else if(d < 15.0)
{
*delta_deg2rad = 10.0;
q = 5.0;
}
else if(d < 25.0)
{
*delta_deg2rad = 20.0;
q = 5.0;
}
else if(d < 35.0)
{
*delta_deg2rad = 30.0;
q = 6.0;
}
else if(d < 45.0)
{
*delta_deg2rad = 40.0;
q = 360.0 / 56.0;
}
else if(d < 55.0)
{
*delta_deg2rad = 50.0;
q = 8.0;
}
else if(d < 65.0)
{
*delta_deg2rad = 60.0;
q = 10.0;
}
else if(d < 75.0)
{
*delta_deg2rad = 70.0;
q = 15.0;
}
else if(d < 85.0)
{
*delta_deg2rad = 80.0;
q = 30.0;
}
else
{
*delta_deg2rad = 90.0;
q = 360.0;
}
p /= q;
i = (int)(p + 0.499999);
p = (double)i * q;
i = (int)(p + 0.499999);
while(i >= 360)
i -= 360;
p = (double)i;
*phi_deg2rad = p;
x->x_delta[index] = (int)(*delta_deg2rad);
x->x_phi[index] = i;
*delta_deg2rad *= x->x_pi_over_180;
*phi_deg2rad *= x->x_pi_over_180;
}
static void bin_ambi_calc_HRTF_do_2d(t_bin_ambi_calc_HRTF *x, int argc, t_atom *argv)
{
double delta=0.0, phi;
int index;
int n_ls = x->x_n_ls;
if(argc < 2)
{
post("bin_ambi_calc_HRTF ERROR: ls-input needs 1 index and 1 angle: ls_index + phi [degree]");
return;
}
index = (int)atom_getint(argv++) - 1;
phi = (double)atom_getfloat(argv);
if(index < 0)
index = 0;
if(index >= n_ls)
index = n_ls - 1;
bin_ambi_calc_HRTF_quant(x, &delta, &phi, index);
}
static void bin_ambi_calc_HRTF_do_3d(t_bin_ambi_calc_HRTF *x, int argc, t_atom *argv)
{
double delta, phi;
int index;
int n_ls=x->x_n_ls;
if(argc < 3)
{
post("bin_ambi_calc_HRTF ERROR: ls-input needs 1 index and 2 angles: ls index + delta [degree] + phi [degree]");
return;
}
index = (int)atom_getint(argv++) - 1;
delta = atom_getfloat(argv++);
phi = atom_getfloat(argv);
if(index < 0)
index = 0;
if(index >= n_ls)
index = n_ls - 1;
bin_ambi_calc_HRTF_quant(x, &delta, &phi, index);
}
static void bin_ambi_calc_HRTF_ls(t_bin_ambi_calc_HRTF *x, t_symbol *s, int argc, t_atom *argv)
{
if(argc >= 3)
bin_ambi_calc_HRTF_do_3d(x, argc, argv);
else
bin_ambi_calc_HRTF_do_2d(x, argc, argv);
}
static void bin_ambi_calc_HRTF_check_fade_out(t_bin_ambi_calc_HRTF *x)
{
t_garray *a;
int npoints;
t_float *fadevec;
if((int)(x->x_beg_fade_out_hrir) == 0)
{
if (!(a = (t_garray *)pd_findbyclass(x->x_s_fade_out_hrir, garray_class)))
error("%s: no such array", x->x_s_fade_out_hrir->s_name);
else if (!garray_getfloatarray(a, &npoints, &fadevec))
error("%s: bad template for bin_ambi_calc_HRTF", x->x_s_fade_out_hrir->s_name);
else if (npoints < x->x_fftsize)
error("%s: bad array-size: %d", x->x_s_fade_out_hrir->s_name, npoints);
else
x->x_beg_fade_out_hrir = fadevec;
}
}
/*
x_prod:
n_ambi columns;
n_ambi rows;
*/
static void bin_ambi_calc_HRTF_load_HRIR(t_bin_ambi_calc_HRTF *x, float findex)
{
int index=(int)findex - 1;
int p;
char buf[60];
if(index < 0)
index = 0;
if(index >= x->x_n_ls)
index = x->x_n_ls - 1;
p = x->x_phi[index];
if(p)/*change*/
p = 360 - p;
if(p < 10)
sprintf(buf, "L%de00%da.wav", x->x_delta[index], p);
else if(p < 100)
sprintf(buf, "L%de0%da.wav", x->x_delta[index], p);
else
sprintf(buf, "L%de%da.wav", x->x_delta[index], p);
x->x_hrir_filename[index] = gensym(buf);
SETSYMBOL(x->x_at, x->x_hrir_filename[index]);
SETSYMBOL(x->x_at+1, x->x_s_hrir[index]);
outlet_list(x->x_obj.ob_outlet, &s_list, 2, x->x_at);
}
static void bin_ambi_calc_HRTF_check_arrays(t_bin_ambi_calc_HRTF *x, float findex)
{
int index=(int)findex - 1;
int j, k, n;
int fftsize = x->x_fftsize;
int fs2=fftsize/2;
t_garray *a;
int npoints;
t_float *vec_hrir, *vec, *vec_fade_out_hrir;
t_float *vec_hrtf_re, *vec_hrtf_im;
t_symbol *hrir, *hrtf_re, *hrtf_im;
float decr, sum;
if(index < 0)
index = 0;
if(index >= x->x_n_ls)
index = x->x_n_ls - 1;
hrir = x->x_s_hrir[index];
hrtf_re = x->x_s_hrtf_re[index];
hrtf_im = x->x_s_hrtf_im[index];
if (!(a = (t_garray *)pd_findbyclass(hrtf_re, garray_class)))
error("%s: no such array", hrtf_re->s_name);
else if (!garray_getfloatarray(a, &npoints, &vec_hrtf_re))
error("%s: bad template for bin_ambi_calc_HRTF", hrtf_re->s_name);
else if (npoints < fftsize)
error("%s: bad array-size: %d", hrtf_re->s_name, npoints);
else if (!(a = (t_garray *)pd_findbyclass(hrtf_im, garray_class)))
error("%s: no such array", hrtf_im->s_name);
else if (!garray_getfloatarray(a, &npoints, &vec_hrtf_im))
error("%s: bad template for bin_ambi_calc_HRTF", hrtf_im->s_name);
else if (npoints < fftsize)
error("%s: bad array-size: %d", hrtf_im->s_name, npoints);
else if (!(a = (t_garray *)pd_findbyclass(hrir, garray_class)))
error("%s: no such array", hrir->s_name);
else if (!garray_getfloatarray(a, &npoints, &vec_hrir))
error("%s: bad template for bin_ambi_calc_HRTF", hrir->s_name);
else
{
x->x_beg_hrtf_re[index] = vec_hrtf_re;
x->x_beg_hrtf_im[index] = vec_hrtf_im;
if(npoints < fftsize)
{
post("warning: %s-array-size: %d", hrir->s_name, npoints);
}
vec = x->x_beg_hrir;
vec += index * fftsize;
if((int)(x->x_beg_fade_out_hrir))
{
vec_fade_out_hrir = x->x_beg_fade_out_hrir;
for(j=0; j<fs2; j++)
vec[j] = vec_hrir[j]*vec_fade_out_hrir[j];
}
else
{
post("no HRIR-fade-out-window found");
n = fs2 * 3;
n /= 4;
for(j=0; j<n; j++)
vec[j] = vec_hrir[j];
sum = 1.0f;
decr = 4.0f / (float)fs2;
for(j=n, k=0; j<fs2; j++, k++)
{
sum -= decr;
vec[j] = vec_hrir[j] * sum;
}
}
}
}
static void bin_ambi_calc_HRTF_calc_fft(t_bin_ambi_calc_HRTF *x, float findex)
{
int index=(int)findex - 1;
int i, j, k, w_index, w_inc, i_inc, v_index, fs1, fs2;
int fftsize = x->x_fftsize;
BIN_AMBI_COMPLEX old1, old2, w;
BIN_AMBI_COMPLEX *sincos = x->x_sin_cos;
BIN_AMBI_COMPLEX *val = x->x_spec;
t_float *vec_hrir, *vec_hrtf_re, *vec_hrtf_im;
int n_ls = x->x_n_ls;
if(index < 0)
index = 0;
if(index >= n_ls)
index = n_ls - 1;
vec_hrtf_re = x->x_beg_hrtf_re[index];
vec_hrtf_im = x->x_beg_hrtf_im[index];
vec_hrir = x->x_beg_hrir;
vec_hrir += index * fftsize;
for(k=0; k<fftsize; k++)
{
val[k].real = vec_hrir[k];
val[k].imag = 0.0f;
}
fs1 = fftsize - 1;
fs2 = fftsize / 2;
i_inc = fs2;
w_inc = 1;
for(i=1; i<fftsize; i<<=1)
{
v_index = 0;
for(j=0; j<i; j++)
{
w_index = 0;
for(k=0; k<i_inc; k++)
{
old1 = val[v_index];
old2 = val[v_index+i_inc];
w = sincos[w_index];
val[v_index+i_inc].real = (old1.real-old2.real)*w.real - (old1.imag-old2.imag)*w.imag;
val[v_index+i_inc].imag = (old1.imag - old2.imag)*w.real + (old1.real-old2.real)*w.imag;
val[v_index].real = old1.real+old2.real;
val[v_index].imag = old1.imag+old2.imag;
w_index += w_inc;
v_index++;
}
v_index += i_inc;
}
w_inc <<= 1;
i_inc >>= 1;
}
j = 0;
for(i=1;i<fs1;i++)
{
k = fs2;
while(k <= j)
{
j = j - k;
k >>= 1;
}
j = j + k;
if(i < j)
{
old1 = val[j];
val[j] = val[i];
val[i] = old1;
}
}
vec_hrtf_re[0] = val[0].real;
vec_hrtf_im[0] = 0.0f;
for( i = 1; i < fs2; i++ )
{
vec_hrtf_re[i] = 2.0f*val[i].real;
vec_hrtf_im[i] = 2.0f*val[i].imag;
}
vec_hrtf_re[fs2] = 0.0f;
vec_hrtf_im[fs2] = 0.0f;
}
static void bin_ambi_calc_HRTF_free(t_bin_ambi_calc_HRTF *x)
{
freebytes(x->x_hrir_filename, x->x_n_ls * sizeof(t_symbol *));
freebytes(x->x_s_hrir, x->x_n_ls * sizeof(t_symbol *));
freebytes(x->x_s_hrtf_re, x->x_n_ls * sizeof(t_symbol *));
freebytes(x->x_s_hrtf_im, x->x_n_ls * sizeof(t_symbol *));
freebytes(x->x_delta, x->x_n_ls * sizeof(int));
freebytes(x->x_phi, x->x_n_ls * sizeof(int));
freebytes(x->x_spec, x->x_fftsize * sizeof(BIN_AMBI_COMPLEX));
freebytes(x->x_sin_cos, x->x_fftsize * sizeof(BIN_AMBI_COMPLEX));
freebytes(x->x_beg_hrir, x->x_fftsize * x->x_n_ls * sizeof(float));
freebytes(x->x_beg_hrtf_re, x->x_n_ls * sizeof(float *));
freebytes(x->x_beg_hrtf_im, x->x_n_ls * sizeof(float *));
}
/*
1.arg: t_symbol *hrir_name;
2.arg: t_symbol *hrtf_re_name;
3.arg: t_symbol *hrtf_im_name;
4.arg: t_symbol *hrir_fade_out_name;
5.arg: int ambi_order;
6.arg: int dim;
7.arg: int n_ambi;
8.arg: int fftsize;
*/
static void *bin_ambi_calc_HRTF_new(t_symbol *s, int argc, t_atom *argv)
{
t_bin_ambi_calc_HRTF *x = (t_bin_ambi_calc_HRTF *)pd_new(bin_ambi_calc_HRTF_class);
char buf[400];
int i, j, fftok;
int n_ls, fftsize;
t_symbol *s_hrir;
t_symbol *s_hrtf_re;
t_symbol *s_hrtf_im;
if((argc >= 6) &&
IS_A_SYMBOL(argv,0) &&
IS_A_SYMBOL(argv,1) &&
IS_A_SYMBOL(argv,2) &&
IS_A_SYMBOL(argv,3) &&
IS_A_FLOAT(argv,4) &&
IS_A_FLOAT(argv,5))
{
s_hrir = (t_symbol *)atom_getsymbolarg(0, argc, argv);
s_hrtf_re = (t_symbol *)atom_getsymbolarg(1, argc, argv);
s_hrtf_im = (t_symbol *)atom_getsymbolarg(2, argc, argv);
x->x_s_fade_out_hrir = (t_symbol *)atom_getsymbolarg(3, argc, argv);
n_ls = (int)atom_getintarg(4, argc, argv);
fftsize = (int)atom_getintarg(5, argc, argv);
if(n_ls < 1)
n_ls = 1;
j = 2;
fftok = 0;
for(i=0; i<21; i++)
{
if(j == fftsize)
{
fftok = 1;
i = 22;
}
j *= 2;
}
if(!fftok)
{
post("bin_ambi_calc_HRTF-ERROR: fftsize not equal to 2 ^ n !!!");
return(0);
}
x->x_n_ls = n_ls;
x->x_fftsize = fftsize;
x->x_hrir_filename = (t_symbol **)getbytes(x->x_n_ls * sizeof(t_symbol *));
x->x_s_hrir = (t_symbol **)getbytes(x->x_n_ls * sizeof(t_symbol *));
x->x_s_hrtf_re = (t_symbol **)getbytes(x->x_n_ls * sizeof(t_symbol *));
x->x_s_hrtf_im = (t_symbol **)getbytes(x->x_n_ls * sizeof(t_symbol *));
for(i=0; i<n_ls; i++)
{
sprintf(buf, "%d%s", i+1, s_hrir->s_name);
x->x_s_hrir[i] = gensym(buf);
sprintf(buf, "%d%s", i+1, s_hrtf_re->s_name);
x->x_s_hrtf_re[i] = gensym(buf);
sprintf(buf, "%d%s", i+1, s_hrtf_im->s_name);
x->x_s_hrtf_im[i] = gensym(buf);
}
x->x_delta = (int *)getbytes(x->x_n_ls * sizeof(int));
x->x_phi = (int *)getbytes(x->x_n_ls * sizeof(int));
x->x_spec = (BIN_AMBI_COMPLEX *)getbytes(x->x_fftsize * sizeof(BIN_AMBI_COMPLEX));
x->x_sin_cos = (BIN_AMBI_COMPLEX *)getbytes(x->x_fftsize * sizeof(BIN_AMBI_COMPLEX));
x->x_beg_fade_out_hrir = (float *)0;
x->x_beg_hrir = (float *)getbytes(x->x_fftsize * x->x_n_ls * sizeof(float));
x->x_beg_hrtf_re = (float **)getbytes(x->x_n_ls * sizeof(float *));
x->x_beg_hrtf_im = (float **)getbytes(x->x_n_ls * sizeof(float *));
x->x_pi_over_180 = 4.0 * atan(1.0) / 180.0;
bin_ambi_calc_HRTF_init_cos(x);
outlet_new(&x->x_obj, &s_list);
return(x);
}
else
{
post("bin_ambi_calc_HRTF-ERROR: need 4 symbols + 2 floats arguments:");
post(" hrir_name + hrtf_re_name + hrtf_im_name + hrir_fade_out_name +");
post(" number_of_loudspeakers + fftsize");
return(0);
}
}
void bin_ambi_calc_HRTF_setup(void)
{
bin_ambi_calc_HRTF_class = class_new(gensym("bin_ambi_calc_HRTF"), (t_newmethod)bin_ambi_calc_HRTF_new, (t_method)bin_ambi_calc_HRTF_free,
sizeof(t_bin_ambi_calc_HRTF), 0, A_GIMME, 0);
class_addmethod(bin_ambi_calc_HRTF_class, (t_method)bin_ambi_calc_HRTF_ls, gensym("ls"), A_GIMME, 0);
class_addmethod(bin_ambi_calc_HRTF_class, (t_method)bin_ambi_calc_HRTF_check_fade_out, gensym("check_fade_out"), 0);
class_addmethod(bin_ambi_calc_HRTF_class, (t_method)bin_ambi_calc_HRTF_load_HRIR, gensym("load_HRIR"), A_FLOAT, 0);
class_addmethod(bin_ambi_calc_HRTF_class, (t_method)bin_ambi_calc_HRTF_check_arrays, gensym("check_arrays"), A_FLOAT, 0);
class_addmethod(bin_ambi_calc_HRTF_class, (t_method)bin_ambi_calc_HRTF_calc_fft, gensym("calc_fft"), A_FLOAT, 0);
class_sethelpsymbol(bin_ambi_calc_HRTF_class, gensym("iemhelp2/help-bin_ambi_calc_HRTF"));
}
/*
Reihenfolge:
n_ls x bin_ambi_calc_HRTF_ls
bin_ambi_calc_HRTF_check_fade_out
n_ls x bin_ambi_calc_HRTF_load_HRIR
n_ls x bin_ambi_calc_HRTF_check_arrays
n_ls x bin_ambi_calc_HRTF_calc_fft
*/
