62 	int	order,			// LP FIR filter order  in sinc_coef()  argument
69 	int half = order >> 1;  in sinc_coef()
70 	if (order & 1) {		// order is odd, center = half + 0.5  in sinc_coef()
76 	} else {	// order is even, center = half  in sinc_coef()
83 	for (; i <= order; i++)		// symmetric FIR  in sinc_coef()
84 		coef[i] = coef[order - i];  in sinc_coef()
93  * convolution of coef[order + 1] and data[length] up-sampled by a factor
100 	    int		order,		// filter order  in poly_conv()  argument
105 	if ((order < 0) || (inc_coef < 1) || (length < 1))  in poly_conv()
109 		double *coef_end = coef + order;  in poly_conv()
155 	order = (fold << 4) - 2;		// filter order = fold * 16 - 2  in ResampleFilter()
156 	coef = new double[order + 1];  in ResampleFilter()
157 	sinc_coef(fold, order, coef);		// required bandwidth = PI/fold  in ResampleFilter()
159 	if (up > 1) {				// need (order/up) states  in ResampleFilter()
160 		num_state = (order + up  - 1) / up;  in ResampleFilter()
165 		num_state = order;  in ResampleFilter()
166 		state = new double[order];	// need order states  in ResampleFilter()
169 	delay = (order + 1) >> 1;	// assuming symmetric FIR  in ResampleFilter()
197 	int init_size = (size <= order)? size : order;  in decimate_noadjust()
200 		    convolve(coef + i + 1, state + i, order - i));  in decimate_noadjust()
202 		*out_ptr++ = double2short(convolve(coef, in_buf + i - order,  in decimate_noadjust()
203 		    order + 1));  in decimate_noadjust()
265 		    order - coef_offset, up, state, num_state));  in interpolate_noadjust()
272 			    order - j, up, in_buf, i) + poly_conv(  in interpolate_noadjust()
273 			    coef + coef_offset, order - coef_offset, up, state,  in interpolate_noadjust()
282 		    order - j, up, in_buf, size) + poly_conv(  in interpolate_noadjust()
283 		    coef + coef_offset, order - coef_offset, up, state,  in interpolate_noadjust()
356 		    order - coef_offset, up, state, num_state));  in filter_noadjust()
364 			    order - j, up, in_buf, i) + poly_conv(  in filter_noadjust()
365 			    coef + coef_offset, order - coef_offset, up, state,  in filter_noadjust()
374 		*out++ = double2short(up * (poly_conv(coef + j, order - j, up,  in filter_noadjust()
376 		    order - coef_offset, up, state, num_state)));  in filter_noadjust()