// 132 tap linear phase fir decimation filter
//     decimation ratio:     MM = 2
//     stopband attenuation: R = -80 dB
//     transition band:      d = 10 % of passband
//     corrects droop of 7-th order comb filter
// usage in midas deck where y is output of modulator or digital recombining
// network in MASH:
// Y=linfir(x=comb(x=y,M=M/MM, k=7),h=hh,M=MM,skip=150,isodd=isodd);



MM = 2;  isodd = 0;
hh = [  0.0004029, 0.0026467, 0.0071896, 0.0112980, 0.0102331, 0.0032328,
       -0.0035446,-0.0034787, 0.0016088, 0.0035158,-0.0006136,-0.0041240,
       -0.0013475, 0.0034686, 0.0033146,-0.0002643,-0.0008125, 0.0012455,
        0.0000936,-0.0033561,-0.0019749, 0.0036525, 0.0039873,-0.0028069,
       -0.0059187,-0.0004777, 0.0026494,-0.0036939,-0.0078212, 0.0008294,
        0.0096870, 0.0023063,-0.0106535,-0.0066344, 0.0084776, 0.0076155,
       -0.0093035,-0.0118142, 0.0071394, 0.0144330,-0.0056618,-0.0188866,
        0.0008812, 0.0206897, 0.0028539,-0.0239463,-0.0098294, 0.0242239,
        0.0160393,-0.0251117,-0.0255312, 0.0213668, 0.0313602,-0.0229360,
       -0.0473136, 0.0154864, 0.0656271, 0.0036599,-0.0767571,-0.0217951,
        0.0999542, 0.0634430,-0.1283444,-0.1589623, 0.1579680, 0.5341726 ];