/**************************************************************************
*
* ROUTINE
*               confg_cb
*
* FUNCTION
*               computes initial states (direct form filters) 
*               for update filter's memory 
*
**************************************************************************/

#include "ccsub.h"
extern int idb, no;
extern float bb[MAXNP+1], e0[MAXLP], fc[MAXNO+1], gamma2;
extern float fc_3[MAXNO+1];
confg_cb(s,l,d1,d2,d3,d4)
int l;
float s[], d1[],d2[],d3[],d4[];
{
  float fctemp[MAXNO+1];
  int i;

  setr(MAXNO+1, 0.0, fctemp);

     for (i = idb-l; i < idb; i++)
            d1[i] = e0[i-149] += d1[i];

    polefilt(fc, no, d2, e0, l);

  for (i = 0; i < l; i++)
    e0[i] = s[i] - e0[i];

    zerofilt(fc,no,d3,e0,l);
    bwexp(gamma2, fc, fctemp, no);
    polefilt(fctemp, no, d4, e0, l);
/*
    bwexp(gamma2, fc_3, fctemp, no);
    polefilt_3rd(fctemp, 3, d4, e0, l);
*/
}