memset(coef, 0, sizeof(coef));  if(dbGain > max_dbGain) {dbGain = max_dbGain;}  S = 1; 
void calc_highshelf_coefs(int32* coef,int32 cutoff_freq,FLOAT_T dbGain,int32 rate){	FLOAT_T a0,a1,a2,b0,b1,b2,omega,sn,cs,A,beta,S;	FLOAT_T max_dbGain = 12;	memset(coef, 0, sizeof(coef));	if(dbGain > max_dbGain) {dbGain = max_dbGain;}	S = 1;	A = pow(10, dbGain / 40);	omega = 2.0 * M_PI * (FLOAT_T)cutoff_freq / (FLOAT_T)rate;	sn = sin(omega);	cs = cos(omega);	beta = sqrt((A * A + 1.0) / S - (A - 1) * (A - 1));	a0 = 1.0 /((A + 1) - (A - 1) * cs + beta * sn);	a1 = -2 * ((A - 1) - (A + 1) * cs);	a2 = -((A + 1) - (A - 1) * cs - beta * sn);	b0 = A * ((A + 1) + (A - 1) * cs + beta * sn);	b1 = -2 * A * ((A - 1) + (A + 1) * cs);	b2 = A * ((A + 1) + (A - 1) * cs - beta * sn);	a1 *= a0;	a2 *= a0;	b1 *= a0;	b2 *= a0;	b0 *= a0;	coef[0] = TIM_FSCALE(a1, 24);	coef[1] = TIM_FSCALE(a2, 24);	coef[2] = TIM_FSCALE(b0, 24);	coef[3] = TIM_FSCALE(b1, 24);	coef[4] = TIM_FSCALE(b2, 24);}
4357	/local/tlutelli/issta_data/temp/c/2005_temp/2005/4357/982c069ceae585dd86c633cbadac26bf8e4dffb6/reverb.c/clean/timidity/reverb.c