aboutsummaryrefslogtreecommitdiff
path: root/src/ambi_decode.c
diff options
context:
space:
mode:
Diffstat (limited to 'src/ambi_decode.c')
-rw-r--r--src/ambi_decode.c724
1 files changed, 724 insertions, 0 deletions
diff --git a/src/ambi_decode.c b/src/ambi_decode.c
new file mode 100644
index 0000000..db030cf
--- /dev/null
+++ b/src/ambi_decode.c
@@ -0,0 +1,724 @@
+/* For information on usage and redistribution, and for a DISCLAIMER OF ALL
+* WARRANTIES, see the file, "LICENSE.txt," in this distribution.
+
+iem_ambi written by Thomas Musil, Copyright (c) IEM KUG Graz Austria 2000 - 2005 */
+
+#ifdef NT
+#pragma warning( disable : 4244 )
+#pragma warning( disable : 4305 )
+#endif
+
+
+#include "m_pd.h"
+#include "iemlib.h"
+#include <math.h>
+#include <stdio.h>
+#include <string.h>
+
+
+
+/* -------------------------- ambi_decode ------------------------------ */
+/*
+ ** berechnet ein reduziertes Ambisonic-Decoder-Set in die HRTF-Spektren **
+ ** Inputs: ls + Liste von 3 floats: Index [1 .. 25] + Elevation [-90 .. +90 degree] + Azimut [0 .. 360 degree] **
+ ** Inputs: calc_inv **
+ ** Inputs: load_HRIR + float index1..25 **
+ ** 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 .. LS25_L_HRIR heissen und existieren **
+ ** setzt voraus , dass die HRTF-tabele-names von LS1_HRTF_re .. LS25_HRTF_re heissen und existieren **
+ ** setzt voraus , dass die HRTF-tabele-names von LS1_HRTF_im .. LS25_HRTF_im heissen und existieren **
+ */
+
+typedef struct _ambi_decode
+{
+ t_object x_obj;
+ t_atom *x_at;
+ double *x_inv_work1;
+ double *x_inv_work2;
+ double *x_inv_buf2;
+ double *x_transp;
+ double *x_ls_encode;
+ double *x_prod;
+ double *x_ambi_channel_weight;
+ double x_sing_range;
+ int x_n_ambi;
+ int x_n_order;
+ int x_n_ls;
+ int x_n_phls;
+ int x_n_dim;
+ t_symbol *x_s_matrix;
+ double x_sqrt3;
+ double x_sqrt10_4;
+ double x_sqrt15_2;
+ double x_sqrt6_4;
+ double x_sqrt35_8;
+ double x_sqrt70_4;
+ double x_sqrt5_2;
+ double x_sqrt126_16;
+ double x_sqrt315_8;
+ double x_sqrt105_4;
+ double x_pi_over_180;
+} t_ambi_decode;
+
+static t_class *ambi_decode_class;
+
+static void ambi_decode_copy_row2buf(t_ambi_decode *x, int row)
+{
+ int n_ambi2 = 2*x->x_n_ambi;
+ int i;
+ double *dw=x->x_inv_work2;
+ double *db=x->x_inv_buf2;
+
+ dw += row*n_ambi2;
+ for(i=0; i<n_ambi2; i++)
+ *db++ = *dw++;
+}
+
+static void ambi_decode_copy_buf2row(t_ambi_decode *x, int row)
+{
+ int n_ambi2 = 2*x->x_n_ambi;
+ int i;
+ double *dw=x->x_inv_work2;
+ double *db=x->x_inv_buf2;
+
+ dw += row*n_ambi2;
+ for(i=0; i<n_ambi2; i++)
+ *dw++ = *db++;
+}
+
+static void ambi_decode_copy_row2row(t_ambi_decode *x, int src_row, int dst_row)
+{
+ int n_ambi2 = 2*x->x_n_ambi;
+ int i;
+ double *dw_src=x->x_inv_work2;
+ double *dw_dst=x->x_inv_work2;
+
+ dw_src += src_row*n_ambi2;
+ dw_dst += dst_row*n_ambi2;
+ for(i=0; i<n_ambi2; i++)
+ *dw_dst++ = *dw_src++;
+}
+
+static void ambi_decode_xch_rows(t_ambi_decode *x, int row1, int row2)
+{
+ ambi_decode_copy_row2buf(x, row1);
+ ambi_decode_copy_row2row(x, row2, row1);
+ ambi_decode_copy_buf2row(x, row2);
+}
+
+static void ambi_decode_mul_row(t_ambi_decode *x, int row, double mul)
+{
+ int n_ambi2 = 2*x->x_n_ambi;
+ int i;
+ double *dw=x->x_inv_work2;
+
+ dw += row*n_ambi2;
+ for(i=0; i<n_ambi2; i++)
+ {
+ (*dw) *= mul;
+ dw++;
+ }
+}
+
+static void ambi_decode_mul_buf_and_add2row(t_ambi_decode *x, int row, double mul)
+{
+ int n_ambi2 = 2*x->x_n_ambi;
+ int i;
+ double *dw=x->x_inv_work2;
+ double *db=x->x_inv_buf2;
+
+ dw += row*n_ambi2;
+ for(i=0; i<n_ambi2; i++)
+ {
+ *dw += (*db)*mul;
+ dw++;
+ db++;
+ }
+}
+
+static int ambi_decode_eval_which_element_of_col_not_zero(t_ambi_decode *x, int col, int start_row)
+{
+ int n_ambi = x->x_n_ambi;
+ int n_ambi2 = 2*n_ambi;
+ int i, j;
+ double *dw=x->x_inv_work2;
+ double singrange=x->x_sing_range;
+ int ret=-1;
+
+ dw += start_row*n_ambi2 + col;
+ j = 0;
+ for(i=start_row; i<n_ambi; i++)
+ {
+ if((*dw > singrange) || (*dw < -singrange))
+ {
+ ret = i;
+ i = n_ambi+1;
+ }
+ dw += n_ambi2;
+ }
+ return(ret);
+}
+
+static void ambi_decode_mul1(t_ambi_decode *x)
+{
+ double *vec1, *beg1=x->x_ls_encode;
+ double *vec2, *beg2=x->x_ls_encode;
+ double *inv=x->x_inv_work1;
+ double sum;
+ int n_ls=x->x_n_ls+x->x_n_phls;
+ int n_ambi=x->x_n_ambi;
+ int i, j, k;
+
+ for(k=0; k<n_ambi; k++)
+ {
+ beg2=x->x_ls_encode;
+ for(j=0; j<n_ambi; j++)
+ {
+ sum = 0.0;
+ vec1 = beg1;
+ vec2 = beg2;
+ for(i=0; i<n_ls; i++)
+ {
+ sum += *vec1++ * *vec2++;
+ }
+ beg2 += n_ls;
+ *inv++ = sum;
+ }
+ beg1 += n_ls;
+ }
+}
+
+static void ambi_decode_mul2(t_ambi_decode *x)
+{
+ int n_ls=x->x_n_ls+x->x_n_phls;
+ int n_ambi=x->x_n_ambi;
+ int n_ambi2=2*n_ambi;
+ int i, j, k;
+ double *vec1, *beg1=x->x_transp;
+ double *vec2, *beg2=x->x_inv_work2+n_ambi;
+ double *vec3=x->x_prod;
+ double *acw_vec=x->x_ambi_channel_weight;
+ double sum;
+
+ for(k=0; k<n_ls; k++)
+ {
+ beg2=x->x_inv_work2+n_ambi;
+ for(j=0; j<n_ambi; j++)
+ {
+ sum = 0.0;
+ vec1 = beg1;
+ vec2 = beg2;
+ for(i=0; i<n_ambi; i++)
+ {
+ sum += *vec1++ * *vec2;
+ vec2 += n_ambi2;
+ }
+ beg2++;
+ *vec3++ = sum * acw_vec[j];
+ }
+ beg1 += n_ambi;
+ }
+}
+
+static void ambi_decode_transp_back(t_ambi_decode *x)
+{
+ double *vec, *transp=x->x_transp;
+ double *straight=x->x_ls_encode;
+ int n_ls=x->x_n_ls+x->x_n_phls;
+ int n_ambi=x->x_n_ambi;
+ int i, j;
+
+ for(j=0; j<n_ambi; j++)
+ {
+ vec = transp;
+ for(i=0; i<n_ls; i++)
+ {
+ *straight++ = *vec;
+ vec += n_ambi;
+ }
+ transp++;
+ }
+}
+
+static void ambi_decode_inverse(t_ambi_decode *x)
+{
+ int n_ambi = x->x_n_ambi;
+ int n_ambi2 = 2*n_ambi;
+ int i, j, nz;
+ int r,c;
+ double *src=x->x_inv_work1;
+ double *db=x->x_inv_work2;
+ double rcp, *dv;
+
+ dv = db;
+ for(i=0; i<n_ambi; i++) /* init */
+ {
+ for(j=0; j<n_ambi; j++)
+ {
+ *dv++ = *src++;
+ }
+ for(j=0; j<n_ambi; j++)
+ {
+ if(j == i)
+ *dv++ = 1.0;
+ else
+ *dv++ = 0.0;
+ }
+ }
+
+ /* make 1 in main-diagonale, and 0 below */
+ for(i=0; i<n_ambi; i++)
+ {
+ nz = ambi_decode_eval_which_element_of_col_not_zero(x, i, i);
+ if(nz < 0)
+ {
+ post("ambi_decode ERROR: matrix not regular !!!!");
+ return;
+ }
+ else
+ {
+ if(nz != i)
+ ambi_decode_xch_rows(x, i, nz);
+ dv = db + i*n_ambi2 + i;
+ rcp = 1.0 /(*dv);
+ ambi_decode_mul_row(x, i, rcp);
+ ambi_decode_copy_row2buf(x, i);
+ for(j=i+1; j<n_ambi; j++)
+ {
+ dv += n_ambi2;
+ rcp = -(*dv);
+ ambi_decode_mul_buf_and_add2row(x, j, rcp);
+ }
+ }
+ }
+
+ /* make 0 above the main diagonale */
+ for(i=n_ambi-1; i>=0; i--)
+ {
+ dv = db + i*n_ambi2 + i;
+ ambi_decode_copy_row2buf(x, i);
+ for(j=i-1; j>=0; j--)
+ {
+ dv -= n_ambi2;
+ rcp = -(*dv);
+ ambi_decode_mul_buf_and_add2row(x, j, rcp);
+ }
+ }
+
+ post("matrix_inverse regular");
+}
+
+static void ambi_decode_pinv(t_ambi_decode *x)
+{
+ t_atom *at=x->x_at;
+ int i, n=x->x_n_ls*x->x_n_ambi;
+ double *dv=x->x_prod;
+
+ ambi_decode_transp_back(x);
+ ambi_decode_mul1(x);
+ ambi_decode_inverse(x);
+ ambi_decode_mul2(x);
+ at += 2;
+ for(i=0; i<n; i++)
+ {
+ SETFLOAT(at, (float)(*dv));
+ dv++;
+ at++;
+ }
+ outlet_anything(x->x_obj.ob_outlet, x->x_s_matrix, n+2, x->x_at);
+}
+
+static void ambi_decode_encode_ls_2d(t_ambi_decode *x, int argc, t_atom *argv, int ls0_ph1)
+{
+ double phi;
+ double *dw = x->x_transp;
+ int index;
+ int n_ls=x->x_n_ls;
+ int n_phls=x->x_n_phls;
+ int order=x->x_n_order;
+
+ if(argc < 2)
+ {
+ post("ambi_decode 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(ls0_ph1)
+ {
+ if(n_phls)
+ {
+ if(index >= n_phls)
+ index = n_phls - 1;
+ index += n_ls;
+ }
+ else
+ return;
+ }
+ else
+ {
+ if(index >= n_ls)
+ index = n_ls - 1;
+ }
+
+ phi *= x->x_pi_over_180;
+
+ dw += index * x->x_n_ambi;
+
+ *dw++ = 1.0;
+ *dw++ = cos(phi);
+ *dw++ = sin(phi);
+
+ if(order >= 2)
+ {
+ *dw++ = cos(2.0*phi);
+ *dw++ = sin(2.0*phi);
+
+ if(order >= 3)
+ {
+ *dw++ = cos(3.0*phi);
+ *dw++ = sin(3.0*phi);
+ if(order >= 4)
+ {
+ *dw++ = cos(4.0*phi);
+ *dw++ = sin(4.0*phi);
+
+ if(order >= 5)
+ {
+ *dw++ = cos(5.0*phi);
+ *dw++ = sin(5.0*phi);
+
+ if(order >= 6)
+ {
+ *dw++ = cos(6.0*phi);
+ *dw++ = sin(6.0*phi);
+
+ if(order >= 7)
+ {
+ *dw++ = cos(7.0*phi);
+ *dw++ = sin(7.0*phi);
+
+ if(order >= 8)
+ {
+ *dw++ = cos(8.0*phi);
+ *dw++ = sin(8.0*phi);
+
+ if(order >= 9)
+ {
+ *dw++ = cos(9.0*phi);
+ *dw++ = sin(9.0*phi);
+
+ if(order >= 10)
+ {
+ *dw++ = cos(10.0*phi);
+ *dw++ = sin(10.0*phi);
+
+ if(order >= 11)
+ {
+ *dw++ = cos(11.0*phi);
+ *dw++ = sin(11.0*phi);
+
+ if(order >= 12)
+ {
+ *dw++ = cos(12.0*phi);
+ *dw++ = sin(12.0*phi);
+ }
+ }
+ }
+ }
+ }
+ }
+ }
+ }
+ }
+ }
+ }
+}
+
+static void ambi_decode_encode_ls_3d(t_ambi_decode *x, int argc, t_atom *argv, int ls0_ph1)
+{
+ double delta, phi;
+ double cd, sd, cd2, cd3, sd2, csd, cp, sp, cp2, sp2, cp3, sp3, cp4, sp4;
+ double *dw = x->x_transp;
+ int index;
+ int n_ls=x->x_n_ls;
+ int n_phls=x->x_n_phls;
+ int order=x->x_n_order;
+
+ if(argc < 3)
+ {
+ post("ambi_decode 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(ls0_ph1)
+ {
+ if(n_phls)
+ {
+ if(index >= n_phls)
+ index = n_phls - 1;
+ index += n_ls;
+ }
+ else
+ return;
+ }
+ else
+ {
+ if(index >= n_ls)
+ index = n_ls - 1;
+ }
+
+ delta *= x->x_pi_over_180;
+ phi *= x->x_pi_over_180;
+
+ dw += index * x->x_n_ambi;
+
+ cd = cos(delta);
+ sd = sin(delta);
+ cp = cos(phi);
+ sp = sin(phi);
+
+
+ *dw++ = 1.0;
+ *dw++ = cd * cp;
+ *dw++ = cd * sp;
+ *dw++ = sd;
+
+ if(order >= 2)
+ {
+ cp2 = cos(2.0*phi);
+ sp2 = sin(2.0*phi);
+ cd2 = cd * cd;
+ sd2 = sd * sd;
+ csd = cd * sd;
+ *dw++ = 0.5 * x->x_sqrt3 * cd2 * cp2;
+ *dw++ = 0.5 * x->x_sqrt3 * cd2 * sp2;
+ *dw++ = x->x_sqrt3 * csd * cp;
+ *dw++ = x->x_sqrt3 * csd * sp;
+ *dw++ = 0.5 * (3.0 * sd2 - 1.0);
+
+ if(order >= 3)
+ {
+ cp3 = cos(3.0*phi);
+ sp3 = sin(3.0*phi);
+ cd3 = cd2 * cd;
+ *dw++ = x->x_sqrt10_4 * cd3 * cp3;
+ *dw++ = x->x_sqrt10_4 * cd3 * sp3;
+ *dw++ = x->x_sqrt15_2 * cd * csd * cp2;
+ *dw++ = x->x_sqrt15_2 * cd * csd * sp2;
+ *dw++ = x->x_sqrt6_4 * cd * (5.0 * sd2 - 1.0) * cp;
+ *dw++ = x->x_sqrt6_4 * cd * (5.0 * sd2 - 1.0) * sp;
+ *dw++ = 0.5 * sd * (5.0 * sd2 - 3.0);
+
+ if(order >= 4)
+ {
+ cp4 = cos(4.0*phi);
+ sp4 = sin(4.0*phi);
+ *dw++ = x->x_sqrt35_8 * cd2 * cd2 * cp4;
+ *dw++ = x->x_sqrt35_8 * cd2 * cd2 * sp4;
+ *dw++ = x->x_sqrt70_4 * cd2 * csd * cp3;
+ *dw++ = x->x_sqrt70_4 * cd2 * csd * sp3;
+ *dw++ = 0.5 * x->x_sqrt5_2 * cd2 * (7.0 * sd2 - 1.0) * cp2;
+ *dw++ = 0.5 * x->x_sqrt5_2 * cd2 * (7.0 * sd2 - 1.0) * sp2;
+ *dw++ = x->x_sqrt10_4 * csd * (7.0 * sd2 - 3.0) * cp;
+ *dw++ = x->x_sqrt10_4 * csd * (7.0 * sd2 - 3.0) * sp;
+ *dw++ = 0.125 * (sd2 * (35.0 * sd2 - 30.0) + 3.0);
+
+ if(order >= 5)
+ {
+ *dw++ = x->x_sqrt126_16 * cd3 * cd2 * cos(5.0*phi);
+ *dw++ = x->x_sqrt126_16 * cd3 * cd2 * sin(5.0*phi);
+ *dw++ = x->x_sqrt315_8 * cd3 * csd * cp4;
+ *dw++ = x->x_sqrt315_8 * cd3 * csd * sp4;
+ *dw++ = 0.25 * x->x_sqrt70_4 * cd3 * (9.0 * sd2 - 1.0) * cp3;
+ *dw++ = 0.25 * x->x_sqrt70_4 * cd3 * (9.0 * sd2 - 1.0) * sp3;
+ *dw++ = x->x_sqrt105_4 * cd * csd * (3.0 * sd2 - 1.0) * cp2;
+ *dw++ = x->x_sqrt105_4 * cd * csd * (3.0 * sd2 - 1.0) * sp2;
+ *dw++ = 0.25 * x->x_sqrt15_2 * cd * (sd2 * (21.0 * sd2 - 14.0) + 1.0) * cp;
+ *dw++ = 0.25 * x->x_sqrt15_2 * cd * (sd2 * (21.0 * sd2 - 14.0) + 1.0) * sp;
+ *dw = 0.125 * sd * (sd2 * (63.0 * sd2 - 70.0) + 15.0);
+ }
+ }
+ }
+ }
+}
+
+static void ambi_decode_ls(t_ambi_decode *x, t_symbol *s, int argc, t_atom *argv)
+{
+ if(x->x_n_dim == 2)
+ ambi_decode_encode_ls_2d(x, argc, argv, 0);
+ else
+ ambi_decode_encode_ls_3d(x, argc, argv, 0);
+}
+
+static void ambi_decode_phls(t_ambi_decode *x, t_symbol *s, int argc, t_atom *argv)
+{
+ if(x->x_n_dim == 2)
+ ambi_decode_encode_ls_2d(x, argc, argv, 1);
+ else
+ ambi_decode_encode_ls_3d(x, argc, argv, 1);
+}
+
+static void ambi_decode_ambi_weight(t_ambi_decode *x, t_symbol *s, int argc, t_atom *argv)
+{
+ if(argc > x->x_n_order)
+ {
+ int i, k=0, n=x->x_n_order;
+ double d;
+
+ x->x_ambi_channel_weight[k] = atom_getfloat(argv++);
+ k++;
+ if(x->x_n_dim == 2)
+ {
+ for(i=1; i<=n; i++)
+ {
+ d = atom_getfloat(argv++);
+ x->x_ambi_channel_weight[k] = d;
+ k++;
+ x->x_ambi_channel_weight[k] = d;
+ k++;
+ }
+ }
+ else
+ {
+ int j, m;
+
+ for(i=1; i<=n; i++)
+ {
+ d = atom_getfloat(argv++);
+ m = 2*i + 1;
+ for(j=0; j<m; j++)
+ {
+ x->x_ambi_channel_weight[k] = d;
+ k++;
+ }
+ }
+ }
+ }
+ else
+ post("ambi_decode-ERROR: ambi_weight needs %d float weights", x->x_n_order+1);
+}
+
+static void ambi_decode_sing_range(t_ambi_decode *x, t_floatarg f)
+{
+ if(f < 0.0f)
+ x->x_sing_range = -(double)f;
+ else
+ x->x_sing_range = (double)f;
+}
+
+static void ambi_decode_free(t_ambi_decode *x)
+{
+ freebytes(x->x_inv_work1, x->x_n_ambi * x->x_n_ambi * sizeof(double));
+ freebytes(x->x_inv_work2, 2 * x->x_n_ambi * x->x_n_ambi * sizeof(double));
+ freebytes(x->x_inv_buf2, 2 * x->x_n_ambi * sizeof(double));
+ freebytes(x->x_transp, (x->x_n_ls+x->x_n_phls) * x->x_n_ambi * sizeof(double));
+ freebytes(x->x_ls_encode, (x->x_n_ls+x->x_n_phls) * x->x_n_ambi * sizeof(double));
+ freebytes(x->x_prod, (x->x_n_ls+x->x_n_phls) * x->x_n_ambi * sizeof(double));
+ freebytes(x->x_ambi_channel_weight, x->x_n_ambi * sizeof(double));
+ freebytes(x->x_at, (x->x_n_ls * x->x_n_ambi + 2) * sizeof(t_atom));
+}
+
+static void *ambi_decode_new(t_symbol *s, int argc, t_atom *argv)
+{
+ t_ambi_decode *x = (t_ambi_decode *)pd_new(ambi_decode_class);
+ int nls, order, dim, i;
+ int nphls=0;/* phantom_loudspeaker */
+
+ if(argc < 3)
+ {
+ post("ambi_decode-ERROR: need following arguments: ambi_order dimension number_of_loudspeakers (number_of_phantom_speakers)");
+ return(0);
+ }
+ else
+ {
+ order = (int)atom_getint(argv++);
+ dim = (int)atom_getint(argv++);
+ nls = (int)atom_getint(argv++);
+ if((argc > 3)&&IS_A_FLOAT(argv,0))
+ nphls=(int)atom_getint(argv);
+
+ if(order < 1)
+ order = 1;
+ if(dim != 3)
+ {
+ dim = 2;
+ if(order > 12)
+ order = 12;
+ x->x_n_ambi = 2*order + 1;
+ }
+ else
+ {
+ if(order > 5)
+ order = 5;
+ x->x_n_ambi = (order + 1)*(order + 1);
+ }
+ x->x_n_dim = dim;
+ x->x_n_order = order;
+ if(nls < 1)
+ nls = 1;
+ if(nphls < 0)
+ nphls = 0;
+ if(nls < x->x_n_ambi)
+ post("ambi_decode-WARNING: Number of Loudspeakers < Number of Ambisonic-Channels !!!!");
+ if(nphls > nls)
+ {
+ post("ambi_decode-WARNING: Number of Phantom-Loudspeakers > Number of Loudspeakers !!!!");
+ nphls = nls;
+ }
+ x->x_n_ls = nls;
+ x->x_n_phls = nphls;
+ x->x_inv_work1 = (double *)getbytes(x->x_n_ambi * x->x_n_ambi * sizeof(double));
+ x->x_inv_work2 = (double *)getbytes(2 * x->x_n_ambi * x->x_n_ambi * sizeof(double));
+ x->x_inv_buf2 = (double *)getbytes(2 * x->x_n_ambi * sizeof(double));
+ x->x_transp = (double *)getbytes((x->x_n_ls+x->x_n_phls) * x->x_n_ambi * sizeof(double));
+ x->x_ls_encode = (double *)getbytes((x->x_n_ls+x->x_n_phls) * x->x_n_ambi * sizeof(double));
+ x->x_prod = (double *)getbytes((x->x_n_ls+x->x_n_phls) * x->x_n_ambi * sizeof(double));
+ x->x_ambi_channel_weight = (double *)getbytes(x->x_n_ambi * sizeof(double));
+ x->x_at = (t_atom *)getbytes((x->x_n_ls * x->x_n_ambi + 2) * sizeof(t_atom));
+ x->x_s_matrix = gensym("matrix");
+ /*change*/
+ SETFLOAT(x->x_at, (float)x->x_n_ls);
+ SETFLOAT(x->x_at+1, (float)x->x_n_ambi);
+ x->x_sqrt3 = sqrt(3.0);
+ x->x_sqrt5_2 = sqrt(5.0) / 2.0;
+ x->x_sqrt6_4 = sqrt(6.0) / 4.0;
+ x->x_sqrt10_4 = sqrt(10.0) / 4.0;
+ x->x_sqrt15_2 = sqrt(15.0) / 2.0;
+ x->x_sqrt35_8 = sqrt(35.0) / 8.0;
+ x->x_sqrt70_4 = sqrt(70.0) / 4.0;
+ x->x_sqrt126_16 = sqrt(126.0) / 16.0;
+ x->x_sqrt315_8 = sqrt(315.0) / 8.0;
+ x->x_sqrt105_4 = sqrt(105.0) / 4.0;
+ x->x_pi_over_180 = 4.0 * atan(1.0) / 180.0;
+ x->x_sing_range = 1.0e-10;
+ for(i=0; i<x->x_n_ambi; i++)
+ x->x_ambi_channel_weight[i] = 1.0;
+ outlet_new(&x->x_obj, &s_list);
+ return (x);
+ }
+}
+
+void ambi_decode_setup(void)
+{
+ ambi_decode_class = class_new(gensym("ambi_decode"), (t_newmethod)ambi_decode_new, (t_method)ambi_decode_free,
+ sizeof(t_ambi_decode), 0, A_GIMME, 0);
+ class_addmethod(ambi_decode_class, (t_method)ambi_decode_ls, gensym("ls"), A_GIMME, 0);
+ class_addmethod(ambi_decode_class, (t_method)ambi_decode_phls, gensym("phls"), A_GIMME, 0);
+ class_addmethod(ambi_decode_class, (t_method)ambi_decode_ambi_weight, gensym("ambi_weight"), A_GIMME, 0);
+ class_addmethod(ambi_decode_class, (t_method)ambi_decode_sing_range, gensym("sing_range"), A_DEFFLOAT, 0);
+ class_addmethod(ambi_decode_class, (t_method)ambi_decode_pinv, gensym("pinv"), 0);
+ class_sethelpsymbol(ambi_decode_class, gensym("iemhelp2/help-ambi_decode"));
+}
generated by cgit v1.2.1 (git 2.18.0) at 2024-04-29 09:59:12 +0000