Lines Matching refs:st
122 void speex_echo_get_residual(SpeexEchoState *st, spx_word32_t *Yout, int len);
155 spx_word32_t Davg1; /* 1st recursive average of the residual power difference */
157 spx_float_t Dvar1; /* Estimated variance of 1st estimator */
404 SpeexEchoState *st = (SpeexEchoState *)speex_alloc(sizeof(SpeexEchoState));
406 st->K = nb_speakers;
407 st->C = nb_mic;
408 C=st->C;
409 K=st->K;
418 st->frame_size = frame_size;
419 st->window_size = 2*frame_size;
420 N = st->window_size;
421 M = st->M = (filter_length+st->frame_size-1)/frame_size;
422 st->cancel_count=0;
423 st->sum_adapt = 0;
424 st->saturated = 0;
425 st->screwed_up = 0;
427 st->sampling_rate = 8000;
428 st->spec_average = DIV32_16(SHL32(EXTEND32(st->frame_size), 15), st->sampling_rate);
430 st->beta0 = DIV32_16(SHL32(EXTEND32(st->frame_size), 16), st->sampling_rate);
431 st->beta_max = DIV32_16(SHL32(EXTEND32(st->frame_size), 14), st->sampling_rate);
433 st->beta0 = (2.0f*st->frame_size)/st->sampling_rate;
434 st->beta_max = (.5f*st->frame_size)/st->sampling_rate;
436 st->leak_estimate = 0;
438 st->fft_table = spx_fft_init(N);
440 st->e = (spx_word16_t*)speex_alloc(C*N*sizeof(spx_word16_t));
441 st->x = (spx_word16_t*)speex_alloc(K*N*sizeof(spx_word16_t));
442 st->input = (spx_word16_t*)speex_alloc(C*st->frame_size*sizeof(spx_word16_t));
443 st->y = (spx_word16_t*)speex_alloc(C*N*sizeof(spx_word16_t));
444 st->last_y = (spx_word16_t*)speex_alloc(C*N*sizeof(spx_word16_t));
445 st->Yf = (spx_word32_t*)speex_alloc((st->frame_size+1)*sizeof(spx_word32_t));
446 st->Rf = (spx_word32_t*)speex_alloc((st->frame_size+1)*sizeof(spx_word32_t));
447 st->Xf = (spx_word32_t*)speex_alloc((st->frame_size+1)*sizeof(spx_word32_t));
448 st->Yh = (spx_word32_t*)speex_alloc((st->frame_size+1)*sizeof(spx_word32_t));
449 st->Eh = (spx_word32_t*)speex_alloc((st->frame_size+1)*sizeof(spx_word32_t));
451 st->X = (spx_word16_t*)speex_alloc(K*(M+1)*N*sizeof(spx_word16_t));
452 st->Y = (spx_word16_t*)speex_alloc(C*N*sizeof(spx_word16_t));
453 st->E = (spx_word16_t*)speex_alloc(C*N*sizeof(spx_word16_t));
454 st->W = (spx_word32_t*)speex_alloc(C*K*M*N*sizeof(spx_word32_t));
456 st->foreground = (spx_word16_t*)speex_alloc(M*N*C*K*sizeof(spx_word16_t));
458 st->PHI = (spx_word32_t*)speex_alloc(N*sizeof(spx_word32_t));
459 st->power = (spx_word32_t*)speex_alloc((frame_size+1)*sizeof(spx_word32_t));
460 st->power_1 = (spx_float_t*)speex_alloc((frame_size+1)*sizeof(spx_float_t));
461 st->window = (spx_word16_t*)speex_alloc(N*sizeof(spx_word16_t));
462 st->prop = (spx_word16_t*)speex_alloc(M*sizeof(spx_word16_t));
463 st->wtmp = (spx_word16_t*)speex_alloc(N*sizeof(spx_word16_t));
465 st->wtmp2 = (spx_word16_t*)speex_alloc(N*sizeof(spx_word16_t));
468 st->window[i] = (16383-SHL16(spx_cos(DIV32_16(MULT16_16(25736,i<<1),N)),1));
469 st->window[N-i-1] = st->window[i];
473 st->window[i] = .5-.5*cos(2*M_PI*i/N);
475 for (i=0;i<=st->frame_size;i++)
476 st->power_1[i] = FLOAT_ONE;
478 st->W[i] = 0;
483 st->prop[0] = QCONST16(.7, 15);
484 sum = EXTEND32(st->prop[0]);
487 st->prop[i] = MULT16_16_Q15(st->prop[i-1], decay);
488 sum = ADD32(sum, EXTEND32(st->prop[i]));
492 st->prop[i] = DIV32(MULT16_16(QCONST16(.8f,15), st->prop[i]),sum);
496 st->memX = (spx_word16_t*)speex_alloc(K*sizeof(spx_word16_t));
497 st->memD = (spx_word16_t*)speex_alloc(C*sizeof(spx_word16_t));
498 st->memE = (spx_word16_t*)speex_alloc(C*sizeof(spx_word16_t));
499 st->preemph = QCONST16(.9,15);
500 if (st->sampling_rate<12000)
501 st->notch_radius = QCONST16(.9, 15);
502 else if (st->sampling_rate<24000)
503 st->notch_radius = QCONST16(.982, 15);
505 st->notch_radius = QCONST16(.992, 15);
507 st->notch_mem = (spx_mem_t*)speex_alloc(2*C*sizeof(spx_mem_t));
508 st->adapted = 0;
509 st->Pey = st->Pyy = FLOAT_ONE;
512 st->Davg1 = st->Davg2 = 0;
513 st->Dvar1 = st->Dvar2 = FLOAT_ZERO;
516 st->play_buf = (spx_int16_t*)speex_alloc(K*(PLAYBACK_DELAY+1)*st->frame_size*sizeof(spx_int16_t));
517 st->play_buf_pos = PLAYBACK_DELAY*st->frame_size;
518 st->play_buf_started = 0;
520 return st;
524 EXPORT void speex_echo_state_reset(SpeexEchoState *st)
527 st->cancel_count=0;
528 st->screwed_up = 0;
529 N = st->window_size;
530 M = st->M;
531 C=st->C;
532 K=st->K;
534 st->W[i] = 0;
537 st->foreground[i] = 0;
540 st->X[i] = 0;
541 for (i=0;i<=st->frame_size;i++)
543 st->power[i] = 0;
544 st->power_1[i] = FLOAT_ONE;
545 st->Eh[i] = 0;
546 st->Yh[i] = 0;
548 for (i=0;i<st->frame_size;i++)
550 st->last_y[i] = 0;
554 st->E[i] = 0;
558 st->x[i] = 0;
561 st->notch_mem[i] = 0;
563 st->memD[i]=st->memE[i]=0;
565 st->memX[i]=0;
567 st->saturated = 0;
568 st->adapted = 0;
569 st->sum_adapt = 0;
570 st->Pey = st->Pyy = FLOAT_ONE;
572 st->Davg1 = st->Davg2 = 0;
573 st->Dvar1 = st->Dvar2 = FLOAT_ZERO;
575 for (i=0;i<3*st->frame_size;i++)
576 st->play_buf[i] = 0;
577 st->play_buf_pos = PLAYBACK_DELAY*st->frame_size;
578 st->play_buf_started = 0;
583 EXPORT void speex_echo_state_destroy(SpeexEchoState *st)
585 spx_fft_destroy(st->fft_table);
587 speex_free(st->e);
588 speex_free(st->x);
589 speex_free(st->input);
590 speex_free(st->y);
591 speex_free(st->last_y);
592 speex_free(st->Yf);
593 speex_free(st->Rf);
594 speex_free(st->Xf);
595 speex_free(st->Yh);
596 speex_free(st->Eh);
598 speex_free(st->X);
599 speex_free(st->Y);
600 speex_free(st->E);
601 speex_free(st->W);
603 speex_free(st->foreground);
605 speex_free(st->PHI);
606 speex_free(st->power);
607 speex_free(st->power_1);
608 speex_free(st->window);
609 speex_free(st->prop);
610 speex_free(st->wtmp);
612 speex_free(st->wtmp2);
614 speex_free(st->memX);
615 speex_free(st->memD);
616 speex_free(st->memE);
617 speex_free(st->notch_mem);
619 speex_free(st->play_buf);
620 speex_free(st);
630 EXPORT void speex_echo_capture(SpeexEchoState *st, const spx_int16_t *rec, spx_int16_t *out)
633 /*speex_warning_int("capture with fill level ", st->play_buf_pos/st->frame_size);*/
634 st->play_buf_started = 1;
635 if (st->play_buf_pos>=st->frame_size)
637 speex_echo_cancellation(st, rec, st->play_buf, out);
638 st->play_buf_pos -= st->frame_size;
639 for (i=0;i<st->play_buf_pos;i++)
640 st->play_buf[i] = st->play_buf[i+st->frame_size];
643 if (st->play_buf_pos!=0)
646 st->play_buf_pos = 0;
648 for (i=0;i<st->frame_size;i++)
653 EXPORT void speex_echo_playback(SpeexEchoState *st, const spx_int16_t *play)
655 /*speex_warning_int("playback with fill level ", st->play_buf_pos/st->frame_size);*/
656 if (!st->play_buf_started)
661 if (st->play_buf_pos<=PLAYBACK_DELAY*st->frame_size)
664 for (i=0;i<st->frame_size;i++)
665 st->play_buf[st->play_buf_pos+i] = play[i];
666 st->play_buf_pos += st->frame_size;
667 if (st->play_buf_pos <= (PLAYBACK_DELAY-1)*st->frame_size)
670 for (i=0;i<st->frame_size;i++)
671 st->play_buf[st->play_buf_pos+i] = play[i];
672 st->play_buf_pos += st->frame_size;
680 EXPORT void speex_echo_cancel(SpeexEchoState *st, const spx_int16_t *in, const spx_int16_t *far_end, spx_int16_t *out, spx_int32_t *Yout)
682 speex_echo_cancellation(st, in, far_end, out);
686 EXPORT void speex_echo_cancellation(SpeexEchoState *st, const spx_int16_t *in, const spx_int16_t *far_end, spx_int16_t *out)
702 N = st->window_size;
703 M = st->M;
704 C = st->C;
705 K = st->K;
707 st->cancel_count++;
719 filter_dc_notch16(in+chan, st->notch_radius, st->input+chan*st->frame_size, st->frame_size, st->notch_mem+2*chan, C);
722 for (i=0;i<st->frame_size;i++)
726 tmp32 = SUB32(EXTEND32(st->input[chan*st->frame_size+i]), EXTEND32(MULT16_16_P15(st->preemph, st->memD[chan])));
731 if (st->saturated == 0)
732 st->saturated = 1;
737 if (st->saturated == 0)
738 st->saturated = 1;
741 st->memD[chan] = st->input[chan*st->frame_size+i];
742 st->input[chan*st->frame_size+i] = EXTRACT16(tmp32);
748 for (i=0;i<st->frame_size;i++)
751 st->x[speak*N+i] = st->x[speak*N+i+st->frame_size];
752 tmp32 = SUB32(EXTEND32(far_end[i*K+speak]), EXTEND32(MULT16_16_P15(st->preemph, st->memX[speak])));
758 st->saturated = M+1;
763 st->saturated = M+1;
766 st->x[speak*N+i+st->frame_size] = EXTRACT16(tmp32);
767 st->memX[speak] = far_end[i*K+speak];
777 st->X[(j+1)*N*K+speak*N+i] = st->X[j*N*K+speak*N+i];
780 spx_fft(st->fft_table, st->x+speak*N, &st->X[speak*N]);
786 Sxx += mdf_inner_prod(st->x+speak*N+st->frame_size, st->x+speak*N+st->frame_size, st->frame_size);
787 power_spectrum_accum(st->X+speak*N, st->Xf, N);
795 spectral_mul_accum16(st->X, st->foreground+chan*N*K*M, st->Y+chan*N, N, M*K);
796 spx_ifft(st->fft_table, st->Y+chan*N, st->e+chan*N);
797 for (i=0;i<st->frame_size;i++)
798 st->e[chan*N+i] = SUB16(st->input[chan*st->frame_size+i], st->e[chan*N+i+st->frame_size]);
799 Sff += mdf_inner_prod(st->e+chan*N, st->e+chan*N, st->frame_size);
805 if (st->adapted)
806 mdf_adjust_prop (st->W, N, M, C*K, st->prop);
808 if (st->saturated == 0)
816 weighted_spectral_mul_conj(st->power_1, FLOAT_SHL(PSEUDOFLOAT(st->prop[j]),-15), &st->X[(j+1)*N*K+speak*N], st->E+chan*N, st->PHI, N);
818 st->W[chan*N*K*M + j*N*K + speak*N + i] += st->PHI[i];
823 st->saturated--;
836 if (j==0 || st->cancel_count%(M-1) == j-1)
840 st->wtmp2[i] = EXTRACT16(PSHR32(st->W[chan*N*K*M + j*N*K + speak*N + i],NORMALIZE_SCALEDOWN+16));
841 spx_ifft(st->fft_table, st->wtmp2, st->wtmp);
842 for (i=0;i<st->frame_size;i++)
844 st->wtmp[i]=0;
846 for (i=st->frame_size;i<N;i++)
848 st->wtmp[i]=SHL16(st->wtmp[i],NORMALIZE_SCALEUP);
850 spx_fft(st->fft_table, st->wtmp, st->wtmp2);
853 st->W[chan*N*K*M + j*N*K + speak*N + i] -= SHL32(EXTEND32(st->wtmp2[i]),16+NORMALIZE_SCALEDOWN-NORMALIZE_SCALEUP-1);
855 spx_ifft(st->fft_table, &st->W[chan*N*K*M + j*N*K + speak*N], st->wtmp);
856 for (i=st->frame_size;i<N;i++)
858 st->wtmp[i]=0;
860 spx_fft(st->fft_table, st->wtmp, &st->W[chan*N*K*M + j*N*K + speak*N]);
868 for (i=0;i<=st->frame_size;i++)
869 st->Rf[i] = st->Yf[i] = st->Xf[i] = 0;
877 spectral_mul_accum(st->X, st->W+chan*N*K*M, st->Y+chan*N, N, M*K);
878 spx_ifft(st->fft_table, st->Y+chan*N, st->y+chan*N);
879 for (i=0;i<st->frame_size;i++)
880 st->e[chan*N+i] = SUB16(st->e[chan*N+i+st->frame_size], st->y[chan*N+i+st->frame_size]);
881 Dbf += 10+mdf_inner_prod(st->e+chan*N, st->e+chan*N, st->frame_size);
882 for (i=0;i<st->frame_size;i++)
883 st->e[chan*N+i] = SUB16(st->input[chan*st->frame_size+i], st->y[chan*N+i+st->frame_size]);
884 See += mdf_inner_prod(st->e+chan*N, st->e+chan*N, st->frame_size);
896 st->Davg1 = ADD32(MULT16_32_Q15(QCONST16(.6f,15),st->Davg1), MULT16_32_Q15(QCONST16(.4f,15),SUB32(Sff,See)));
897 st->Davg2 = ADD32(MULT16_32_Q15(QCONST16(.85f,15),st->Davg2), MULT16_32_Q15(QCONST16(.15f,15),SUB32(Sff,See)));
898 st->Dvar1 = FLOAT_ADD(FLOAT_MULT(VAR1_SMOOTH, st->Dvar1), FLOAT_MUL32U(MULT16_32_Q15(QCONST16(.4f,15),Sff), MULT16_32_Q15(QCONST16(.4f,15),Dbf)));
899 st->Dvar2 = FLOAT_ADD(FLOAT_MULT(VAR2_SMOOTH, st->Dvar2), FLOAT_MUL32U(MULT16_32_Q15(QCONST16(.15f,15),Sff), MULT16_32_Q15(QCONST16(.15f,15),Dbf)));
902 st->Davg1 = .6*st->Davg1 + .4*(Sff-See);
903 st->Davg2 = .85*st->Davg2 + .15*(Sff-See);
904 st->Dvar1 = .36*st->Dvar1 + .16*Sff*Dbf;
905 st->Dvar2 = .7225*st->Dvar2 + .0225*Sff*Dbf;
913 else if (FLOAT_GT(FLOAT_MUL32U(st->Davg1, ABS32(st->Davg1)), FLOAT_MULT(VAR1_UPDATE,(st->Dvar1))))
915 else if (FLOAT_GT(FLOAT_MUL32U(st->Davg2, ABS32(st->Davg2)), FLOAT_MULT(VAR2_UPDATE,(st->Dvar2))))
921 st->Davg1 = st->Davg2 = 0;
922 st->Dvar1 = st->Dvar2 = FLOAT_ZERO;
925 st->foreground[i] = EXTRACT16(PSHR32(st->W[i],16));
928 for (i=0;i<st->frame_size;i++)
929 st->e[chan*N+i+st->frame_size] = MULT16_16_Q15(st->window[i+st->frame_size],st->e[chan*N+i+st->frame_size]) + MULT16_16_Q15(st->window[i],st->y[chan*N+i+st->frame_size]);
935 if (FLOAT_GT(FLOAT_MUL32U(NEG32(st->Davg1), ABS32(st->Davg1)), FLOAT_MULT(VAR_BACKTRACK,st->Dvar1)))
937 if (FLOAT_GT(FLOAT_MUL32U(NEG32(st->Davg2), ABS32(st->Davg2)), FLOAT_MULT(VAR_BACKTRACK,st->Dvar2)))
943 st->W[i] = SHL32(EXTEND32(st->foreground[i]),16);
947 for (i=0;i<st->frame_size;i++)
948 st->y[chan*N+i+st->frame_size] = st->e[chan*N+i+st->frame_size];
949 for (i=0;i<st->frame_size;i++)
950 st->e[chan*N+i] = SUB16(st->input[chan*st->frame_size+i], st->y[chan*N+i+st->frame_size]);
953 st->Davg1 = st->Davg2 = 0;
954 st->Dvar1 = st->Dvar2 = FLOAT_ZERO;
963 for (i=0;i<st->frame_size;i++)
967 tmp_out = SUB32(EXTEND32(st->input[chan*st->frame_size+i]), EXTEND32(st->e[chan*N+i+st->frame_size]));
969 tmp_out = SUB32(EXTEND32(st->input[chan*st->frame_size+i]), EXTEND32(st->y[chan*N+i+st->frame_size]));
971 tmp_out = ADD32(tmp_out, EXTEND32(MULT16_16_P15(st->preemph, st->memE[chan])));
975 if (st->saturated == 0)
976 st->saturated = 1;
979 st->memE[chan] = tmp_out;
983 dump_audio(in, far_end, out, st->frame_size);
987 for (i=0;i<st->frame_size;i++)
989 st->e[chan*N+i+st->frame_size] = st->e[chan*N+i];
990 st->e[chan*N+i] = 0;
995 Sey += mdf_inner_prod(st->e+chan*N+st->frame_size, st->y+chan*N+st->frame_size, st->frame_size);
996 Syy += mdf_inner_prod(st->y+chan*N+st->frame_size, st->y+chan*N+st->frame_size, st->frame_size);
997 Sdd += mdf_inner_prod(st->input+chan*st->frame_size, st->input+chan*st->frame_size, st->frame_size);
1000 spx_fft(st->fft_table, st->e+chan*N, st->E+chan*N);
1001 for (i=0;i<st->frame_size;i++)
1002 st->y[i+chan*N] = 0;
1003 spx_fft(st->fft_table, st->y+chan*N, st->Y+chan*N);
1006 power_spectrum_accum(st->E+chan*N, st->Rf, N);
1007 power_spectrum_accum(st->Y+chan*N, st->Yf, N);
1011 /*printf ("%f %f %f %f\n", Sff, See, Syy, Sdd, st->update_cond);*/
1021 st->screwed_up += 50;
1022 for (i=0;i<st->frame_size*C;i++)
1027 st->screwed_up++;
1030 st->screwed_up=0;
1032 if (st->screwed_up>=50)
1035 speex_echo_state_reset(st);
1044 Sxx += mdf_inner_prod(st->x+speak*N+st->frame_size, st->x+speak*N+st->frame_size, st->frame_size);
1045 power_spectrum_accum(st->X+speak*N, st->Xf, N);
1050 for (j=0;j<=st->frame_size;j++)
1051 st->power[j] = MULT16_32_Q15(ss_1,st->power[j]) + 1 + MULT16_32_Q15(ss,st->Xf[j]);
1054 for (j=st->frame_size;j>=0;j--)
1057 Eh = PSEUDOFLOAT(st->Rf[j] - st->Eh[j]);
1058 Yh = PSEUDOFLOAT(st->Yf[j] - st->Yh[j]);
1062 st->Eh[j] = MAC16_32_Q15(MULT16_32_Q15(SUB16(32767,st->spec_average),st->Eh[j]), st->spec_average, st->Rf[j]);
1063 st->Yh[j] = MAC16_32_Q15(MULT16_32_Q15(SUB16(32767,st->spec_average),st->Yh[j]), st->spec_average, st->Yf[j]);
1065 st->Eh[j] = (1-st->spec_average)*st->Eh[j] + st->spec_average*st->Rf[j];
1066 st->Yh[j] = (1-st->spec_average)*st->Yh[j] + st->spec_average*st->Yf[j];
1074 tmp32 = MULT16_32_Q15(st->beta0,Syy);
1075 if (tmp32 > MULT16_32_Q15(st->beta_max,See))
1076 tmp32 = MULT16_32_Q15(st->beta_max,See);
1080 st->Pey = FLOAT_ADD(FLOAT_MULT(alpha_1,st->Pey) , FLOAT_MULT(alpha,Pey));
1081 st->Pyy = FLOAT_ADD(FLOAT_MULT(alpha_1,st->Pyy) , FLOAT_MULT(alpha,Pyy));
1082 if (FLOAT_LT(st->Pyy, FLOAT_ONE))
1083 st->Pyy = FLOAT_ONE;
1085 if (FLOAT_LT(st->Pey, FLOAT_MULT(MIN_LEAK,st->Pyy)))
1086 st->Pey = FLOAT_MULT(MIN_LEAK,st->Pyy);
1087 if (FLOAT_GT(st->Pey, st->Pyy))
1088 st->Pey = st->Pyy;
1090 st->leak_estimate = FLOAT_EXTRACT16(FLOAT_SHL(FLOAT_DIVU(st->Pey, st->Pyy),14));
1092 if (st->leak_estimate > 16383)
1093 st->leak_estimate = 32767;
1095 st->leak_estimate = SHL16(st->leak_estimate,1);
1096 /*printf ("%f\n", st->leak_estimate);*/
1100 tmp32 = MULT16_32_Q15(st->leak_estimate,Syy);
1115 RER = (.0001*Sxx + 3.*MULT16_32_Q15(st->leak_estimate,Syy)) / See;
1124 if (!st->adapted && st->sum_adapt > SHL32(EXTEND32(M),15) && MULT16_32_Q15(st->leak_estimate,Syy) > MULT16_32_Q15(QCONST16(.03f,15),Syy))
1126 st->adapted = 1;
1129 if (st->adapted)
1132 for (i=0;i<=st->frame_size;i++)
1136 r = MULT16_32_Q15(st->leak_estimate,SHL32(st->Yf[i],3));
1137 e = SHL32(st->Rf[i],3)+1;
1146 /*st->power_1[i] = adapt_rate*r/(e*(1+st->power[i]));*/
1147 st->power_1[i] = FLOAT_SHL(FLOAT_DIV32_FLOAT(r,FLOAT_MUL32U(e,st->power[i]+10)),WEIGHT_SHIFT+16);
1165 for (i=0;i<=st->frame_size;i++)
1166 st->power_1[i] = FLOAT_SHL(FLOAT_DIV32(EXTEND32(adapt_rate),ADD32(st->power[i],10)),WEIGHT_SHIFT+1);
1170 st->sum_adapt = ADD32(st->sum_adapt,adapt_rate);
1174 for (i=0;i<st->frame_size;i++)
1175 st->last_y[i] = st->last_y[st->frame_size+i];
1176 if (st->adapted)
1179 for (i=0;i<st->frame_size;i++)
1180 st->last_y[st->frame_size+i] = in[i]-out[i];
1184 st->last_y[i] = st->x[i];*/
1190 void speex_echo_get_residual(SpeexEchoState *st, spx_word32_t *residual_echo, int len)
1196 N = st->window_size;
1200 st->y[i] = MULT16_16_Q15(st->window[i],st->last_y[i]);
1203 spx_fft(st->fft_table, st->y, st->Y);
1204 power_spectrum(st->Y, residual_echo, N);
1207 if (st->leak_estimate > 16383)
1210 leak2 = SHL16(st->leak_estimate, 1);
1212 if (st->leak_estimate>.5)
1215 leak2 = 2*st->leak_estimate;
1218 for (i=0;i<=st->frame_size;i++)
1223 EXPORT int speex_echo_ctl(SpeexEchoState *st, int request, void *ptr)
1229 (*(int*)ptr) = st->frame_size;
1232 st->sampling_rate = (*(int*)ptr);
1233 st->spec_average = DIV32_16(SHL32(EXTEND32(st->frame_size), 15), st->sampling_rate);
1235 st->beta0 = DIV32_16(SHL32(EXTEND32(st->frame_size), 16), st->sampling_rate);
1236 st->beta_max = DIV32_16(SHL32(EXTEND32(st->frame_size), 14), st->sampling_rate);
1238 st->beta0 = (2.0f*st->frame_size)/st->sampling_rate;
1239 st->beta_max = (.5f*st->frame_size)/st->sampling_rate;
1241 if (st->sampling_rate<12000)
1242 st->notch_radius = QCONST16(.9, 15);
1243 else if (st->sampling_rate<24000)
1244 st->notch_radius = QCONST16(.982, 15);
1246 st->notch_radius = QCONST16(.992, 15);
1249 (*(int*)ptr) = st->sampling_rate;
1253 *((spx_int32_t *)ptr) = st->M * st->frame_size;
1257 int M = st->M, N = st->window_size, n = st->frame_size, i, j;
1264 st->wtmp2[i] = EXTRACT16(PSHR32(st->W[j*N+i],16+NORMALIZE_SCALEDOWN));
1265 spx_ifft(st->fft_table, st->wtmp2, st->wtmp);
1267 spx_ifft(st->fft_table, &st->W[j*N], st->wtmp);
1270 filt[j*n+i] = PSHR32(MULT16_16(32767,st->wtmp[i]), WEIGHT_SHIFT-NORMALIZE_SCALEDOWN);