void
pwf_zresp(real input[], real output[])
{
   int j,k;
   real tmp;

   for (k=0; k<IDIM; k++) {
        tmp = input[k];
        for (j=LPCW-1; j>=1; j--) {
           input[k] += zirwfir[j] * pwf_z_coeff[j+1];
           zirwfir[j] = zirwfir[j-1];
        }
        input[k] += zirwfir[0] * pwf_z_coeff[1];
        zirwfir[0] = tmp;
        for (j=LPCW-1; j>=1; j--) {
            input[k] -= zirwiir[j] * pwf_p_coeff[j+1];
            zirwiir[j] = zirwiir[j-1];
        }
        output[k] = input[k] - zirwiir[0] * pwf_p_coeff[1];
        zirwiir[0] = output[k];
   }
}